=========== distance-east-japan-earthquake-cumulative-human-damage ================== ==== 東日本大震災の主要被災市町村の ==== 【震源からの距離】と【人的被害の累和】との関連グラフ ============================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0043.00,aaa // ケース始点,終点番号 ,★第１系列名　緯度 ,bbb // 空白で同一ケース範囲,★第２系列名　経度 ,ccc // 空白で同一ケース範囲,★第３系列名　浸水域人口 ,ddd // 空白で同一ケース範囲,★第４系列名　人的被害 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 40.5122 141.4883 5229 2 40.6833 141.3688 1924 2 40.5992 141.3978 3820 0 40.4525 141.6211 1189 0 39.6414 141.9572 18378 901 39.0823 141.7083 19073 463 40.1903 141.7753 7171 4 39.0150 141.6295 16640 2203 39.2759 141.8856 13164 1355 39.3599 141.9066 11915 1711 39.4675 141.9489 11418 938 39.8431 141.7968 1137 7 39.9303 141.8889 1582 36 40.0053 141.8861 1115 1 40.1103 141.8179 3177 38 40.4086 141.7187 2733 0 38.2681 140.8697 29962 860 38.4342 141.3028 112276 5734 38.3145 141.0222 18718 22 38.9083 141.5700 40331 1533 38.1717 140.8919 12155 1052 38.2940 141.0042 17144 188 38.1045 140.8699 8051 186 38.4264 141.2104 34014 1769 38.0378 140.8527 14080 273 37.9625 140.8778 8990 749 38.3814 141.0692 4053 4 38.3047 141.0594 9149 76 38.3303 140.9756 542 4 38.4455 141.4444 8048 1156 38.6778 141.4464 14389 1173 37.0506 140.8878 32520 383 37.7967 140.9197 10436 469 37.6422 140.9571 13377 938 37.2145 140.9947 1385 3 37.2825 140.9935 1746 13 37.3456 141.0086 1401 11 37.4044 140.9835 1127 44 37.4492 141.0125 1278 35 37.4947 141.0008 3356 186 37.8763 140.9195 4666 115 36.8020 140.7511 7212 6 35.7203 140.6467 8303 15 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 c,ccc // 浸水域人口 d,ddd // 人的被害 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 r=(d/c)*100 // 人的被害浸水域人口比 P=:ci(a) // 個体識別文字列P作成 U=(38.1) // ★震源の緯度のスカラー　 ■震源の替わりに特定市の緯度の指示も可 V=(142.9) // ★震源の経度のスカラー　 ■震源の替わりに特定市の経度の指示も可 A=abs(a-U) // 各市の緯度の震源からの乖離 B=abs(b-V) // 各市の経度の震源からの乖離 =pr*(a,b,c,d,r,A,B,P) // 数値プリント ..................... 起点（震源）から各終点(各市町村庁）までの距離の計測処理 n=cos(U/180*3.14159) // 起点緯度（ラジアン）のcos l=cos(a/180*3.14159) // 終点緯度（ラジアン）のcos m=cos((U+a)/2/180*3.14159) // 起点終点の平均緯度（ラジアン）のcos　 M=(1+((6378.137**2-6356.752**2)/6356.752**2)*(m**2))**0.5 M=(6378.137**2/6356.752)/(M**3) // 起点終点の平均緯度の子午線曲率半径 N=(1+((6378.137**2-6356.752**2)/6356.752**2)*(n**2))**0.5 N=(6378.137**2/6356.752)/N // 起点緯度の卯酉線曲率半径 L=(1+((6378.137**2-6356.752**2)/6356.752**2)*(l**2))**0.5 L=(6378.137**2/6356.752)/L // 終点緯度の卯酉線曲率半径 =pr*(M,N,L) // 数値プリント h=(A/180*3.14159*M) // 南北の緯度差（ラジアン）の子午線弧長 i=(B/180*3.14159*N*n) // 起点緯度の東西の経度差（ラジアン）弧長 j=(B/180*3.14159*L*l) // 終点緯度の東西の経度差（ラジアン）弧長 t=(h**2+i*j)**0.5 // 起点（震源）から各終点までの距離　km T=(t*1000) // 同上 m 単位 =pr*(h,i,j,t,P,T) // 数値プリント ..................... j=r.l(t)blank // tの小さい順(定数項blankで欠測値にも末尾の順位)の順位変量j ....j=r.g(t)blank // 大きい順の場合は先頭の....を取る c=pmt(c,j) // 並び替え（順序数変量jによる） 浸水域人口 d=pmt(d,j) // 並び替え（順序数変量jによる） 人的被害 r=pmt(r,j) // 並び替え（順序数変量jによる） 人的被害浸水域人口比 t=pmt(t,j) // 並び替え（順序数変量jによる） 震源から市町村庁までの距離 P=pmt(P,j) // 並び替え（順序数変量jによる） 個体識別（印字）文字列 P,nam,:ci,P=pmt(P,j) // 印字系列Pが文字系列であることを示す変量名に変更 ..................... C=cum(c) // 浸水域人口の累和C D=cum(d) // 人的被害の累和D R=(D/C*100)) // 累和の人的被害浸水域人口比R =pr*(C,D,R,t,P) // 数値プリント ---------------- $r // 回帰分析 F,@"*,logR=(logt) // 累和の人的被害浸水域人口比Rを【震源からの距離】で説明する対数回帰係数F .... ■【震源からの距離】を【対象市町村の特定市からの距離】 とする場合は前行を F,@"*,R=(t) に変更 =============== $$g // グラフセクション $g // スケールの目盛り指示コマンド（標準10ポイント） C,002 // 変量Cの目盛りを細かく2ポイントごとに D,002 R,002 $z // ゼロ軸表示 ■【震源】を【特定市】 とする場合 先頭に // を挿入してコメントに R // 変量Rのゼロ軸表示 ■【震源】を【特定市】 とする場合 先頭に // を挿入してコメントに $3 // ３次元図　 　横軸に【震源からの距離】t C,t, ,P=c // 縦軸C,横軸t,奥行軸なし,個体識別P=バブル変量c　【浸水域人口の累和】 D,t, ,P=d // 縦軸D,横軸t,奥行軸なし,個体識別P=バブル変量d　【人的被害の累和】 ...... 変量直後の*は対数関数表示を意味する R*t* ,P=d,F // 縦軸R* 横軸t* 奥行軸なし,個体識別P=バブル変量d,関数F【累和の人的被害浸水域人口比】 .... ■【震源からの距離】を【対象市町村の特定市からの距離】 とする場合は前行を R,t, ,P=d,F に変更 ======================= $$ // 終了セクション