Dear all, I need help to get different 68 plots specifc fitted curves in one plot with respective field data observations (age vs dominant height).
aspdomH2<-groupedData(domH2~age|plotno,data=aspdomH2) names(aspdomH2) plotno, age, origin, soilcharacter, domH2, plotno-different plot no. I have 68 plots age- every plot have from age 5 to 30 years origin- two, native aspen and hybrid aspen domH2<-dominant height soilcharacter-3, clay, silt and mold. both origin have different soil charcter #### then I fit model fm2cham.nlme<-nlme(domH2~cham(age,b0,b1,b2), data=aspdomH2, fixed = list(b0~1+origin+soilcharacter,b1~ 1,b2 ~ 1+origin+soilcharacter), random = b0+b2~1|plotno, start=c(b0=26.3387,0,0,0,b1=0.1065,b2=1.9453,0,0,0), weights=varPower(form = ~age, 0.5), correlation=corAR1()) #### parameter values Fixed effects: list(b0 ~ 1 + origin + soil character, b1 ~ 1, b2 ~ 1 + origin + soil character) Value b0.(Intercept) 21.081124 b0.origin1 7.735064 b0.soilcharactermold 10.689051 b0.soilcharactersilt 3.906585 b1 0.079035 b2.(Intercept) 1.616360 b2.origin1 -0.384421 b2.soilcharactermold 0.612285 b2.soilcharactersilt 0.527462 ##### I can easily get the augmented plot. I got different 68 plots specific curves. ####### Any body can help me how can I will get all 68 plots specific fitted curves in one plot with respective plot specific age and dominant height obervations? (x axis will represent age, y axis will represent dominant height and fitted curves of all 68 plots) [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.