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

## ----eval = FALSE-------------------------------------------------------------
# devtools::install_github("chaochengstat/mrPStrata")

## ----import, message=FALSE, warning=FALSE-------------------------------------
library("mrPStrata")

## -----------------------------------------------------------------------------
attach(sim_data)
head(sim_data)

## -----------------------------------------------------------------------------
res = mrPStrata(times=c(1,2,3,4,5,6,7,8),
                data = sim_data,
                Xpi_names = c("X1","X2","X3","X4","X5"),
                Xe_names = c("X1","X2","X3","X4","X5"),
                Xc_names = c("X1","X2","X3","X4","X5"),
                Xt_names = c("X1","X2","X3","X4","X5"),
                Z_name = "z",
                S_name = "s",
                U_name ="U",
                delta_name = "delta",
                B=20)

## -----------------------------------------------------------------------------
print(res$Compliers)

## ----out.width="90%",dpi=300,fig.width=10, fig.height=7-----------------------
plot.psce(res)

## -----------------------------------------------------------------------------
res = mrPStrata_PI_SA(times=c(1,2,3,4,5,6,7,8),
                      data = sim_data,
                      Xpi_names = c("X1","X2","X3","X4","X5"),
                      Xe_names = c("X1","X2","X3","X4","X5"),
                      Xc_names = c("X1","X2","X3","X4","X5"),
                      Xt_names = c("X1","X2","X3","X4","X5"),
                      Z_name = "z",
                      S_name = "s",
                      U_name ="U",
                      delta_name = "delta",
                      xi0 = 0.05,
                      xi1 = 0.05,
                      eta0=1,
                      eta1=1,
                      B=20)

## ----out.width="90%",dpi=300,fig.width=10, fig.height=7-----------------------
plot.psce(res)

## -----------------------------------------------------------------------------
res = mrPStrata_MO_SA(times=c(1,2,3,4,5,6,7,8),
                      data = sim_data,
                      Xpi_names = c("X1","X2","X3","X4","X5"),
                      Xe_names = c("X1","X2","X3","X4","X5"),
                      Xc_names = c("X1","X2","X3","X4","X5"),
                      Xt_names = c("X1","X2","X3","X4","X5"),
                      Z_name = "z",
                      S_name = "s",
                      U_name ="U",
                      delta_name = "delta",
                      zeta=0.02,
                      B=20)

## ----out.width="90%",dpi=300,fig.width=12, fig.height=7-----------------------
plot.psce(res)