第五章 結論與建議
第二節 後續研究建議
國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
48
最後本文應用先前的模擬結果，檢測民國 98 年癌症死亡資料中，核一廠、
核三廠及麥寮六輕周圍是否有較高的癌症死亡率，檢定結果顯示男性在核一廠及 麥寮六輕周圍有較高的癌症死亡率，女性則是在麥寮六輕周圍有較高的癌症死亡 率。模擬過程中發現在研究區域半徑的變動範圍內，不論半徑的大小，核一廠及 麥寮六輕的檢定顯著，說明這兩個位置附近的癌症發生率特別高，需要趕快找出 可能的原因。其中相同的核能設施，但卻有不同的檢定結果，差異的原因可能來 自於檢測位置造成的邊緣效應、鄰近人口的多寡及運轉時間差異等造成的影響。
需要注意的是檢定結果顯著並非說明檢測位置上的設施造成較高的癌症死亡 率，只能說明檢測位置周圍有較高癌症死亡率。
第二節 後續研究建議
本文比較不同焦點檢定方法的優劣，藉由比較結果選擇合適的檢定方法檢測 臺灣癌症死亡資料，在模擬的過程中，發現雖然焦點檢定在檢測已知位置時，具 有較好的檢定力，但由於統計量中是以各區塊之間的距離作為權重的參數，檢測 不同位置具有不同的檢定統計量，如果檢測研究區域的所有位置，可能需要大量 的電腦模擬，除此之外，焦點檢定通常假設群聚可能的型態，如果檢測所有可能 的地點時，預先假設群聚型態是否合適可能是檢定時會產生的問題，故焦點檢定 方法不適合去檢測研究區域中群聚發生的位置。
由以上的模擬結果，本文建議未來在檢測群聚是否發生的問題方面，可以結 合局部檢定方法與焦點檢定方法，舉例來說，通常不知道哪些事件是否有聚集的 狀況，可以先藉由局部檢定方法，由較大的研究區域做一個概觀的分析，從中找 出可能發生群聚的事件及群聚的位置，再分析出影響此事件的因素及影響的範 圍，最後由於焦點檢定具有較好的檢定力，藉由焦點檢定方法來判定群聚是否發 生。如果已知一些可能發生群聚的位置則可以直接使用焦點檢定方法檢定群聚是 否發生，例如高壓電塔、焚化爐、行動電話基地台等等都是被懷疑對周圍環境可
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
49
能產生影響，可以藉由焦點檢定方法來檢定這些已知位置的周圍是否出現群聚。
在焦點檢定方法的研究上，本文建議未來可以嘗試下列幾個方向：
(1) 如何設定合適的曝露函數
在本文的模擬過程中，可以發現檢定群聚是否發生很重要的因素在於風險分 佈型態的假設，模擬結果顯示合適的曝露函數設定會具有相當好的檢定力，在未 來的研究中可以進一步討論如何由資料中估計風險分佈型態，進而設定合適的曝 露函數。
(2) 如何檢測區域範圍的污染或是多個污染源
應用檢測方法時，通常將污染源以一個檢測位置表示，實際上污染源可能無 法以單一位置表示，例如污染源可能是一個區域範圍，這此情況之下該如何設定 檢測位置。檢定方法通常檢測單一的污染源，實際應用時可能發生在一個區域中 有兩個以上同類型的污染源，該如何檢定此情況下是否發生群聚。
(3) 在檢定方法中考量時間的影響
本文在實證分析中使用的是民國 98 年的癌症死亡資料，雖然檢定結果顯示 核一廠及麥竂六輕周圍有較高的癌症死亡率，但是單年資料的檢定結果可能只是 這段時間內的一個特例，如果檢測群聚可以加入時間的影響，用來判斷檢測位置 周圍是否在這一段時期皆有較高的癌症死亡率，檢測結果可能會較為可信。
(4) 考慮區域與區域之間相互影響的情況
在本文的假設中，皆設定每個區塊的事件數相互獨立，實際上對於某些事件 可能有相互影響的情況，例如已經絕跡的天花等傳染性疾病。當區塊之間相互影 響時，本文所介紹的方法將不適用此情形，如何檢測群聚的發生。
(5) 區塊人口分佈的影響
由於臺灣各鄉鎮市人口分佈不均，如果以各鄉鎮市的地理中心做為區塊中心 可能會造成錯誤的檢定結果，未來在資訊充足的情況下，可以考慮各區塊的人口 分佈，並由人口分佈的分析中設定合適的區塊中心。
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
50
參考文獻
Besag, J. and Newell, J. “The detection of clusters in rare diseases”, Journal of Royal Statistical Society, Series A 154:143-155 (1991).
Bithell, J.F. “The choice of test for detecting raised disease risk near a point source”, Statistical in Medicine 14:2309-2322 (1995).
Diggle, P.J. “A point process modelling approach to raised incidence of a rare phenomenon in the vicinity of a prespecified point. Journal of the Royal Statistical Society, Series A 153:349-362 (1990).
Kulldorff, M. and Nagarwalla N. “Spatial disease clusters: detection and inference”, Statistical in Medicine 14:799-810 (1995).
Kulldorff, M. “A spatial scan statistic”, Communications in Statistics–Theory and Methods, 26(6):
1481-96 (1997).
Kulldorff , M. ”Tests for spatial randomness adjusting for an underlying inhomogeneity: a general Framework”, Journal of the American Statistical Association 101:1289-1305 (2006).
Knox, G. “Detection of low intensity epidemicity”, British Journal of Preventive and Social Medicine 17:121-127 (1963).
Lawson AB. “On the analysis of mortality events associated with a prespecified fixed point”, Journal of Royal Statistical Society, Series A 156:363-377 (1993) .
Lawson, A.B. and Waller, L. “A review of point pattern methods for spatial modelling of events around sources of pollution”, Environmetrics 7:471-488 (1996).
Lawson, A.B. and Denison, D. “Spatial Cluster Modelling”, CRC Press, New York (2002).
Lawson, A.B.” Statistical methods in spatial epidemiology”, Wiley, New York (2006).
Lawson, A.B, Williams F.L.R. “Armadale: a case-study in environmental epidemiology”, Journal of Royal Statistical Society, Series A 157:285-298 (1994).
Marshal, R. C.”A review of the statistical analysis of spatial patterns of disease”, Journal of Royal Statistical Society, Series A 154:421-441 (1991).
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
51
Puett, R.C., Lawson, A.B., Clark, A.B., Herbert, J.R. and Kulldorff, M. “Power evaluation of focused cluster tests”, Environmental and Ecological Statistics 7:303-316 (2010).
Stone, R.A. “Investigations of excess environmental risks around putative sources: statistical problems and a proposed test”, Statistical in Medicine 7:649-660(1988).
Sun, Y. “Determining the size of spatial clusters in focused tests: Comparing two methods by means of simulation in a GIS”, Journal of Geographical Systems 4:359-370(2002).
Tango, T. “A class of tests for detecting general and focused clustering of rare diseases”, Statistical in Medicine 14:2323-2334(1995).
Tango, T. “Score tests for detecting excess risks around putative sources”, Statistical in Medicine 21:497-514 (2002).
Tango, T. “Statistical methods for disease clustering”, Springer, New York (2010).
Waller, L.A. “Book Review of Methods for investigating localized clustering of disease”, Statistics in Medicine 17:1531-1535 (1998).
Waller, L.A. and Turnbull, B.W. “The effects of scale on tests for disease clustering”, Statistical in Medicine 12:1869-1884 (1993).
Waller, L.A. “Statistical power and design of focused clustering studies”, Statistical in Medicine 15:765-782 (1996).
Waller, L.A and Lawson, A.B. “The power of focused tests to detect disease clustering”, Statistical in Medicine 14:2291-2308 (1995).
陳楓玲. “A Simulation Study of Detecting Cluster: An Application to Taiwan Cancer Mortality”, 臺 灣智慧科技與應用統計學術研究會 (2003)。
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
52
附錄圖一、圓形群聚下－不同曝露函數設定的 ST 檢定力比較
研究區域設定為 20×20 個單位正方形所組成的方形研究區域。群聚中心置於 座標(10.4999, 10.4998)。事件發生率為萬分之五。各區人口數皆為一萬人。在模 擬過程中，設定熱點群聚模型參數 r 為 2，線性模型參數 r 為 2，指數模型參數 s 為 1。變動模型參數 RR 由 1.1 到 3，間隔 0.2。模擬次數 1000 次。顯著水準 0.05。
1.5 2.0 2.5 3.0
0.00.20.40.60.81.0
Exponential Decline Clustering
RR
Power
true 1/rank 1/(1+d) 1/(1+d)^2
1.5 2.0 2.5 3.0
0.00.20.40.60.81.0
Linear Decline Clustering
RR
Power
true 1/rank 1/(1+d) 1/(1+d)^2
1.5 2.0 2.5 3.0
0.00.20.40.60.81.0
Hot Spot Clustering
RR
Power
true 1/rank 1/(1+d) 1/(1+d)^2
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
53
附錄圖二、民國 98 年臺灣女性及男性人口數
‧
Radius of study region
power
Radius of study region
power
Radius of study region
power
Radius of study region
power
‧
Radius of study region
power
Radius of study region
power 置於座標(10.4999, 10.4998)；事件發生率為萬分之五，各區人口數皆為一萬人。
模擬次數 1000 次。顯著水準 0.05。
ST EST Pmax MLR KN
檢定統計量 1.590045 2.169933 1.822785 6.528192 30.441103
附錄表二、實證模擬的檢定統計量
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
56
檢測位置：麥寮六輕
ST EST Pmax MLR KN 半徑 50 公里 1.612470 1.965101 2.949467 1.415314 1.346399 半徑 60 公里 1.605934 2.076570 3.042521 1.443530 1.342839 半徑 70 公里 1.797012 2.197090 3.657484 1.451191 1.407876 半徑 80 公里 1.685300 2.174427 3.511849 1.442107 1.393328 半徑 90 公里 1.679711 2.158487 3.426398 1.418047 1.371508 半徑 100 公里 1.763065 2.233654 3.520118 1.428906 1.385263
檢測位置：核一廠
ST EST Pmax MLR KN 半徑 50 公里 1.550707 2.174641 3.239298 1.753100 1.333276 半徑 60 公里 1.683931 2.206187 3.374693 1.751556 1.359253 半徑 70 公里 1.522641 2.125355 3.072313 1.766015 1.326297 半徑 80 公里 1.622108 2.264008 3.545846 1.765939 1.369256 半徑 90 公里 1.719430 2.210415 3.304712 1.733786 1.367128 半徑 100 公里 1.689190 2.241932 3.401418 1.764130 1.363203
檢測位置：核三廠
ST EST Pmax MLR KN 半徑 50 公里 1.657118 1.866266 1.188509 1.362767 1.178981 半徑 60 公里 1.658458 1.888742 1.528353 1.415302 1.203870 半徑 70 公里 1.708616 2.041720 2.593681 1.483291 1.349681
‧ 國
立 政 治 大 學
‧
N a tio na
l C h engchi U ni ve rs it y
57
半徑 80 公里 1.856726 2.183806 2.830849 1.514279 1.353792 半徑 90 公里 1.796789 2.006460 2.766682 1.490371 1.313969 半徑 100 公里 1.689720 2.036051 3.242649 1.486590 1.373939