Date created: 20:44:49 21 Oct 2022.
Install commands that are unavailable or out of date.
. *<< dd_do : noout > >
. clear all
.
. cap noi which predictms
C:\Users\CISS Fondecyt\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:\Users\CISS Fondecyt\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:\Users\CISS Fondecyt\ado\plus\s\sumat.ado
*! Part of package matrixtools v. 0.25
*! 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:\Users\CISS Fondecyt\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:\Users\CISS Fondecyt\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"
. 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)'"
. }
. else display "This file does not exist"
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)"
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)'"
Location= C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date: 21 Oct 2022, considering an OS Windows for
> the user: CISS Fondecyt
. }
. else display "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"
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)'"
Location= C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date: 21 Oct 2022, considering an OS Windows for
> the user: CISS Fondecyt
. }
. else display "This file does not exist"
.
Path data= C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)_mult_state_ags;
Timestamp: 21 Oct 2022, considering that is a Windows OS for the username: CISS Fondecyt
The file is located and named as: C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/archivo_multiestado0_oct_22_corr3_3y.dta
timevar0 not found
Defined a transition matrix.
. *https://twitter.com/RDAnalyticsAB/status/1480550270773911556/photo/1
.
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3.sters"
.
. local labs `" "Est" "Cons" "Par1" "Par2" "Par3" "Time in days""'
. esttab m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn m2_6_ggam m2_7_ggam m2_8_ggam using tab1_oct_22_corr3.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.html)
The selected models and their parameters can be summarised here:
| 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.45 | [0.39,0.51] | -0.20 | [-0.24,-0.17] | -1.24 | [-1.48,-1.00] | -1.19 | [-1.34,-1.05] | -0.17 | [-0.43,0.10] | 0.02 | [-0.13,0.17] | -0.01 | [-0.49,0.47] | -0.26 | [-0.49,-0.03] | |
| Cons | ||||||||||||||||
| -7.19 | [-7.23,-7.14] | -5.42 | [-5.45,-5.39] | 10.03 | [9.70,10.36] | 9.58 | [9.40,9.75] | 8.06 | [7.75,8.37] | 7.48 | [7.28,7.69] | 6.78 | [5.82,7.75] | 7.30 | [6.98,7.63] | |
| Par1 | ||||||||||||||||
| -0.00 | [-0.00,-0.00] | -0.00 | [-0.00,-0.00] | 1.15 | [1.11,1.19] | 1.23 | [1.18,1.28] | 0.58 | [0.51,0.65] | 0.77 | [0.71,0.82] | 0.72 | [0.59,0.85] | 0.58 | [0.48,0.68] | |
| 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.17 | [-0.33,-0.02] | -0.64 | [-0.87,-0.42] | -1.56 | [-2.48,-0.64] | -0.39 | [-0.73,-0.06] | |||||||||
| Dist | gompertz | gompertz | lognormal | ggamma | lognormal | ggamma | ggamma | ggamma | ||||||||
| Log-likelihood | -44369.46 | -92016.08 | -13892.34 | -36434.4 | -4088.633 | -11852.53 | -1140.921 | -3962.794 | ||||||||
| Cases | 22452 | 22452 | 5356 | 13067 | 1518 | 4046 | 447 | 1305 | ||||||||
| Note. | ||||||||||||||||
. *To check parameters
. *cap qui noi estimates restore m2_1_rp9
. *ereturn display
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3.sters"
.
.
.
. local labs `" "Est" "Cons" "Knots" "Time in days""'
. esttab m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2 using tab1_oct_22_corr3_rp.html, ///
> varlabels(none) ///
> eqlabels(`labs', lhs("Parameters")) /// transform(#1: exp(@) exp(@)) ///
> label cells ("b(fmt(2) label(Coef)) ci(fmt(2)par label(CI95%))") compress nogap /// /* _sign */
> stats(family1 ndistap1 ll N, fmt(2 %9.0g %9.0g %9.0g 3 2 %9.0g) labels(Dist df Log-likelihood "Cases")) style(html) ///
> title(Selected models: Royston-Parmar) ///
> mtitles("1st tr" "2nd tr" "3rd tr" "4th tr" "5th tr" "6th tr" "7th tr" "8th tr") ///
> hlinechar(`=char(151)') ///
> addnotes("Note.") sfmt(%15,3g) plain wide replace
(output written to tab1_oct_22_corr3_rp.html)
The selected models and their parameters can be summarised here:
| 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.43 | [0.38,0.49] | -0.21 | [-0.24,-0.17] | 0.51 | [0.40,0.63] | 0.45 | [0.39,0.51] | 0.14 | [-0.09,0.36] | 0.00 | [-0.11,0.11] | -0.09 | [-0.49,0.32] | 0.14 | [-0.05,0.33] | |
| Cons | ||||||||||||||||
| -3.00 | [-3.07,-2.93] | -0.80 | [-0.82,-0.77] | -1.88 | [-2.03,-1.73] | -1.67 | [-1.73,-1.61] | -1.56 | [-1.81,-1.31] | -1.29 | [-1.40,-1.18] | -1.79 | [-2.33,-1.24] | -1.35 | [-1.58,-1.11] | |
| Knots | ||||||||||||||||
| 1.88 | [1.82,1.95] | 1.01 | [1.00,1.03] | 0.82 | [0.78,0.86] | 0.81 | [0.79,0.84] | 1.03 | [0.94,1.12] | 0.97 | [0.91,1.02] | 1.14 | [0.94,1.33] | 1.12 | [1.02,1.22] | |
| -0.02 | [-0.05,0.01] | 0.32 | [0.30,0.33] | 0.05 | [0.02,0.08] | 0.10 | [0.08,0.12] | 0.20 | [0.11,0.28] | 0.24 | [0.19,0.28] | 0.37 | [0.18,0.56] | 0.31 | [0.23,0.39] | |
| 0.59 | [0.56,0.62] | 0.24 | [0.23,0.25] | -0.01 | [-0.02,0.01] | 0.02 | [0.01,0.03] | 0.09 | [0.05,0.14] | 0.09 | [0.06,0.12] | -0.02 | [-0.12,0.09] | |||
| 0.16 | [0.14,0.18] | 0.11 | [0.10,0.11] | 0.04 | [0.03,0.06] | 0.04 | [0.03,0.05] | -0.01 | [-0.03,0.01] | 0.07 | [0.02,0.12] | |||||
| -0.01 | [-0.02,0.00] | 0.04 | [0.04,0.05] | 0.01 | [0.00,0.03] | 0.00 | [-0.00,0.01] | |||||||||
| 0.01 | [0.00,0.01] | 0.02 | [0.02,0.02] | 0.01 | [-0.00,0.02] | 0.00 | [-0.00,0.01] | |||||||||
| -0.01 | [-0.01,-0.00] | 0.01 | [0.01,0.01] | 0.01 | [-0.00,0.01] | 0.00 | [-0.00,0.01] | |||||||||
| -0.02 | [-0.02,-0.01] | 0.01 | [0.00,0.01] | 0.01 | [-0.00,0.01] | |||||||||||
| -0.01 | [-0.01,-0.00] | 0.01 | [0.00,0.01] | 0.00 | [-0.00,0.01] | |||||||||||
| 0.00 | [0.00,0.01] | 0.01 | [0.00,0.01] | |||||||||||||
| Time in days | ||||||||||||||||
| 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] | |||||
| Dist | rp | rp | rp | rp | rp | rp | rp | rp | ||||||||
| df | 9 | 10 | 10 | 7 | 3 | 4 | 4 | 2 | ||||||||
| Log-likelihood | -41618.14 | -91559.78 | -13877.87 | -36441.73 | -4075.904 | -11819.03 | -1137.579 | -3961.152 | ||||||||
| Cases | 22452 | 22452 | 5356 | 13067 | 1518 | 4046 | 447 | 1305 | ||||||||
| Note. | ||||||||||||||||
=============================================================================
=============================================================================
(70,463 missing values generated)
.
. qui estread "${pathdata2}parmodels_m2_jun_22_corr3.sters"
.
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
.
. 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.sters"
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
. 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.sters"
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
.
. 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.sters"
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
.
. 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.sters"
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
.
. 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.sters"
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
. 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.sters"
.
. predictms, transmatrix(mat_nine_states) models(m2_1_gom m2_2_gom m2_3_logn m2_4_ggam m2_5_logn 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 for
> ma de hacer bootstrap
Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 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
Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 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
.
. 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= 01:36:30 23 Oct 2022