plot xaxp issue

together <- rbind(boy, girl)
together$sex <- factor(rep(c("boy", "girl"), c(8,7)))

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-bounces at r-project.org [mailto:r-help-bounces at r-
project.org] On Behalf Of Elaine Kuo
Sent: Monday, January 07, 2013 9:27 AM
To: r-help at 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)