運(yùn)用PowerShell計(jì)算行×列表χ2檢驗(yàn)

1、運(yùn)用PowerShell盤算行列表2查驗(yàn)【摘要】微軟已經(jīng)推出成效強(qiáng)盛的pershell1.0，并決定在將來的新操縱體系中內(nèi)置pershell，pershell自己具有美滿的數(shù)學(xué)盤算和多維數(shù)組成效，只需通過簡樸足本編程就可以便利地盤算醫(yī)學(xué)統(tǒng)計(jì)中常常碰到的行列表2查驗(yàn)題目?！娟P(guān)鍵詞】pershelkeyrdsntingenytable;hi-squaretest;pershell1有關(guān)pershell的先容pershell是一款基于工具并創(chuàng)立在.net框架之上，可以或許同時支持i、ad、.net，辦理員將可以或許在下令行的情況中利用和.net類庫中的全部成效，可以說成效非常強(qiáng)盛。用戶不必要別的安裝

2、軟件，僅通過pershell舉行簡樸的足本編程就可以完成尋常碰到的種種醫(yī)學(xué)統(tǒng)計(jì)使命，一樣平常來說沒有需要安裝宏大專業(yè)統(tǒng)計(jì)軟件，也沒有需要別的安裝其他編程語言的情況。2運(yùn)用pershell足本編程盤算行列表卡方筆者曾經(jīng)討論過運(yùn)用批處置懲罰的要領(lǐng)完成四格表卡方查驗(yàn)的盤算1，由于批處置懲罰沒有數(shù)組的成效，要想實(shí)現(xiàn)多行多列表格2查驗(yàn)是很困難的。pershell具有美滿的數(shù)學(xué)盤算的本領(lǐng)，支持多維數(shù)組，因此要實(shí)現(xiàn)行列表2查驗(yàn)是相稱簡樸的。我們實(shí)驗(yàn)運(yùn)用pershell足本編程要領(lǐng)盤算醫(yī)學(xué)統(tǒng)計(jì)中常常碰到的行列表2查驗(yàn)。起首我們利用pershell對數(shù)組賦值的要領(lǐng)在操縱臺把表格數(shù)據(jù)直接賦值給數(shù)組$r，然后我們編

3、寫盤算足本kf.ps1，對數(shù)組$r舉行闡發(fā)并盤算出2查驗(yàn)效果。kf.ps1足本的源代碼如下：#對輸入的數(shù)據(jù)斷定行數(shù)和列數(shù)$r=$r.unt$=$r0.untrite-hst您輸入的是$r行*$列表!#盤算各行之和并賦值到一維數(shù)組$rt中，并盤算表格全部數(shù)值之和n$rt=)fr$i=0;$i-lt$r;$i+)fr$j=0;$j-lt$;$j+)$tep+=$r$i$j$rt=$rt+$tep$n+=$tep$tep=0#盤算各列之和并賦值到一維數(shù)組$t中$t=)fr$j=0;$j-lt$;$j+)fr$i=0;$i-lt$r;$i+)$tep+=$r$i$j$t=$t+$tep$tep=0#盤

4、算卡方fr$i=0;$i-lt$r;$i+)fr$j=0;$j-lt$;$j+)$tep+=$r$i$j*$r$i$j/$rt$i*$t$j)$x2=$n*$tep-1)rite-hstn=$n卡方=$x2在操縱臺下令行輸入:kf.ps1即可運(yùn)行該足本。3測試步伐3.1例1對648例消化道潰瘍病人住院天數(shù)闡發(fā)2某病案統(tǒng)計(jì)室對病院收治的648例消化道潰瘍病人術(shù)前住院天數(shù)舉行了闡發(fā)，探究術(shù)前住院天數(shù)對住院總天數(shù)的影響，病案統(tǒng)計(jì)室利用第四軍醫(yī)大學(xué)統(tǒng)計(jì)教研室體例的spl軟件舉行統(tǒng)計(jì)。表1術(shù)前住院天數(shù)對住院總天數(shù)的影響略此例原文中盤算效果為2=97.54,p0.01，差異術(shù)前住院天數(shù)間、住院總天數(shù)有明顯

5、性差異，說明術(shù)前住院天數(shù)直接影響住院。如今我們在pershell足原來盤算，起首在操縱臺界面輸入這個行列表的數(shù)據(jù)：ps:$r=115,220,105,72),3,52,35,3),0,3,20,20)這里簡樸說明一下，這種情勢的輸入是pershell創(chuàng)立二維數(shù)組的一種方法，行列表中每一行的數(shù)據(jù)放在一對括號內(nèi)，一行內(nèi)的數(shù)據(jù)用逗號分開，行與行之間的數(shù)據(jù)也用逗號離隔。通過這一行的輸入，我們已經(jīng)把行列表全部數(shù)據(jù)賦值到變量$r中。然后運(yùn)行我們已經(jīng)編寫好的盤算足本kf.ps1：ps:kf.ps1您輸入的是3行*4列表!n=648卡方=97.5386448154744此例的盤算效果說明pershell足本編

6、程的盤算效果正確。3.2例2扭挫傷患者708人和腰肌勞損患者347人，都用磁場療法治療，有用顯效與好轉(zhuǎn))人數(shù)別離為673人和312人，無效例數(shù)別離為35人和35人，有服從別離為95106%和89191%。問兩組有服從有無差異?3表2兩種疾病用磁場療法的療效略查驗(yàn)假設(shè)h0:兩組總體有服從相稱。備擇假設(shè)h1:兩組總體有服從不相稱。查驗(yàn)水準(zhǔn):雙側(cè)=0.05，盤算效果為2=9.94，p0.01，回絕h0,以為磁場療法對兩種疾病的有服從不雷同。我們通過pershell足本驗(yàn)證此例的盤算效果，起首在操縱臺界面輸入數(shù)據(jù)：ps:$r=673,35),312,35)運(yùn)行盤算足本ps:kf.ps1您輸入的是2行*

7、2列表!n=1055卡方=9.9427369384943本例是四格表的卡方統(tǒng)計(jì)題目，說明本研究的pershell步伐同樣實(shí)用于四格表卡方統(tǒng)計(jì)。通過以上二例的盤算，我們以為通過pershell足本盤算行列表卡方是一種便利易用的要領(lǐng)。4討論以下幾點(diǎn)是對本研究中足本步伐的總結(jié)說明：假設(shè)編寫的kf.ps1足本不克不及運(yùn)行，那么應(yīng)該開啟ps1足本的運(yùn)行權(quán)限，在操縱臺輸入set-exeutinpliyunrestrited即可。微軟從寧靜角度思量封閉了足本運(yùn)行權(quán)限。微軟操縱體系要求運(yùn)行ps1足本時要鍵入完備的途徑，在本文中的完備途徑是:kf.ps1，而不克不及只鍵入文件名kf.ps1。不少醫(yī)學(xué)事情者利用exel做卡方盤算，缺點(diǎn)很顯著，在exel表格中對行列數(shù)的限定是預(yù)先設(shè)定的，而本研究的足本步伐沒有行列數(shù)的限定，可以機(jī)動盤算多行多列的表格卡方。本研究的足本步