That's perfect! Many thanks.


On 12 August 2014 14:32, Richard M. Heiberger <rmh@temple.edu> wrote:

> Yes, use xlim=c(0, 30) in your definition of P1
>
> On Tue, Aug 12, 2014 at 7:26 AM, Naser Jamil <jamilnaser79@gmail.com>
> wrote:
> > 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 <dulcalma@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: mackay@northnet.com.au
> >>
> >> -----Original Message-----
> >> From: r-help-bounces@r-project.org [mailto:r-help-bounces@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.
>

	[[alternative HTML version deleted]]

