Thanks a lot.
Please kindly indicate the meaning of the c(8,7).
> together <- rbind(boy, girl)
> together$sex <- factor(rep(c("boy", "girl"), c(8,7)))
Elaine
On Mon, Jan 7, 2013 at 4:55 PM, PIKAL Petr <petr.pi...@precheza.cz> wrote:
> 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.
