Unique contribution of SSS to each model (semi-partial correlations)

load(file="~/Desktop/Research/HANDLS Data/MeasurementData1.rdata")
library(ppcor)

## Loading required package: MASS

Make Variables Numeric

##Make Sex Numeric
MeasurementData1$SexNum <- ifelse(MeasurementData1$Sex == "Men",1,0)

##Make Race Numeric
MeasurementData1$RaceNum<- ifelse(MeasurementData1$Race == "White",1,0)

Unique Contribution of self-rated health to the Model

spcor.test(MeasurementData1$sHealthNum,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","RaceNum","Educ","Employment01","acasiIncomx01","Neighborhood02","CES")])

##        estimate         p.value   statistic    n gp  Method
## 1 0.07933045288 0.0003028704894 3.618979278 2077  7 pearson

Unique Contribution of Race to the Model

spcor.test(MeasurementData1$RaceNum,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","sHealthNum","Educ","Employment01","acasiIncomx01","Neighborhood02","CES")])

##         estimate         p.value    statistic    n gp  Method
## 1 -0.08103773942 0.0002235156953 -3.697373455 2077  7 pearson

Unique Contribution of Education to the Model

spcor.test(MeasurementData1$Educ,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","RaceNum","sHealthNum","Employment01","acasiIncomx01","Neighborhood02","CES")])

##        estimate         p.value   statistic    n gp  Method
## 1 0.09445947212 1.672285848e-05 4.314862823 2077  7 pearson

Unique Contribution of Employment to the Model

spcor.test(MeasurementData1$Employment01,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","RaceNum","sHealthNum","Educ","acasiIncomx01","Neighborhood02","CES")])

##         estimate      p.value     statistic    n gp  Method
## 1 -0.02035707631 0.3545876032 -0.9259353797 2077  7 pearson

Unique Contribution of Income to the Model

spcor.test(MeasurementData1$acasiIncomx01,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","RaceNum","sHealthNum","Educ","Employment01","Neighborhood02","CES")])

##         estimate        p.value    statistic    n gp  Method
## 1 -0.06001803008 0.006305275628 -2.734265081 2077  7 pearson

Unique Contribution of Neighborhood to the Model

spcor.test(MeasurementData1$Neighborhood02,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","RaceNum","sHealthNum","Educ","Employment01","acasiIncomx01","CES")])

##         estimate        p.value    statistic    n gp  Method
## 1 -0.06965289726 0.001519397879 -3.175195807 2077  7 pearson

Unique Contribution of CES to the Model

spcor.test(MeasurementData1$CES,MeasurementData1$MCLcountry,MeasurementData1[,c("SexNum","RaceNum","sHealthNum","Educ","Employment01","acasiIncomx01","Neighborhood02")])

##        estimate         p.value    statistic    n gp  Method
## 1 -0.1806406928 1.212151649e-16 -8.352082734 2077  7 pearson

Unique Contribution of Sex to the Model

spcor.test(MeasurementData1$SexNum,MeasurementData1$MCLcountry,MeasurementData1[,c("CES","RaceNum","sHealthNum","Educ","Employment01","acasiIncomx01","Neighborhood02")])

##          estimate      p.value    statistic    n gp  Method
## 1 0.0003061665441 0.9888927501 0.0139230063 2077  7 pearson