=========== distance-east-japan-earthquake-human-damage ================== ==== 東日本大震災の主要被災市町村の ==== 人的被害浸水域人口比と【震源からの距離】【老年人口比】との関連グラフ ============================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0043.00,aaa // ケース始点,終点番号 ,★第１系列名　緯度 ,bbb // 空白で同一ケース範囲,★第２系列名　経度 ,eee // 空白で同一ケース範囲,★第３系列名　老年人口比2005 ,ccc // 空白で同一ケース範囲,★第４系列名　浸水域人口 ,ddd // 空白で同一ケース範囲,★第５系列名　人的被害 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 40.5122 141.4883 19.58 5229 2 40.6833 141.3688 18.13 1924 2 40.5992 141.3978 18.50 3820 0 40.4525 141.6211 18.00 1189 0 39.6414 141.9572 27.25 18378 901 39.0823 141.7083 27.04 19073 463 40.1903 141.7753 23.73 7171 4 39.0150 141.6295 30.47 16640 2203 39.2759 141.8856 31.20 13164 1355 39.3599 141.9066 28.49 11915 1711 39.4675 141.9489 28.06 11418 938 39.8431 141.7968 34.28 1137 7 39.9303 141.8889 30.04 1582 36 40.0053 141.8861 27.87 1115 1 40.1103 141.8179 26.78 3177 38 40.4086 141.7187 26.66 2733 0 38.2681 140.8697 15.78 29962 860 38.4342 141.3028 24.17 112276 5734 38.3145 141.0222 23.49 18718 22 38.9083 141.5700 26.83 40331 1533 38.1717 140.8919 17.30 12155 1052 38.2940 141.0042 15.31 17144 188 38.1045 140.8699 17.54 8051 186 38.4264 141.2104 20.47 34014 1769 38.0378 140.8527 20.57 14080 273 37.9625 140.8778 27.76 8990 749 38.3814 141.0692 26.68 4053 4 38.3047 141.0594 18.26 9149 76 38.3303 140.9756 13.81 542 4 38.4455 141.4444 30.00 8048 1156 38.6778 141.4464 27.61 14389 1173 37.0506 140.8878 22.42 32520 383 37.7967 140.9197 23.87 10436 469 37.6422 140.9571 24.49 13377 938 37.2145 140.9947 22.48 1385 3 37.2825 140.9935 24.13 1746 13 37.3456 141.0086 19.84 1401 11 37.4044 140.9835 20.31 1127 44 37.4492 141.0125 24.46 1278 35 37.4947 141.0008 24.07 3356 186 37.8763 140.9195 25.49 4666 115 36.8020 140.7511 23.04 7212 6 35.7203 140.6467 22.15 8303 15 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 e,eee // 老年人口比2005 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,e,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) // 数値プリント ---------------- $r // 回帰分析 F,@"*,r=(t,e) // 人的被害浸水域人口比を【震源からの距離】【老年人口比】で説明する回帰係数F ,run,d=(c,t,e) // 人的被害を【浸水域人口】【震源からの距離】【老年人口比】で説明する回帰結果 =============== $$g // グラフセクション $c // 散布図　 d,t,*,P // 縦軸d,横軸t,回帰線描画*,個体識別P 人的被害 <-（距離） r,t,*,P // 縦軸r,横軸t,回帰線描画*,個体識別P 人的被害浸水域人口比 <-（距離） r,e,*,P // 縦軸r,横軸e,回帰線描画*,個体識別P 人的被害浸水域人口比 <-（老年人口比） $3 // ３次元図　縦軸:人的被害浸水域人口比　横軸:震源からの距離，奥行軸:老年人口比 r,t,e,P=d,F // 縦軸r,横軸t,奥行軸e,個体識別P=バブル変量d,関数F ======================= $$ // 終了セクション