Displaying 1 result from an estimated 1 matches for "body_status".
2013 Jan 12
0
two-group comparison in boxplot
...=2, body weight51-100
the fifth boxplot: sex =1, body weight >101
the sixth boxplot: sex =2, body weight >101)
Please kindly advise how to modify the code below for the above purpose.
Thank you.
Elaine
Code
library(HH)
# data input 1 (all data)
dataN<-read.csv("H:/body_status.csv",header=T, row.names=1)
dataN$sex <-factor(dataN$sex)
dim(dataN)
dataN[1,]
str(dataN)
sex.colors <- c("darkblue","red2")
sex.code <- c("Boy", "Girl")
levels(dataN$sex) <- sex.code
# graph
par(mai=...