运用PowerShell计算行×列表χ2检验的论文.docVIP

　　运用PowerShell计算行×列表χ2检验的论文 【摘要】 微软已经推出功能强大的poicrosoft have designed po. poathematical putation ability and can manage multidimension array expediently. so simple po can be used to calculate contingency table chi-square test in medical statistics. 　　key ei、、ado、.，管理员将能够在命令行的环境中使用和.类库中的全部功能，可以说功能十分强大。用户不需要另外安装软件，仅通过pop += $rc［$i］［$j］ } $rt = $rt + $temp $n += $temp $temp = 0 } 　　# 计算各列之和并赋值到一维数组$ct中 　　$ct = @（) 　　for （$j=0; $j -lt $c; $j++) { for （$i=0; $i -lt $r; $i++) { $temp += $rc［$i］［$j］ } $ct = $ct + $temp $temp = 0 } 　　# 计算卡方 　　for （$i=0; $i -lt $r; $i++) { for （$j=0; $j -lt $c; $j++) { $temp += $rc［$i］［$j］*$rc［$i］［$j］ / （$rt［$i］*$ct［$j］) } } 　　$x2 = $n * （$temp - 1) 　　 软件进行统计。表1 术前住院天数对住院总天数的影响（略） 此例原文中计算结果为χ2 = 97.54 , plt;0.01 ，不同术前住院天数间、 住院总天数有显著性差异，说明术前住院天数直接影响住院。 现在我们在powershell脚本来计算，首先在控制台界面输入这个行列表的数据： ps c:＼gt; $rc = （115,220,105,72),（3,52,35,3),（0,3,20,20) 这里简单解释一下，这种形式的输入是powershell创建二维数组的一种方式，行列表中每一行的数据放在一对括号内，一行内的数据用逗号分隔，行与行之间的数据也用逗号隔开。通过这一行的输入，我们已经把行列表所有数据赋值到变量$rc中。 然后运行我们已经编写好的计算脚本kf.ps1： ps c:＼gt; c:＼kf.ps1 您输入的是 3行* 4列表! n=648 卡方=97.5386448154744 此例的计算结果说明powershell脚本编程的计算结果准确。 3.2 例2 扭挫伤患者708 人和腰肌劳损患者347人，都用磁场疗法治疗，有效（显效与好转) 人数分别为673 人和312 人，无效例数分别为35人和35人，有效率分别为95106 % 和89191 %。问两组有效率有无差别? ［3］表2 两种疾病用磁场疗法的疗效（略） 　　检验假设h0 :两组总体有效率相等。 备择假设h1 :两组总体有效率不相等。 检验水准:双侧α= 0.05，计算结果为χ2 = 9.94，plt;0.01，拒绝h0 ,认为磁场疗法对两种疾病的有效率不相同。 我们通过powershell脚本验证此例的计算结果，首先在控制台界面输入数据： ps c:＼gt; $rc = （673,35),（312,35) 运行计算脚本 ps c:＼gt; c:＼kf.ps1 您输入的是 2行* 2列表! n=1055 卡方=9.9427369384943 本例是四格表的卡方统计问题，说明本研究的powershell程序同样适用于四格表卡方统计。通过以上二例的计算，我们认为通过powershell脚本计算行列表卡方是一种方便易用的方法。 　　4 讨论 以下几点是对本研究中脚本程序的总结说明： ① 如果编写的kf.ps1脚本不能运行，那么应该开启ps1脚本的运行权限，在控制台输入set-executionpolicy unrestricted 即可。微软从安全角度考虑关闭了脚本运行权限。 ② 微软操作系统要求运行ps1脚本时要键入完整的路径，在本文中的完整路径是c:＼kf.ps1，而不能只键入文件名kf.ps1。 ③ 不少医学工作者使用excel做卡方计算，缺点很明显，在excel表格中对行列数的限制是预先设定的，而本研究的脚本程序没有行列数的限制，可以灵活计算多行多列的表格卡方。 ④ 本研究的脚本程序只是根据公式完成了计算行列表卡方最基本的功能，在输入数据的时候请保持数据格式的正确性??因为程序对录入数据的格式不作检查。另外如果需要校正检验，可以对脚本适