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)

pdf("fig61.pdf",width=7.5,height=2.8)

par(mfrow=c(1,3), mar=c(4,4,1,1))

plot(data_stroke$age_ctrl, data_stroke$systolic_ctrl, 
  type="p",pch=".",xlab="Age",ylab="Systolic Blood Pressures",
  mgp=c(2,1,0),xlim=range(data_stroke$age_ctrl), 
  ylim=range(data_stroke$systolic_ctrl), cex=0.2)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,1],lty=1,cex=6,lwd=1.5)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,1]+
        qnorm(0.975)*sqrt(Estimated_Regular_Pattern$var_est[,1,1]),
        lty=2,cex=6,lwd=1.5)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,1]-
        qnorm(0.975)*sqrt(Estimated_Regular_Pattern$var_est[,1,1]),
        lty=2,cex=6,lwd=1.5)

plot(data_stroke$age_ctrl, data_stroke$diastolic_ctrl, 
  type="p",pch=".",xlab="Age",ylab="Diastolic Blood Pressures",
  mgp=c(2,1,0),xlim=range(data_stroke$age_ctrl), 
  ylim=c(30,140), cex=0.2)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,2],lty=1,cex=6,lwd=1.5)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,2]+
        qnorm(0.975)*sqrt(Estimated_Regular_Pattern$var_est[,2,2]),
        lty=2,cex=6,lwd=1.5)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,2]-
        qnorm(0.975)*sqrt(Estimated_Regular_Pattern$var_est[,2,2]),
        lty=2,cex=6,lwd=1.5)

plot(data_stroke$age_ctrl, data_stroke$cholesterol_ctrl, 
  type="p",pch=".",xlab="Age",ylab="Cholesterol Levels",
  mgp=c(2,1,0),xlim=range(data_stroke$age_ctrl), 
  ylim=c(100,400), cex=0.2)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,3],lty=1,cex=6,lwd=1.5)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,3]+
        qnorm(0.975)*sqrt(Estimated_Regular_Pattern$var_est[,3,3]),
        lty=2,cex=6,lwd=1.5)
  lines(Estimated_Regular_Pattern$grid,
        Estimated_Regular_Pattern$mean_est[,3]-
        qnorm(0.975)*sqrt(Estimated_Regular_Pattern$var_est[,3,3]),
        lty=2,cex=6,lwd=1.5)

graphics.off()