我已经为以下data set编写了两个不同患者群体之间的Anova计算代码
require(reshape2)
require(ggplot2)
require(ggpubr)
df.m <- melt(HC_BLVP_WM,id.var= "Subjects")
p <- ggboxplot(df.m , x= "Subjects" , y= "value",
color = "Subjects", palette = "jco",
add = "jitter",facet.by = "variable",
scales="free",short.panel.labs = TRUE)+
xlab("Subjects") + ylab("Volume mm³")+labs(colour = "Subjects")
p + stat_compare_means(label="p.signif",method = "anova",paired = FALSE)
当我运行代码并尝试使用ggboxplot
绘图时,我遇到以下警告消息
" Warning message:
In stringUnit(units) : reached elapsed time limit "
我无法在图像中看到任何Anova p值或统计显着性标记(*),如here所示,我在Ubuntu 16.04上使用Rstudio 1.0.153。有人可以告诉我问题在哪里,它是在我的代码中还是R或其中一个软件包中的错误?我以前从未遇到过这个问题;我能够成功绘制每个箱图上方的p值。
sessionInfo()的结果是
R version 3.4.2 (2017-09-28)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.3 LTS
Matrix products: default
BLAS: /usr/lib/libblas/libblas.so.3.6.0
LAPACK: /usr/lib/lapack/liblapack.so.3.6.0
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=de_DE.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=de_DE.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=de_DE.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=de_DE.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] reshape2_1.4.2 ggpubr_0.1.5 magrittr_1.5 ggplot2_2.2.1
loaded via a namespace (and not attached):
[1] Rcpp_0.12.13 dplyr_0.7.4 assertthat_0.2.0 grid_3.4.2
[5] plyr_1.8.4 R6_2.2.2 gtable_0.2.0 scales_0.5.0
[9] stringi_1.1.5 rlang_0.1.2 lazyeval_0.2.0 bindrcpp_0.2
[13] tools_3.4.2 stringr_1.2.0 glue_1.1.1 munsell_0.4.3
[17] compiler_3.4.2 pkgconfig_2.0.1 colorspace_1.3-2 bindr_0.1
[21] tibble_1.3.4