R语言绘制箱线图

1.数据处理,把数据处理成箱图所需要的形式:

a1<-data.frame(Average_Degree=c(13,4,2,12.15,6.667,10.67,2.286,1,1,5,11.2,11.5,11.3,11.5,11.5),Method=rep("DYNMOGA",times=15))
a2<-data.frame(Average_Degree=c(11.57,8.2,11.88,3.5,5.111,9.579,4.75,2,8,2.5,12.53,14.88,6.75,18,9,8.833,12.53,5.75,16.35,3,17.2,5,2,11.47,14.12,10,11.13,12.43,6.75,4.222,2,7.4,5,14.12,17,15.79,8.667,11.16,6,8.211,1,4.118,17.7,14.47,14.17,11.69,6.2,9.667,10.15,18,9.077,4.667,7.2,16.18,10.83,7.636,1.333,6,1,11.33,16.22,4,12.86,11.07,11.54,12.75,10.75,13.47,13.07,1,1,7.273,15.65,13.91,3,8,1,3,12,4.909,9.444,6,1,9.667,10.44,13.13,2,9,4.444,5,3.2,1,12.86,10.29,14.73,3.75,6.25,10.92,16.9,4.857,1,11.54,6,6.889,2,8,12.27,4.667,16.65,13.56,2,1,6.222,12.75,1.333,3,5.143,8,3.5,11.62,1,2.571,13.88,9.833,10.67,2,14.38,8,1,12.93,5.429,8.4,12.52,9.538,13.26,10.07,18,3,14,7,11,13.64,4.667,1,8.6,9.818,11.81,4,8,5.429,1,16.19,4,9.538,5.714,13.07,15.26,5,9.429,7.538,6.8,11.38,10.71,9.429,3.556,16.63,4.667,7,18,7.4,10.27,3.333,1,8,10.67,12.63,2,11.8,17.04,8.857,6.5,9,13,14.5,16.92,5.143,10,13,10.5,8.545,4,2,7.4,7.8,3.6,3,8.909),Method=rep("DrNMF",times=197))
a3<-data.frame(Average_Degree=c(5.536,3.25,6.341,0.6667,2.222,8.727,3.448,1.333,3.533,1.4,6.214,5.787,2.296,2.75,1.333,0.6667,3.071,4.784,9.567,4.196,2.4,3,4,2,3.385,6.714,8.104,0.4706,1.882,1.538,5.625,2.941,6.542,1.917,6.361,7.27,2.833,5.579,0.5,3.238,3.3,0.6,1.375,2.077,4.952,4.296,3.529,5.136,6.885,2,1.1,2,2,6.541,1.556,4.813,2.182,8.381,0.6667,3.4,1.765,6.48,4.731,0.8889,4.69,1.286,3,3.471,1.8,3.316,6.791,8.133,2.133,1.714,1.5,4.545,2.75,6.387,7.843,2.353,6.935,5.02,3.667,6.186,3.556,5.606,1.889,4.118,6.667,0.3333,2.824,0.6154,3,10.24,0.5714,3.6,4.5,4.488,4,2,2.667,6.935,2.857,4.375,0.75,1.5,4.969,3.875,5.333,1.455,1,1.714,6.05,3.818,5.458,4.167,0.6667,9.275,3.056,5.375,1.2,9.25,2.359,2.471,5.647,5.304,1.429,4.237,4.323,4.596,1.286,5.462,2.429,2.522,0.4444,0.5,4.208,1.091,2.632,5.429,8.333,7.6,1.25,1.6,2.667,1.176,2.313,8.621,3,5.333,6.9,2.6,6.973,1.083,5.269,1.833,3.867,4.839,3.538,6.258,1.429,6.1,5.514,3.375,4.118,2,1.667,2.444,6.8,3.806,5.182,2.88,2,0.4,1.091,0.5,1.8,1.455,3.067,9.05,1.867,5.059,1.143,2.4,3.152,4.364,1.25,6.571,2.778,2.2,7.729,3.125,4.444,4,6.491,3.385,1,3.25,5.463,4.222,3.023,3.067,9.091,4.941,3.688,1,6.643,1.125,5.45,7.194,5.294,0.9231,1.6,1.6,0.6667,1.176,2.933,2.182,6.806,2.08,10.28,4.182,2,4.194,4.489,4.816,8.442,6.875,10,2.452,3.4,1.412,2.667,3.784,4.8,3.037,6.421,6.923,6.52,7,2.167,4.75,5.424,7.026,3.75,7.121,2.909,2.2,0.6154,2.091,4.977,4.947,0.7273,7.532,0.7692,5.815,8.69,1.667),Method=rep("TANMF",times=258))
a4<-data.frame(Average_Degree=c(27.8,15.56,31.13,20.14,4,30.67,2,0.6667,33,0.6667,12.75,30.1,6.444,24.31,31.39,15,1,27.08,26.71,2.5,3,1,12.8,13,16,5.81,19.09,11,16.62,8,1,18.83,22.9,6,27.56,6,14.11,14.24,2.286,17,6.5,11.22,33.69,25.07,1,13,16.78,11,5.833,1,22.85,11.69,5,4,24.63,32.39,30.86,30.68,32.17,0.3333,3.2,14.88,6.5,3.556,19.69,12.71,26.67,6.154,24.15,25.08,15.13,7.111,17.9,5,32.07,12.13,3,1,8.6,7,12.22,1,23.68,2,24.38,2,32.45,8.286,14.71,34.8,28.05,2,15.89,16.62,14.5,18.43,0.5,1,2.5,6.333,10.62,30.1),Method=rep("PisCES",times=102))

a5<-data.frame(Average_Degree=c(4.25,8.476,13.89,6.8,14.09,11.8,5,1.2,15.26,14.11,7.667,6.667,7.143,7.714,6.833,16.57,12.92,11.47,4.2,9.429,2,4.25,17.17,11.33,13.31,1,6,0,1,9.875,8.182,5.6,3.429,15,4,11.2,13.54,0.6667,7.368,18.86,12.76,3.429,2.571,8.8,3,0,13,13.53,3,4,18.58,7,13.22,6.533,11.78,6.909,9.846,9.636,6,17.68,5.5,4.75,4.25,1.333,1.333,7.867,16.07,6.923,6.286,2.8,2.444,16.55,4.909,0,9,1.333,9.467,0.6667,15.79,12.25,7.444,18.11,12.71,17,15.05,7.714,1.5,3,17.2,17.24,4.286,8.615,7,6,9,1,8,14,3.25,3,12.11,16.07,3.333,5.25,0,7.429,13.92,16.26,2,15.6,6.444,0,7.6,2,0,6,4,14,6.25,1,9.789,16.67,9.273,4.667,14.33,9.333,12.82,9.375,13.29,14.93,8.706,16.4,11.4,5.8,6.308,6.286,1,5.111,0,8.381,8.308,17.85,10.42,15.33,2,2.333,7.846,7.273,2,12.74,12.42,13.13,2.5,0,14,3,7.6,14.87,1.333,17.8,7.875,17.16,5,15.92,5.556,5.6,13.71,13.67,12.38,13.89,4.25,15.74,3.667,2.8,5.667),Method=rep("Semi_NMF",times=175))

输入完数据,把数据按照行合并为一个数据集:

a<-rbind(a3,a1,a2,a4,a5)

最后简单的画图:

my_comparisons <- list(c("TANMF","DYNMOGA"), c("TANMF", "DrNMF"),c("TANMF", "PisCES"),c("TANMF","Semi_NMF"))

ggboxplot(a, x="Method", y="Average_Degree",color="Method", palette = "jco")+ stat_compare_means(method="t.test",comparisons = my_comparisons)

涉及复杂参数的调整:

最后如果涉及多个图片的时候需要使用Rmisc包的multiplot函数:

p1<- ggboxplot

p2<- ggboxplot

p3<- ggboxplot

p4<- ggboxplot

multiplot(p1,p2,p3,p4,cols=2)