群体间相关性分析及箱线图制作

数据准备与作图

getwd() #查看工作路径
read.csv("tiller analyse.csv",header = T) #读取数据表格
df<-read.csv("tiller analyse.csv",header = T) #普通赋值
library(reshape2)
df1<-melt(df)
head(df1)  #数据类型转换
library(ggplot2)
ggplot(df1,aes(x=variable,y=value))+geom_boxplot(aes(fill=variable))+theme_bw()+labs(x="variable",y="tiller number")+
  annotate("text",x=1,y=29,label="ab",size=5)+annotate("text",x=2,y=18,label="abc",size=5)+annotate("text",x=3,y=24,label="bc",size=5)+
  annotate("text",x=4,y=25,label="bc",size=5)+annotate("text",x=5,y=20,label="c",size=5)+annotate("text",x=6,y=23,label="abc",size=5)+
  annotate("text",x=7,y=25,label="a",size=5) # ggplot作图加标记等
<