Script to test simex functionality using parallel processing and sgpData_LONG

simex_simulation.R

##############################################################################################
###
### Test script of Simex functionality
###
##############################################################################################

### Load SGP package

require(SGP)
options(error=recover)


### prepareSGP

Demonstration_SGP_Simex <- prepareSGP(sgpData_LONG)


### analyzeSGP

Demonstration_SGP_Simex <- analyzeSGP(Demonstration_SGP_Simex,
                                  sgp.percentiles=TRUE,
                                  sgp.projections = FALSE,
                                  sgp.projections.lagged = FALSE,
                                  sgp.percentiles.baseline = FALSE,
                                  sgp.projections.baseline = FALSE,
                                  sgp.projections.lagged.baseline = FALSE,
                                  simulate.sgps=FALSE,
                                  calculate.simex=list(state="DEMO",lambda=seq(0,2,0.5), simulation.iterations=25, extrapolation="linear"),
                                  goodness.of.fit.print=FALSE,
                                  parallel.config=list(BACKEND="PARALLEL", WORKERS=list(PERCENTILES=20)))


### combineSGP

Demonstration_SGP_Simex <- combineSGP(Demonstration_SGP_Simex)


### Save results

save(Demonstration_SGP_Simex, file="Demonstration_SGP_Simex.Rdata")