[求助]stata单因素方差分析两两比较的统计量问题

希望会编程解决问题的朋友帮忙

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次以上）时，则由于其检验水准选择得过低，结论偏于保守。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：单因素方差分析 Stata tata 计量问题 方差分析 Stata 单因素 方差分析 统计量

oneway命令不符合你的要求么？

oneway能够进行两两比较，但是完不成统计学老师提出的要求。贴图如下

[此贴子已经被作者于2008-4-12 9:43:25编辑过]

有没有高手相助。

愿意使用自己所有的金币的一半感谢帮助我的人。

凡是给出解决方案的朋友，证实确实是后，可以上传个文件，卖价是我所拥有现金的量的一半即可。第二个给出解决方案的，将是我当时现金量的一半

其实根据上面的公式很容易求出来啊，只不过你给出的结果，很难看出你的每一个处理中有多少个数据，因此没有办法给你算了，呵呵……

[此贴子已经被作者于2008-4-13 21:22:59编辑过]

其实，你用-ttest-命令不就行了？

. ttest lfr* if group==1 | group==2, by(group)

whgyu大哥，你好，用你的提示进行，算出的p值真的和两两比较式算出的p值相同，现在的问题是：

1：ttest是Two-sample t test with equal variances，这个t值和贝氏校正的t值是否相等？

2：你的ttest这样是用，和 Two-group mean-comparison test

        ttest varname [if] [in] , by(groupvar) [options1]有什么区别吗，应该没有吧？

其实Bonferroni correction的本质是改变显著性的threshold, 对test本身不会有什么影响。比如判断一个test的显著性用是否小于0.05，那么同时做两个test, Bonferroni就要求把阈值降到0.05/2=0.025. 而test statistic（就是t值)是不会改变的。