重复测量方差分析事后检验

Question

在spss中进行单因素重复测量方差分析后，得到事后检验的表格中时候均值差，请问怎么求 benferroni corrected的 pairwise comparison的t value？
或者在r语言进行pairwise.t.test后只能输出p value，有什么方法能输出t value 吗？

Answer 1

在SPSS中进行单因素重复测量方差分析后，得到事后检验表格中的均值差，可以通过如下步骤求得Benferroni校正后的pairwise comparison的t值：

执行单因素重复测量方差分析。
在SPSS中，打开"Analyze"菜单，选择"General Linear Model"，然后选择"Repeated Measures"，将变量添加到"Within-Subjects Factor"中，并选择适当的错误类型和校正方法，然后单击"OK"。
在"Repeated Measures: Options"对话框中勾选"Descriptive Statistics"和"Estimates of Effect Size"，然后单击"Continue"。
单击"Post Hoc"按钮，然后选择Benferroni校正方法。
单击"OK"，然后在结果表格中找到校正后的均值差和标准误。
计算t值：将均值差除以标准误，然后将结果与自由度和显著性水平一起输入到t分布表中，以获取临界值。然后，将t值与临界值进行比较以确定是否存在显著差异。
在R语言中执行pairwise.t.test后，会输出p值而不是t值。如果您想获得t值，可以使用以下公式将p值转换为t值：

t = q_t(1 - p/2, df)

其中，q_t是t分布的临界值函数，p是pairwise.t.test输出的p值，df是自由度。您可以使用qt函数来获取t分布的临界值。例如，以下代码将计算第一组和第二组之间的t值：

pvalue <- pairwise.t.test(x, g, p.adjust.method = "bonferroni")$p.value[1,2]
df <- n - 1
tvalue <- qt(1 - pvalue/2, df)

其中，x是数据向量，g是分组向量，n是每组的样本量。您可以根据需要修改代码以计算其他组之间的t值。