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.