统计学里Non-parametric Statistics课或应用生物统计学课都会讲到类似的例子,出现这样分组A都高但在总样本中却B高的原因是不平衡的样本分布, 在小肾结石组, 虽然A的成功率很高, 可是总样本数少, 那么他的成功率在总样本中的权重则低了。你把小肾结石的A跟大肾结石的B成功率比, 把大肾结石的A跟小肾结石的B比, 他们样本数接近了。但你可以看到小肾结石的B成功率87%》大肾结石的A的成功率,以为样本大, 所以这87%占了很大分量。 利用Cochran-Mantel-Haenszel Test, 你可以ACCOUNT这样的情况。 我做的结果是B的效果好于A(总结论), 但是这个差别在统计中并不显著(p=0.055),或者说刚刚搭上显著的边。 老新有SAS的话, 运行下面的代码就可以得到结果了。在SAS的结果中找到(Col1 Risk), 就是你要的结果。 *************************************************************************; data kidneyStone;www.ddhw.cominput drug $ Effect Stone count; datalines ; A 1 1 83 A 0 1 4 B 1 1 234 B 0 1 36 A 1 2 192 A 0 2 71 B 1 2 55 B 0 2 25 ;www.ddhw.com run ;www.ddhw.comproc freq data = kidneyStone;www.ddhw.comtable stone* Drug * Effect /cmh ;www.ddhw.comweight count ;www.ddhw.comrun ;
本贴由[NuBHH]最后编辑于:2012-12-7 8:21:8 |