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@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.

Reply via email to