=============== skyline-nutrients-uc ============ ============ ユーザデータセクション $$u ----------------- $c // クロスセクションデータ属性コマンド ---------- クロスセクションでは県名や企業名等の文字データを扱うことも多い． ---------- 各文字変量には漢字２文字（英字４文字）のみ入力される． ---------- それを超える文字は無視される． ---------- 文字系列変量名の先頭は「:n1,」「:n2,」…「:n6,」のいずれかを用いる ---------- 識別文字系列変量名の先頭は「:ci,」を用いる 0001.00,0012.00,aa // ケース始点,終点番号 ,数値系列変量名;単位　蛋白質 ,bb // 空白で同一ケース範囲,数値系列変量名;単位　脂質 ,cc // 空白で同一ケース範囲,数値系列変量名;単位　炭水化物 ----------------- データ入力指示コマンド $d ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストされたい． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 8.3 13.7 54.7 5.4 15.2 23.7 12.9 5.6 63.3 12.4 1.5 84.3 13.4 2.3 87.1 5.7 13.1 44.8 13.2 0.4 72 8.7 1.1 72.7 4.9 3.1 31 19.4 19.7 85.5 17.8 10.4 78.3 13.7 9.3 73.4 ================== 変量分析セクション $$v ----------------- 変量記号割当 $a a,aa //　蛋白質 b,bb // 脂質 c,cc // 炭水化物 ----------------- $d // 数値出力範囲 all // 全範囲 ------------------ $t // 変数変換コマンド ------------------- ★ 分母xとしてa,b,cの項目のいずれか２項目以内を選ぶ ....x=(a) // 蛋白質 これを分母xに選ぶ場合には，先頭....を取る x=(b) // 脂質 　　　　他を分母xに選ぶ場合には，先頭に....を付ける ....x=(c) // 炭水化物 以下　同様 ....x=(a+b) // 蛋白質+脂質 ....x=(a+c) // 蛋白質+炭水化物 ....x=(b+c) // 脂質+炭水化物 ------------------- y=(a+b+c) // ★ ３大栄養素合計 y s=(y)/x // ★ 比率　この場合　 ３大栄養素合計y / 分母要素 x P=:ci(y) // 個体識別文字列P作成 =pr*(y,x,s,a,b,c,P) // 数値プリント ............ ★次行の７カラム目の３大栄養素合計yをxやsへの変更で別変量での並び替え可 j=r.g(y)blank // ３大栄養素合計yの大きい順(定数項blankで欠測値にも末尾の順位)の順位変量j ....j=r.l(y)blank // 小さい順の場合は先頭の....を取る x=pmt(x,j) // 並び替え（順序数変量jによる） y=pmt(y,j) s=pmt(s,j) a=pmt(a,j) b=pmt(b,j) c=pmt(c,j) P=pmt(P,j) P,nam,:ci,P=pmt(P,j) // 個体識別文字系列の変量名の先頭は :ci が必須なので，変量名再設定 =pr*(y,x,s,a,b,c,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) // ★ h 比率 = ３大栄養素合計y / 分母要素 x の参考値として １ .=(0,h) // スカイライン図上の比率 h の横線 y=0*x+h の右辺係数 [0,h] の関数「.」 +=(h,0) // 扇形散布図上の比率 h の斜線 y=h*x+0 の右辺係数 [h,0] の関数「+」 z=(0*y) // すべてゼロの数値の変量zを作成（ 扇形散布図の原点に利用） k=(-1,-1,+20) // ３次元関数k b= -a -c +20 （つまり a+b+c = 20）★ +20　変更対象 l=(-1,-1,+40) // ３次元関数l b= -a -c +40 （つまり a+b+c = 40）★ +40　変更対象 m=(-1,-1,+60) // ３次元関数m b= -a -c +60 （つまり a+b+c = 60）★ +60　変更対象 n=(-1,-1,+80) // ３次元関数n b= -a -c +80 （つまり a+b+c = 80）★ +80　変更対象 o=(-1,-1,+100) // ３次元関数o b= -a -c +100 （つまり a+b+c = 100）★ +100　変更対象 $r // 回帰分析 ,run,y=(x) // yを被説明（従属）変数とし，xを説明（独立）変数とする回帰 ,run,y=(a,b,c) // 被説明変数y，説明変数a,b,cによる重回帰の計測 =============== $$g // グラフセクション $d // 表示範囲 all // 全範囲 $g // スケールの目盛り指示コマンド（標準10ポイント） s,002 // 変量sの目盛りを細かく2ポイントごとに y,002 x,002 $z // ゼロ軸表示 syx // 変量s,y,xのゼロ軸表示 $p // プロット x,y,s // 変量x,y,sをを別スケール ------------ 比率　------------------------------------- $3 // ３次元図　スカイライン図 s,q, ,P,.,* // 縦軸s,横軸q,奥行軸なし,個体識別P,関数.,合成用保存* s,r, ,P,* // 縦軸s,横軸r,奥行軸なし,個体識別P,合成用保存* // 合成 比率スカイライン図（リンク面描画，３次元図圧縮） ........ $3 // ３次元図　扇形散布図 y,x, ,P,+,* // 縦軸y,横軸x,奥行軸なし,個体識別P,関数+,合成用保存* z,z, ,P,* // 縦軸z,横軸z,奥行軸なし,個体識別P,合成用保存【原点】 // 合成（２次元図上の散布点と原点のリンク，３次元図圧縮を利用） ------------ ３要素立体図　------------------------------------- $3 // 　３次元三色虫ピングラフ b,a,c,P=y,k,l,m,n,o,* // 縦軸b,横軸a,奥行軸c,印字P=バブル変量y,関数k,l,m,n,o,合成用保存* z,z,z,P,* // 縦軸z,横軸z,奥行軸z,個体識別P,合成用保存【原点】 // 合成（３次元図上の散布点の垂線,バブル,塗りつぶし色,原点とのリンクを利用） ======================= $$ // 終了セクション