Dear Med venlig hilsen and all Thanks a lot for your attention and valuable codes. I produced the graph using the follwing code. I fitted mixed effects dominant height model. The grpah is showing plot-specific dominant height against Age.
aspdomH2<-groupedData(domH2~age|plotno,data=aspdomH2)
fm3cham.nlme<-nlme(domH2~cham(age,b0,b1,b2),
data=aspdomH22,
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())
windows()
plot(augpred(fm3cham.nlme)
Note: if you look the attched graph, the highlighted heading of the text of the
each plot (i.e. ASS/3, HASS/3.....) are too big front size. I want to reduce
the front size. Actually, this is my plot id. To overcome the problem I used
your suggested code, but I got the folllowing error
plot(augPred(fm3cham.nlme),strip = strip.custom(par.strip.text = list(cex =
0.5)))
Error in plot.augPred(augPred(fm3cham.nlme), strip =
strip.custom(par.strip.text = list(cex = 0.5))) :
could not find function "strip.custom"
how I will fix the problem?
Thanks for your help. Actually I am new in R
Best regards
Salam
________________________________________
From: Frede Aakmann Tøgersen <[email protected]>
Sent: Thursday, May 15, 2014 3:09 PM
To: Mir Salam; [email protected]
Subject: RE: needs help in augented plot
Dear Salam
So you apparently did not read the posting guide as you were supposed to when
joining this list????
Every email from the list has this footer:
> ______________________________________________
> [email protected] 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.
So why didn't you provide commented, minimal, self-contained, reproducible
code????
Well the man page for augPred has this example:
fm1 <- lme(Orthodont, random = ~1)
augFM1 <- augPred(fm1, length.out = 2, level = c(0,1))
>From this you can
plot(augFM1)
Is it something like this you want:
plot(augFM1, par.settings = list(par.main.text = list(cex = 0.5)), main =
"Testing")
????
Or this:
plot(augFM1, strip = strip.custom(par.strip.text = list(cex = 0.5)))
????
If the latter I will leave it as an exercise to you to change the default strip
height. Google for r lattice strip height
Have a nice day.
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance & Modeling
Technology & Service Solutions
T +45 9730 5135
M +45 2547 6050
[email protected]
http://www.vestas.com
Company reg. name: Vestas Wind Systems A/S
This e-mail is subject to our e-mail disclaimer statement.
Please refer to www.vestas.com/legal/notice
If you have received this e-mail in error please contact the sender.
> -----Original Message-----
> From: [email protected] [mailto:[email protected]]
> On Behalf Of Mir Salam
> Sent: 15. maj 2014 13:22
> To: [email protected]
> Subject: [R] needs help in augented plot
>
> Dear all,
>
> I have problem in augmented plot to reduce the text (font ) size in heading
> of each plot. I have 68 plots and in augmented plot, I have heading of 68
> plots
> in a panel. The headline of each plot coded with plot characteristic, but the
> text seems to be big. How I will reduce the text size?
>
> R code
>
> plot(augPred(fm3cham.nlme),cex=0.5))
>
>
> Best regards
> Salam
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> [email protected] 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.
________________________________________
From: Mir Salam
Sent: Thursday, May 15, 2014 9:37 PM
To: Frede Aakmann Tøgersen
Cc: [email protected]
Subject: RE: needs help in augented plot
Dear Med venlig hilsen
Thanks a lot for your attention and valuable codes. I produced the graph using
the follwing code.
aspdomH2<-groupedData(domH2~age|plotno,data=aspdomH2)
fm3cham.nlme<-nlme(domH2~cham(age,b0,b1,b2),
data=aspdomH22,
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())
windows()
plot(augpred(fm3cham.nlme)
Note: if you look the attched graph, the highlighted heading of the text of the
each plot (i.e. ASS/3, HASS/3.....) are too big front size. I want to reduce
the front size. Actually, this is my plot id. To overcome the problem I used
your suggested code, but I got the folllowing error
plot(augPred(fm3cham.nlme),strip = strip.custom(par.strip.text = list(cex =
0.5)))
Error in plot.augPred(augPred(fm3cham.nlme), strip =
strip.custom(par.strip.text = list(cex = 0.5))) :
could not find function "strip.custom"
how I will fix the problem?
Thanks for your help. Actually I am new in R
Best regards
Salam
________________________________________
From: Frede Aakmann Tøgersen <[email protected]>
Sent: Thursday, May 15, 2014 3:09 PM
To: Mir Salam; [email protected]
Subject: RE: needs help in augented plot
Dear Salam
So you apparently did not read the posting guide as you were supposed to when
joining this list????
Every email from the list has this footer:
> ______________________________________________
> [email protected] 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.
So why didn't you provide commented, minimal, self-contained, reproducible
code????
Well the man page for augPred has this example:
fm1 <- lme(Orthodont, random = ~1)
augFM1 <- augPred(fm1, length.out = 2, level = c(0,1))
>From this you can
plot(augFM1)
Is it something like this you want:
plot(augFM1, par.settings = list(par.main.text = list(cex = 0.5)), main =
"Testing")
????
Or this:
plot(augFM1, strip = strip.custom(par.strip.text = list(cex = 0.5)))
????
If the latter I will leave it as an exercise to you to change the default strip
height. Google for r lattice strip height
Have a nice day.
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance & Modeling
Technology & Service Solutions
T +45 9730 5135
M +45 2547 6050
[email protected]
http://www.vestas.com
Company reg. name: Vestas Wind Systems A/S
This e-mail is subject to our e-mail disclaimer statement.
Please refer to www.vestas.com/legal/notice
If you have received this e-mail in error please contact the sender.
> -----Original Message-----
> From: [email protected] [mailto:[email protected]]
> On Behalf Of Mir Salam
> Sent: 15. maj 2014 13:22
> To: [email protected]
> Subject: [R] needs help in augented plot
>
> Dear all,
>
> I have problem in augmented plot to reduce the text (font ) size in heading
> of each plot. I have 68 plots and in augmented plot, I have heading of 68
> plots
> in a panel. The headline of each plot coded with plot characteristic, but the
> text seems to be big. How I will reduce the text size?
>
> R code
>
> plot(augPred(fm3cham.nlme),cex=0.5))
>
>
> Best regards
> Salam
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> [email protected] 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.
Augmented plot.pdf
Description: Augmented plot.pdf
______________________________________________ [email protected] 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.
