# NOT RUN {
my_population <- declare_population(N = 500, noise = rnorm(N))
my_potential_outcomes <- declare_potential_outcomes(
Y_Z_0 = noise, Y_Z_1 = noise +
rnorm(N, mean = 2, sd = 2))
my_assignment <- declare_assignment()
my_estimand <- declare_estimand(ATE = mean(Y_Z_1 - Y_Z_0))
my_estimator <- declare_estimator(Y ~ Z, estimand = my_estimand)
my_reveal <- declare_reveal()
design <- my_population +
my_potential_outcomes +
my_estimand +
my_assignment +
my_reveal +
my_estimator
# }
# NOT RUN {
# using built-in defaults:
diagnosis <- diagnose_design(design)
diagnosis
# }
# NOT RUN {
# using a user-defined diagnosand
my_diagnosand <- declare_diagnosands(absolute_error = mean(abs(estimate - estimand)))
# }
# NOT RUN {
diagnosis <- diagnose_design(design, diagnosands = my_diagnosand)
diagnosis
get_diagnosands(diagnosis)
get_simulations(diagnosis)
reshape_diagnosis(diagnosis)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab