I just realized that I was actually using a different random number
generator, could that be a valid reason for the discrepancy?
The code should be:
RNGkind("L'Ecuyer")
set.seed(883)
x <- rnorm(100)
x %*% x - sum(x^2) # equal to 1.421085e-14
Regards,
Alexis Sarda.
On Tue, Sep 20, 2016 at
> Alexis Sarda
> on Tue, 20 Sep 2016 17:33:49 +0200 writes:
> I just realized that I was actually using a different random number
> generator, could that be a valid reason for the discrepancy?
> The code should be:
> RNGkind("L'Ecuyer")
>
> peter dalgaard
> on Fri, 16 Sep 2016 13:33:11 +0200 writes:
> On 16 Sep 2016, at 12:41 , Alexis Sarda wrote:
>> Hello,
>>
>> while testing the crossprod() function under Linux, I noticed the
following:
>>
>>
On 16 Sep 2016, at 12:41 , Alexis Sarda wrote:
> Hello,
>
> while testing the crossprod() function under Linux, I noticed the following:
>
> set.seed(883)
> x <- rnorm(100)
> x %*% x - sum(x^2) # equal to 1.421085e-14
>
> Is this difference normal? It seems to be
Hello,
while testing the crossprod() function under Linux, I noticed the following:
set.seed(883)
x <- rnorm(100)
x %*% x - sum(x^2) # equal to 1.421085e-14
Is this difference normal? It seems to be rather large for double precision.
Regards,
Alexis.
[[alternative HTML version