咨询QQ：3025393450
周长符号
欢迎登陆官网：y0/datablog
r语言ggplot2误差棒图快速指南数据分析报告
来源：
给直方图和线图添加误差棒
准备数据
这里使用ToothGrowth 数据集。
风险代理律师library(ggplot2) df <- ToothGrowth df$do <- as.factor(df$do) head(df) ## len supp do ## 1 4.2 VC 0.5 ## 2 11.5 VC 0.5 ## 3 7.3 VC 0.5 ## 4 5.8 VC 0.5 ## 5 6.4 VC 0.5 ## 6 10.0 VC 0.5
len ：牙齿长度
绵羊用英语怎么说do : 剂量(0.5, 1, 2) 单位是毫克
supp : 支持类型(VC or OJ)
竭尽全力的拼音
在下面的例子中，我们将绘制每组中牙齿长度的均值。标准差用来绘制图形中的误差棒。
首先，下面的帮助函数会用来计算每组中兴趣变量的均值和标准差：高中生个人简历
怎么申请退款咨询QQ：3025393450
欢迎登陆官网：y0/datablog
#+++++++++++++++++++++++++ # Function to calculate the mean and the standard deviation # for each group #+++++++++++++++++++++++++ # data : a data frame # varname : the name of a column containing the variable #to be summariezed # groupnames : vector of column names to be ud as # grouping variables data_summary <- function(data, varname, groupnames){ require(plyr) summary_func <- function(x, col){ c(mean = mean(x[[col]],
<=TRUE), sd = sd(x[[col]], na.rm=TRUE)) } data_sum<-ddply(data,
香积寺
groupnames, .fun=summary_func, varname) data_sum <- rename(data_sum, c("mean" = varname)) return(data_sum) }
统计数据
df2 <- data_summary(ToothGrowth, varname="len", groupnames=c("supp", "do")) # 把剂量转换为因子变量df2$do=as.factor(df2$do) head(df2) ## supp do len sd ## 1 OJ 0.5 13.23 4.459709 ## 2 OJ 1 22.70 3.910953 ## 3 OJ 2 26.06 2.655058 ## 4 VC 0.5 7.98 2.746634 ## 5 VC 1 16.77 2.515309 ## 6 VC 2 26.14 4.797731 有误差棒的直方图
函数 geom_errorbar()可以用来生成误差棒：
library(ggplot2) # Default bar plot p<- ggplot(df2, aes(x=do, y=len, fill=supp)) +
geom_bar(stat="identity", color="black", position=position_dodge()) +
geom_errorbar(aes(ymin=len-sd, ymax=len+sd), width=.2, position=position_dodge(.9))
print(p) # Finished bar plot p+labs(title="Tooth length per do", x="Do (mg)", y = "Length")+ theme_classic() + scale_fill_manual(values=c('#999999','#E69F00'))
咨询QQ：3025393450
欢迎登陆官网：y0/datablog
咨询QQ：3025393450
欢迎登陆官网：y0/datablog
美国一号
