=============== skyline-fatty-acids-uc ============ ============ ユーザデータセクション $$u ----------------- $c // クロスセクションデータ属性コマンド ---------- クロスセクションでは県名や企業名等の文字データを扱うことも多い． ---------- 各文字変量には漢字２文字（英字４文字）のみ入力される． ---------- それを超える文字は無視される． ---------- 文字系列変量名の先頭は「:n1,」「:n2,」…「:n6,」のいずれかを用いる ---------- 識別文字系列変量名の先頭は「:ci,」を用いる 0001.00,0012.00,aa // ケース始点,終点番号 ,数値系列変量名;単位　飽和脂肪酸 ,bb // 空白で同一ケース範囲,数値系列変量名;単位　一価不飽和脂肪酸 ,cc // 空白で同一ケース範囲,数値系列変量名;単位　多価不飽和脂肪酸 ----------------- データ入力指示コマンド $d ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストされたい． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 1.33 1.50 1.04 0.14 0.05 0.31 0.28 0.11 0.61 1.26 1.86 1.55 0.51 0.20 1.12 0.62 0.82 0.90 0.29 0.21 0.31 0.55 0.52 0.69 0.40 0.42 0.80 0.49 0.70 1.37 0.25 0.10 0.56 0.25 0.19 0.28 ================== 変量分析セクション $$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,+1) // ３次元関数k b= -a -c +1 （つまり a+b+c = 1）★ +1　変更対象 l=(-1,-1,+2) // ３次元関数l b= -a -c +2 （つまり a+b+c = 2）★ +2　変更対象 m=(-1,-1,+3) // ３次元関数m b= -a -c +3 （つまり a+b+c = 3）★ +3 変更対象 n=(-1,-1,+4) // ３次元関数n b= -a -c +4 （つまり a+b+c = 4）★ +4　変更対象 o=(-1,-1,+5) // ３次元関数o b= -a -c +5 （つまり a+b+c = 5）★ +5　変更対象 $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,合成用保存【原点】 // 合成（３次元図上の散布点の垂線,バブル,塗りつぶし色,原点とのリンクを利用） ======================= $$ // 終了セクション