if (FALSE) {
set.seed(0)
dat_phaseI <- simulate_data_RoMFCC(p_cellwise = 0.05,
p_casewise = 0.05,
outlier = "outlier_E",
M_outlier_cell = 0.03,
M_outlier_case = 0.01,
max_n_cellwise = 10)
dat_phaseII <- simulate_data_RoMFCC(OC = "OC_E",
M_OC = 0.01,
which_OC = 5)
mfdobj_phaseI <- get_mfd_list(dat_phaseI$X_mat_list, n_basis = 5)
mfdobj_phaseII <- get_mfd_list(dat_phaseII$X_mat_list, n_basis = 5)
mfdobj_training_phaseI <- mfdobj_phaseI[1:333, ]
mfdobj_tuning_phaseI <- mfdobj_phaseI[334:1000, ]
out_phaseI <- RoAMFEWMA_PhaseI(mfdobj = mfdobj_training_phaseI,
mfdobj_tuning = mfdobj_tuning_phaseI)
out_phaseII <- RoAMFEWMA_PhaseII(mfdobj_2 = mfdobj_phaseII,
mod_1 = out_phaseI)
plot_control_charts(out_phaseII$cc)
}
Run the code above in your browser using DataLab