On Tue, 2010-10-05 at 13:56 -0700, 1Rnwb wrote: > Hello List, > I have a longitudinal samples on multiple individuals, for initial analysis, > i am doing some plots to see how the variable changes with the age. In some > of the individuals one or two data points are missing, > IL1Ra Age.at.Sample.Collection Subject.ID > 6.575466 2.004106 00709-0 > NA 2.162902 00709-0 > 5.115421 3.312799 00709-0 > 6.119626 3.515400 00709-0 > 5.705648 4.156057 00709-0 > 7.298887 4.405201 00709-0 > 6.195063 4.577686 00709-0 > 6.574864 5.541409 00709-0 > When I run > cap<-paste(sub.id,cases2$group[1],sep=" ") > plot(cases2$IL1Ra~cases2$Age.at.Sample.Collection, ylab="eo_conc",xlab="Age > at Sample",pch=1, ylim=c(3,10), xlim=c(0,22),main=cap) > points(cases2$MIP1b~cases2$Age.at.Sample.Collection,pch=2,ylim=c(1.5,10),col="red") > points(cases2$MCP1~cases2$Age.at.Sample.Collection,,pch=3, > ylim=c(1.5,10),col="brown") > lines(predict(loess(cases2$IL1Ra~cases2$Age.at.Sample.Collection))~cases2$Age.at.Sample.Collection, > lwd=1, lty=1) > lines(predict(loess(cases2$MIP1b~cases2$Age.at.Sample.Collection))~cases2$Age.at.Sample.Collection, > col="red", lwd=1, lty=1) > lines(predict(loess(cases2$MCP1~cases2$Age.at.Sample.Collection))~cases2$Age.at.Sample.Collection,col="brown", > lwd=1, lty=1) > > I get error > Error in model.frame.default(formula = predict(loess(cases2$IL1Ra ~ > cases2$Age.at.Sample.Collection)) ~ : > variable lengths differ (found for 'cases2$Age.at.Sample.Collection') > In addition: There were 50 or more warnings (use warnings() to see the first > 50) > > I tried using na.action=na.omit, which is not working at all. > > I would appreciate some pointers in this regards to get these charts. > Thanks
That's the second time this has come up today. Try na.action = na.exclude instead of na.omit in you loess() call. na.exclude will remove the NA from the data used to fit the model, but ensure theye are put back in to the results extracted from the loess model using the predict() extractor function. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/ UK. WC1E 6BT. [w] http://www.freshwaters.org.uk %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% ______________________________________________ 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.
