時空掃描統(tǒng)計量-空間統(tǒng)計

1、時空掃描統(tǒng)計量-空間統(tǒng)計1.1時空掃描統(tǒng)計量空間自相關分析具有一定的局限性，其沒有納入時間因素故而不能判斷地理空間上聚集區(qū)域隨時間變化的趨勢。而時空掃描很好的解決了這一難題，不僅能夠揭示地理空間上聚集區(qū)域隨時間變化的規(guī)律，同時還能得出聚集區(qū)域的相對危險度，更加精確的對空間聚集區(qū)域位置進行了定位。圖錯誤！文檔中沒有指定樣式的文字。-1時空掃描原理示意圖時空聚集性分析同時考慮了時間和空間兩個因素，主要采用移動窗口法（movingwindows）,如圖錯誤！文檔中沒有指定樣式的文字。-1所示，在地理空間上創(chuàng)建掃描窗口（圓柱體），圓柱體的底面對應研究的地理區(qū)域，圓柱體的高對應掃描時間問隔，圓柱體的半徑

2、對應掃描的風險人口。掃描窗口的大小和位置在地理空間上不斷的變化，其過程為：隨機選取地理空間上任一地理位置作為掃描窗口的底面中心，相應的地理區(qū)域（掃描窗口的底面積）不斷變化，對應時問問隔（掃描窗口的高）也不斷變化，直到達到預先規(guī)定的上限為止。然后對每一個掃描窗口，根據(jù)實際屬性值（如發(fā)病數(shù)）和人口數(shù)可計算出預期屬性值（發(fā)病數(shù)），然后利用掃描窗口內(nèi)外的實際值和預期值構(gòu)造對數(shù)似然比（loglikelihoodratio,LLR）檢驗統(tǒng)計量，最終選取LLR值最大的窗口為高聚集窗口。針對每個掃描窗口，在泊松分布假設前提下，有:Uz=亮（式錯誤！文檔中沒有指定樣式的文字。-1）其中Uz表示隨機假設下時空窗口

3、Z中的理論事件數(shù)；mz表示時空窗口Z中的總案例數(shù)；mG表示研究區(qū)域內(nèi)的總案例數(shù)；nG表示研究時空范圍內(nèi)的總事件數(shù)。接下來，根據(jù)窗口內(nèi)外的理論與實際事件數(shù)構(gòu)造對數(shù)似然比（loglikelihoodratioLLR）統(tǒng)計量來描述窗口內(nèi)事件異常聚集的程度，表達為：、nZ廣nG_nZnzInG-nZ._LZLUZJ,UGUZjz.、LLR=上=工_（式錯誤！文檔中沒有指定樣式的文字。-2）L。n）、UgJ其中nZ表示時空窗口Z中的是事件數(shù)，nG是研究區(qū)域總事件數(shù)，Uz是時空窗口Z中的理論事件數(shù)，UG是研究區(qū)域總的理論事件數(shù)，Ug=ZUz。LLR是單調(diào)遞減的，P值是根據(jù)蒙特卡羅假設檢驗得出。當LLR取最

4、大值且P<0.05時，可以認為該區(qū)域存在聚集特征?！景咐?-41使用廣東某疾病數(shù)據(jù)（GUangdong_case_popU.xls包含病例和人口數(shù)據(jù)）和地理坐標數(shù)據(jù)（GUangdong_geo.xlS數(shù)據(jù)，以廣東省區(qū)縣地理位置為基本單元、4。風險人口為半徑、月為時間問隔，運用SaTsca詆件對廣東省某疾病進行回顧性時空掃描統(tǒng)計，計算RR值和時空聚類的LLR值。部分數(shù)據(jù)分別展示如表錯誤！文檔中沒有指定樣式的文字。-1表錯誤！文檔中沒有指定樣式的文字。-2所示：表錯誤！文檔中沒有指定樣式的文字。-1廣東某疾病數(shù)據(jù)（部分數(shù)據(jù)）qUxiancasepopUtime4401036793026Jan

5、-094401031793026Feb-0944010336793026Mar-09440103208793026Apr-09440103180793026May-09440103222793026Jun-09440103142793026Jul-0944010362793026Aug-0944010385793026Sep-0944010395793026Oct-09表錯誤！文檔中沒有指定樣式的文字。-2廣東地理坐標數(shù)據(jù)（部分數(shù)據(jù)）quxianlonlat440781112.7122.03440785112.2822.24440981110.9722.03440983111.1222.434

6、41223112.4323.67441224112.1723.97441225111.7123.56441226111.9823.27441402116.1224.28441421116.1924.3600注意：數(shù)據(jù)集中白行需要按地點-時間的順序規(guī)范排列。如果某地某個時間段沒有發(fā)生事件，則不能略去相應的行，而要將該行對應的事件數(shù)記為步驟1:打開文件打開SaTScan,單擊【File->NewSessiori,彈出對話框，如圖錯誤!文檔中沒有指定樣式的文字。-2所示，圖中顯示有三個界面，分別為：Input、Analysis、Output。圖錯誤！文檔中沒有指定樣式的文字。-2時空掃描界面步

7、驟2:Input界面在Input界面中，有多個文件選項，分別為:CaseFile、PopulationFile、CoordinateFile。(1)CaseFile,操作步驟如圖錯誤！文檔中沒有指定樣式的文字。-3(ae)所示：國ImportFileWizardCaseFile:匚.ImA.ESqati<51Pmm心&1tsi.Sangdengue叁占e自中，lTheexpectedformaofthecaefile,usingtiiePoissonprobabibtymodel1阮鞍力二dumberpfCos«><Daie/Timi>C0variaseJ>,.<CoriaieN>Cl&arImjortIftheselectetifileisnotSaTScanformattedfwhilfspacedelimitedorfieldsarenotin±eexpectedai