地震前兆の目撃報告が集まる掲示板

誠さん

おまたせしました。

ご紹介頂いたサイトを眺めてみましたが、数値だけが並んでいるのではなく、文章として表現されているため、拾い上げるのが大変そうですね。

(1) 間隔はまちまちになってしまうと思いますが、まずは拾えるだけ
「日時」と「溶岩底の高さ」
を拾って、片っ端からExcelに入れて見やすいリストの形に整理しておきます。

(2) 勤務表みたいにカレンダーを「１日１行」で縦一列に並べた表を用意して、そこに(1)のリストの中から目星をつけた「前兆現象」、つまり「溶岩底の高さが△△m以上上昇（または下降）」の発生日に印がわりの「1」を入れてゆきます。

(3) 同じ表に「対応地震」と自分で適当に決めた「震源地域」と「規模（マグニチュードでも震度でも）」の範囲におさまる地震の発生日に印がわりの「1」を入れてゆきます。

ここまでで、貼付けた表の左上の部分が出来上がります。ここが出発点です。

(4) 目星をつけた「前兆現象」と「対応地震」との関係が以下のいずれだったか、を日ごとに拾ってゆきます。
A=前兆現象の発生後、X日以内に対応地震が発生した《的中》
B=前兆現象の発生後、X日以内に対応地震が発生しなかった《ハズレ》
C=前兆現象の発生がなかったが、X日以内に対応地震が発生した《検出モレ》
D=前兆現象の発生がなく、X日以内に対応地震も発生しなかった《安全確認成功》
この作業に手と頭で拾うと疲れたり間違えますし、後で出てくる(7)の繰り返し作業が大変ですから、Excelに自動的に拾ってもらいましょう。
例えば、この表で5月1日の「前兆現象」がC5セル、「対応地震」がD5セルとすると、

5月1日の「１日以内」のABCD判定の計算式はこうなります。（当日は前後関係が不明のため除くこと、としました）
=IF(C5*D6>=1,"A",IF(C5=1,(IF(D6=0,"B")),IF(C5=0,(IF(D6>=1,"C","D")))))

「２日以内」ならこう
=IF(C5*SUM(D6:D7)>=1,"A",IF(C5=1,(IF(SUM(D6:D7)=0,"B")),IF(C5=0,(IF(SUM(D6:D7)>=1,"C","D")))))

「３日以内」ならこうですね。
=IF(C5*SUM(D6:D8)>=1,"A",IF(C5=1,(IF(SUM(D6:D8)=0,"B")),IF(C5=0,(IF(SUM(D6:D8)>=1,"C","D")))))

(5) Excelのcountif関数を使って、ABCDの各出現回数を求めます。

(6) A〜Dの回数を計算サイト
http:/
/www.kisnet.or.jp/nappa/software/star/freq/2x2.htm
に入力してpの値を求めます。
0.10<p：傾向なし
0.05<p<0.10：有為傾向あり
p<0.05：95％の確率で有為
となります。

(7) pが0.1以上の場合は有為傾向なしということですので、条件を変えてみます。つまり(2)と(3)の条件を様々に振ってみるとABCDの値が変わってきますので、その中でpが最も小さな値となるパターンを試行錯誤的に探してみる、ということになります。

私が過去、これを試してみた例を>>2457と>>12889に紹介していますので、そちらもご覧になってみてください。
試行錯誤をする際、まずは上記計算サイトの４つのハコにいろんな数字を試しに入れてみて、それぞれがどんな大小関係になるとpの値が小さくなるか、を理解しておくと、(2)や(3)の作業でムダが省けて作業効率が上がると思います。

> 溶岩底の反応とＤｓｔ値が同じようになっているようです。

については、「対応地震」の部分を「Dst値の△△以上の増加（または減少）」の発生日に置き換えて同じように拾ってみると、溶岩底の高さとDst値との関係の有無を判定できます。