set.seed(123) drug <- as.integer(rnorm(200, 2, 2)) drug[drug<0] <- 0 placebo <- as.integer(rnorm(200, 4, 3)) placebo[placebo<0] <- 0 duration <- c(drug, placebo) trt <- factor(rep(c('Drug','Placebo'), each=200)) ddrug <- density(drug) dplac <- density(placebo) #png('C:/Research/Graphics/Graphs4Display/webpages/classes/pages/images/boxplotbasic.png', # width=500, height=500) boxplot(duration~trt, col=c('blue','grey60'), horizontal=TRUE, main='horizontal=TRUE', varwidth=TRUE, notch=TRUE) #dev.off()