add output parameter posterior

lgnakano · Feb 16, 2020 · 3e101cd · 3e101cd
1 parent 107273b
commit 3e101cd
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/R/qh.eSIR.R b/R/qh.eSIR.R
@@ -427,7 +427,7 @@ qh.eSIR<-function (Y,R, phi=NULL,change_time=NULL,begin_str="01/23/2020",T_fin=2
 
   if(save_files) write.csv(out_table,file=paste0(file_add,casename,"_summary.csv"))
 
-  return(list(casename=casename,incidence_mean=incidence_mean,incidence_ci=incidence_ci,out_table=out_table,forecast_infection=plot,forecast_removed=plot2,first_stat_mean=as.character(first_order_change_date),first_stat_ci=as.character(first_order_change_date_ci),second_stat_mean=as.character(second_order_change_date),second_stat_ci=as.character(second_order_change_date_ci),dic_val=dic_val,plot_list=plot_list,plot2_list=plot2_list))
+  return(list(casename=casename,incidence_mean=incidence_mean,incidence_ci=incidence_ci,out_table=out_table,forecast_infection=plot,forecast_removed=plot2,first_stat_mean=as.character(first_order_change_date),first_stat_ci=as.character(first_order_change_date_ci),second_stat_mean=as.character(second_order_change_date),second_stat_ci=as.character(second_order_change_date_ci),dic_val=dic_val,plot_list=plot_list,plot2_list=plot2_list,theta_p=theta_p,theta_pp=theta_pp,Y_pp=Y_pp,R_pp=R_pp))
 }
 
 if(F){

diff --git a/R/tvt.eSIR.R b/R/tvt.eSIR.R
@@ -384,7 +384,7 @@ out_table<-matrix(c(theta_p_mean,theta_p_ci,R0_p_mean,R0_p_ci,gamma_p_mean,gamma
 
   if(save_files) write.csv(out_table,file=paste0(file_add,casename,"_summary.csv"))
 
-  return(list(casename=casename,incidence_mean=incidence_mean,incidence_ci=incidence_ci,out_table=out_table,forecast_infection=plot,forecast_removed=plot2,first_stat_mean=as.character(first_order_change_date),first_stat_ci=as.character(first_order_change_date_ci),second_stat_mean=as.character(second_order_change_date),second_stat_ci=as.character(second_order_change_date_ci),dic_val=dic_val,plot_list=plot_list,plot2_list=plot2_list))
+  return(list(casename=casename,incidence_mean=incidence_mean,incidence_ci=incidence_ci,out_table=out_table,forecast_infection=plot,forecast_removed=plot2,first_stat_mean=as.character(first_order_change_date),first_stat_ci=as.character(first_order_change_date_ci),second_stat_mean=as.character(second_order_change_date),second_stat_ci=as.character(second_order_change_date_ci),dic_val=dic_val,plot_list=plot_list,plot2_list=plot2_list,theta_p=theta_p,theta_pp=theta_pp,Y_pp=Y_pp,R_pp=R_pp))
 }
 
 if(F){
-Original file line number
+Diff line change
@@ Expand Up @@
       if(save_files) write.csv(out_table,file=paste0(file_add,casename,"_summary.csv"))
-      return(list(casename=casename,incidence_mean=incidence_mean,incidence_ci=incidence_ci,out_table=out_table,forecast_infection=plot,forecast_removed=plot2,first_stat_mean=as.character(first_order_change_date),first_stat_ci=as.character(first_order_change_date_ci),second_stat_mean=as.character(second_order_change_date),second_stat_ci=as.character(second_order_change_date_ci),dic_val=dic_val,plot_list=plot_list,plot2_list=plot2_list))
+      return(list(casename=casename,incidence_mean=incidence_mean,incidence_ci=incidence_ci,out_table=out_table,forecast_infection=plot,forecast_removed=plot2,first_stat_mean=as.character(first_order_change_date),first_stat_ci=as.character(first_order_change_date_ci),second_stat_mean=as.character(second_order_change_date),second_stat_ci=as.character(second_order_change_date_ci),dic_val=dic_val,plot_list=plot_list,plot2_list=plot2_list,theta_p=theta_p,theta_pp=theta_pp,Y_pp=Y_pp,R_pp=R_pp))
     }
     if(F){
@@ Expand Down @@