doc/Data.R

In wal615/prime.total.effect: Total effect estimation of the environmental mixtures

## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = ""
)

## ----hemoglobin.PCB, echo = TRUE----------------------------------------------
library(prime.total.effect)
dim(hemoglobin.PCB)

# there are still many NAs in the data set, to remove the all the NAs, you could use na.omit

PCB.no.missing <- na.omit(hemoglobin.PCB)

# response is LBXGH
head(PCB.no.missing)

# covariates are PCBs
head(PCB.no.missing)

dim(hemoglobin.PCB)

dim(PCB.no.missing)