AI代码解释 stat.test2<-ToothGrowth%>%mutate(dose=as.factor(dose))%>%group_by(dose)%>%t_test(len~supp)%>%adjust_pvalue()%>%add_significance("p.adj")%>%add_xy_position(x="dose")stat.test3<-ToothGrowth%>%t_test(len~dose,p.adjust.method="bonferroni")%>%adjust_pvalue()%>%add_s...
remotes::install_github("csdaw/ggprism") 添加P值 GraphPad Prism 9 已经支持自动添加P值了,因此ggprism也增加了add_pvalue函数,添加P值和显著性标记。 这个函数其实是对stat_pvalue_manual函数的重写,stat_pvalue_manual又是基于ggsignif包写的,所以掌握了ggsignif,你就掌握了R语言添加显著性标记的技能。 基本...
#Addp-valuescomparinggroups#Specifythecomparisonsyouwantmy_comparisons<-list(c("0.5","1"),c("1","2"),c("0.5","2"))p+stat_compare_means(comparisons=my_comparisons)+#Addpairwisecomparisonsp-valuestat_compare_means(label.y=50)#Addglobalp-value 添加分布p值的带有抖动点的箱线图 里面有箱线...
library(survival) head(colon[, c(1:4)])#Fit survival curvesfit <- survfit(Surv(time, status)~adhere, data = colon) library(survminer) ggsurv <- ggsurvplot(fit, data = colon, palette = "jco", #jco palette pval = TRUE, pval.coord=c(500, 0.4), #Add p-value risk.table = TRUE...
adj y.position 1 0.5 len OJ VC 3.1697328 0.0064 T-test 0.0128 35 2 1.0 len OJ VC 4.0327696 0.0010 T-test 0.0030 35 3 2.0 len OJ VC -0.0461361 0.9600 T-test 0.9600 35 # Create a box plot and add the p-value p <- ggboxplot( ToothGrowth, x = "supp", y = "len", color = ...
ggdensity(df,x="weight",add = "mean",rug = TRUE,color = "sex",fill = "sex",palette= c("#00AFBB", "#E7B800")) 柱状图 } 1 gghistogram(df,x="weight",add = "mean",rug = TRUE,color = "sex",fill = "sex",palette ...
#数据重铸 upper_tri_melt <- melt(upper_tri, na.rm = TRUE) ggplot(data=upper_tri_melt, aes(Var1, y=Var2, fill=value))+ geom_tile(color="white")+ scale_fill_gradient2(low = "blue", high = "red", mid = "white", midpoint = 0, limit=c(-1, 1), space = "Lab", name=...
#Add p-values comparing groups#Specify the comparisons you wantmy_comparisons<-list( c("0.5","1"), c("1","2"), c("0.5","2") )p+stat_compare_means(comparisons=my_comparisons)+#Add pairwise comparisons p-valuestat_compare_means(label.y=50)#Add global p-value ...
P.Value <-c(a$P.Value) FC <-c(a$FC) df <-data.frame(P.Value, FC) df.G <-subset(df, log2(FC) <-1& P.Value <0.05)#define Green df.G <-cbind(df.G, rep(1, nrow(df.G))) colnames(df.G)[3] <-"Color" df.B <-subset(df, (log2(FC) >=-1& log2(FC) <=1) ...
gg1 <- gg + theme( # 设置标题大小,face="bold"字体加粗 plot.title=element_text(size=30, face="bold"), axis.text.x=element_text(size=15), axis.text.y=element_text(size=15), axis.title.x=element_text(size=25), axis.title.y=element_text(size=25)) + # add title and axis text...