[R] na.action in model.tables and TukeyHSD
ripley at stats.ox.ac.uk
Thu Feb 27 19:34:02 CET 2003
This ia already fixed in R-devel. The answer is the same: don't use
na.omit implicitly: use it explicitly.
On Wed, 26 Feb 2003, CG Pettersson wrote:
> Hello everybody!
> I use R 1.6.2 in Windows, and have a problem controlling the na.action.
> In a dataset with twelve trials, one of the trials lack any readings of the variable "STS.SH" (standing power at harvest)
> Fitting an aov() object with the call:
> led1t7sts.aov <- aov(STS.SH ~ Trial/Block + Treatment + Treatment:Trial, data = led1t7, na.action=na.exclude)
> seems to work as it produces an object with 10 df for the factor "Trial".
> But when I use model.tables or TukeyHSD on the object I get this:
> > model.tables(led1t7sts.aov, "means")
> Error in replications(paste("~", paste(names(tables), collapse = "+")), :
> na.action must be a function
> I have tried to use "na.action=na.exclude" inside the model.tables call as well, without any bettering.
> I can naturally cope with the problem by taking the whole trial away from the dataset, but it doesn´t feel very sophisticated...;-)
> (Prof. Ripley answered a similar question from me two weeks ago. The answer was good but didn´t work as the reason of the error was the same as this time: a whole
> trial with only na:s in it).
> CG Pettersson
> cg.pettersson at evp.slu.se
> R-help at stat.math.ethz.ch mailing list
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
More information about the R-help