SAS探討慢性B型肝炎病毒感染與糖尿病及代謝症候群之關聯性HW

一、研究背景
過去已有研究指出,慢性C型肝炎病毒感染與糖代謝混亂、胰島素阻抗甚至是第二型糖尿病有關。由於胰島素阻抗是導致代謝症候群的主要原因,因此也有結果指出,慢性C型肝炎病毒感染會增加其發生代謝症候群的危險性,但這些結果,在慢性B型肝炎病毒感染者身上尚未看到一致性。
本研究希望藉由台灣B型肝炎高盛行率的優勢,利用健檢資料(包括身高體重、血壓、血液生化值等),使用橫斷研究以探討慢性B型肝炎病毒感染與糖尿病及代謝症候群之關聯性。本資料庫包含1,000位民眾以進行分析討論。

健檢資料庫內容包括:
1.  基本人口學資料:性別、年齡。
2.  體位測量:身高、體重、收縮壓、舒張壓、腰圍、臀圍。
3.  生化測量:空腹血糖、總膽紅素、GOTGPT、肌酸酐、尿酸、三酸甘油脂、總膽固醇、高密度脂蛋白膽固醇、低密度脂蛋白膽固醇。
4.  慢性發炎指標:C反應蛋白。
5.  肝炎指標:B肝表面抗原。
6.  生活習慣:抽菸、飲酒。
7.  個人用藥史:是否長期服用高血壓藥物、是否長期服用糖尿病藥物、是否長期服用高血脂藥物。
注意:以下5項危險因子中,若包含3項或以上者,即可判定為代謝症候群(Metabolic Syndrome)
危 險 因 子
異 常 值
腹部肥胖(Central obesity)
腰圍(Waist Circumference, WC:
男性 90 cm ;女性 80 cm
血壓(BP)
收縮壓130 mm/Hg或舒張壓85 mm/Hg
或已服用降高血壓藥物
高密度脂蛋白膽固醇(HDL-C)
男性 <40 mg/dl ; 女性 <50 mg/dl
空腹血糖值(Fasting glucose)
100 mg/dl或已服用治療糖尿病的藥物
三酸甘油酯(Triglyceride)
150 mg/dl或已服用治療高血脂的藥物
 (國民健康局,2007)
二、資料庫變項說明
Variable
Variable description
Codes
sex
性別
0=女生
1=男生
age
年齡

height
身高

weight
體重

sbp
收縮壓

dbp
舒張壓

glucose
空腹葡萄糖值

tftbil
總膽紅素

tfsgot
肝功能指數(GOT)

tfsgpt
肝功能指數(GPT)

ufcre
腎功能指標(肌酸酐)

ufua
尿酸值(UA)

bgtg
三酸甘油脂

bgchol
總膽固醇值

crp
C反應蛋白

bghdlc
高密度脂蛋白膽固醇

bgldlc
密度脂蛋白膽固醇

wawidth
腰圍

hipwidth
臀圍

smokeornot
抽菸
1=從未           2=二手菸
3=戒菸           4=每週1-3
5=每週4-5      6=每天
drinkornot
飲酒
1=從未           2=戒酒
3=每週1-2      4=每週3-4      5=每天
mdrug04
使用高血壓藥物
T:有服用        F:無服用
mdrug05
使用糖尿病藥物
T:有服用        F:無服用
mdrug07
使用高血脂藥物
T:有服用        F:無服用
hbv
B型肝炎表面抗原(HBsAg)
0=無資料         1=陽性
2=陰性




































1.          本研究的母群體與研究樣本分別為何?
2.          本研究中的主要「自變項」是什麼?主要「依變項」是什麼?
3.          請描述有無代謝症候群者,基本人口學資料體位測量生化測量肝炎指標生活習慣及個人用藥史的分佈情形。
4.          本研究中B型肝炎與代謝症候群兩者的分佈為何?請畫圖並描述之。
5.          承上,請問代謝症候群中B型肝炎的比例為何?
6.          承上,請問B型肝炎患者中得到代謝症候群的比例為何?
7.          本研究中不同年齡層與代謝症候群兩者的分佈為何?請畫圖並描述之。
8.          請利用身高體重計算出BMI,並依照台灣對於肥胖的定義來描述分佈情況。(可以參考衛生福利部國民健康署網站)
9.          本研究中肥胖程度與代謝症候群兩者的分佈為何?請畫圖並描述之。
10.     請利用圖表方式呈現年齡與膽固醇之間的關係。
---------------------------------------------------------------------
/*讀資料,把資料異常的弄成missing*/
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";
/* MetabolicSyndrome=1有代謝症候群/=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;
/*計算出BMI BMI = 體重(公斤) / 身高2(公尺2)*/
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 print data= aa ;
run;

/*基本人口學資料(性別、年齡)*/
proc freq data=aa ; 
table (sex age)*MetabolicSyndrome  ;
RUN;

/*體位測量(身高、體重、收縮壓、舒張壓、腰圍、臀圍)*/
proc means data=aa;
class MetabolicSyndrome;
var  height weight  sbp  dbp  wawidth  hipwidth  ;
RUN;

/*生化測量(空腹血糖、總膽紅素、GOT、GPT、肌酸酐、
尿酸、三酸甘油脂、總膽固醇、高密度脂蛋白膽固醇、低密度脂蛋白膽固醇)*/
proc means data=aa;
class MetabolicSyndrome;
var glucose tftbil tfsgot tfsgpt  ufcre  ufua bgtg  bgchol bghdlc bgldlc;
RUN;

/*肝炎指標(B肝表面抗原)*/
proc freq data=aa ; 
table (hbv)*MetabolicSyndrome  ;
RUN;

/*生活習慣(抽菸、飲酒)*/
proc freq data=aa ; 
table (smokeornot drinkornot )*MetabolicSyndrome  ;
RUN;

/*個人用藥史(是否長期服用高血壓藥物、*/
/*是否長期服用糖尿病藥物、是否長期服用高血脂藥物)*/
proc freq data=aa ; 
table (mdrug04 mdrug05 mdrug07)*MetabolicSyndrome  ;
RUN;

/*B型肝炎與代謝症候群兩者的分佈*/
proc freq data=aa ; 
table (hbv)*MetabolicSyndrome;
RUN;

PROC SGPANEL DATA=aa NOAUTOLEGEND;
PANELBY MetabolicSyndrome;
HISTOGRAM hbv;
DENSITY hbv;
RUN;

/*不同年齡層與代謝症候群兩者的分佈*/
proc freq data=aa ; 
table (age)*MetabolicSyndrome;
RUN;

PROC SGPANEL DATA=aa NOAUTOLEGEND;
PANELBY MetabolicSyndrome;
HISTOGRAM age;
DENSITY age;
RUN;

/*肥胖的定義來描述分佈情況*/
proc freq data=aa ; 
table BB;
RUN;

/*肥胖程度與代謝症候群*/
proc freq data=aa ; 
table BB*MetabolicSyndrome;
RUN;

PROC SGPANEL DATA=aa NOAUTOLEGEND;
PANELBY MetabolicSyndrome;
HISTOGRAM BB;
DENSITY BB;
RUN;

/*年齡與膽固醇之間的關係*/
PROC GPLOT;
PLOT bghdlc*age;
RUN;

PROC GPLOT;
PLOT bgldlc*age;
RUN;

留言