[R] How to remove NAs and lme function

[Andrew Robinson]

Jen,

try

na.action = na.exclude

Andrew


On Wed, May 28, 2008 9:26 pm, Jen_mp3 wrote:
>
> I am working on a project to find a model for the concentration of
> dissolved
> oxygen in the river clyde. Ive fitted a linear mixed model as
> lme(DOW~Temperature+Salinity+Year+factor(Station)*factor(Depth),
> random~1|id), where id is an identifier of the day over 20 years defined
> as
> Day*10000 + Month*100 + (1900 - Year).
> Anyway, there are some NAs for the concentration of dissolved oxygen in
> the
> water so I know you add in na.action = na.omit and that omits the NAs so
> there are 9008 observations in the model, but it doesnt do it for the
> whole
> data set where there are 10965 including observations with NAs. I would
> like
> to plot the residuals from the model against the Salinity, Temperature and
> Year, but when I try, it seems to want to take the observations of these
> variables from the full data set and the residuals from the model which of
> course doesnt work. I have tried using
> data1 <- data[data$DOW != "NA",] on the whole data set but it doesnt work.
> How can I remove the NAs from a data set?
>
> --
> View this message in context:
> http://www.nabble.com/How-to-remove-NAs-and-lme-function-tp17510564p17510564.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>


Andrew Robinson
Senior Lecturer in Statistics                       Tel: +61-3-8344-6410
Department of Mathematics and Statistics            Fax: +61-3-8344 4599
University of Melbourne, VIC 3010 Australia
Email: a.robinson at ms.unimelb.edu.au    Website: http://www.ms.unimelb.edu.au