最近有朋友问我给他们黄芪眼贴回复视力的数据画一个图,我就给他们画了一个boxplot和barplot的图,把代码分享到这里,有需要的小伙伴可以自取
数据整理
高度近视改善情况
H<-c(0.025,0.125,0.275,0.05,0.375,0.1,0.45,0.325,0.35,0.1,0.3,0.2, 0.18, 0.05, 0.225, 0.5, 0.065, 0.05, 0.15, 0.175, 0.125, 0.125)
中度近视改善
M<-c(0.2,0.1,0.4,0.4,0.1,0.25,0.05)
低度近视改善
L<-c(0.4,0.1,0.4,0.3,0.7,0.3,0.3,0.2,0.1)
数据表格整理
eye.data<-data.frame(c(H,M,L),c(rep(“1H”,length(H)),rep(“2M”,length(M)),rep(“3L”,length(L))))
names(eye.data)<-data.frame(“value”,“myopia”)
table(round(H,1))
画图成barplot
plot(table(round(H,1)+0.0),type=“h”,col=“red3”,xlim=c(0,0.5),ylim = c(0,8),xlab=“近视度数改善”,ylab=“人数”)
par(new=TRUE)
plot(table(round(M,1)+0.01),type=“h”,col=“green3”,xlim=c(0,0.5),ylim = c(0,8),xlab=“”,ylab=“”,axes = F)
par(new=TRUE)
plot(table(round(L,1)++0.02),type=“h”,col=“blue”,xlim=c(0,0.5),ylim = c(0,8),xlab=“”,ylab=“”,axes = F)
legend(“topright”,c(“高”,“中”,“低”),col=c(“red3”,“green3”,“blue3”),pch=c(1,2,3))
这里其实还有别的办法,我是好久不写代码了,手生又没空去找函数,大家可以根据实际情况改
library(ggplot2)
ggplot(data = eye.data) + geom_boxplot(aes(x = myopia, y = value ,color = factor(myopia))) +
geom_jitter(aes(x = myopia, y = value, color = factor(myopia)),position = position_jitterdodge())
这样就可以了,这里用Rstudio来操作的,所以没有保存成pdf或者png,可以参考以前的或者别的博文去保存