My usual advice on getting nonstandard F tests out of anova() is to fit the
models explicitly and compare.
So how about this?
fit1 <- lm(diff(extra,10) ~ 1, sleep)
fit0 <- update(fit1, ~ -1)
anova(fit0, fit1)
-pd
> On 26 Dec 2022, at 13:49 , Gabor Grothendieck wrote:
>
> Suppose we want to
FWIW, I suspect the problem might have something to do with:
> 1/1e-308
[1] 1e+308
> 1/1e-309
[1] Inf
> 1e-309 > 0
[1] TRUE
> 1e-324 > 0
[1] FALSE
That is 1e-309 starts being treated as zero by division, but not by comparison
(which happens at 1e-324). This is not quite the range of values you
Suppose we want to perform a paired test using the sleep data frame
with anova in R. Then this works and gives the same p value as
t.test(extra ~ group, sleep, paired = TRUE, var.equal = TRUE)
ones <- rep(1, 10)
anova(lm(diff(extra, 10) ~ ones + 0, sleep)
This gives output but does not
Dear R core team:
Sys.timezone() return wrong result on windows 10 (simple chinese version), it
return Asia/Taipei ,but The time zone is " China Standard Time (UTC+08:00)
Beijing, Chongqing, Hong Kong, Urumqi" .
after some digging , the code related this issue at
Thanks for all of your replies.
I would like to clarify that I am not expecting optim to return a good
result in these cases, as I am aware of the difficulties that come with the
numerical precision for such small numbers. I would be happy for optim to
take zero steps and return the initial value
