=============== ternary-trans-fatty-acids-nutrition-jp ============ ==== 「国民健康・栄養調査報告」2007年の一人一日当たり主要食品摂取の ==== トランス脂肪酸を含む脂肪酸構成の三色三角バブルグラフ ========================================================================== ================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0018.00,aa // ケース始点,終点番号 ,第１系列変量名;単位　飽和脂肪酸 ,bb // 空白で同一ケース範囲,第２系列変量名;単位　トランス脂肪酸 ,cc // 空白で同一ケース範囲,第３系列変量名;単位　その他不飽和脂肪酸 ----------------- データ入力指示コマンド $d ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストされたい． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 0.4309 0.0528 0.7701 0.1995 0.0118 0.3205 0.3662 0.0054 0.4020 0.4345 0.0095 1.7208 1.9016 0.0813 3.2290 0.0560 0.0079 0.0429 2.1180 0.0827 0.8172 0.3840 0.0198 0.1575 0.4081 0.0096 0.1710 0.5777 0.0400 0.2928 0.5045 0.0195 0.1816 0.2405 0.0770 0.5253 0.8666 0.1102 6.3710 0.0393 0.0014 0.0520 0.2851 0.0467 0.4345 0.1863 0.0269 0.1286 0.2548 0.0323 1.8737 0.2275 0.0421 2.0976 ================== 変量分析セクション $$v ----------------- 変量記号割当 $a a,aa //　飽和脂肪酸 b,bb // トランス脂肪酸 c,cc // その他不飽和脂肪酸 ----------------- $d // 数値出力範囲 all // 全範囲 ------------------ $t // 変数変換コマンド ------------------- ★ X,Y,Zの各変量と上記のa,b,cの入力変量とを対応させる Y=(a) //　飽和脂肪酸 X=(b) // トランス脂肪酸 Z=(c) // その他不飽和脂肪酸 ....... S=(X+Y+Z) // 脂肪酸総量 S ....... L=(X) // ★バブル変量として選択 トランス脂肪酸 ....... x=(X/S)*100 // トランス脂肪酸構成比 x y=(Y/S)*100 // 飽和脂肪酸構成比％ y z=(Z/S)*100 // その他不飽和脂肪酸構成比％ z p=:ci(x) // データの散布点印字用の文字系列p ............. =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 // 変数変換コマンド f=(-1,-1,+100) // 関数f y= -x -z +100 （つまり x+y+z = 100） ............. i=(100,50,0,0,0,50) // 三角形の頂点と中点の座標 j=(0,50,100,50,0,0) k=(0,0,0,50,100,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) // 数値プリントしてチェックするには先頭....を取る ------------------------　三角グラフ平面用に変換 Y=(y) // データの三角グラフ平面への縦軸変換 X=(2*x+y)/1.7320508 // データの三角グラフ平面への横軸変換 U=(1.732,0) // 関数U Y=1.732X+0 V=(-1.732,200) // 関数V Y=-1.732X+200 ............. 三角形の頂点と中点の座標 J=(j) // 三角形の頂点と中点の三角グラフ平面への縦軸変換 I=(2*i+j)/1.7320508 // 三角形の頂点と中点の三角グラフ平面への横軸変換 ------------- 中三角形の頂点の座標 a=(0,0,75) // ★小三角形の場合は a=(0,0,50) b=(0,75,0) // ★ b=(0,50,0) c=(100,25,25) // ★ c=(100,50,50) .... // @は @=(0*x) として定義済みで，原点の変量（ケースの数はデータ分） a=(@,a) // 原点の変量と中三角形の頂点を連結した変量 b=(@,b) c=(@,c) --------------------- B=(b) // 中三角形の頂点の三角グラフ平面への縦軸変換 A=(2*a+b)/1.7320508 // 中三角形の頂点の三角グラフ平面への横軸変換 v=(-1.732,150) // 中三角形の右辺関数v Y=-1.732X+ (75*2) ★小三角形の場合 v=(-1.732,100) ========================= グラフセクション $$g ----------------- ゼロ軸表示 $z xyzXY // 変量xyzXY についてゼロ軸表示 ----------------- 目盛 $g X,001 // X変量の目盛 1間隔（標準は10間隔） Y,001 // X変量の目盛 1間隔（標準は10間隔） ----------------- ３次元図 $3 // 　三角グラフ立体 j,i,k,Q,* // 縦軸j，横軸i，奥行軸k，散布点印字Q，合成用保存* y,x,z,p=L,f,* // 縦軸y，横軸x，奥行軸z，印字p=バブル変量L，関数f，合成用保存* // 合成 $3 // 　三角グラフ平面 J,I, ,Q,* // 縦軸J，横軸I，奥行軸なし，印字Q，合成用保存* Y,X, ,p=L,U,V,* // 縦軸Y,横軸X,奥行軸なし,印字p=バブル変量L,関数U,V,合成用保存* // 合成 $3 // 　中三角グラフ平面 B,A, ,Q,* // 縦軸B，横軸A，奥行軸なし，印字Q，合成用保存* Y,X, ,p=L,U,v,* // 縦軸Y，横軸X，奥行軸なし,印字p=バブル変量L,関数U,v,合成用保存* // 合成 ========================= 終了セクション $$