<> tag dd_include failed to find file {bf:G:/Mi unidad/Alvacast/SISTRAT 2019 (git... r(601); <> tag dd_include failed to find file {bf:G:\Mi unidad\Alvacast\SISTRAT 2019 (git... r(601); <> tag dd_include failed to find file {bf:C:\Users\CISS Fondecyt\Mi unidad\Alvaca... r(601);

Database (step 2)

Date created: 18:13:33 8 Nov 2022.

Install commands that are unavailable or out of date.

. *<< dd_do : noout > >
. clear all

. 
. cap noi which predictms
c:\ado\plus\p\predictms.ado
*! version 4.4.0 22dec2021 MJC

. if _rc==111 {
.         cap noi net install multistate, from("https://www.mjcrowther.co.uk/code/multistate") 
.         }

. cap noi which merlin
c:\ado\plus\m\merlin.ado
*! version 2.1.5 18mar2022 MJC

. if _rc==111 {
.         cap noi net install merlin, from("https://www.mjcrowther.co.uk/code/merlin/") 
.         }

. cap noi which sumat
c:\ado\plus\s\sumat.ado
*! Part of package matrixtools v. 0.28
*! Support: Niels Henrik Bruun, niels.henrik.bruun@gmail.com
*! 2021-01-03 toxl added

. if _rc==111 {
.         cap noi scc install matrixtools
.         }

. cap noi which estwrite
c:\ado\plus\e\estwrite.ado
*! version 1.2.4 04sep2009
*! version 1.0.1 15may2007 (renamed from -eststo- to -estwrite-; -append- added)
*! version 1.0.0 29apr2005 Ben Jann (ETH Zurich)

. if _rc==111 {
.         cap noi ssc install estwrite
.         }

. cap noi which winsor2
c:\ado\plus\w\winsor2.ado
*! Inspirit of -winsor-(NJ Cox) and -winsorizeJ-(J Caskey)
*! Lian Yujun, arlionn@163.com, 2013-12-25
*! 1.1 2014.12.16

. if _rc==111 {
.         cap noi ssc install winsor2
.         }       

.

We need to obtain the file and the work folder.

. mata : st_numscalar("OK", direxists("/volumes/sdrive/data//"))

. if scalar(OK) == 1 {
.         cap noi cd "/volumes/sdrive/data//"
.         global pathdata "/volumes/sdrive/data//"
.         di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }

. else display "This file does not exist"
This file does not exist

. 
. mata : st_numscalar("OK", direxists("E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\"))

. if scalar(OK) == 1 {
.         cap noi cd "E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags
.         global pathdata "E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
.         global pathdata2 "E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
.         di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
Location= E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date:  8 Nov 2022, considering an OS Windows for the user: andre
. }

. else display "This file does not exist"

. 
. mata : st_numscalar("OK", direxists("G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\"))

. if scalar(OK) == 1 {
.         cap noi cd "G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
.         global pathdata "G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
.         global pathdata2 "G:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
.         di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }

. else display "This file does not exist"
This file does not exist

.                 
. mata : st_numscalar("OK", direxists("C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\"))

. if scalar(OK) == 1 {
.         cap noi cd "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)"
.         global pathdata "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
.         global pathdata2 "C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
.         di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }

. else display "This file does not exist"
This file does not exist

. 
. mata : st_numscalar("OK", direxists("C:\Users\andre\Desktop\_mult_state_ags\"))

. if scalar(OK) == 1 {
.         cap noi cd "C:\Users\andre\Desktop\_mult_state_ags"
.         global pathdata "C:\Users\andre\Desktop\_mult_state_ags"
.         global pathdata2 "C:/Users/andre/Desktop/_mult_state_ags/"
.         di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }

. else display "This file does not exist"
This file does not exist

. 
. mata : st_numscalar("OK", direxists("C:\Users\CISS Fondecyt\OneDrive\Documentos\"))

. if scalar(OK) == 1 {
.         cap noi cd "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
.         global pathdata "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
.         global pathdata2 "C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
.         di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }

. else display "This file does not exist"
This file does not exist

.

Path data= E:\Mi unidad\Alvacast\SISTRAT 2019 (github)_mult_state_ags;

Timestamp: 8 Nov 2022, considering that is a Windows OS for the username: andre

The file is located and named as: E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/archivo_multiestado0_oct_22_corr3_cc_3y.dta

timevar0 not found

Defined a transition matrix.

. *https://twitter.com/RDAnalyticsAB/status/1480550270773911556/photo/1
. //Cambia el 3 desde logn por ggam
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. local labs `" "Est" "Cons" "Par1" "Par2" "Par3" "Time in days""' 

. esttab m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam using tab1_oct_22_corr3_cc.html, /// transform(1 exp(@) @ @ @) ///
>         label cells ("b(fmt(2) label(Coef)) ci(fmt(2)par label(CI95%))") compress nogap /// /* _sign */
>         stats(family1 ll N, fmt(2 %9.0g %9.0g 3 2 %9.0g) labels(Dist Log-likelihood "Cases")) style(html)  ///
>         title(Selected models: Parametric) ///
>         mtitles("1st tr" "2nd tr" "3rd tr" "4th tr" "5th tr" "6th tr" "7th tr" "8th tr") ///
>         hlinechar(`=char(151)') ///
>         varlabels(none) ///
>         eqlabels(`labs', lhs("Parameters")) ///         //* eqlabels(none, lhs("Parameters")) ///
>         addnotes("Note.") sfmt(%15,3g) plain wide replace
(output written to tab1_oct_22_corr3_cc.html)

The selected models and their parameters can be summarised here:

Selected models: Parametric
1st tr 2nd tr 3rd tr 4th tr 5th tr 6th tr 7th tr 8th tr
Coef CI95% Coef CI95% Coef CI95% Coef CI95% Coef CI95% Coef CI95% Coef CI95% Coef CI95%
Est
0.40[0.34,0.47] -0.21[-0.26,-0.17] -0.92[-1.18,-0.67] -1.17[-1.32,-1.02] -0.22[-0.51,0.07] 0.09[-0.08,0.25] 0.16[-0.39,0.70] -0.32[-0.57,-0.06]
Cons
-7.19[-7.25,-7.13] -5.42[-5.45,-5.39] 9.68[9.36,10.01] 9.29[9.13,9.45] 7.63[7.20,8.07] 7.41[7.19,7.63] 6.60[5.45,7.75] 7.33[7.00,7.67]
Par1
-0.00[-0.00,-0.00] -0.00[-0.00,-0.00] 0.72[0.54,0.90] 1.09[1.06,1.12] 0.69[0.58,0.80] 0.72[0.66,0.79] 0.72[0.57,0.86] 0.54[0.43,0.65]
Par2
-0.00[-0.00,-0.00] -0.00[-0.00,-0.00] 0.00[-0.00,0.00] 0.00[0.00,0.00] 0.00[-0.00,0.00] 0.00[0.00,0.00]
Par3
0.69[0.40,0.98] -0.73[-1.26,-0.21] -0.50[-0.75,-0.26] -1.40[-2.45,-0.34] -0.35[-0.70,0.01]
Dist gompertz gompertz ggamma lognormal ggamma ggamma ggamma ggamma
Log-likelihood -30076.72 -62806.16 -10311.65 -26946.33 -3101.883 -9113.612 -877.0382 -3196.825
Cases 15494 15494 3617 8903 1153 3049 341 1011
Note.

=============================================================================

Transition probabilities & Lengths of Stay, Parametric Models

=============================================================================

(48,882 missing values generated)

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(1) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

. 

. rename (_prob_at1_1_*lci) (fprob_from1a_*lci_3y)

. rename (_prob_at1_1_*uci) (fprob_from1a_*uci_3y)

. rename (_prob_at1_1_*) (fprob_from1a_*_3y)

. 
. rename (_los_at1_1_*lci) (flos_from1a_*lci_3y)

. rename (_los_at1_1_*uci) (flos_from1a_*uci_3y)

. rename (_los_at1_1_*) (flos_from1a_*_3y)

. 
. rename (_prob_at2_1_*lci) (fprob_from1b_*lci_3y)

. rename (_prob_at2_1_*uci) (fprob_from1b_*uci_3y)

. rename (_prob_at2_1_*) (fprob_from1b_*_3y)

. 
. rename (_los_at2_1_*lci) (flos_from1b_*lci_3y)

. rename (_los_at2_1_*uci) (flos_from1b_*uci_3y)

. rename (_los_at2_1_*) (flos_from1b_*_3y)

. 
. rename (_ratio_prob_at2_1_*_lci) (fratio_prob_from1_*_lci_3y)

. rename (_ratio_prob_at2_1_*_uci) (fratio_prob_from1_*_uci_3y)

. rename (_ratio_prob_at2_1_*) (fratio_prob_from1_*_3y)

. rename (_diff_prob_at2_1_*_lci) (fdiff_prob_from1_*_lci_3y)

. rename (_diff_prob_at2_1_*_uci) (fdiff_prob_from1_*_uci_3y)

. rename (_diff_prob_at2_1_*) (fdiff_prob_from1_*_3y)  

. 
. rename (_ratio_los_at2_1_*_lci) (fratio_los_from1_*_lci_3y)

. rename (_ratio_los_at2_1_*_uci) (fratio_los_from1_*_uci_3y)

. rename (_ratio_los_at2_1_*) (fratio_los_from1_*_3y)

. rename (_diff_los_at2_1_*_lci) (fdiff_los_from1_*_lci_3y)

. rename (_diff_los_at2_1_*_uci) (fdiff_los_from1_*_uci_3y)

. rename (_diff_los_at2_1_*) (fdiff_los_from1_*_3y)  

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(2) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000


. rename (_prob_at1_2_*lci) (fprob_from2a_*lci_3y)

. rename (_prob_at1_2_*uci) (fprob_from2a_*uci_3y)

. rename (_prob_at1_2_*) (fprob_from2a_*_3y)

. 
. rename (_los_at1_2_*lci) (flos_from2a_*lci_3y)

. rename (_los_at1_2_*uci) (flos_from2a_*uci_3y)

. rename (_los_at1_2_*) (flos_from2a_*_3y)

. 
. rename (_prob_at2_2_*lci) (fprob_from2b_*lci_3y)

. rename (_prob_at2_2_*uci) (fprob_from2b_*uci_3y)

. rename (_prob_at2_2_*) (fprob_from2b_*_3y)

. 
. rename (_los_at2_2_*lci) (flos_from2b_*lci_3y)

. rename (_los_at2_2_*uci) (flos_from2b_*uci_3y)

. rename (_los_at2_2_*) (flos_from2b_*_3y)

. 
. rename (_ratio_prob_at2_2_*_lci) (fratio_prob_from2_*_lci_3y)

. rename (_ratio_prob_at2_2_*_uci) (fratio_prob_from2_*_uci_3y)

. rename (_ratio_prob_at2_2_*) (fratio_prob_from2_*_3y)

. rename (_diff_prob_at2_2_*_lci) (fdiff_prob_from2_*_lci_3y)

. rename (_diff_prob_at2_2_*_uci) (fdiff_prob_from2_*_uci_3y)

. rename (_diff_prob_at2_2_*) (fdiff_prob_from2_*_3y)  

. 
. rename (_ratio_los_at2_2_*_lci) (fratio_los_from2_*_lci_3y)

. rename (_ratio_los_at2_2_*_uci) (fratio_los_from2_*_uci_3y)

. rename (_ratio_los_at2_2_*) (fratio_los_from2_*_3y)

. rename (_diff_los_at2_2_*_lci) (fdiff_los_from2_*_lci_3y)

. rename (_diff_los_at2_2_*_uci) (fdiff_los_from2_*_uci_3y)

. rename (_diff_los_at2_2_*) (fdiff_los_from2_*_3y) 

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(3) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

. 

. rename (_prob_at1_3_*lci) (fprob_from3a_*lci_3y)

. rename (_prob_at1_3_*uci) (fprob_from3a_*uci_3y)

. rename (_prob_at1_3_*) (fprob_from3a_*_3y)

. 
. rename (_los_at1_3_*lci) (flos_from3a_*lci_3y)

. rename (_los_at1_3_*uci) (flos_from3a_*uci_3y)

. rename (_los_at1_3_*) (flos_from3a_*_3y)

. 
. rename (_prob_at2_3_*lci) (fprob_from3b_*lci_3y)

. rename (_prob_at2_3_*uci) (fprob_from3b_*uci_3y)

. rename (_prob_at2_3_*) (fprob_from3b_*_3y)

. 
. rename (_los_at2_3_*lci) (flos_from3b_*lci_3y)

. rename (_los_at2_3_*uci) (flos_from3b_*uci_3y)

. rename (_los_at2_3_*) (flos_from3b_*_3y)

. 
. rename (_ratio_prob_at2_3_*_lci) (fratio_prob_from3_*_lci_3y)

. rename (_ratio_prob_at2_3_*_uci) (fratio_prob_from3_*_uci_3y)

. rename (_ratio_prob_at2_3_*) (fratio_prob_from3_*_3y)

. rename (_diff_prob_at2_3_*_lci) (fdiff_prob_from3_*_lci_3y)

. rename (_diff_prob_at2_3_*_uci) (fdiff_prob_from3_*_uci_3y)

. rename (_diff_prob_at2_3_*) (fdiff_prob_from3_*_3y)  

. 
. rename (_ratio_los_at2_3_*_lci) (fratio_los_from3_*_lci_3y)

. rename (_ratio_los_at2_3_*_uci) (fratio_los_from3_*_uci_3y)

. rename (_ratio_los_at2_3_*) (fratio_los_from3_*_3y)

. rename (_diff_los_at2_3_*_lci) (fdiff_los_from3_*_lci_3y)

. rename (_diff_los_at2_3_*_uci) (fdiff_los_from3_*_uci_3y)

. rename (_diff_los_at2_3_*) (fdiff_los_from3_*_3y) 

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(4) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

. 

. rename (_prob_at1_4_*lci) (fprob_from4a_*lci_3y)

. rename (_prob_at1_4_*uci) (fprob_from4a_*uci_3y)

. rename (_prob_at1_4_*) (fprob_from4a_*_3y)

. 
. rename (_los_at1_4_*lci) (flos_from4a_*lci_3y)

. rename (_los_at1_4_*uci) (flos_from4a_*uci_3y)

. rename (_los_at1_4_*) (flos_from4a_*_3y)

. 
. rename (_prob_at2_4_*lci) (fprob_from4b_*lci_3y)

. rename (_prob_at2_4_*uci) (fprob_from4b_*uci_3y)

. rename (_prob_at2_4_*) (fprob_from4b_*_3y)

. 
. rename (_los_at2_4_*lci) (flos_from4b_*lci_3y)

. rename (_los_at2_4_*uci) (flos_from4b_*uci_3y)

. rename (_los_at2_4_*) (flos_from4b_*_3y)

. 
. rename (_ratio_prob_at2_4_*_lci) (fratio_prob_from4_*_lci_3y)

. rename (_ratio_prob_at2_4_*_uci) (fratio_prob_from4_*_uci_3y)

. rename (_ratio_prob_at2_4_*) (fratio_prob_from4_*_3y)

. rename (_diff_prob_at2_4_*_lci) (fdiff_prob_from4_*_lci_3y)

. rename (_diff_prob_at2_4_*_uci) (fdiff_prob_from4_*_uci_3y)

. rename (_diff_prob_at2_4_*) (fdiff_prob_from4_*_3y)  

. 
. rename (_ratio_los_at2_4_*_lci) (fratio_los_from4_*_lci_3y)

. rename (_ratio_los_at2_4_*_uci) (fratio_los_from4_*_uci_3y)

. rename (_ratio_los_at2_4_*) (fratio_los_from4_*_3y)

. rename (_diff_los_at2_4_*_lci) (fdiff_los_from4_*_lci_3y)

. rename (_diff_los_at2_4_*_uci) (fdiff_los_from4_*_uci_3y)

. rename (_diff_los_at2_4_*) (fdiff_los_from4_*_3y) 

. 

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(5) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

. 

. rename (_prob_at1_5_*lci) (fprob_from5a_*lci_3y)

. rename (_prob_at1_5_*uci) (fprob_from5a_*uci_3y)

. rename (_prob_at1_5_*) (fprob_from5a_*_3y)

. 
. rename (_los_at1_5_*lci) (flos_from5a_*lci_3y)

. rename (_los_at1_5_*uci) (flos_from5a_*uci_3y)

. rename (_los_at1_5_*) (flos_from5a_*_3y)

. 
. rename (_prob_at2_5_*lci) (fprob_from5b_*lci_3y)

. rename (_prob_at2_5_*uci) (fprob_from5b_*uci_3y)

. rename (_prob_at2_5_*) (fprob_from5b_*_3y)

. 
. rename (_los_at2_5_*lci) (flos_from5b_*lci_3y)

. rename (_los_at2_5_*uci) (flos_from5b_*uci_3y)

. rename (_los_at2_5_*) (flos_from5b_*_3y)

. 
. rename (_ratio_prob_at2_5_*_lci) (fratio_prob_from5_*_lci_3y)

. rename (_ratio_prob_at2_5_*_uci) (fratio_prob_from5_*_uci_3y)

. rename (_ratio_prob_at2_5_*) (fratio_prob_from5_*_3y)

. rename (_diff_prob_at2_5_*_lci) (fdiff_prob_from5_*_lci_3y)

. rename (_diff_prob_at2_5_*_uci) (fdiff_prob_from5_*_uci_3y)

. rename (_diff_prob_at2_5_*) (fdiff_prob_from5_*_3y)  

. 
. rename (_ratio_los_at2_5_*_lci) (fratio_los_from5_*_lci_3y)

. rename (_ratio_los_at2_5_*_uci) (fratio_los_from5_*_uci_3y)

. rename (_ratio_los_at2_5_*) (fratio_los_from5_*_3y)

. rename (_diff_los_at2_5_*_lci) (fdiff_los_from5_*_lci_3y)

. rename (_diff_los_at2_5_*_uci) (fdiff_los_from5_*_uci_3y)

. rename (_diff_los_at2_5_*) (fdiff_los_from5_*_3y)  

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(6) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000


. rename (_prob_at1_6_*lci) (fprob_from6a_*lci_3y)

. rename (_prob_at1_6_*uci) (fprob_from6a_*uci_3y)

. rename (_prob_at1_6_*) (fprob_from6a_*_3y)

. 
. rename (_los_at1_6_*lci) (flos_from6a_*lci_3y)

. rename (_los_at1_6_*uci) (flos_from6a_*uci_3y)

. rename (_los_at1_6_*) (flos_from6a_*_3y)

. 
. rename (_prob_at2_6_*lci) (fprob_from6b_*lci_3y)

. rename (_prob_at2_6_*uci) (fprob_from6b_*uci_3y)

. rename (_prob_at2_6_*) (fprob_from6b_*_3y)

. 
. rename (_los_at2_6_*lci) (flos_from6b_*lci_3y)

. rename (_los_at2_6_*uci) (flos_from6b_*uci_3y)

. rename (_los_at2_6_*) (flos_from6b_*_3y)

. 
. rename (_ratio_prob_at2_6_*_lci) (fratio_prob_from6_*_lci_3y)

. rename (_ratio_prob_at2_6_*_uci) (fratio_prob_from6_*_uci_3y)

. rename (_ratio_prob_at2_6_*) (fratio_prob_from6_*_3y)

. rename (_diff_prob_at2_6_*_lci) (fdiff_prob_from6_*_lci_3y)

. rename (_diff_prob_at2_6_*_uci) (fdiff_prob_from6_*_uci_3y)

. rename (_diff_prob_at2_6_*) (fdiff_prob_from6_*_3y)  

. 
. rename (_ratio_los_at2_6_*_lci) (fratio_los_from6_*_lci_3y)

. rename (_ratio_los_at2_6_*_uci) (fratio_los_from6_*_uci_3y)

. rename (_ratio_los_at2_6_*) (fratio_los_from6_*_3y)

. rename (_diff_los_at2_6_*_lci) (fdiff_los_from6_*_lci_3y)

. rename (_diff_los_at2_6_*_uci) (fdiff_los_from6_*_uci_3y)

. rename (_diff_los_at2_6_*) (fdiff_los_from6_*_3y) 

. 
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3_cc.sters"

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_ggam m2_4_logn m2_5_ggam m2_6_ggam m2_7_ggam m2_8_ggam) ///
>         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
>         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
>         from(7) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una forma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

Calculating CIs via parametric bootstrap for at2() (2000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
..................................................  1050
..................................................  1100
..................................................  1150
..................................................  1200
..................................................  1250
..................................................  1300
..................................................  1350
..................................................  1400
..................................................  1450
..................................................  1500
..................................................  1550
..................................................  1600
..................................................  1650
..................................................  1700
..................................................  1750
..................................................  1800
..................................................  1850
..................................................  1900
..................................................  1950
..................................................  2000

. 

. rename (_prob_at1_7_*lci) (fprob_from7a_*lci_3y)

. rename (_prob_at1_7_*uci) (fprob_from7a_*uci_3y)

. rename (_prob_at1_7_*) (fprob_from7a_*_3y)

. 
. rename (_los_at1_7_*lci) (flos_from7a_*lci_3y)

. rename (_los_at1_7_*uci) (flos_from7a_*uci_3y)

. rename (_los_at1_7_*) (flos_from7a_*_3y)

. 
. rename (_prob_at2_7_*lci) (fprob_from7b_*lci_3y)

. rename (_prob_at2_7_*uci) (fprob_from7b_*uci_3y)

. rename (_prob_at2_7_*) (fprob_from7b_*_3y)

. 
. rename (_los_at2_7_*lci) (flos_from7b_*lci_3y)

. rename (_los_at2_7_*uci) (flos_from7b_*uci_3y)

. rename (_los_at2_7_*) (flos_from7b_*_3y)

. 
. rename (_ratio_prob_at2_7_*_lci) (fratio_prob_from7_*_lci_3y)

. rename (_ratio_prob_at2_7_*_uci) (fratio_prob_from7_*_uci_3y)

. rename (_ratio_prob_at2_7_*) (fratio_prob_from7_*_3y)

. rename (_diff_prob_at2_7_*_lci) (fdiff_prob_from7_*_lci_3y)

. rename (_diff_prob_at2_7_*_uci) (fdiff_prob_from7_*_uci_3y)

. rename (_diff_prob_at2_7_*) (fdiff_prob_from7_*_3y)  

. 
. rename (_ratio_los_at2_7_*_lci) (fratio_los_from7_*_lci_3y)

. rename (_ratio_los_at2_7_*_uci) (fratio_los_from7_*_uci_3y)

. rename (_ratio_los_at2_7_*) (fratio_los_from7_*_3y)

. rename (_diff_los_at2_7_*_lci) (fdiff_los_from7_*_lci_3y)

. rename (_diff_los_at2_7_*_uci) (fdiff_los_from7_*_uci_3y)

. rename (_diff_los_at2_7_*) (fdiff_los_from7_*_3y) 

.

Saved at= 18:13:04 10 Nov 2022