1.何謂type I error、type II error和power。
若 H0（虛無假說）為真，但結論卻否決 H0，則犯了第一型錯誤(Type I Error )。
若 H1（對立假說）為真，但結論卻接受 H0，則犯了第二型錯誤(Type II Error)。
發生type II error的機率稱β，而當1-β時，就是我們拒絕錯誤虛無假設，而做出正確統計推論的機率，也稱為檢定的強度power，以此定義，若β越小的話，power就越大。
討論代謝症候群與Ｂ型肝炎的關係，請問可以利用何種假設檢定？
卡方檢定(因為兩者皆為獨立的類別變項)
(1). 承上，虛無假設與對立假設為何？
虛無假設是代謝症候群與B型肝炎沒有關係；對立假設是代謝症候群與B型肝炎有關係。
(1). 統計結果為何？並解釋之。
卡方檢定的P值0.0016小於0.05 ，表示
代謝症候群與B型肝炎有顯著關係。
(2). 若此樣本未達常態分佈，應該改用何種統計方法?
Kruskal-Wallis tset
1. 欲探討年齡與代謝症候群之關係，請問可以利用何種假設檢定？T test
(1). 承上，虛無假設與對立假設為何？
虛無假設是年齡與代謝症候群無關係；年齡與代謝症候群有相關。
(2). 統計結果為何？並解釋之。
小於0.05 有顯著，所以是拒絕虛無假設，代表年齡與代謝症候群有相關。
(3). 若此樣本未達常態分佈，應該改用何種統計方法?
無母數分析Mann-Whitney U-test
2. 欲探討肥胖對於空腹葡萄糖之關係，請問可以利用何種假設檢定？ANOVA檢定。
(1). 承上，虛無假設與對立假設為何？虛無假設肥胖對於空腹葡萄糖無關係；肥胖對於空腹葡萄糖有關係。
(2). 統計結果為何？並解釋之。
小於0.05 ，拒絕虛無假設，所以虛無假設肥胖對於空腹葡萄糖有關係。
(3). 若此樣本未達常態分佈，應該改用何種統計方法?
Post-hoc: LSD
3. 請問本研究之樣本數是否足以說明代謝症候群與Ｂ型肝炎的關係。
是，樣本數足以說明代謝症候群與Ｂ型肝炎的關係，因為每個細格小於五的，不到全部的 20%。
(如有錯誤，請多指正謝謝)
---------------------------------------------------------------
data aa ;
set work.b;
if height =0 then height = .;
if weight =0 then weight=.;
if sbp = 0 then sbp =.;
if dbp = 0 then dbp=.;
if glucose = 0 then glucose=.;
if tftbil = 0 then tftbil=.;
if tfsgot = 0 then tfsgot=.;
if tfsgpt = 0 then tfsgpt=.;
if ufcre = 0 then ufcre=.;
if ufua = 0 then ufua=.;
if bgtg = 0 then bgtg=.;
if bgchol = 0 then bgchol=.;
if bgtg = 0 then bgtg=.;
if crp = 0 then crp=.;
if bghdlc = 0 then bghdlc=.;
if bgldlc = 0 then bgldlc=.;
if wawidth = 0 then wawidth=.;
if hipwidth = 0 then hipwidth=.;
if (smokeornot=0 or smokeornot >6 )then smokeornot=.;
if (drinkornot=0 or drinkornot>5) then drinkornot=.;
if (hbv=0 or hbv>2 ) then hbv=. ;
MS=0;
MetabolicSyndrome="0";
if(sex=0 and wawidth >=80 ) then MS=MS+1;
if(sex=1 and wawidth >=90 ) then MS=MS+1;
if(sbp >=130 or dbp >=85 or mdrug04='T' ) then MS=MS+1;
if(sex=1 and bghdlc<40 ) then MS=MS+1;
if(sex=0 and bghdlc<50) then MS=MS+1;
if( glucose>=100 or mdrug05='T') then MS=MS+1;
if(bgtg >=150 or mdrug07='T') then MS=MS+1;
if(MS>=3) then MetabolicSyndrome="1";
if(MS<3) then MetabolicSyndrome="0";
drop MS;
if(weight=. or height=.) then BMI=.;
else BMI=weight/((height/100)**2);
if (BMI=. ) then BB=".";
else if (BMI<18.5 ) then BB=0; /*??*/
else if(18.5 <=BMI< 24 ) then BB=1;/*??*/
else if( 24 <=BMI) then BB=2;/*??*/
RUN;
proc freq data=aa;
table MetabolicSyndrome *hbv / EXPECTED CHISQ;
RUN;
Proc TTEST data=aa;
class MetabolicSyndrome;
VAR age;
RUN;
proc anova data=aa;
class BB;
model glucose=BB;
RUN;
--------------------------------------------------------------
留言
張貼留言