error message in plot(aov-object) -- repost
Pascal A. Niklaus wrote:
Hi all, I posted this question several days ago, but did not get any answer until now. Since I still have no clue about the source of this error message, I repost a description of the problem including some code: A student at our institute fitted an aov model, and got the following error message:
> plot(p.aov)
Hit <Return> to see next plot:
Hit <Return> to see next plot:
Error in plot.window(xlim, ylim, log, asp, ...) :
need finite ylim values
The second plot was not produced. I then tried to produced the plot
manually, which worked, and checked the range of the residuals:
> qqnorm(resid(p.aov))
> range(resid(p.aov))
[1] -0.2428688 0.2020649 Her data set and the model are as follows: Data: http://www.bot.unibas.ch/~pascal/ParzJan.csv Model (I don't understand the model, but that's not the issue here) p <- read.csv("ParzJan.csv",header=T) p.aov <-aov(terms(JI~SP+Ex+Neig+pH+K/G+H+D+tN+br+aN:br+tN:D +br:D+aN:br:D+br:tN+aN:br:D + aN:br:tN+ K:D+ K:tN+ K:br +K:aN:br+ G:K:D+ G:K:tN+ G:K:br+G:K:aN:br+ D:H+ tN:H+ br:H + aN:br:H, keep.order=T), data=p) Any help regarding the source of this plot() error message is appreciated.
Without the data nobody was able to reproduce your problem ..
plot.lm() produces the error message, in particular the lines
hii <- lm.influence(x, do.coef = FALSE)$hat
which appears to be 1 in some cases, and
rs <- r.w/(s * sqrt(1 - hii))
(you get a division by zero and the outcome is Inf, hence it's not
possible to scale the plot any more).
I'd call it a bug, but I don't know what the bug really is. I'll think
about it...
Uwe Ligges
Pascal P.S.: She used R 1.7.1 on a Mac, but the problem was reproducible on R 1.8.1 under Linux.
______________________________________________ R-help at stat.math.ethz.ch mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html