stata单因素方差分析两两比较的统计量问题
师妹毕业，做得单因素方差分析两两比较，用的spss的snk两两比较，但是把关的统计老师要q值，软件没有给出，我原来用stata，但是stata也没有给出统计量，如果用贝氏校正的话，能否自己编程，计算贝氏的t值呀。

关于贝氏两两比较的原理，从网上搜到了一点文字，贴到这里来：

1．Bonferroni法

 Bonferroni提出，如果在α水准上进行c次假设检验，当无效假设为真时，至少有一次拒绝无效假设的累积I类错误概率α’不超过c×a，即有不等式α’< c×a。因此可以重新选择I类错误概率水准α，以便使累积I类错误概率α’＝0.5。根据Bonferroni不等式可得到要重新选择的α水准为α＝（0.05/3）＝0.016。只有当t检验的I类错误概率等于或小于0.016时才能拒绝无效假设。这样当无效假设为真时，其累积I类错误概率不超过α’＝0.05。

 用Bonferroni法进行多个平均值之间的两两比较时所用的t检验公式为：

 t＝（Yi－Yh）/Se

 式中的分母Se为两平均值之差的标准误，计算公式为：

 Se＝[ MS组内×（1/ni＋1/nh）]1/2

 Yi、Yh及ni、nh分别是两个比较组的平均值及观察例数。

 当比较次数不多时，Bonferroni法的效果较好。但当比较次数较多（例如在10次以上）时，则由于其检验水准选择得过低，结论偏于保守。