Hi
Instead of two plots with par(new = TRUE) try to put boys and girls together
(quite natural thing, they will be pleased 8-)
together <- rbind(boy, girl)
together$sex <- factor(rep(c("boy", "girl"), c(8,7)))
plot(together$body_length, together$body_weight, ....,
col=c("firebrick3","saddlebrown")[as.numeric(together$sex)], ....)
Regards
Petr
> -----Original Message-----
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of Elaine Kuo
> Sent: Monday, January 07, 2013 9:27 AM
> To: r-help@r-project.org
> Subject: [R] plot xaxp issue
>
> Hello,
>
> I have data of Body length and Body weight of 8 boys and 7 girls.
>
> I want to draw the plot of Body length (for X) and Body weight (for Y)
> based on sex.
> Then the two plots want to be overlapped for comparison.
>
> I used the code below but found the unit length of X axis of boy and
> girl plot are not the same.
> For instance, the length between 0 and 1 of boy plot is larger than
> that in girl plot.
> The same thing happened to Y axis as well.
> (In other words, though axap and yaxp were set to be the same, the
> display were not the same.)
>
> Please kindly advise correction of the code.
> Thank you.
>
> Elaine
>
> # plot code
> boy<-read.csv("H:/boy_data.csv",header=T)
> girl<-read.csv("H:/girl_data.csv",header=T)
> par(mai=c(1.03,1.03,0.4,0.4))
>
> plot(boy$body_length, boy$body_weight,
> xlab=" body_length (cm)",
> ylab=" body_weight ( kg )",
> xaxp=c(0,200,4),
> yaxp=c(0,100,4),
> type="p",
> pch=1,lwd=1.0,
> cex.lab=1.4, cex.axis=1.2,
> font.axis=2,
> cex=1.5,
> las=1,
> bty="l",col="firebrick3")
>
> boyline<-lm(body_weight ~ body_length, boy)
> summary(boyline)
> abline(boyline,col="firebrick3",lwd=2)
>
> #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> # graph
> par(mai=c(1.03,1.03,0.4,0.4))
>
> par(new=T)
>
> plot(girl$body_length, girl$body_weight,
> xlab=" body_length (cm)",
> ylab=" body_weight ( kg )",
> xaxp=c(0,200,4),
> yaxp=c(0,100,4),
> type="p",
> pch=1,lwd=1.0,
> cex.lab=1.4, cex.axis=1.2,
> font.axis=2,
> cex=1.5,
> las=1,
> bty="l",col="saddlebrown")
>
>
> girlline<-lm(body_weight~ body_length, girl)
> summary(girlline)
> abline(girlline,col="saddlebrown",lwd=2)
>
> ______________________________________________
> 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.
______________________________________________
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.
