manuscript/justify.R

In mac-theobio/effects: isolated confidence intervals and bias-corrected predictor effects

library(shellpipes)
library(vareffects); varefftheme()
library(dplyr)

loadEnvironments()
startGraphics()
 
set.seed(9991)

############################################################################################
# Global variables
nHH_obs <- 50
perHH <- 1

############################################################################################
# Figure1 : Simulate data for justifying proposed approach
############################################################################################

#### Simulation

## No interaction
justify_sim <- linearsim(nHH=nHH_obs
	, perHH=perHH
	, form=~1+x1+x2+x3
	, betas=c(5, -3, 1, 2)
	, pgausian=list(p=3)
	, pcat=list(p=0)
	, link_scale=TRUE
)
justify_sim_df <- justify_sim$data
justify_sim_df$x2 <- justify_sim_df$x2 - min(justify_sim_df$x2)
justify_sim_betas <- justify_sim$betas
justify_sim_betas
head(justify_sim_df)

true_prop_df <- (justify_sim_df
	%>% summarize_all(mean)
)
true_prop_df

## With interaction
justify_inter_sim <- linearsim(nHH=nHH_obs
	, perHH=perHH
	, form=~1+x1+x2*x3
	, betas=c(5, -3, 1, 2, 5)
	, pgausian=list(p=3)
	, pcat=list(p=0)
	, link_scale=TRUE
)
justify_inter_sim_df <- justify_inter_sim$data 
justify_inter_sim_betas <- justify_inter_sim$betas
head(justify_inter_sim_df)
true_prop_inter_df <- (justify_inter_sim_df
	%>% summarize_all(mean)
)
true_prop_inter_df

saveEnvironment()