Task 3: How to Check Distributions and Describe the Impact of Influential Outliers in NHANES II Analyses

In this task, you will check for outliers and their potential impact using the following steps:

 

Step 1 Check distributions by running a univariate analysis

Before you analyze your data, it is very important that you check the distribution and normality of the data and identify outliers for continuous variables.

 

Program to Plot Distribution of Continuous Variables

Statements Explanation
proc univariate data =demo3_nh 2 normal plot;

Use the proc univariate procedure to get all default descriptive statistics, such as mean, minimum and maximum values, standard deviation, and skewness, etc. Use the normal plot option to obtain a plot of normality.

where n2ah0047>=20 ;

Use the where statement to select the participants who were age 20 years and older.

id seqn;

Use the id statement to list the sequence numbers associated with extreme values in the output.

var n2lb0421;

Use the var statement to list the variables of interest.

Highlighted items from the univariate analysis output:

 

Step 2 Plot graph of survey weight against the distribution of the variable

In this example, you will plot the examined sample survey weight (n2ah0282) against the distribution of the cholesterol variable to determine whether the extreme observations are influential outliers.

 

Plot Exam Weight Against Cholesterol

Statements Explanation
symbol1 value =dot height = .2;

Use the option statements, symbol and height, to format the output of the plot.

proc gplot data =demo3_nh2;
where n2ah0047>=20 ;

      plot n2ah0282*n2lb0421/ frame ;

      title 'NHANES II, adults age 20 years and older' ;

run ;

Use the proc gplot procedure to plot the total serum cholesterol (n2lb0421) by the corresponding sample weight for each observation in the dataset.  Use the where statement to select the participants who were age 20 years and older.

 

Highlighted items from plotting the survey weight against the distribution of the cholesterol variable:

 

Step 3 Identify outliers and compare estimates with outliers deleted against the original estimates with outliers included

In this step you will:

Program to Create Dataset Without Outliers and Output Means of Both Datasets

Statements Explanation
data temp4_nh2;

set demo3_nh2;

Use the data and set statements to refer to your analytic dataset.

if seqn in (2474, 18622) then delete;

Use the if, then statements to delete the outliers using their SEQN previously identified in the plot of survey weight versus distribution of the variable.  The SEQNs associated with these outliers are listed in the proc univariate output under extreme observations.

proc means data =temp4_nh2 mean stderr maxdec = 1 ;

Use the proc means procedure to determine the mean and standard error for the dataset without the outliers.

where n2ah0047>=20 ;

Use the where statement to select the participants who were age 20 years and older.

var n2lb0421;

Use the var statement to indicate the variable of interest.

class n2ah0056;

Use the class statement to group the variable of interest by race categories.

weight n2ah0282;

Use the weight statement to account for the unequal probability of sampling and non-response. In this example, the examined sample weight is used.

proc means data=demo3_nh2 mean stderr maxdec=1;

Use the proc means procedure to determine the mean and standard error for the dataset with the outliers.

where n2ah0047>=20 ;

Use the where statement to select the participants who were age 20 years and older.

var n2lb0421;

Use the var statement to indicate the variable of interest.

class n2ah0056;

Use the class statement to group the variable of interest by race categories.

weight n2ah0282;

Use the weight statement to account for the unequal probability of sampling and non-response. In this example, the examined sample weight for all six years of data is used.

 

Highlighted items from comparison of the results with and without outliers:

 

close window icon Close Window