=========== distance-hyogo-kobe-population-from-osaka ================== ==== 兵庫県内市区の人口変化2005年比と ==== 【中心地の大阪駅からの距離】【老年人口比】との関連グラフ ============================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0037.00,bbb // ケース始点,終点番号 ,★第１系列名　経度 ,aaa // 空白で同一ケース範囲,★第２系列名　緯度 ,yyy // 空白で同一ケース範囲,★第３系列名　2010年人口 ,xxx // 空白で同一ケース範囲,★第３系列名　2005年人口 ,ccc // 空白で同一ケース範囲,★第５系列名　2010年年少人口 ,ddd // 空白で同一ケース範囲,★第６系列名　2010年生産年齢人口 ,eee // 空白で同一ケース範囲,★第７系列名　2010年老年人口 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 135.1955 34.6900 210408 206037 29617 138228 41526 135.2394 34.7122 133451 128050 16396 85667 29860 135.1653 34.6805 108304 106985 10535 65135 30002 135.1509 34.6657 101624 103791 10839 60393 29901 135.1303 34.6504 167475 171628 20147 104546 42060 135.0569 34.6306 220411 222729 28711 136506 54608 135.1444 34.7273 226836 225945 31772 142101 52492 135.1978 34.6950 126393 116591 10829 82184 28696 134.9816 34.6832 249298 243637 36117 166199 45073 134.6855 34.8153 536270 536232 80093 338884 115703 135.4064 34.7333 453748 462647 53922 289125 106070 134.9975 34.6431 290959 291027 40266 184936 61866 135.3419 34.7378 482640 465337 71847 313110 92399 134.8956 34.3425 47254 50030 6109 27608 13484 135.3042 34.7269 93238 90590 12635 58820 21570 135.4007 34.7842 196127 192250 29644 127087 39275 134.4681 34.8036 31158 32475 3617 18337 9149 134.8200 35.5445 85592 89208 11893 49523 24144 134.8408 34.7567 266937 267100 39271 172184 55078 134.3902 34.7550 50523 51794 6970 30603 12888 134.9694 34.9935 42802 43953 5977 25061 11734 135.3602 34.7999 225700 219862 32168 142599 50453 134.9901 34.7968 81009 84361 10197 49512 21262 134.7906 34.7660 93901 94813 13450 60280 20059 135.4171 34.8300 156423 157668 21727 94176 40398 134.9316 34.8532 49680 49761 7638 30893 11125 135.2254 34.8895 114216 113572 16051 78588 19286 134.8419 34.9279 47993 49396 6162 29457 12364 135.2193 35.0758 43263 45245 5398 25493 12346 134.7675 35.4047 26501 28306 3316 14419 8759 135.0360 35.1772 67757 70810 9427 38768 19512 134.7792 34.2957 49834 52283 6387 28791 14616 134.8531 35.3397 32814 34791 4320 18524 9948 134.9147 34.4399 46459 49078 5377 26020 15062 134.5495 35.0045 40938 43302 5726 23842 11369 134.9735 34.9174 40181 39970 5805 25515 8861 134.5455 34.8579 80518 81561 11253 49966 19223 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て b,bbb // 経度 a,aaa // 緯度 y,yyy // 2010年人口 x,xxx // 2005年人口 c,ccc // 2010年年少人口 d,ddd // 2010年生産年齢人口 e,eee // 2010年老年人口 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 r=(y/x)*100 // 人口変化2005年比 g=(c+d+e) // 2010年人口（年齢不詳をの除く） s=(e/g)*100 // 老年人口比2005 P=:ci(a) // 個体識別文字列P作成 U=(34.7023) // ★中心地の緯度のスカラー　ここでは大阪駅の緯度 V=(135.4951) // ★中心地の経度のスカラー 大阪駅の経度 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 // 散布図　 y,t,*,P // 縦軸y,横軸t,回帰線描画*,個体識別P 人口2010 <-（距離） r,t,*,P // 縦軸r,横軸t,回帰線描画*,個体識別P 人口変化2005年比 <-（距離） r,s,*,P // 縦軸r,横軸s,回帰線描画*,個体識別P 人口変化2005年比 <-（老年人口比） $3 // ３次元図　縦軸:人口変化2005年比　横軸:中心地からの距離，奥行軸:老年人口比 r,t,s,P=y,F,* // 縦軸r,横軸t,奥行軸s,個体識別P=バブル変量y,関数F,合成用保存* 【観測値】 R,t,s,P,* // 縦軸R,横軸t,奥行軸s,個体識別P,合成用保存* 【推定値】 // 合成 ======================= $$ // 終了セクション