============= ternary2-population-pref-2000-2008-uc ================ =============　三角三色バブルグラフの異時点合成 ============================================== ==================　2000年 =================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド ---------- クロスセクションでは県名や企業名等の文字データを扱うことも多い． ---------- 各文字変量には漢字２文字（英字４文字）のみ入力される． ---------- それを超える文字は無視される． ---------- 文字系列変量名の先頭は「:n1,」「:n2,」…「:n6,」のいずれかを用いる ---------- 識別文字系列変量名の先頭は「:ci,」を用いる 0001.00,0047.00,aa // ケース始点,終点番号 ,★第１系列名 年少人口 ,bb // 空白で同一ケース範囲,★第２系列名 生産年齢人口 ,cc // 空白で同一ケース範囲,★第３系列名 老年人口 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 792352 3832902 1031552 223141 964661 287099 212470 899177 303988 353516 1601826 409156 163095 746252 279764 186182 772100 285590 341038 1353500 431797 458501 2030360 495693 306905 1352311 344506 306895 1346441 367117 1024787 5011202 889243 842534 4235925 837017 1420919 8685878 1910456 1184231 6121470 1169528 365667 1581186 526112 157179 730541 232733 175569 781137 219666 130143 529017 169489 137594 576767 173580 334306 1404575 475127 322769 1401064 383168 568986 2532256 665574 1081280 4914857 1019999 283081 1222594 350959 220072 906629 215552 360531 1810233 459273 1249955 6224186 1315213 830112 3776483 939950 213822 987435 239432 159496 683805 226323 93584 383921 134984 111982 460103 189031 291346 1265122 393658 428035 1916796 531537 213578 974131 339836 117217 525724 180637 148215 659881 214242 219340 953189 320078 111740 509050 191729 742740 3393080 870290 144028 553351 179132 243046 956692 315871 288654 1173790 396020 179439 774403 265901 187431 740401 241754 280717 1101401 403239 264279 861826 182557 ====================== $$v // 変量分析セクション ((B(( // 繰り返しルーティン B の始点 $a // 変量記号の割り当てコマンド X,aa //　年少人口 Z,bb // 生産年齢人口 Y,cc // 老年人口 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換コマンド S=(X+Y+Z) // 人口総数 x=(X/S*100) // 年少人口構成比％ y=(Y/S*100) // 生産年齢人口構成比％ z=(Z/S*100) // 老年人口構成比％ p=:ci(x) // ケース識別文字系列作成 =pr*(X,Y,Z,S,x,y,z,p) // 数値プリント -------------------- $r // 回帰コマンド ,run,y=(x,z) // 被説明変数y，説明変数x,zによる重回帰の計測 ,run,Y=(X,Z) // 被説明変数Y，説明変数X,Zによる重回帰の計測 $t // 変数変換コマンド Y=(y) // データのy変量の三角グラフ平面への縦軸変換 X=(2*x+y)/1.7320508 // データのx変量の三角グラフ平面への横軸変換 ............. ))B)) // 繰り返しルーティン B の終点 $t // 変数変換コマンド ?x=(x) // 異時点グラフ作成用に古い時点の変量の記号を変更 ?y=(y) ?z=(z) ?Y=(Y) ?X=(X) ?S=(S) ================================================ ====================== 2008年 ================ $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0047.00,aa // ケース始点,終点番号 ,★第１系列名 年少人口 ,bb // 空白で同一ケース範囲,★第２系列名 生産年齢人口 ,cc // 空白で同一ケース範囲,★第３系列名 老年人口 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 675000 3555000 1305000 179000 873000 340000 174000 823000 355000 317000 1520000 504000 128000 665000 315000 152000 720000 317000 288000 1268000 496000 404000 1928000 632000 276000 1311000 424000 279000 1281000 452000 975000 4777000 1361000 820000 4069000 1233000 1517000 8723000 2599000 1202000 6000000 1715000 308000 1473000 611000 143000 681000 278000 161000 739000 267000 116000 499000 197000 120000 545000 206000 297000 1319000 554000 294000 1325000 481000 516000 2424000 860000 1088000 4896000 1419000 258000 1183000 434000 212000 914000 276000 347000 1694000 588000 1213000 5725000 1868000 778000 3576000 1233000 188000 898000 318000 132000 615000 264000 78000 366000 152000 93000 424000 207000 270000 1206000 473000 394000 1815000 660000 186000 883000 394000 100000 487000 207000 137000 617000 249000 188000 886000 370000 95000 464000 215000 702000 3269000 1084000 124000 528000 205000 195000 882000 363000 254000 1110000 457000 160000 730000 311000 157000 693000 286000 239000 1031000 447000 246000 893000 237000 ====================== $$v // 変量分析セクション ((B)) // 繰り返しルーティン B の実行 ================================================ $t // 変数変換コマンド f=(-1,-1,+100) // 関数f y= -x -z +100 （つまり x+y+z = 100） ............. i=(100,0,0,50,50,0) // 三角形の頂点と中点の座標 j=(0,100,0,0,50,50) k=(0,0,100,50,0,50) Q=:ci(i)******** // 三角形の頂点と中点の３次元図印字用の文字系列Q ............. @=(0*x) // 原点の変量（ケースの数はデータ分） i=(@,i) // 原点の変量と三角形の頂点と中点を連結した変量 j=(@,j) k=(@,k) Q=(p,Q) // データの散布点印字変量pと頂点と中点の印字変量Qの連結 Q,nam,:ci,Q=(p,Q) // 印字変量Qが文字系列であることを示す変量名に変更 .... =pr*(i,j,k,Q) // 数値プリントしてチェックするには先頭....を取る ------------------------　三角グラフ平面用に変換 U=(1.732,0) // 関数U Y=1.732X+0 三角形の左辺 V=(-1.732,200) // 関数V Y=-1.732X+200 三角形の右辺 W=(-1.732,0) // 関数W W=-1.732X+0 原点を通り三角形右辺に平行(負値の拡張三角グラフ用) ............. 三角形の頂点と中点の座標 J=(j) // 三角形の頂点と中点の三角グラフ平面への縦軸変換 I=(2*i+j)/1.7320508 // 三角形の頂点と中点の三角グラフ平面への横軸変換 --------------------- a=(0,0,0,25,50) // 小さい三角形の頂点と中点のの座標 b=(0,50,25,0,0) c=(100,50,75,75,50) .... // @は @=(0*x) として定義済みで，原点の変量（ケースの数はデータ分） a=(@,a) // 原点の変量と小さい三角形の頂点と中点を連結した変量 b=(@,b) c=(@,c) ............. B=(b) // 小さい三角形の頂点と中点の三角グラフ平面への縦軸変換 A=(2*a+b)/1.7320508 // 小さい三角形の頂点と中点の三角グラフ平面への横軸変換 v=(-1.732,100) // 関数v Y=-1.732X+ (50*2) 小さい三角形の右辺 ========================= グラフセクション $$g ----------------- 目盛 $g X,001 // X変量の目盛 1間隔（標準は10間隔） Y,001 // Y変量の目盛 1間隔（標準は10間隔） y,001 // y変量の目盛 1間隔（標準は10間隔） ----------------- ３次元図 $3 // 三角グラフ立体 j,i,k,Q,* // 縦軸j，横軸i，奥行軸k，散布点印字Q，合成用保存* y,x,z,p=S,f,* // 縦軸y，横軸x，奥行軸z，印字p=バブル変量S，関数f,合成用保存* ?y,?x,?z,p=?S,f,* // 縦軸?y横軸?x奥行軸?z,印字p=バブル?S,関数f合成用保存* // 合成 $3 // 　三角グラフ平面 J,I, ,Q,* // 縦軸J，横軸I，奥行軸なし，印字Q，合成用保存* Y,X, ,p=S,U,V,W,* // 縦軸Y横軸X奥行軸なし印字p=バブル変量S,関数U,V,W合成保存* ?Y,?X, ,p=?S,U,V,W,* // 縦軸?Y横軸?X奥行軸なし印字p=バブルS,関数U,V,W合成保存* // 合成 ----------------- 小さいグラフ $3 // 　小さい三角グラフ立体 b,a,c,Q,* // 縦軸b，横軸a，奥行軸c，散布点印字Q，合成用保存* y,x,z,p=S,f,* // 縦軸y，横軸x，奥行軸z，印字p=バブル変量S，関数f,合成用保存* ?y,?x,?z,p=?S,f,* // 縦軸?y横軸?x奥行軸?z,印字p=バブル?S,関数f,合成用保存* // 合成 $3 // 　小さい三角グラフ平面 B,A, ,Q,* // 縦軸B，横軸A，奥行軸なし，印字Q，合成用保存* Y,X, ,p=S,U,v,W,* // 縦軸Y横軸X奥行軸なし印字p=バブル変量S,関数U,v,W合成保存* ?Y,?X, ,p=?S,U,v,W,* // 縦軸?Y横軸?X奥行軸なし印字p=バブルS,関数U,v,W合成保存* // 合成 ========================= 終了セクション $$