data(data_stroke)

nrow=length(data_stroke$systolic_ctrl[,1])
ncol=length(data_stroke$systolic_ctrl[1,])

Estimated_Regular_Pattern = estimate_pattern_long_md(
  data_array=array(c(data_stroke$systolic_ctrl,
                     data_stroke$diastolic_ctrl,
                     data_stroke$cholesterol_ctrl),
                   dim=c(nrow,ncol,3)),
  time_matrix=data_stroke$age_ctrl,
  nobs=rep(ncol,nrow),
  time_unit=1,
  estimation_method="meanvar",
  bw_mean=10,
  bw_var=10)

Monitoring_Results_IC = monitor_long_md(
  data_array_new=array(c(data_stroke$systolic_ctrl,
                         data_stroke$diastolic_ctrl,
                         data_stroke$cholesterol_ctrl),
                       dim=c(nrow,ncol,3)),
  time_matrix_new=data_stroke$age_ctrl,
  nobs_new=rep(ncol,nrow),
  pattern=Estimated_Regular_Pattern,
  side="upward",
  method="multivariate EWMA",
  parameter=0.2)

CL=search_CL(
  chart_matrix=Monitoring_Results_IC$chart_matrix,
  time_matrix=data_stroke$age_ctrl,
  nobs=rep(ncol,nrow),
  starttime=data_stroke$age_ctrl[,1],
  endtime=data_stroke$age_ctrl[,ncol],
  time_unit=1,
  ATS_nominal=15)

nrow_case=dim(data_stroke$systolic_case)[1]
ncol_case=dim(data_stroke$systolic_case)[2]

Monitoring_Results_OC=monitor_long_md(
  data_array_new=array(c(data_stroke$systolic_case,
                         data_stroke$diastolic_case,
                         data_stroke$cholesterol_case),
                       dim=c(nrow_case,ncol_case,3)),
  time_matrix_new=data_stroke$age_case,
  nobs_new=rep(ncol_case,nrow_case),
  pattern=Estimated_Regular_Pattern,
  side="upward",
  method="multivariate EWMA",
  parameter=0.2)

Evaluation=evaluate_control_chart_two_groups(
  chart_matrix_IC=Monitoring_Results_IC$chart_matrix,
  time_matrix_IC=data_stroke$age_ctrl,
  nobs_IC=rep(ncol,nrow),
  starttime_IC=data_stroke$age_ctrl[,1],
  endtime_IC=data_stroke$age_ctrl[,ncol],
  chart_matrix_OC=Monitoring_Results_OC$chart_matrix,
  time_matrix_OC=data_stroke$age_case,
  nobs_OC=rep(ncol_case,nrow_case),
  starttime_OC=data_stroke$age_case[,1],
  endtime_OC=data_stroke$age_case[,ncol_case],
  time_unit=1,
  no_signal_action="maxtime")

pdf("fig63.pdf",width=6.5,height=3.3)

par(mfrow=c(1,2), mar=c(4,4,1,1), mgp=c(2,1,0))

plot_evaluation(Evaluation)

graphics.off()