MEPS data

Description

2008 MEPS data.

Usage

data(meps)

Format

meps is a 18592 row data frame with the following columns

bmi

body mass index.

age

age in years.

gender

equal to 1 if male.

race

levels: 2 white, 3 black, 4 native American, 5 others.

education

years of education.

health

levels: 5 excellent, 6 very good, 7 good, 8 fair, 9 poor.

limitation

equal to 1 if health limits physical activity.

region

levels: 2 northeast, 3 mid-west, 4 south, 5 west.

private

equal to 1 if individual has private health insurance.

visits.hosp

equal to 1 if at least one visit to hospital outpatient departments.

diabetes

equal to 1 if diabetic.

hypertension

equal to 1 if hypertensive.

hyperlipidemia

equal to 1 if hyperlipidemic.

income

income (000's).

Source

The data have been obtained from http://www.meps.ahrq.gov/.

Examples

## Not run:  

###################################################
###################################################

library("SemiParBIVProbit")
data("meps", package = "SemiParBIVProbit") 

###################################################
# Bivariate brobit models with endogenous treatment
###################################################

treat.eq <- private ~ s(bmi) + s(income) + s(age) + s(education) +
                      as.factor(health) + as.factor(race) +
                      as.factor(limitation) + as.factor(region) + 
                      gender  + hypertension + hyperlipidemia + diabetes
out.eq <- visits.hosp ~ private + s(bmi) + s(income) + s(age) + 
                        s(education) + as.factor(health) + 
                        as.factor(race) + as.factor(limitation) + 
                        as.factor(region) + gender + hypertension + 
                        hyperlipidemia + diabetes

f.list <- list(treat.eq, out.eq) 
bpN    <- SemiParBIVProbit(f.list, data = meps)
bpF    <- SemiParBIVProbit(f.list, data = meps, BivD = "F")
bpC0   <- SemiParBIVProbit(f.list, data = meps, BivD = "C0")
bpC180 <- SemiParBIVProbit(f.list, data = meps, BivD = "C180")
bpJ0   <- SemiParBIVProbit(f.list, data = meps, BivD = "J0")
bpJ180 <- SemiParBIVProbit(f.list, data = meps, BivD = "J180")
bpG0   <- SemiParBIVProbit(f.list, data = meps, BivD = "G0")
bpG180 <- SemiParBIVProbit(f.list, data = meps, BivD = "G180")

conv.check(bpJ0)

AIC(bpN, bpF, bpC0, bpC180, bpJ0, bpJ180, bpG0, bpG180) 

set.seed(1)
summary(bpJ0, cm.plot = TRUE, cex.axis = 1.6, 
        cex.lab = 1.6, cex.main = 1.7)

#dev.copy(postscript, "contplot.eps")
#dev.off()

par(mfrow = c(2, 2), mar = c(4.5, 4.5, 2, 2), 
    cex.axis = 1.6, cex.lab = 1.6)
plot(bpJ0, eq = 1, seWithMean = TRUE, scale = 0, shade = TRUE, 
     pages = 1, jit = TRUE)

#dev.copy(postscript, "spline1.eps")
#dev.off() 

par(mfrow = c(2, 2), mar = c(4.5, 4.5, 2, 2), 
    cex.axis = 1.6, cex.lab = 1.6)
plot(bpJ0, eq = 2, seWithMean = TRUE, scale = 0, shade = TRUE, 
     pages = 1, jit = TRUE)

#dev.copy(postscript, "spline2.eps")
#dev.off() 

set.seed(1)
AT(bpJ0, nm.end = "private", hd.plot = TRUE, cex.axis = 1.5, 
   cex.lab = 1.5, cex.main = 1.6)

#dev.copy(postscript, "hd.plotAT.eps")
#dev.off()

AT(bpJ0, nm.end = "private", type = "univariate")

AT(bpJ0, nm.end = "private", type = "naive")


## End(Not run)

#

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(SemiParBIVProbit)
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-12. For overview type 'help("mgcv-package")'.

This is SemiParBIVProbit 3.7-1.
For overview type 'help("SemiParBIVProbit-package")'.

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/SemiParBIVProbit/meps.Rd_%03d_medium.png", width=480, height=480)
> ### Name: meps
> ### Title: MEPS data
> ### Aliases: meps
> 
> ### ** Examples
> 
> 
> ## Not run: 
> ##D  
> ##D 
> ##D ###################################################
> ##D ###################################################
> ##D 
> ##D library("SemiParBIVProbit")
> ##D data("meps", package = "SemiParBIVProbit") 
> ##D 
> ##D ###################################################
> ##D # Bivariate brobit models with endogenous treatment
> ##D ###################################################
> ##D 
> ##D treat.eq <- private ~ s(bmi) + s(income) + s(age) + s(education) +
> ##D                       as.factor(health) + as.factor(race) +
> ##D                       as.factor(limitation) + as.factor(region) + 
> ##D                       gender  + hypertension + hyperlipidemia + diabetes
> ##D out.eq <- visits.hosp ~ private + s(bmi) + s(income) + s(age) + 
> ##D                         s(education) + as.factor(health) + 
> ##D                         as.factor(race) + as.factor(limitation) + 
> ##D                         as.factor(region) + gender + hypertension + 
> ##D                         hyperlipidemia + diabetes
> ##D 
> ##D f.list <- list(treat.eq, out.eq) 
> ##D bpN    <- SemiParBIVProbit(f.list, data = meps)
> ##D bpF    <- SemiParBIVProbit(f.list, data = meps, BivD = "F")
> ##D bpC0   <- SemiParBIVProbit(f.list, data = meps, BivD = "C0")
> ##D bpC180 <- SemiParBIVProbit(f.list, data = meps, BivD = "C180")
> ##D bpJ0   <- SemiParBIVProbit(f.list, data = meps, BivD = "J0")
> ##D bpJ180 <- SemiParBIVProbit(f.list, data = meps, BivD = "J180")
> ##D bpG0   <- SemiParBIVProbit(f.list, data = meps, BivD = "G0")
> ##D bpG180 <- SemiParBIVProbit(f.list, data = meps, BivD = "G180")
> ##D 
> ##D conv.check(bpJ0)
> ##D 
> ##D AIC(bpN, bpF, bpC0, bpC180, bpJ0, bpJ180, bpG0, bpG180) 
> ##D 
> ##D set.seed(1)
> ##D summary(bpJ0, cm.plot = TRUE, cex.axis = 1.6, 
> ##D         cex.lab = 1.6, cex.main = 1.7)
> ##D 
> ##D #dev.copy(postscript, "contplot.eps")
> ##D #dev.off()
> ##D 
> ##D par(mfrow = c(2, 2), mar = c(4.5, 4.5, 2, 2), 
> ##D     cex.axis = 1.6, cex.lab = 1.6)
> ##D plot(bpJ0, eq = 1, seWithMean = TRUE, scale = 0, shade = TRUE, 
> ##D      pages = 1, jit = TRUE)
> ##D 
> ##D #dev.copy(postscript, "spline1.eps")
> ##D #dev.off() 
> ##D 
> ##D par(mfrow = c(2, 2), mar = c(4.5, 4.5, 2, 2), 
> ##D     cex.axis = 1.6, cex.lab = 1.6)
> ##D plot(bpJ0, eq = 2, seWithMean = TRUE, scale = 0, shade = TRUE, 
> ##D      pages = 1, jit = TRUE)
> ##D 
> ##D #dev.copy(postscript, "spline2.eps")
> ##D #dev.off() 
> ##D 
> ##D set.seed(1)
> ##D AT(bpJ0, nm.end = "private", hd.plot = TRUE, cex.axis = 1.5, 
> ##D    cex.lab = 1.5, cex.main = 1.6)
> ##D 
> ##D #dev.copy(postscript, "hd.plotAT.eps")
> ##D #dev.off()
> ##D 
> ##D AT(bpJ0, nm.end = "private", type = "univariate")
> ##D 
> ##D AT(bpJ0, nm.end = "private", type = "naive")
> ##D 
> ## End(Not run)
> 
> #
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>