本文選自美國NCGR占星研究協會NCGR占星調查研究年鑑 – 2018-19年冬季號 第41-48頁 ***特別感謝 Astrid Lian 翻譯好文與會員共享*** 正文 預測結果 原文 正文 地震有季節性嗎？對92年期間全球規模7以上地震分布之分析 摘要 本文假設數據的分類和呈現方式，會對這些數據得出的結論產生重大影響。 序言 突然間，沒有警告，沒有憐憫，只有混亂、破壞和死亡。一場我們稱之為地震的大規模震盪襲來，幸運的人立即身亡，而受困的倖存者在痛苦、黑暗和絕望中度過數日，直到他們因負傷及乾渴慢慢死去。救援人員在瓦礫堆中瘋狂地搜索，希望逐漸消失，空氣瀰漫著死亡的惡臭。 在人類忍受許多所謂的「神的行為」中，地震是目前為止最殘酷的一種。如果上帝的旨意「行在地上，如同行在天上」，我們是否可以察看天空，得知即將發生的事？地震可以預測嗎？地震有季節性嗎？關於美國地質調查局（United States Geological Survey / USGS）的一份聲明指出： 「不存在所謂的『地震天氣』。據統計，在寒冷的天氣、炎熱的天氣、多雨的天氣等，地震的分布大致相同。已知有關主要風暴系統（颱風、颶風等）的極大低壓變化，會引發地殼斷層的滑動（慢地震／slow earthquakes），也可能在引發一些破壞性地震上發揮作用。然而，這些數字很小，沒有統計學上的意義。」（Malcolm Johnston博士提供）[1] 。 而2010年1月不幸的海地地震，正是向我提出調查這一現象的挑戰。多年來，我一直觀察熱帶氣旋和其他惡劣天氣，在我看來地震和颶風常常同時出現在新聞頭條上。從那時的觀察過程中，我相信自己的學習有些收獲，現在將與大家分享其中的一個結果。 研究計劃 本文假設，以曆法月份來組織和記錄科學數據的傳統方式，可能會掩蓋這些數據原本存在的自然模式。 我的做法是，從美國地質調查局的地震項目中下載了1205次規模7.0以上的地震，這些地震發生在1925年至2016年、共橫跨了92年期間；然後，按照三種處理方法來組織這些數據：首先，我以傳統方法、即按曆法月份進行分類後，把數據全加總起來，這是科學機構的通用方法。再來，我根據黃道星座（Signs of the Zodiac）、大約每個月的21日到下個月21日，整理相同的數據。最後，作為對照組，我把黃道拆開來、從每個月的5日到下個月5日左右，將數據進行分類。我接著採用基本的統計程式，用圖形說明所獲得的分布，再透過卡方適合度（Chi-square Goodness of Fit）[2]和二項分布（Binomial Distribution）[3]檢定對其進行分析。變量和檢定的具體細節列於附錄中。 顯著水準(p)（Significance Level） 研究人員選擇顯著水準(p)作為衡量統計檢定成功／失敗（即存在「證據」）的臨界值（Critical Value），是假定「虛無假設」（null hypothesis）為真時獲得樣本結果的可能性；即他／她願意接受正確拒絕「虛無假設」的機率水準，那麼他／她的研究假說可能就是正確的。另外，顯著水準是受眾群體和其他研究人員在審查檢定時考量的因素，可依此決定這項資訊是否可信、有效或有用。 顯著水準以機率(p)表示，為小數形式從0.0（無機會）到1.0（肯定有機會）不等，這是一個事件因正常抽樣變化或僅僅是「隨機」（random chance）發生的可能性。 顯著水準不是固定不變的，視情況而定。最常見的是： (p) = 0.10，即10%或十分之一的機會（中等顯著）(p) = 0.05，即5%或二十分之一的機會（最普遍接受的水準）(p) = 0.01，即1%或百分之一的機會（非常高度的顯著水準） 物理科學量化檢定的顯著水準通常相當嚴格：(p)為0.05或(p) 0.01。社會科學質性檢定（如占星學）的水準比較寬鬆，認為0.10的(p)便具有中等意義[4]。 由於我們的實驗是對質性／分類變量（qualitative/ categorical variables）的次數分布進行比較，我選擇小於或等於(p)=0.10的臨界值，這是一個中等顯著水準，以證明占星學對這份地震活動樣本的影響。 處理方法一：按曆法月份 圖一以圖例形式顯示了1925年至2016年期間所發生的1205次規模7.0以上地震（EQ）之曆法月份分布。卡方適合度檢定得出的機率為(p)=0.30，表示類似的模式可能在100個隨機樣本中出現30次；或者，換句話說，這種模式是隨機的結果為30%，遠遠超過普遍接受的顯著水準。請注意，黑色的分布線大致沿著灰色的預期平均線，只有9月和11月兩處明顯偏離。 表格1是通過二項分布檢定法所得到的（單側；one-tailed）機率[5]總結。如表格所示，9月的地震低於預期，其顯著水準為(p)=0.04，11月的地震高於預期，其顯著水準為(p)=0.01，兩者都認為是高度顯著。 然而，我們必須記住，這種地震分布模式之於整體卡方值為隨機抽樣變化發生的可能性是30%，因此，在相同的1205次地震資料進行第二系列12個二項式檢定的結果，確實增加了12個類別（星座，譯註：應為月份）檢定中，結果至少有一個類別(p)小於0.05的機會。這就是所謂的「邦佛洛尼」（Bonferroni）效應。打個比方，這就像對著一個距離很近的大靶子開了12槍，你射中「靶心」（Bulls@Eyes）的機會就會增加。 以這種傳統方式來看，似乎支持了美國地質調查局（「不存在所謂的地震天氣」）堅持的「虛無假設」立場。 曆法月份分類：左為【圖1】卡方適合度檢定分布；右為【表格1】二項分布檢定法的機率總結 ＊ ＊ ＊ ＊ ＊ ＊ ＊ 處理方法二：按黃道星座 當相同的1,205次地震按黃道星座分類時，出現了一個非常不同且獨特的模式。這種分布顯示了對稱性和內聚性。 卡方檢定（The Chi-square test）返回的機率為(p)=0.10，這是一個中等但仍然顯著的結果，表示類似的模式在100個隨機樣本中只出現10次，或者說，這個結果只有10%隨機發生的可能性，與月份分類明顯不同。還注意到，黑色的分布線在12個星座中只有5個星座沿著灰色的預期平均線，這一圈的其餘部分有5個星座顯而易見地偏離，與月份分布截然不同。 此外，在對12個類別（星座）逐一進行二項分布檢定時，黃道星座明顯偏離預期平均值的結果有五個，而月份分布中只有兩個。表格2的二項式總結包含了五個中等到高度顯著的結果，而考量到這種地震分布模式的整體卡方值(p)0.10較低，這些結果獲得了某種可信度，而月份分布卻沒有。 再用一次剛剛的比喻，這個「黃道靶子」更小，更遠，而我們一樣開12槍卻命中了5次靶心。 我認為，當這個樣本按黃道星座分類時，其結果的可信度確實挑戰了美國地質調查局的立場，即「這些數字很小，沒有統計學上的意義」。 黃道星座分類：左為【圖2】卡方適合度檢定分布；右為【表格2】二項分布檢定法的機率總結 ＊ ＊ ＊ ＊ ＊ ＊ ＊ 處理方法三：分割黃道 最後，我們得到圖三，以圖形展示當我們在對照組抵消黃道的假設影響時會發生什麼。 請記住，我們觀察到月份分類都包含了一個星座的2/3和下一個星座的1/3，可能會使結果偏向於擁有較大比例天數的星座。 如果我們把黃道在月份分類比例大於2/3-1/3的影響削弱，會得到一個非常不同的結果。當透過「分割」一半黃道的影響，創建另一12個類別的對照組，並將地震的分布均勻地分配到這些中和類別裡，卡方檢定返回一個非常高的機率值，(p)=0.56，意味著100個隨機樣本中會產生56個像這樣的模式，或者說，這種分布有56%的可能性是隨機的結果。 最後，由於這種分布為隨機的機率如此之高（0.56），所得到的兩個顯著的二項式結果（表格3）幾乎毫無意義，完全是隨機抽樣變化的結果。還是用那個比喻：一個伸手就摸得到的特大靶子。 分割黃道分類：左為【圖3】卡方適合度檢定分布；右為【表格3】二項分布檢定法的機率總結 ＊ ＊ ＊ ＊ ＊ ＊ ＊ 小結 綜上所述，我們考慮了三種分布，按其顯著水準的順序排列： 黃道星座(p) = 0.10 （有五個顯著的類別）曆法月份(p) = 0.30 （有兩個顯著的類別）分割黃道的檢驗(p) = 0.56 （有兩個顯著的類別） 顯然，黃道星座的分布是最不可能偶然發生的，而且12個類別（星座）中顯示出五個顯著的（單側）二項式機率（即五次靶心），本研究的證據支持這樣的假設：按黃道星座的分類，呈現了這種自然現象最清晰的畫面，我們可以據此進行調查。 然而，我相信這項研究最重要的發現是，實際的地震在分類和重新分類下對每種分布的「影響」，第一是整個黃道的分布，然後是曆法月份（2/3-1/3黃道）的分布，最後是分割成一半的黃道，以便盡可能地減少其影響，同時保持地震資料在其自然時間序列中發生的完整性。 而地震模式的機率從黃道星座劃分的10%低點，跳到無黃道星座之對照組56%的高點（跟擲硬幣的機率差不多），成為黃道星座對地震活動具影響力的證據。 【圖5】地震模式於黃道與無黃道分布之機率；中和黃道後，分布降低至「虛無假設」之水準 ＊ ＊ ＊ ＊ ＊ ＊ ＊ 那麼，這是否意味著我們可以用太陽星座學來預測地震？不，當然不是。但它確實意味著我們在研究因果時，可以知道要去何時何地找關係。例如，我們的樣本在這92年中，速度較快的行星已經繞行了許多圈；冥王星和海王星還沒繞完，而天王星只繞了一圈。冥王星在那些星座花費了最多的時間，也顯示出最大的意義，所以這些慢速者也許應該成為我們研究的基礎。 哦，還有一件事：似乎確實有一個地震的季節。在這項研究中，巨蟹座和處女座的重大地震次數為83次，但在獅子座期間這數字躍升了50%，達123次，所以在獅子座期間遇上重大地震的機率大幅提高。 因此，如果你有業務往來於全球地震活躍的地區，可能你的行程會想跳過雙子座和獅子座期間，改在巨蟹座或處女座的月份出差！這就是我現在的全部內容，謝謝你的關注。 結語：地震預測 我確實相信，地震，或者說，地震發生的時間，可以有一定的預測性，因為過去我在這方面有些成功的經驗。截至本文交稿給NCGR的日期，即2018年6月27日，而自2018年1、2月這上一個交食週期以來，還沒有發生規模7以上的地震。就在地平線上，我看到一些極端的天文事件將在今年夏天發生：一連串與近地點／遠地點週期重合的食相，以及火星，在南赤緯極端地出界，正對齊了月交點，並將與在金牛座和摩羯座的土星／天王星形成半六分相／四分相。因此，根據這些訊息，我要做一個有憑據地推測，即規模7以上的地震，至少有一次，可能還有另一次，將在下列任何日期的前後一或兩天發生： 2018年7月13日 2018年7月19日 2018年7月27日 2018年8月1日 2018年8月10日 2018年8月20日 二項分布計算器[6] 6個(n) 正中1個(K) 的機率是p=0.179，所以，看起來我有大約18%的機率預測成功。時間會證明一切。 ＊ ＊ ＊ ＊ ＊ ＊ ＊ 附錄 為了避免無聊的重複，我只納入變量的具體細節，以及按黃道星座劃分地震分布的統計檢定。 地震（EQ）頻率：M7 + (1925-2016) 92年 x 365.25= 33,603天為母體1205次eq ÷ 33,603天 = 每天(p) 0.03585次M7+EQ M8+ (1901-2016) (用於增補)116年 x 365.25 = 42,369天93年eq ÷ 42,369天 = 每天(p) 0.00219次M8+EQ 黃道變量／分割黃道的檢驗365.25 ÷ 12星座 = 每個星座30.4天92年 x 30.4 = 每個星座2800.25天(p) 0.03585 x 每個星座2800.25天 = 預期100.4天 曆法月份變量／30天的月份92年 x 30 = 每月2760天(p) 0.03585 x 每月2760天 =…
歡迎加入NCGR會員！