=========== distance-east-japan-earthquake-human-damage-2012-01 ============ ==== 東日本大震災の主要被災市町村の ==== 人的被害と【浸水域人口】【震源からの距離】【老年人口比】との関連グラフ ================================================================================ $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0043.00,bbb // ケース始点,終点番号 ,★第１系列名　経度 ,aaa // 空白で同一ケース範囲,★第２系列名　緯度 ,eee // 空白で同一ケース範囲,★第３系列名　老年人口比2010 ,ccc // 空白で同一ケース範囲,★第４系列名　浸水域人口 ,ddd // 空白で同一ケース範囲,★第５系列名　人的被害 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 141.4883 40.5122 23.24 5229 2 141.3688 40.6833 20.42 1924 2 141.3978 40.5992 20.91 3820 0 141.6211 40.4525 21.79 1189 0 141.9572 39.6414 30.92 18378 534 141.7083 39.0823 30.88 19073 427 141.7753 40.1903 26.36 7171 4 141.6295 39.0150 34.87 16640 1852 141.8856 39.2759 34.81 13164 1056 141.9066 39.3599 32.39 11915 1307 141.9489 39.4675 31.84 11418 769 141.7968 39.8431 37.80 1137 7 141.8889 39.9303 33.85 1582 29 141.8861 40.0053 31.51 1115 1 141.8179 40.1103 30.07 3177 38 141.7187 40.4086 30.52 2733 0 140.8697 38.2681 18.58 29962 730 141.3028 38.4342 27.30 112276 3777 141.0222 38.3145 27.47 18718 32 141.5700 38.9083 30.81 40331 1373 140.8919 38.1717 19.13 12155 967 141.0042 38.2940 18.38 17144 189 140.8699 38.1045 19.76 8051 183 141.2104 38.4264 23.17 34014 1113 140.8527 38.0378 23.41 14080 269 140.8778 37.9625 31.63 8990 690 141.0692 38.3814 30.85 4053 2 141.0594 38.3047 21.56 9149 75 140.9756 38.3303 16.28 542 46 141.4444 38.4455 33.50 8048 936 141.4464 38.6778 30.06 14389 876 140.8878 37.0506 25.08 32520 348 140.9197 37.7967 25.47 10436 458 140.9571 37.6422 26.62 13377 646 140.9947 37.2145 23.77 1385 3 140.9935 37.2825 25.91 1746 25 141.0086 37.3456 21.11 1401 73 140.9835 37.4044 20.96 1127 41 141.0125 37.4492 27.13 1278 54 141.0008 37.4947 26.65 3356 184 140.9195 37.8763 26.94 4666 114 140.7511 36.8020 25.72 7212 6 140.6467 35.7203 24.09 8303 15 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 e,eee // 老年人口比2010 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 M,@"*,d=(c,t,e) // 人的被害を【浸水域人口】【震源からの距離】【老年人口比】で説明する回帰係数M ---------------- $t // 変数変換 i=@."(M)1. // 回帰係数ベクトルMの第１説明変数の係数のみのスカラー　定数指示1.の.は必須 j=@."(M)4. // 回帰係数ベクトルMの定数項のみのスカラー　定数指示4.の.は必須 m=@."(M)2.,3. // 回帰係数ベクトルMの第２・第３説明変数の係数ベクトル抽出　定数数字後の.は必須 =pr*(M,i,j,m) // 数値プリント Z=(d-i*c-j) // 【浸水域人口】と定数による説明部分からの人的被害データの乖離 =pr*(d,i*c+j,Z,P) // 数値プリント =============== $$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 ..........　縦軸:人的被害データの【浸水域人口】による説明部分からの乖離 .......... 　横軸:【震源からの距離】，奥行軸:【老年人口比】 Z,t,e,P=c,m // 縦軸Z,横軸t,奥行軸e,個体識別P=バブル変量c,関数m ======================= $$ // 終了セクション