Try something like plot(x,y) box.dat <- boxplot(x=split(counts, dial), plot=FALSE) bxp(box.dat, add=TRUE, at=c(-1, 0, 1), show.names=FALSE)
-----Original Message----- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Karla Meurk Sent: Tuesday, June 29, 2004 21:46 PM To: [EMAIL PROTECTED] Subject: [R] funny plotting Hi, I just wanted to plot a boxplot with a nice curve going through it, I thought this would be a simple task but for some reason I can't get the two graphs on the same page accurately. Enclosed is the code showing the two plots seperately and together. I would have thought it should work if I could use boxplot() then plot() overlayed but it won't allow the command add=TRUE (which has worked for me in the past). Thanks Carla P.S. please excuse the clumsy code! #Section 2 Data Set particle dial<-rbind(-1, -1, 0, 0, 0, 0, 1, 1, 1) counts<-rbind(2, 3, 6, 7, 8, 9, 10, 12, 15) particle<-as.data.frame(cbind(dial,counts),row.names=NULL) names(particle)<-c("dial","counts") attach(particle) pois.particle<-glm(counts~dial,family=poisson) x<-seq(-2,2,length=20) y<-predict(pois.particle,data.frame(dial=x),type="response") #Overlaying plots x11() boxplot(counts~dial,main="Boxplot of counts for dial setting and poisson fit",ylim=c(0,25)) lines(x,y) #The seperate plots x11() boxplot(counts~dial,ylim=c(0,25)) x11() plot(x,y,ylim=c(0,25),type="l") ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
