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.
