[R] graphs

Mir Salam mir.salam at uef.fi
Sun Mar 8 22:56:53 CET 2015


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]]



More information about the R-help mailing list