R语言画棒状图(bar chart)和误差棒(error bar)

假设我们现在有CC,CG,GG三种基因型及三种基因型对应的表型,我们现在想要画出不同的基因型对应表型的棒状图及误差棒。整个命令最重要的就是最后一句了,用arrows函数画误差棒。用到的R语言如下:

data<-read.csv("E:/model/data.csv",sep=" ",header=T)#导入数据data
mean_CC<-mean(data[,1])#计算CC基因型对应的表型的平均值
mean_GG<-mean(data[,2])
mean_GC<-mean(data[,3])
install.packages("plotrix")#安装plotrix包计算标准误
library("plotrix")
SE_CC<-std.error(data[,1])#计算CC基因型对应表型平均值的标准误
SE_GG<-std.error(data[,2])
SE_GC<-std.error(data[,3])
density<-c(mean_CC,mean_GC,mean_GG)#画分布图
snp_name<-c("CC","GC","GG")
barCenters<-barplot(beside=T,density,snp_name,border="black",col=c("purple","green3","blue"),density=c(5,5,5),width=c(2,2,2),space=c(1.5,1.5,1.5),ylim=c(0,10))#beside:逻辑参数。如果FALSE,那么将绘画堆叠式的条形;如果是T
  • 0
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值