=============== ternary-trans-fatty-acids-consume-jp ============ ==== 全国2009年家計世帯一人当たり主要食品のトランス脂肪酸を含む ==== 脂肪酸構成の三色三角バブルグラフ ========================================================================== ================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0012.00,aa // ケース始点,終点番号 ,第１系列変量名;単位　飽和脂肪酸 ,bb // 空白で同一ケース範囲,第２系列変量名;単位　トランス脂肪酸 ,cc // 空白で同一ケース範囲,第３系列変量名;単位　その他不飽和脂肪酸 ----------------- データ入力指示コマンド $d ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストされたい． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 85.10 10.43 152.10 238.47 14.14 383.07 84.06 1.23 92.27 68.58 1.04 74.75 275.63 11.78 468.02 657.12 25.66 253.54 0.63 0.03 0.27 78.51 3.04 28.26 123.16 6.36 50.53 293.12 37.27 2154.86 97.49 31.24 212.98 105.19 19.45 970.03 ================== 変量分析セクション $$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,合成用保存* // 合成 ========================= 終了セクション $$