=========== distance-east-japan-earthquake-population-sendai ================== ==== 東日本大震災の主要被災市町村の ==== 人口変化2005年比と【中心地の仙台からの距離】【老年人口比】との関連グラフ ============================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0043.00,aaa // ケース始点,終点番号 ,★第１系列名　緯度 ,bbb // 空白で同一ケース範囲,★第２系列名　経度 ,ccc // 空白で同一ケース範囲,★第３系列名　人口2010 ,ddd // 空白で同一ケース範囲,★第４系列名　人口2005 ,eee // 空白で同一ケース範囲,★第５系列名　老年人口2005 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 40.5122 141.4883 237473 244700 47912 40.6833 141.3688 41260 42425 7692 40.5992 141.3978 24188 24172 4472 40.4525 141.6211 14702 15356 2764 39.6414 141.9572 59442 63588 17325 39.0823 141.7083 40738 43331 11716 40.1903 141.7753 36875 39141 9290 39.0150 141.6295 23302 24709 7528 39.2759 141.8856 39578 42987 13411 39.3599 141.9066 15277 16516 4705 39.4675 141.9489 18625 20142 5651 39.8431 141.7968 10804 11914 4084 39.9303 141.8889 3843 4241 1274 40.0053 141.8861 3088 3358 936 40.1103 141.8179 4632 5019 1344 40.4086 141.7187 17910 19524 5205 38.2681 140.8697 1045903 1025126 161795 38.4342 141.3028 160704 167324 40435 38.3145 141.0222 56490 59357 13943 38.9083 141.5700 73494 78011 20931 38.1717 140.8919 73140 68662 11876 38.2940 141.0042 62979 62717 9600 38.1045 140.8699 44198 43921 7703 38.4264 141.2104 42908 43235 8849 38.0378 140.8527 34846 35132 7226 37.9625 140.8778 16711 17713 4918 38.3814 141.0692 15089 16193 4320 38.3047 141.0594 20419 21068 3847 38.3303 140.9756 34000 32257 4456 38.4455 141.4444 10051 10723 3217 38.6778 141.4464 17431 18645 5147 37.0506 140.8878 342198 354492 79472 37.7967 140.9197 37796 38630 9221 37.6422 140.9571 70895 72837 17836 37.2145 140.9947 5418 5533 1244 37.2825 140.9935 7701 8188 1976 37.3456 141.0086 15996 15910 3156 37.4044 140.9835 11511 10992 2232 37.4492 141.0125 6932 7170 1754 37.4947 141.0008 20908 21615 5203 37.8763 140.9195 8218 8584 2188 36.8020 140.7511 47026 49645 11439 35.7203 140.6467 69074 70643 15645 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 c,ccc // 人口2010 d,ddd // 人口2005 e,eee // 老年人口2005 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 r=(c/d)*100 // 人口変化2005年比 s=(e/d)*100 // 老年人口比2005 P=:ci(a) // 個体識別文字列P作成 U=(38.2681) // ★中心地の緯度のスカラー　ここでは仙台の緯度 V=(140.8697) // ★中心地の経度のスカラー 仙台の経度 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,s) // 人口変化2005年比を【中心地からの距離】【老年人口比】で説明する回帰係数F R,est,r=(t,s) // 人口変化2005年比を【中心地からの距離】【老年人口比】で説明する推定値 ,pr*,r,R,t,s,c,P // 数値プリント =============== $$g // グラフセクション $c // 散布図　 c,t,*,P // 縦軸d,横軸t,回帰線描画*,個体識別P 人口2010 <-（距離） r,t,*,P // 縦軸r,横軸t,回帰線描画*,個体識別P 人口変化2005年比 <-（距離） r,s,*,P // 縦軸r,横軸e,回帰線描画*,個体識別P 人口変化2005年比 <-（老年人口比） $3 // ３次元図　縦軸:人口変化2005年比　横軸:中心地からの距離，奥行軸:老年人口比 r,t,s,P=c,F,* // 縦軸r,横軸t,奥行軸e,個体識別P=バブル変量c,関数F,合成用保存* 【観測値】 R,t,s,P,* // 縦軸R,横軸t,奥行軸e,個体識別P,合成用保存* 【推定値】 // 合成 ======================= $$ // 終了セクション