=============== skyline-shortage-of-doctors ============= ==== ｢必要医師数実態調査｣の医師不足度のスカイライン図 ==== 医師不足率の扇形散布図 ========================================================================== ================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0047.00,aaa // ケース始点,終点番号 ,★第１系列名 ,bbb // 空白で同一ケース範囲,★第２系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 7567 1007 1520 484 1600 640 2408 360 1482 302 1513 360 2397 555 3292 492 2836 486 2490 469 6757 705 6812 803 20161 1656 7527 716 2698 591 1736 291 2119 235 1233 224 1047 302 2718 485 2314 559 4149 861 8267 928 1982 400 1892 407 4260 490 13008 1219 7393 986 2115 347 1812 278 1037 199 1133 320 3358 419 3971 607 2132 307 1268 280 1637 313 2128 370 1501 361 7976 843 1378 167 1944 229 2839 356 1812 473 1566 283 2483 526 1776 323 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て x,aaa // 分母　現員医師数 b,bbb // 　必要医師数 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 y=(x+b) // 分子 　 現員医師数＋必要医師数 s=(y/x) // 比率 医師不足度 　（現員医師数＋必要医師数）／現員医師数 ................... P=:ci(y) // 個体識別文字列P作成 =pr*(y,x,s,P) // 数値プリント .................... q=cum(x) // xの累和 q<i>=x<1>+x<2>+...+x<i-1>+x<i> r=(q-x) // 直前までの累和 r<i>=x<1>+x<2>+...+x<i-1> =q<i>-x<i> .................... h=(1.0) // ★ h 目安となる医師不足度1.0　つまり１倍　　 .=(0,h) // スカイライン図上の医師不足度1.0％の横線 y=0*x+h の右辺係数の関数「.」 ................... B=@.s(b) // 必要医師数の総和　つまり全国必要医師数 スカラー X=@.s(x) // 現員医師数の総和　つまり全国現員医師数 スカラー H=(B/X) // 全国の医師不足率 　　全国必要医師数／全国現員医師数 +=(H,0) // 扇形散布図上の比率(H)斜線 b=H*x+0 の右辺係数の関数「+」 z=(0*x) // すべてゼロの数値の変量zを作成（扇形散布図の原点に利用） =pr*(X,B,H) // 数値プリント --------------------- $r // 回帰分析 ,run,b=(x) // bを被説明（従属）変数とし，xを説明（独立）変数とする回帰 =============== $$g // グラフセクション $d // 表示範囲 all // 全範囲 $g // スケールの目盛り指示コマンド（標準10ポイント） s,001 // 変量sの目盛りを細かく1ポイントごとに y,001 x,001 $z // ゼロ軸表示 syx // 変量s,y,xのゼロ軸表示 $p // プロット xy,s // 変量x,yを同一スケール，変量sを別スケール --------------------- $3 // ３次元図　医師不足度の スカイライン図 s,q, ,P,.,* // 縦軸s,横軸q,奥行軸なし,個体識別P,関数.,合成用保存* s,r, ,P,* // 縦軸s,横軸r,奥行軸なし,個体識別P,合成用保存* // 合成 比率スカイライン図（リンク面描画，３次元図圧縮） ........ $3 // ３次元図　 医師不足率の 扇形散布図 b,x, ,P,+,* // 縦軸b,横軸x,奥行軸なし,個体識別P,関数+,合成用保存* z,z, ,P,* // 縦軸z,横軸z,奥行軸なし,個体識別P,合成用保存*【原点】 // 合成（２次元図上の散布点と原点のリンク，３次元図圧縮を利用） ======================= $$ // 終了セクション