文件名称:MyFisher22:一个非常紧凑的例程,用于 Fisher 在 2x2 矩阵上的精确检验; 功效和样本量计算-matlab开发
文件大小:5KB
文件格式:ZIP
更新时间:2024-06-18 09:46:49
matlab
Fisher 对 2x2 列联表的精确检验允许在以下情况下计算精确概率:由于小区频率小,更快速的正态近似和卡方计算可能不准确。 Fisher 精确检验涉及计算多个阶乘以获得观察到的和每个更极端的表的概率。 阶乘增长很快,因此有必要使用阶乘的对数。 这个计算在 Matlab 中非常容易,因为 x!=gamma(x+1) 和 log(x!)=gammaln(x+1)。 我多次重写了这个函数:现在完全矢量化、预分配、使用递归关系对 2x2 矩阵进行 Fisher 精确测试以及使用对数大大加快了执行速度。 它比之前提交的 Fisherextest 更快。 事实上,我执行了这个测试,比较了两个脚本的核心(删除了输入错误检查、显示结果和计算能力的代码)。 X=[70 30; 29 80](要评估 100 个表) 次=零(1,1000); 对于 I=1:1000,抽动; myfisher22(X);
【文件预览】:
myfisher22.zip