运用powershell计算行×列表χ2检验

1/6运用POWERSHELL计算行列表2检验POWERSHELL是微软在06年底设计完成的新一代命令行外壳和脚本语言，用于进行系统管理和自动化处理。LOCALHOST从微软的发展报告中我们可以了解到，POWERSHELL将内置在新版的VISTA和WINDOWSSERVER008操作系统中，这说明POWERSHELL将会随着微软操作系统被广泛地使用，现在已经进入测试阶段。目前电脑用户使用较多的WINDOWSXPSP2也可以在NET的基础上安装使用POWERSHELL。POWERSHELL是一款基于对象并建立在NET框架之上，能够同时支持WMI、COM、ADO、NET，管理员将能够在命令行的环境中使用COM和NET类库中的全部功能，可以说功能十分强大。用户不需要另外安装软件，仅通过POWERSHELL进行简单的脚本编程就可以完成平常遇到的各种医学统计任务，一般来说没有必要安装庞大专业统计软件，也没有必要另外安装其他编程语言的环境。笔者曾经讨论过运用批处理的方法完成四格表卡方检验的计算1，由于批处理没有数组的功能，要想实现多行多列表格2检验是很困难的。POWERSHELL具有完善的数学计算的能力，支持多维数组，因此要实现行列表2检验是相当简单的。我们知道，行列表2检验基本计算公式为2/62N（ANRNC1，从编程计算的角度来说关键是二维数组的恰当运用。对二维数组ARRAY分解并构造出两个一维数组，分别代表行和列的数值，计算出每行的和NR以及每列的和NC，最后由公式计算最后的结果。我们尝试运用POWERSHELL脚本编程方法计算医学统计中经常遇到的行列表2检验。首先我们利用POWERSHELL对数组赋值的方法在控制台把表格数据直接赋值给数组RC，然后我们编写计算脚本，对数组RC进行分析并计算出2检验结果。脚本的源代码如下对输入的数据判断行数和列数RCRC0COUNTWRITEHOST“您输入的是R行C列表“计算各行之和并赋值到一维数组RT中，并计算表格全部数值之和NRT（FOR（I0ILTRIFOR（J0JLTCJ3/6TEMPRCIJRTRTTEMPNTEMPTEMP0计算各列之和并赋值到一维数组CT中CT（FOR（J0JLTCJFOR（I0ILTRITEMPRCIJCTCTTEMPTEMP0计算卡方FOR（I0ILTRIFOR（J0JLTCJ4/6TEMPRCIJRCIJ/（RTICTJXN（TEMP1WRITEHOST“NN卡方X“在控制台命令行输入C即可运行该脚本。测试程序1例1对64例消化道溃疡病人住院天数分析2某病案统计室对医院收治的64例消化道溃疡病人术前住院天数进行了分析，探讨术前住院天数对住院总天数的影响，病案统计室使用第四军医大学统计教研室编制的SPLM软件进行统计。表1术前住院天数对住院总天数的影响此例原文中计算结果为,P现在我们在POWERSHELL脚本来计算，首先在控制台界面输入这个行列表的数据PSCRC检验假设H0两组总体有效率相等。备择假设H1两组总体有效率不相等。5/6检验水准双侧，计算结果为94，P我们通过POWERSHELL脚本验证此例的计算结果，首先在控制台界面输入数据PSCRC（673,35,（312,35运行计算脚本PSCC您输入的是行列表N105卡方本例是四格表的卡方统计问题，说明本研究的POWERSHELL程序同样适用于四格表卡方统计。通过以上二例的计算，我们认为通过POWERSHELL脚本计算行列表卡方是一种方便易用的方法。以下几点是对本研究中脚本程序的总结说明如果编写的脚本不能运行，那么应该开启PS1脚本的运行权限，在控制台输入SETEXECUTIONPOLICYUNRESTRICTED即可。微软从安全角度考虑关闭了脚本运行权限。微软操作系统要求运行PS1脚本时要键入完整的路径，在本文中的完整路径是C，而不能只键入文件名。不少医学工作者使用EXCEL做卡方计算，缺点很明显，在EXCEL表格中对行列数的限制是预先设定的，6/6而本研究的脚本程序没有行列数的限制，可以灵活计算多行多列的表格卡方。本研究的脚本程序只是根据公式完成了计算行列表卡方最基本的功能，在输入数据的时候请保持数据格式的正确性，因为程序对录入数据的格式不作检查。另外如果需要校正检验，可以对脚本适当改动，读者可以通过对POWERSHELL的学习把程序继续完善。本研究的脚本程序在XPSPPOWERSHELL和VISTAPOWERSHELL的环境下运行通过。【参考文献】1王正凯,颜景红运用批处