*2005-2006*
*body measures variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\BMX_D.XPT";
proc copy in=XP out=NH;
run;
*soy consumption variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DBQ_D.XPT";
proc copy in=XP out=NH;
run;
*cancer variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\MCQ_D.XPT";
proc copy in=XP out=NH;
run;
*Urinary phytoestrogens for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PHYTO_D.XPT";
proc copy in=XP out=NH;
run;
*Menopause variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\RHQ_D.XPT";
proc copy in=XP out=NH;
run;
*Creatinine variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\ALB_CR_D.XPT";
proc copy in=XP out=NH;
run;
*Cotinine variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\COT_D.XPT";
proc copy in=XP out=NH;
run;
*Demographic variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DEMO_D.XPT";
proc copy in=XP out=NH;
run;
*Diabetes variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DIQ_D.XPT";
proc copy in=XP out=NH;
run;
*Dietary information variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\FFQRAW_D.XPT";
proc copy in=XP out=NH;
run;
*Kidney variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\KIQ_U_D.XPT";
proc copy in=XP out=NH;
run;
*Physical activity variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PAQ_D.XPT";
proc copy in=XP out=NH;
run;
*Cadmium variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PBCD_D.XPT";
proc copy in=XP out=NH;
run;
*Smoking variable for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\SMQ_D.XPT";
proc copy in=XP out=NH;
run;
*Obesity/weight history for 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\WHQ_D.XPT";
proc copy in=XP out=NH;
run;
*C-reactive Protein data 2005-2006*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\CRP_D";
proc copy in=XP out=NH;
run;
*Attempt to append and merge data. Lists the contents of the 2005-2006 files*;
libname NH "C:\NHANES\DATA";
proc contents data =NH.BMX_D varnum;
title "2005-2006 Body Measures";
proc contents data =NH.DBQ_D varnum;
title "2005-2006 Diet Behavior & Nutrition";
proc contents data =NH.MCQ_D varnum;
title "2005-2006 Medical Conditions";
proc contents data =NH.PHYTO_D varnum;
title "2005-2006 Urinary Phytoestrogens";
proc contents data =NH.RHQ_D varnum;
title "2005-2006 Reproductive Health";
proc contents data =NH.ALB_CR_D varnum;
title "2005-2006 Creatinine, urine (mg/dL)";
proc contents data =NH.COT_D varnum;
title "2005-2006 Cotinine (ng/mL)";
proc contents data =NH.DEMO_D varnum;
title "2005-2006 Demographics (Race/ethnicity, education level, income)";
proc contents data =NH.DIQ_D varnum;
title "2005-2006 Diabetes Status";
proc contents data =NH.FFQRAW_D varnum;
title "2005-2006 Food Frequency Questionnaire (Dietary Data)";
proc contents data =NH.KIQ_U_D varnum;
title "2005-2006 Kidney/Renal Failure";
proc contents data =NH.PAQ_D varnum;
title "2005-2006 Physical Activity";
proc contents data =NH.PBCD_D varnum;
title "2005-2006 Cadmium (nm/L)";
proc contents data =NH.SMQ_D varnum;
title "2005-2006 Smoking Status";
proc contents data =NH.WHQ_D varnum;
title "2005-2006 Weight Status";
proc contents data =NH.CRP_D varnum;
title "2005-2006 C-Reactive Protein";
run;
*END 2005-2006*
*IMPORTING THE DATA FOR 2007-2008*;
*2007-2008*
*body measures variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\BMX_E.XPT";
proc copy in=XP out=NH;
run;
*soy consumption variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DBQ_E.XPT";
proc copy in=XP out=NH;
run;
*cancer variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\MCQ_E.XPT";
proc copy in=XP out=NH;
run;
*Urinary phytoestrogens for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PHYTO_E.XPT";
proc copy in=XP out=NH;
run;
*Menopause variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\RHQ_E.XPT";
proc copy in=XP out=NH;
run;
*Creatinine variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\ALB_CR_E.XPT";
proc copy in=XP out=NH;
run;
*Cotinine variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\COTNAL_E.XPT";
proc copy in=XP out=NH;
run;
*Demographic variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DEMO_E.XPT";
proc copy in=XP out=NH;
run;
*Diabetes variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DIQ_E.XPT";
proc copy in=XP out=NH;
run;
*Dietary information variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\KIQ_U_E.XPT";
proc copy in=XP out=NH;
run;
*Physical activity variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PAQ_E.XPT";
proc copy in=XP out=NH;
run;
*Cadmium variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PBCD_E.XPT";
proc copy in=XP out=NH;
run;
*Smoking variable for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\SMQ_E.XPT";
proc copy in=XP out=NH;
run;
*Obesity/weight history for 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\WHQ_E.XPT";
proc copy in=XP out=NH;
run;
*C-reactive protein data 2007-2008*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\CRP_E.XPT";
proc copy in=XP out=NH;
run;
*Attempt to append and merge data. Lists the contents of the 2007-2008 files*;
libname NH "C:\NHANES\DATA";
proc contents data =NH.BMX_E varnum;
title "2007-2008 Body Measures";
proc contents data =NH.DBQ_E varnum;
title "2007-2008 Diet Behavior & Nutrition";
proc contents data =NH.MCQ_E varnum;
title "2007-2008 Medical Conditions";
proc contents data =NH.PHYTO_E varnum;
title "2007-2008 Urinary Phytoestrogens";
proc contents data =NH.RHQ_E varnum;
title "2007-2008 Reproductive Health";
proc contents data =NH.ALB_CR_E;
title "2007-2008 Creatine, urine";
proc contents data =NH.COTNAL_E varnum;
title "2007-2008 Cotinine (ng/mL)";
proc contents data =NH.DEMO_E varnum;
title "2007-2008 Demographics (Race/ethnicity, education level, income)";
proc contents data =NH.DIQ_E varnum;
title "2007-2008 Diabetes Status";
proc contents data =NH.KIQ_U_E varnum;
title "2007-2008 Kidney/Renal Failure";
proc contents data =NH.PAQ_E varnum;
title "2007-2008 Physical Activity";
proc contents data =NH.PBCD_E varnum;
title "2007-2008 Cadmium (nm/L)";
proc contents data =NH.SMQ_E varnum;
title "2007-2008 Smoking Status";
proc contents data =NH.WHQ_E varnum;
title "2007-2008 Weight Status";
proc contents data =NH.CRP_E varnum;
title "2007-2008 C-Reactive Protein";
run;
*END 2007-2008*
*IMPORTING THE DATA FOR 2009-2010*;
*body measures variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\BMX_F.XPT";
proc copy in=XP out=NH;
run;
*Soy consumption variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "\C:\NHANES\TEMP\DBQ_F.XPT";
proc copy in=XP out=NH;
run;
*Cancer variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\MCQ_F.XPT";
proc copy in=XP out=NH;
run;
*Urinary phytoestrogens variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PHYTO_F.XPT";
proc copy in=XP out=NH;
run;
*Menopause variable, for which 2009-2010 is not public. See the different link where
you were able to obtain this from Martha Singer.*;
libname NH "C:\NHANES\DATA";
libname XP xport "\C:\NHANES\TEMP\RHQ_F.XPT";
proc copy in=XP out=NH;
run;
*Creatinine variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\ALB_CR_F.XPT";
proc copy in=XP out=NH;
run;
*Cotinine variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\COTNAL_F.XPT";
proc copy in=XP out=NH;
run;
*Demographic variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DEMO_F.XPT";
proc copy in=XP out=NH;
run;
*Diabetes variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\DIQ_F.XPT";
proc copy in=XP out=NH;
run;
*Dietary information variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\KIQ_U_F.XPT";
proc copy in=XP out=NH;
run;
*Physical activity variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PAQ_F.XPT";
proc copy in=XP out=NH;
run;
*Cadmium variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\PBCD_F.XPT";
proc copy in=XP out=NH;
run;
*Smoking variable for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\SMQ_F.XPT";
proc copy in=XP out=NH;
run;
*Obesity/weight history for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\WHQ_F.XPT";
proc copy in=XP out=NH;
run;
*C-reactive protein for 2009-2010*;
libname NH "C:\NHANES\DATA";
libname XP xport "C:\NHANES\TEMP\CRP_F.XPT";
proc copy in=XP out=NH;
run
;
*Attempt to append and merge data. Lists the contents of the 2009-2010 files*;
libname NH "C:\NHANES\DATA";
data DBQ_F;
(rename=(DBQ073U=DBD072U));
proc contents data =NH.BMX_F varnum;
title "2009-2010 Body Measures";
proc contents data =NH.DBQ_F varnum;
title "2009-2010 Diet Behavior & Nutrition";
proc contents data =NH.MCQ_F varnum;
title "2009-2010 Medical Conditions";
proc contents data =NH.PHYTO_F varnum;
title "2009-2010 Urinary Phytoestrogens";
proc contents data =NH.RHQ_F varnum;
title "2009-2010 Reproductive Health";
proc contents data =NH.ALB_CR_F;
title "2009-2010 Creatine, urine";
proc contents data =NH.COTNAL_F varnum;
title "2009-2010 Cotinine (ng/mL)";
proc contents data =NH.DEMO_F varnum;
title "2009-2010 Demographics (Race/ethnicity, education level, income)";
proc contents data =NH.DIQ_F varnum;
title "2009-2010 Diabetes Status";
proc contents data =NH.KIQ_U_F varnum;
title "2009-2010 Kidney/Renal Failure";
proc contents data =NH.PAQ_F varnum;
title "2009-2010 Physical Activity";
proc contents data =NH.PBCD_F varnum;
title "2009-2010 Cadmium (nm/L)";
proc contents data =NH.SMQ_F varnum;
title "2009-2010 Smoking Status";
proc contents data =NH.WHQ_F varnum;
title "2009-2010 Weight Status";
proc contents data =NH.CRP_F varnum;
title "2009-2010 C-Reactive Protein";
run;
*END 2009-2010*
*Attempt to fix DBQ data means*;
data DBQ_6yr;
libname NH "C:\NHANES\DATA";
title "Means for Diet Behaviour & Nutrition 2005-2010";
set nh.DBQ_D (KEEP=SEQN DBD222U)
NH.DBQ_E (keep=SEQN DBD222U DBQ915)
NH.DBQ_F (keep=SEQN DBQ223U DBQ915 DBQ223E in=a);
if a then DBD222U = DBQ223U;
drop DBQ223U;
proc sort data=DBQ_6yr;
by SEQN;
proc contents;
proc means N;
run;
*^^^So, running just this section, it is telling you that for Diet Behaviour, and Nutrition,
there are a total of 31,034 observations for 2005-2006, 2007-2008, 2009-2010 combined.
Of those 31,034 observations for DBD223U You drink another type of milk, there were an
average of 1,093 observations. For DBQ223E You drink soy milk, there was an average of
235 observations, and for DBQ915 Self-perceived vegetarian, there were 19,775
observations*;
title "Means for Body Measures 2005-2010";
libname NH "C:\NHANES\DATA";
data BMX_6yr;
set
NH.BMX_D (keep=SEQN BMXWAIST BMXBMI) NH.BMX_E (keep= SEQN BMXWAIST BMXBMI)
NH.BMX_F (keep= SEQN BMXWAIST BMXBMI);
proc sort data=BMX_6yr;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Medical Conditions 2005-2010";
libname NH "C:\NHANES\DATA";
data MCQ_6yr;
set NH.MCQ_D (keep=SEQN MCQ220 MCQ230A MCQ230B MCQ230C MCQ230D MCQ240C
MCQ240CC MCQ240E MCQ240F MCQ240S)
NH.MCQ_E (keep=SEQN MCQ220 MCQ230A MCQ230B MCQ230C MCQ230D MCQ240C
MCQ240CC MCQ240E MCQ240F MCQ240S)
NH.MCQ_F (keep=SEQN MCQ220 MCQ230A MCQ230B MCQ230C MCQ230D MCQ240C
MCQ240CC MCQ240E MCQ240F MCQ240S);
proc sort data=MCQ_6yr;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Urinary Phytoestrogens 2005-2010";
libname NH "C:\NHANES\DATA";
data PHYTO_6yr;
set NH.PHYTO_D (keep=SEQN URXDAZ URXEQU URXGNS)
NH.PHYTO_E (keep=SEQN URXDAZ URXEQU URXGNS)
NH.PHYTO_F (keep=SEQN URXDAZ URXEQU URXGNS);
proc sort data=PHYTO_6yr;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Reproductive Health 2005-2010";
libname NH "C:\NHANES\DATA";
data RHQ_6yr;
set NH.RHQ_D (keep=SEQN RHQ031 RHD042 RHQ070)
NH.RHQ_E (keep=SEQN RHQ031 RHD042 RHQ070)
NH.RHQ_F (keep=SEQN RHQ031 RHD042 RHQ070);
proc sort data=RHQ_6yr;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Creatinine, urine (ng/mL) 2005-2010";
libname NH "C:\NHANES\DATA";
data ALB_CR_6yr;
set NH.ALB_CR_D (keep= SEQN URXUCR) NH.ALB_CR_E (keep= SEQN URXUCR) NH.ALB_CR_F (keep=
SEQN URXUCR);
proc sort data=ALB_CR_6yr;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Cotinine, urine (ng/mL) 2007-2010";
libname NH "C:\NHANES\DATA";
data COTNAL_4yr;
set NH.COTNAL_E (keep=SEQN LBXCOT) NH.COTNAL_F (keep=SEQN LBXCOT);
proc sort data=COTNAL_4yr;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Cotinine, urine (ng/mL) 2005-2006";
libname NH "C:\NHANES\DATA";
data COT_D;
set NH.COT_D (keep= SEQN LBXCOT);
proc sort data=COT_D;
by SEQN;
proc contents;
proc means N;
run;
title "Means for Demographics (Race/ethnicity, education level, income) 2005-2010";
libname NH "C:\NHANES\DATA";
data DEMO_6yr;
set NH.DEMO_D (keep= SEQN DMDEDUC2 RIDRETH1 INDFMINC)
NH.DEMO_E (keep= SEQN DMDEDUC2 RIDRETH1 INDFMIN2)
NH.DEMO_F (keep=SEQN DMDEDUC2 RIDRETH1 INDHHIN2);
proc sort data=DEMO_6yr;
by SEQN;
proc contents;
proc means N;
run;
*Attempt to fix DIQ data means*;
libname NH "C:\NHANES\DATA";
data DIQ_6yr;
title "Means for Diabetes 2005-2010";
set NH.DIQ_D (keep=SEQN DID070)
NH.DIQ_E (keep=SEQN DID070)
NH.DIQ_F (keep=SEQN DIQ070 in=a);
if a then DID070 = DIQ070;
drop DIQ070;
proc sort data=DIQ_6yr;
by SEQN;
proc contents;
proc means N;
run;
*^^^So, this is saying that for Take diabetic pills to lower blood sugar (DID070),
the average was 3,073 observations between 2005-2010.*;
title "Means for Food Frequency Questionnaire 2005-2006";
libname NH "C:\NHANES\DATA";
data FFQRAW_D;
set NH.FFQRAW_D (keep=SEQN FFQ0007A FFQ0096 FFQ0130A);
proc sort data=FFQRAW_D;
by SEQN;
proc contents;
proc means N;
title "Means for Food Frequency Questionnaire 2005-2006";
run;
*For the question "What kind of milk did you drink? (FFQ0007A)" 91 people answered soy milk.
Soy milk is code or value 5.
For the question "Did you eat tofu? (FFQ0096)" I want 5-6 times per week, 1 time per day,
and 2 or more times per day, which correspond with code or value numbers 9, 10, and 11.
9 is for 5-6 times per week, 10 is for 1 time per day, and 11 is for 2 or more times per
day*;
proc means data=FFQRAW_D N Nmiss min max;
where FFQ0007A=5;
var FFQ0007A FFQ0096 FFQ0130A;
title "Check missing, min, max values for numeric variables";
proc freq data=FFQRAW_D;
where FFQ0007A=5;
table FFQ0007A FFQ0096 FFQ0130A/list missing;
title 'Check frequency distribution for categorical variables';
run;
title "Means for Kidney/Renal Failure 2005-2010";
libname NH "C:\NHANES\DATA";
data KIQ_U_6yr;
set NH.KIQ_U_D (keep= SEQN KIQ022) NH.KIQ_U_E (keep=SEQN KIQ022)
NH.KIQ_U_F (keep=SEQN KIQ022);
proc sort data=KIQ_U_6yr;
by SEQN;
proc contents;
proc means N;
title "Means for Kidney/Renal Failure 2005-2010";
run;
title "Means for Physical Activity 2005-2010";
libname NH "C:\NHANES\DATA";
data PAQ_6yr;
set NH.PAQ_D (keep= SEQN PAD590) NH.PAQ_E (keep=SEQN PAD590)
NH.PAQ_F (keep=SEQN PAD590);
proc sort data=PAQ_6yr;
by SEQN;
proc contents;
proc means N;
title "Means for Physical Activity 2005-2010";
run;
title "Means for Blood Cadmium (nm/L) 2005-2010";
libname NH "C:\NHANES\DATA";
data PBCD_6yr;
set NH.PBCD_D (keep= SEQN LBDBCDSI) NH.PBCD_E (keep=SEQN LBDBCDSI)
NH.PBCD_F (keep=SEQN LBDBCDSI);
proc sort data=PBCD_6yr;
by SEQN;
proc contents;
proc means N;
title "Means for Blood Cadmium (nm/L) 2005-2010";
run;
title "Means for Smoking Status 2005-2010";
libname NH "C:\NHANES\DATA";
data SMQ_6yr;
set NH.SMQ_D (keep= SEQN SMQ040) NH.SMQ_E (keep=SEQN SMQ040)
NH.SMQ_F (keep=SEQN SMQ040);
proc sort data=SMQ_6yr;
by SEQN;
proc contents;
proc means N;
title "Means for Smoking Status 2005-2010";
run;
title "Means for Weight Status 2005-2010";
libname NH "C:\NHANES\DATA";
data WHQ_6yr;
set NH.WHQ_D (keep= SEQN WHQ030) NH.WHQ_E (keep=SEQN WHQ030)
NH.WHQ_F (keep=SEQN WHQ030);
proc sort data=WHQ_6yr;
by SEQN;
proc contents;
proc means N;
title "Means for Weight Status 2005-2010";
run;
title "Means for C-Reactive Protein 2005-2010";
libname NH "C:\NHANES\DATA";
data CRP_6yr;
set NH.CRP_D (keep= SEQN LBXCRP) NH.CRP_E (keep=SEQN LBXCRP)
NH.CRP_F (keep=SEQN LBXCRP);
proc sort data=CRP_6yr;
by SEQN;
proc contents;
proc means N;
title "Means for C-Reactive Protein 2005-2010";
run;
*HERE USE PROC UNIVARIATE TO OBSERVE THE DISTRIBUTIONS.
Use proc freq to look at the frequency in categorical variables.
Looking at creatinine first, which is a continuous variable...*;
title "Means for Creatinine, urine (ng/mL) 2005-2010";
libname NH "C:\NHANES\DATA";
data ALB_CR_6yr;
set NH.ALB_CR_D (keep= SEQN URXUCR) NH.ALB_CR_E (keep= SEQN URXUCR) NH.ALB_CR_F (keep=
SEQN URXUCR);
proc sort data=ALB_CR_6yr;
by SEQN;
proc univariate;
Histogram / normal;
run;
*A normal creatinine level is 0.6 to 1.2mg/dL per deciliter in adult males and
0.5 to 1.1 mg/dL in adult females.*;
*Looking at distributions for BMI and waist circumference.
These are continuous variables.
Normal BMI is between 18.5 and 24.9
Normal waist circumference is less than 40 inches/101.6cm around for men,
less than 35 inches/88.9cm for women.*;
title "Means for BMI (kg/m**2) 2005-2010";
libname NH "C:\NHANES\DATA";
data BMX_6yr;
set NH.BMX_D (keep= SEQN BMXBMI BMXWAIST) NH.BMX_E (keep= SEQN BMXBMI BMXWAIST)
NH.BMX_F (keep= SEQN BMXBMI BMXWAIST);
proc sort data=BMX_6yr;
by SEQN;
proc univariate;
Histogram / normal;
run;
*So, the average is 88.6cm (which is good)*;
*Looking at the distribution for cotinine, a biomarker for exposure to tobacco smoke.
Continuous variable.
Cotinine levels in a nonsmoker are generally less than 10 ng/mL. Cotinine levels in
a light smoker or someone exposed to secondhand smoke are 11 ng/mL to 30 ng/mL.
Cotinine levels in a heavy smoker may be more than 500 ng/mL.*;
title "Means for Cotinine (ng/mL) 2005-2006";
libname NH "C:\NHANES\DATA";
data COT_D;
set NH.COT_D (keep= SEQN LBXCOT);
proc sort data=COT_D;
by SEQN;
proc univariate;
Histogram / normal;
run;
*Looking at the distribution for cotinine in 2007-2008*;
title "Means for Cotinine (ng/mL) 2007-2010";
libname NH "C:\NHANES\DATA";
data COTNAL_4yr;
set NH.COTNAL_E (keep=SEQN LBXCOT) NH.COTNAL_F (keep=SEQN LBXCOT);
proc sort data=COTNAL_4yr;
by SEQN;
proc univariate;
Histogram / normal;
run;
*Means for C-Reactive Protein, 2005-2010*;
libname NH "C:\NHANES\DATA";
data CRP_D;
set NH.CRP_D (keep= SEQN LBXCRP);
proc sort data=CRP_D;
by SEQN;
proc univariate;
Histogram / normal;
run;
*Looking at the distribution for C-reactive Protein in 2005-2010*;
title "Means for C-Reactive Protein (mg/dL) 2007-2010";
libname NH "C:\NHANES\DATA";
data CRP_6yr;
set NH.CRP_D (keep=SEQN LBXCRP) NH.CRP_E (keep=SEQN LBXCRP) NH.CRP_F (keep=SEQN LBXCRP);
proc sort data=CRP_6yr;
by SEQN;
proc univariate;
Histogram / normal;
run;
*The analysis I want to run is to see whether there's a correlation between urinary phytoestrogens and cancers.
So here I am looking at Ever told you had cancer or malignancy while accounting for urinary phytoestrogens (URXDAZ).
Urinary phytoestrogens are a strong determinant. So if I'm looking at this correlation and I have a wide urinary
phytoestrogen range, I want to account for urinary phytoestrogens. Look at the correlation between
Ever told you had cancer or malignancy and BMI while accounting for the urinary phytoestrogen.*;
Proc corr data=NH pearson spearman;
libname NH "C:\NHANES\DATA";
var BMXBMI;
with MCQ220;
partial URXDAZ;
title "Correlation Between Ever told you had cancer or malignancy & Urinary Phytoestrogens";
run;
*HOW DO I WRITE THIS CODE?!?!!?*;
Proc corr data=BMX_D data=MCQ_D data=PHYTO_D spearman;
libname NH "C:\NHANES\TEMP";
var BMXBMI;
with MCQ220;
partial URXDAZ;
title "Correlation Between Ever told you had cancer or malignancy & Urinary Phytoestrogens";
run;
*???!?!?*;
romaincrisis
turnnoffyourmind
areudirty
dailybridgerton
andulkaphoto
vaporwave
titsay