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.
