# SAS演習 ## グループ5 * **4618055 時廣夢** * **4618049 田村彰汰** * **4618031 小林俊哉** * **4618068 眞玉大貴** * **4618029 黒木裕介** ### 課題1 > 連続量のデータ「age」,「sbp_b」,「dbp_b」,「sbp_a」,「dbp_a」の基本統計量(データ数,平均,標準偏差,中央値,最小値,最大値)を群ごとに計算した結果を以下に示す。 ```sas= PROC MEANS DATA=yobi n mean std median min max; VAR ; class group; RUN; ``` ![](https://i.imgur.com/2J0pJOz.png) ----- > 2値変数のデータ「sex」の各水準の度数と割合を群ごとに計算した結果を以下に示す。 ```sas= /* 課題1 */ proc freq data = yobi; table group * sex / nocol nopercent; run; ``` ![](https://i.imgur.com/SQKhGtv.png) ### 課題2 ```sas= /* 課題2 */ proc sgscatter data=yobi; plot age * sbp_a / reg; run; proc corr data=yobi; var age sbp_a; run; proc sgscatter data=yobi; plot sbp_b * sbp_a / reg; run; proc corr data=yobi; var sbp_b sbp_a; run; proc sgscatter data=yobi; plot dbp_b * sbp_a / reg; run; proc corr data=yobi; var dbp_b sbp_a; run; proc sgscatter data=yobi; plot dbp_a * sbp_a / reg; run; proc corr data=yobi; var dbp_a sbp_a; run; ``` ![](https://i.imgur.com/cGg4jnu.png) ![](https://i.imgur.com/FeQEoKd.png) ---- ![](https://i.imgur.com/7C0c8Gm.png) ![](https://i.imgur.com/Zs3itsR.png) ---- ![](https://i.imgur.com/8KHh8uP.png) ![](https://i.imgur.com/VCdz06q.png) ---- ![](https://i.imgur.com/fbG0uJs.png) ![](https://i.imgur.com/ZlbvoQG.png) 変数`age`と`dbp_a`は、`sbp_a`変数との相関係数$r > 0.5$なので、この二つの変数は2か月後の収縮期血圧とある程度関連があると考えられる。 ### 課題3 ```sas= proc sgscatter data=yobi; plot sbp_a * group /reg; run; ``` ![](https://i.imgur.com/lOzo0FJ.png) ### 課題4 課題1の結果より、収縮期血圧と年齢のバランスは取れているが、性別のバランスは崩れている。 ### 課題5 ```sas= /* 課題5 */ proc ttest data=yobi; var sbp_a sbp_b; class group; run; ``` ![](https://i.imgur.com/fxw55N9.png) ![](https://i.imgur.com/sgNvI8t.png) ### 課題6 ```sas= /* 課題6 */ data yobi; set yobi; diff=sbp_a-sbp_b; run; PROC Ttest data=yobi; var diff; class group; run; ``` ![](https://i.imgur.com/h7mqf5m.png) ![](https://i.imgur.com/TgRBimW.png) ### 課題7 * 主要変数 : 収縮期血圧の変化量 ```sas= /* 課題7 */ proc power; twosamplemeans test = diff meandiff = -10 stddev = 10 power = 0.8 npergroup = .; run; ``` ![](https://i.imgur.com/Y6VvOOi.png) ### 課題8 ```sas= PROC MEANS DATA=main n mean std median min max; VAR ; class group; RUN; ``` ![](https://i.imgur.com/gFIF9ro.png) ```sas= proc freq data = main; table group * sex / nocol nopercent; run; ``` ![](https://i.imgur.com/m6KYVZN.png) ### 課題9 課題8の結果より、収縮期血圧、年齢、性別ともにバランスが取れていることがわかる。 ### 課題10 ```sas= /* 課題10 */ data main; set main; diff2=sbp_a-sbp_b; run; proc ttest data=main; var diff2; class group; run; ``` ![](https://i.imgur.com/Gm4W15O.png) ![](https://i.imgur.com/AvnQyky.png) 「2群間の主要変数(収縮期血圧の変化量)の母平均に差がない」という帰無仮説のもとt検定を行った結果、P値が0.0001未満であり、有意水準5%のもとで帰無仮説を棄却できる。つまり、2群間の主要変数の母平均に差があると言える。 また、得られた母平均の差$>10$であり十分な差があるため、臨床的に意味があると言える。 ### 課題11 「臨床的に意味がある差があるときには,母平均に差がないという帰無仮説をできるだけ高い確率で棄却する」という標本サイズ を設計し、被験薬の降圧効果を確かめることができた。よって、本試験の目的が達成できたと言える。