=========== distance3-centerkobe-pref-trans-fatty-acids-total-cow-oil ===== ====== 県庁所在全市のトランス脂肪酸摂取量【計】を【牛由来食品】【液体食品】 ====== トランス摂取量と神戸市から距離で説明する重回帰の推定値のグラフ =============================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0051.00,aaa // ケース始点,終点番号 ,★第１系列名 ,bbb // 空白で同一ケース範囲,★第２系列名 ,ccc // 空白で同一ケース範囲,★第３系列名 ,ddd // 空白で同一ケース範囲,★第４系列名 ,eee // 空白で同一ケース範囲,★第５系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 43.064 141.347 162.78 45.98 61.82 40.824 140.740 164.17 38.25 67.77 39.704 141.153 143.72 47.10 50.34 38.269 140.872 149.12 44.81 54.86 39.719 140.102 152.52 35.93 64.07 38.240 140.364 168.62 50.25 63.59 37.750 140.468 147.78 37.27 68.41 36.342 140.447 159.85 43.14 67.71 36.565 139.883 156.94 43.04 66.16 36.391 139.060 156.92 40.22 57.90 35.857 139.649 168.54 51.27 61.51 35.605 140.123 173.47 57.34 54.00 35.689 139.692 160.05 52.44 49.23 35.448 139.642 176.53 55.26 58.36 37.903 139.023 155.10 44.88 48.34 36.695 137.211 141.20 45.28 43.52 36.595 136.626 152.20 45.25 52.75 36.065 136.222 141.48 44.62 49.12 35.664 138.568 172.58 44.56 68.19 36.651 138.181 186.98 46.16 83.25 35.391 136.722 160.18 44.70 53.50 34.977 138.383 154.26 44.94 56.48 35.180 136.907 157.98 47.27 49.77 34.730 136.509 157.21 51.44 56.02 35.005 135.869 174.58 53.25 55.48 35.021 135.756 187.69 58.19 59.21 34.686 135.520 171.01 53.94 48.44 34.691 135.183 190.92 58.37 52.34 34.685 135.833 191.39 60.67 61.38 34.227 135.167 166.42 50.28 50.15 35.504 134.238 172.04 50.20 61.60 35.472 133.050 173.31 42.63 68.34 34.662 133.934 162.49 51.71 51.07 34.397 132.460 179.72 56.89 59.28 34.186 131.471 169.12 47.34 60.57 34.066 134.559 173.35 52.63 60.34 34.340 134.043 162.06 53.84 47.73 33.842 132.766 176.05 51.13 64.37 33.560 133.531 148.66 37.65 58.38 33.607 130.418 152.39 48.99 50.18 33.249 130.300 151.73 46.80 57.28 32.745 129.874 169.22 44.78 63.26 32.790 130.742 157.08 47.52 60.50 33.238 131.613 166.73 49.92 60.65 31.911 131.424 144.44 37.61 63.49 31.560 130.558 175.01 42.25 80.22 26.212 127.681 147.40 39.63 60.80 35.531 139.703 166.13 49.14 56.50 34.711 137.726 142.48 42.91 48.08 34.573 135.483 176.62 49.65 57.03 33.883 130.875 180.38 52.61 63.10 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 c,ccc // トランス脂肪酸摂取量【計】 d,ddd // 【牛由来食品】トランス脂肪酸摂取量 e,eee // 【液体食品】トランス脂肪酸摂取量 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 P=:ci(a) // 個体識別文字列P作成 U=@."(a)28. // ★神戸市（順序数28番）の緯度のスカラー抽出 定数数字後の.は必須 V=@."(b)28. // ★神戸市（順序数28番）の経度のスカラー抽出 定数数字後の.は必須 A=abs(a-U) // 各市の緯度の神戸市からの乖離 B=abs(b-V) // 各市の経度の神戸市からの乖離 =pr*(a,b,c,d,e,A,B,P) // 数値プリント ..................... 起点（神戸市）から各終点(各都市）までの距離の計測処理 n=cos(U/180*3.14159) // 起点緯度（ラジアン）のcos l=cos(a/180*3.14159) // 終点緯度（ラジアン）のcos m=cos((U+a)/2/180*3.14159) // 起点終点の平均緯度（ラジアン）のcos　 M=(1+((6378.137**2-6356.752**2)/6356.752**2)*(m**2))**0.5 M=(6378.137**2/6356.752)/(M**3) // 起点終点の平均緯度の子午線曲率半径 N=(1+((6378.137**2-6356.752**2)/6356.752**2)*(n**2))**0.5 N=(6378.137**2/6356.752)/N // 起点緯度の卯酉線曲率半径 L=(1+((6378.137**2-6356.752**2)/6356.752**2)*(l**2))**0.5 L=(6378.137**2/6356.752)/L // 終点緯度の卯酉線曲率半径 =pr*(M,N,L) // 数値プリント h=(A/180*3.14159*M) // 南北の緯度差（ラジアン）の子午線弧長 i=(B/180*3.14159*N*n) // 起点緯度の東西の経度差（ラジアン）弧長 j=(B/180*3.14159*L*l) // 終点緯度の東西の経度差（ラジアン）の弧長 t=(h**2+i*j)**0.5 // 起点（神戸市）から各終点までの距離　km T=(t*1000) // 同上 m 単位 =pr*(h,i,j,t,P,T) // 数値プリント ---------------- $r // 回帰分析 トランス摂取量【計】cを F,@"*,c=(d,e) // 【牛由来食品】トランスd【液体食品】トランスeで説明する回帰係数F C,est,c=(d,e,t) // 【牛由来食品】d【液体食品】eと神戸市からの距離tで説明する推定値C ,pr*,c,C,c-C,d,e,t,P // 数値プリント =============== $$g // グラフセクション $c // 散布図　横軸に神戸市からの距離 c,t,*,P // 縦軸c,横軸t（距離）,回帰線描画*,個体識別P トランス【計】 d,t,*,P // 縦軸d,横軸t（距離）,回帰線描画*,個体識別P 【牛由来】トランス e,t,*,P // 縦軸e,横軸t（距離）,回帰線描画*,個体識別P 【液体食品】トランス $3 // ３次元図　縦軸にトランス【計】,横軸に【牛由来食品】,奥行軸に【液体食品】 c,d,e,P=t,F,* // 縦軸c,横軸d,奥行軸e,個体識別P=バブル変量t,関数F,合成用保存* C,d,e,P,* // 縦軸C,横軸d,奥行軸e,個体識別P,合成用保存* // 合成 ======================= $$ // 終了セクション