使用多个分组执行配对T检验

Question

我试图在两个叫做n和t的状态之间进行t检验。这些种群通过1至4次测量。每个时间点然后属于它们之上的另外两个组。例如：

组 - ＆gt;子组 - ＆gt;时间点 - ＆gt;状态 - ＆gt; Var_measured

 Var_measured<-sample(248)

 State<-rep(c("t","n"), 124)

 Time<-as.factor(sample(1:4, 62,replace = TRUE))

 Sub_group<-sample(c("Met","Tet","Pep","Cep"),62, replace = TRUE)

 Group<-sample(c("AA","Gly"),62, replace = TRUE)

 df<-data.frame(Var_measured,State,Time,Sub_group,Group)

因此，从这个数据框中，我需要一个来自状态n和t的每个时间点的每个子组之间的t检验的p值列表。这最好用一个例子可视化：

p<-ggplot(df, aes(x= Time, y= Var_measured, fill= State))+
    geom_boxplot()+
    # coord_cartesian(ylim = c(0, 150000))+
    theme(axis.text.x=element_text(angle=90)) +
    facet_grid(Group ~ Sub_group, scales = "free_y")
p

我最终希望能够在箱形图上显示有意义的标记（如下面的“无可挑剔的”油漆作业所示）。任何帮助将不胜感激。

最佳， 中号