感謝您的回覆。我意識到我想要一個barplot而不是一個直方圖。這裏是我想出瞭解決方案:
smoke=read.csv("SmokingEffect.csv",header=TRUE)
smokes=subset(smoke,select=c(Age,Smoke,FEV))
library(plyr)
smokesmeans <- ddply(smokes, c("Age","Smoke"), summarize, mean=mean(FEV),
sem=sd(FEV)/sqrt(length(FEV)))
smokesmeans <- transform(smokesmeans, lower=mean-sem, upper=mean+sem)
smokesmeans[,2] <- sapply(smokesmeans[,2], as.character)
library(ggplot2)
plotation <- qplot(x=Age, y=mean, fill=Smoke, data=smokesmeans,
geom="bar",stat="identity",position="dodge",main="distribution of FEV",
ylab="mean FEV")
plotation <- plotation + geom_errorbar(aes(ymax=upper,
ymin=lower), position=position_dodge(0.9), data=smokesmeans)
png(myplot.png)
plotation
dev.off()
輸出看起來是這樣的:需要在你的問題
更多信息,得到一個很好的答案。你能否描述你的數據 - 是肺還是年齡分類?你可以創建一些模擬數據,並顯示你已經嘗試過的代碼。 – user20650 2014-10-26 21:23:41