Dear Richard and Duncan,
your suggestions are absolutely serving what I need. But I would like to
see x-axis to be up to 30 instead of 20. Do you have any suggestion on that?
Many thanks for your kind help.
Regards,
Jamil.
On 12 August 2014 01:22, Duncan Mackay <dulca...@bigpond.com> wrote:
> Hi
>
> If you want a 1 package and 1 function approach try this
>
> xyplot(conc ~ time | factor(subject, levels = c(2,1,3)), data = data.d,
> par.settings = list(strip.background = list(col = "transparent")),
> layout = c(3,1),
> aspect = 1,
> type = c("b","g"),
> scales = list(alternating = FALSE),
> panel = function(x,y,...){
>
> panel.xyplot(x,y,...)
>
> # f1<-function(x,v,cl,t)
> # (x/v)*exp(-(cl/v)*t) f1(0.5,0.5,0.06,t),
> panel.curve((0.5/0.5)*exp(-(0.06/0.5)*x),0,30)
>
> }
> )
>
> # par.settings ... if you are publishing show text better
> # with factor if you want 1:3 omit the levels
> # has advantage of doing more things than in groupedData as Doug Bates has
> said
>
> Regards
>
> Duncan Mackay
> Department of Agronomy and Soil Science
> University of New England
> Armidale NSW 2351
> Email: home: mac...@northnet.com.au
>
> -----Original Message-----
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On
> Behalf Of Naser Jamil
> Sent: Monday, 11 August 2014 19:06
> To: R help
> Subject: [R] Superimposing graphs
>
> Dear R-user,
> May I seek your help to sort out a little problem. I have the following
> codes
> to draw two graphs. I want to superimpose the second one on each of the
> first one.
>
> ########################################
>
> library(nlme)
> subject<-c(1,1,1,2,2,2,3,3,3)
> time<-c(0.0,5.4,21.0,0.0,5.4,21.0,0.0,5.4,21.0)
> con.cohort<-c(1.10971703,0.54535512,0.07176724,0.75912539,0.47825282,
> 0.10593292,1.20808375,0.47638394,0.02808967)
>
> data.d=data.frame(subject=subject,time=time,conc=con.cohort)
> grouped.data<-groupedData(formula=conc~time | subject, data =data.d)
>
> plot(grouped.data)
>
> ##########################################
>
> f1<-function(x,v,cl,t) {
> (x/v)*exp(-(cl/v)*t)
> }
> t<-seq(0,30, .01)
> plot(t,f1(0.5,0.5,0.06,t),type="l",pch=18, ylim=c(), xlab="time",
> ylab="conc")
>
>
> ###########################################
>
> Any suggestion will really be helpful.
>
>
> Regards,
>
> Jamil.
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.
>
>
[[alternative HTML version deleted]]
______________________________________________
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.
