==== skylineRATIO2-trans-fatty-acids-consume-kobe-tokyo ========= ==== 2009年家計世帯一人当たり主要食品トランス脂肪酸 ==== 神戸市・東京都区部の対全国比の合成スカイライン図・扇形散布図 ========================================================================== ================== $$u // ユーザデータ・セクション ---------------------　神戸市分 ------------------------ $c // クロスセクションデータ属性コマンド 0001.00,0012.00,aaa // ケース始点,終点番号 ,★第１系列名 ,bbb // 空白で同一ケース範囲,★第２系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 8282 0.163 10071 0.204 674 0.128 686 0.128 3492 0.521 31070 0.091 161 0.024 212 1.951 936 0.826 2591 1.395 634 7.004 1309 1.237 ---------------------　東京都区部分 ------------------------ $c // クロスセクションデータ属性コマンド 0001.00,0012.00,ccc // ケース始点,終点番号 ,★第３系列名 ,ddd // 空白で同一ケース範囲,★第４系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 6817 0.163 7507 0.204 888 0.128 657 0.128 2143 0.521 30699 0.091 182 0.024 233 1.951 1059 0.826 2196 1.395 428 7.004 1503 1.237 ---------------------　全国分 ------------------------ $c // クロスセクションデータ属性コマンド 0001.00,0012.00,GGG // ケース始点,終点番号 ,★第３系列名 ,HHH // 空白で同一ケース範囲,★第４系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 6399 0.163 6932 0.204 964 0.128 811 0.128 2261 0.521 28202 0.091 144 0.024 156 1.951 770 0.826 2672 1.395 446 7.004 1572 1.237 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 　神戸市の購入数量 b,bbb // 　　　含有率％ c,ccc // 　東京都区部の購入数量 d,ddd // 　　　含有率％ 【実際には含有率cはbと同じ】 G,GGG // 　全国の購入数量 H,HHH // 　含有率％ 【実際には含有率Hはbと同じ】 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 u=(a*(b/100)) // 分子 　神戸市の摂取量 = 神戸市の購入数量 * （含有率％ / 100） v=(c*(d/100)) // 分子 　東京都区部の摂取量 = 東京都区部の購入数量 * （含有率％ / 100） Y=(G*(H/100)) // 分母 　全国の摂取量 = 全国の購入数量 * （含有率％ / 100） k=(u/Y*100) // 神戸市の対全国比　％ l=(v/Y*100) // 東京都区部の対全国比　％ ................... p=:ci(u) // 個体識別文字列p作成 =pr*(a,c,G,H,u,v,Y,k,l,p) // 数値プリント .................... i=@.a(Y) //　全国摂取量Yの食品別の平均値@.aのスカラーi（区切りに利用） I=:ci(i)* //　スカラーiに文字 "*" の文字列変量I　作成（スカイライン区切りに利用） ?Y=(Y,i,Y) //　全国摂取量　全国分Y＋スカラーi＋全国分Y の連結変量?Y ?K=(k,0,l) //　対全国比 神戸市分k＋数値0＋東京分lの連結変量?K ?P=(p,I,p) //　文字列変量 神戸市分p＋文字列I＋東京分pの連結変量?P q=cum(?Y) // 分母変量?Yの累和 q<i>=?Y<1>+?Y<2>+...+?Y<i-1>+?Y<i> r=(q-?Y) // 直前までの累和 r<i>=?Y<1>+?Y<2>+...+?Y<i-1> =q<i>-?Y<i> .................... h=(100) // ★ h 対全国比 100％　　 .=(0,h) // スカイライン図上の対全国比 100％の横線 y=0*x+h の右辺係数の関数「.」 +=(h/100,0) // 扇形散布図上の対全国比( h％ / 100)斜線 y=(h/100)*x+0 の右辺係数の関数「+」 z=(0*u) // すべてゼロの数値の変量zを作成（扇形散布図の原点に利用） ---------------- $r // 回帰分析 ,run,u=(Y) // uを被説明（従属）変数とし，Yを説明（独立）変数とする回帰 ,run,v=(Y) // vを被説明（従属）変数とし，Yを説明（独立）変数とする回帰 =============== $$g // グラフセクション $d // 表示範囲 all // 全範囲 $g // スケールの目盛り指示コマンド（標準10ポイント） ?K,001 // 変量?Kの目盛りを細かく1ポイントごとに u,001 $z // ゼロ軸表示 ?KuvY // 変量?K,u,v,Yのゼロ軸表示 $p // プロット uvY,kl // 変量u v Y,別スケールでk l ---------------------------------- $3 // ３次元図　スカイライン図【対全国比】分母：全国摂取量 ?K,q, ,?P,.,* // 縦軸?K,横軸q,奥行軸なし,個体識別?P,関数.,合成用保存* ?K,r, ,?P,* // 縦軸?K,横軸r,奥行軸なし,個体識別?P,合成用保存* // 合成 比率スカイライン図（リンク面描画，３次元図圧縮） ........ $3 // ３次元図　扇形散布図 u,Y, ,p,+,* // 縦軸u,横軸Y,奥行軸なし,個体識別p,関数+,合成用保存* 【神戸市・全国摂取量】 z,z, ,p,* // 縦軸z,横軸z,奥行軸なし,個体識別p,合成用保存* 【原点】 v,Y, ,p,+,* // 縦軸v,横軸Y,奥行軸なし,個体識別p,関数+,合成用保存* 【東京都区部・全国摂取量】 // 合成（２次元図上の散布点と原点のリンク，３次元図圧縮を利用） ======================= $$ // 終了セクション