=========== skyline-ascending-age-65over-population-prefecture-2010 ========= ====== 都道府県別2010年老年人口比率の ====== 昇順スカイライン図・昇順累積散布図・昇順累積相対散布図【ローレンツ曲線】 =============================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0047.00,xxx // ケース始点,終点番号 ,★第１系列名　分子 老年人口（65歳以上人口） ,ddd // 空白で同一ケース範囲,★第２系列名　分母 人口総数（年齢不詳除く） $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 1358068 5497549 352768 1368197 360498 1325082 520794 2330633 320450 1084144 321722 1165591 504451 2016978 665065 2956404 438196 1989293 470520 1997353 1464860 7167636 1320120 6128826 2642231 12969827 1819503 8996103 621187 2364157 285102 1089110 275337 1160571 200942 798543 211581 858373 569301 2146726 499399 2071947 891807 3743297 1492085 7348784 447103 1842552 288788 1397124 605709 2593965 1962748 8776018 1281486 5556205 333746 1392819 270846 993424 153614 583663 207398 713769 484718 1928064 676660 2828506 404694 1446699 209926 778310 253245 980366 378591 1422761 218148 758486 1123376 5035432 208096 846749 369290 1420134 463266 1806312 316750 1189703 291301 1130743 449692 1699221 240507 1384780 $l // 入力変量のリスト出力コマンド ====================== $$v // 変量分析セクション $a // 変量記号の割り当て d,ddd // d 分母 x,xxx // x 分子 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 $$v // 変量分析セクション $a // 変量記号の割り当て d,ddd // d 分母 x,xxx // x 分子 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 U=(100) // ★比率の表示単位 百分率100 10万人当たり　100000 q=(x/d)*U // 比率q t=t..(x) // 順序数（第1ケース1，第2ケース2，第3ケース3，…） j=r.l(q)blank // ★選択変量（ここではq）での順位（昇順）j blankで 欠落値も最後の順位をつける m=:ci(d) // 個体識別文字列m作成 =pr*(t,d,x,q,j,m) // 数値プリント ............................... e=pmt(d,j) // 順位jで並び替え後の 分母e y=pmt(x,j) // 順位jで並び替え後の 分子y r=pmt(q,j) // 順位jで並び替え後の 比率r n=pmt(m,j) // 順位jで並び替え後の 個体識別文字n n,nam,:ci,n=pmt(m,j) // 並び替え後の個体識別変量nを文字列を示す変量名（先頭「:ci,」）に変更 =pr*(e,y,r,n) // 数値プリント ................................ f=cum(e) // 累積分母f f<i>=e<1>+e<2>+...+e<i-1>+e<i> g=(f-e) // １つ前までの累積分母 g<i>=e<1>+e<2>+...+e<i-1> =f<i>-e<i> z=cum(y) // 累積分子z z<i>=y<1>+y<2>+...+y<i-1>+y<i> v=(z-y) // １つ前までの累積分子 v<i>=y<1>+y<2>+...+y<i-1> =z<i>-y<i> a=@.s(d) // 分母合計a（スカラー） u=@.s(x) // 分子合計u（スカラー） s=(u/a) // 全体比率s（スカラー） .......................... h=(f/a) // 累積相対分母h w=(z/u) // 累積相対分子w =pr*(f,z,h,w,n) // 数値プリント .......................... ジニ係数の計測 k=(h-h1)*(w+w1)*0.5 // 第2〜終点ケースのローレンツ曲線と水平軸の間の各台形面積 i=@."(h)1. // 累積相対分母hの第１ケースの値（スカラー） l=@."(w)1. // 累積相対分子wの第１ケースの値（スカラー） k=@.s(k) // 第2〜終点ケースのローレンツ曲線と水平軸の間の各台形面積の合計（スカラー） k=(0.5-k-i*l*0.5)*2 // 45度線下の三角形面積0.5とローレンツ曲線と水平軸の間の面積の差を２倍 【ジニ係数】 k,nam,ジニ係数 // 変量kに変量名を付与 =pr*(a,u,s,k) // 数値プリント .................... >=(0,s*U) // 全体比率s（スカラー）（表示単位U）の横線 y=0*x+s*U の右辺係数 [0,s*U] の関数「>」 ^=(s,0) // 累積散布図の比率の斜線 y=s*x+0 の右辺係数 [s,0] の関数「^」 :=(1,0) // 累積相対散布図【ローレンツ曲線】の45度の斜線 y=1*x+0 の右辺係数 [1,0] の関数「:」 .................... o=(0*x) // すべてゼロの数値の変量oを作成（図の原点に利用） =============== $$g // グラフセクション $g // スケールの目盛り指示コマンド（標準10ポイント） r,001 // 変量rの目盛りを細かく1ポイントごとに z,001 w,001 $z // ゼロ軸表示　なおゼロ軸表示を抑止するには，この行および次行の先頭に「//」を挿入 rzw // 変量r z wのゼロ軸表示 ------------------------------------------------- $3 // ３次元図　昇順スカイライン図 r,f, ,n,>,* // 縦軸r,横軸f,奥行軸なし,個体識別n,関数>,合成用保存* r,g, ,n,* // 縦軸r,横軸g,奥行軸なし,個体識別n,合成用保存* // 合成 昇順スカイライン図（リンク縦面描画，３次元図圧縮） ........ $3 // ３次元図　昇順累積散布図 v,g, ,n,^,* // 縦軸v,横軸g,奥行軸なし,個体識別n,関数^,合成用保存* z,f, ,n,* // 縦軸z,横軸f,奥行軸なし,個体識別n,合成用保存* // 合成 昇順累積散布図（リンク縦面描画，３次元図圧縮） ........ $3 // ３次元図　昇順累積相対散布図【ローレンツ曲線】（バブルサイズは分子に比例） w,h, ,n=y,:,* // 縦軸w,横軸h,奥行軸なし,個体識別n=バブルy,関数:,合成用保存* o,o, ,n,* // 縦軸o,横軸o,奥行軸なし,個体識別n,合成用保存*【原点】 // 合成（散布点と原点のリンク直線描画，３次元図圧縮を利用） ======================= $$ // 終了セクション