Database (step 3)
Date created: 14:16:44 26 Jul 2021.
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.3.0 14mar2021 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.0.2 19mar2021 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)\_WO vs MG\"))
. if scalar(OK) == 1 {
. cap noi cd "E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG"
E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG
. global pathdata "E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG"
. global pathdata2 "E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_WO vs MG/"
. 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)\_WO vs MG; Date: 26 Jul 2021, considering an OS Windows for the user: andre
. }
. else display "This file does not exist"
.
. mata : st_numscalar("OK", direxists("C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG\"))
. 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)\_WO vs MG"
. global pathdata2 "C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_WO vs MG/"
. 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\_WO vs MG\"))
. if scalar(OK) == 1 {
. cap noi cd "C:\Users\andre\Desktop\_WO vs MG"
. 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)\_WO vs MG"
. global pathdata "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG"
. global pathdata2 "C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_WO vs MG/"
. 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)\_WO vs MG\"))
. if scalar(OK) == 1 {
. cap noi cd "G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG"
. global pathdata "G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_WO vs MG"
. global pathdata2 "G:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_WO vs MG/"
. 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)_WO vs MG;
Timestamp: 26 Jul 2021, considering that is a Windows OS for the username: andre
First we open the files and drop the variables that would mistakenly amplify the sample, and define labels.
The file is located and named as: E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_WO vs MG/archivo_multiestado0_ago_c.dta
Calculated transition probabilities from a determined states in times 90, 365 (1 year), 1095 (3 years) & 1826 (5 years).
variable trp_ajprob* not found
variable P_AJ* not found
From state is an absorbing state
r(198);
=============================================================================
Figures of Transition Probabilities
=============================================================================
. *Assessing and relaxing the Markov assumption in the illness-death model
. *Jonathan Broomfield1, Caroline E. Weibull2, and Michael J. Crowther
. tw (rarea fprob_from1a_2_lci fprob_from1a_2_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from1b_2_lci fprob_from1b_2_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from1a_2 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from1b_2 timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_12 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_12 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).4, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(11) ring(0) c(1)) ///
> title("Transition Probabilities") name(trans_models12, replace) ///
> subtitle("from Admission to Treatment Completion") /// *text(0.2 4 "Outpatient", size(small)) /// *text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_12.gph, replace
. tw (rarea fprob_from1a_3_lci fprob_from1a_3_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from1b_3_lci fprob_from1b_3_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from1a_3 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from1b_3 timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_13 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_13 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).3, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(5) ring(0) c(1)) ///
> title("Transition Probabilities") name(trans_models13, replace) ///
> subtitle("from Admission to Discharge w/o Completion") /// *text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_23.gph, replace
. tw (rarea fprob_from1a_4_lci fprob_from1a_4_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from1b_4_lci fprob_from1b_4_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from1a_4 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from1b_4 timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_14 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_14 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).3, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(5) ring(0) c(1)) ///
> title("Transition Probabilities") name(trans_models14, replace) ///
> subtitle("from Admission to Readmission") /// *text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_23.gph, replace
. tw (rarea fprob_from2a_4_lci fprob_from2a_4_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from2b_4_lci fprob_from2b_4_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from2a_4 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from2b_4 timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_24 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_24 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).4, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(5) ring(0) c(1)) ///
> title("Transition Probabilities") name(trans_models24, replace) ///
> subtitle("from Treatment Completion to Readmission") /// *text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_34.gph, replace
. tw (rarea fprob_from3a_4_lci fprob_from3a_4_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from3b_4_lci fprob_from3b_4_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from3a_4 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from3b_4 timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_34 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_34 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.2)1, ///
> angle(h) format(%4.1f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(11) ring(0) c(1)) ///
> title("Transition Probabilities") name(trans_models34, replace) ///
> subtitle("from Discharge w/o Completion to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_45.gph, replace
=============================================================================
Figures of Transition Probabilities (RP)
=============================================================================
. *Assessing and relaxing the Markov assumption in the illness-death model
. *Jonathan Broomfield1, Caroline E. Weibull2, and Michael J. Crowther
. tw (rarea fprob_from1a_2_lci_rp fprob_from1a_2_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from1b_2_lci_rp fprob_from1b_2_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from1a_2_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from1b_2_rp timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_12 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_12 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).4, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP (RP)" 4 "WO (RP)" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(11) ring(0) c(1)) ///
> title("Transition Probabilities (RP)") name(trans_models12rp, replace) ///
> subtitle("from Admission to Treatment Completion") /// *text(0.2 4 "Outpatient", size(small)) /// *text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
variable fprob_from1a_2_lci_rp not found
r(111);
. * graph export prob_trans_12_rp.gph, replace
. tw (rarea fprob_from1a_3_lci_rp fprob_from1a_3_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from1b_3_lci_rp fprob_from1b_3_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from1a_3_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from1b_3_rp timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_13 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_13 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).3, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP (RP)" 4 "WO (RP)" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(5) ring(0) c(1)) ///
> title("Transition Probabilities (RP)") name(trans_models13rp, replace) ///
> subtitle("from Admission to Discharge w/o Completion") /// *text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
variable fprob_from1a_3_lci_rp not found
r(111);
. *graph export prob_trans_23_rp.gph, replace
. tw (rarea fprob_from1a_4_lci_rp fprob_from1a_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from1b_4_lci_rp fprob_from1b_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from1a_4_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from1b_4_rp timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_14 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_14 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).3, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP (RP)" 4 "WO (RP)" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(5) ring(0) c(1)) ///
> title("Transition Probabilities (RP)") name(trans_models14rp, replace) ///
> subtitle("from Admission to Readmission") /// *text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
variable fprob_from1a_4_lci_rp not found
r(111);
. *graph export prob_trans_23_rp.gph, replace
. tw (rarea fprob_from2a_4_lci_rp fprob_from2a_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from2b_4_lci_rp fprob_from2b_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from2a_4_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from2b_4_rp timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_24 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_24 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.05).4, ///
> angle(h) format(%4.2f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP (RP)" 4 "WO (RP)" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(5) ring(0) c(1)) ///
> title("Transition Probabilities (RP)") name(trans_models24rp, replace) ///
> subtitle("from Treatment Completion to Readmission") /// *text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_34_rp.gph, replace
. tw (rarea fprob_from3a_4_lci_rp fprob_from3a_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea fprob_from3b_4_lci_rp fprob_from3b_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line fprob_from3a_4_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line fprob_from3b_4_rp timevar01, sort connect(stairstep stairstep) lcolor(black)) ///
> (line trp_ajprob_30_1826_34 _t_scaled if tipo_de_programa_2==1 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50) lpattern("-")) /
> //
> (line trp_ajprob_30_1826_34 _t_scaled if tipo_de_programa_2==2 & _t_scaled<=5, sort connect(stairstep stairstep) lcolor(red%50)), ///
> xtitle(Time (years)) ytitle(Probability) xlab(0(.3)5,labsize(small)) ylab(0(0.2)1, ///
> angle(h) format(%4.1f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP (RP)" 4 "WO (RP)" 5 "GP(LMAJ)" 6 "WO(LMAJ)") pos(11) ring(0) c(1)) ///
> title("Transition Probabilities (RP)") name(trans_models34rp, replace) ///
> subtitle("from Discharge w/o Completion to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export prob_trans_45_rp.gph, replace
=============================================================================
Figures of Expected Lengths of Stay
=============================================================================
. tw (rarea flos_from1a_2_lci flos_from1a_2_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from1b_2_lci flos_from1b_2_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from1a_2 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from1b_2 timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)365, /// *730
> angle(h) format(%4.0f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO") pos(11) ring(0) c(1)) ///
> title("Length of Stay") name(los_models12, replace) ///
> subtitle("from Admission to Treatment Completion") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models12.gph, replace
. tw (rarea flos_from1a_3_lci flos_from1a_3_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from1b_3_lci flos_from1b_3_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from1a_3 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from1b_3 timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)400, /// *730
> angle(h) format(%4.0f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO") pos(11) ring(0) c(1)) ///
> title("Length of Stay") name(los_models13, replace) ///
> subtitle("from Readmission to Treatment w/o Completion") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models23.gph, replace
. tw (rarea flos_from1a_4_lci flos_from1a_4_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from1b_4_lci flos_from1b_4_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from1a_4 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from1b_4 timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)540, /// *730
> angle(h) format(%4.0f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO") pos(11) ring(0) c(1)) ///
> title("Length of Stay") name(los_models14, replace) ///
> subtitle("from Admission to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models34.gph, replace
. tw (rarea flos_from2a_4_lci flos_from2a_4_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from2b_4_lci flos_from2b_4_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from2a_4 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from2b_4 timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)1080, /// *730
> angle(h) format(%4.0f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO") pos(11) ring(0) c(1)) ///
> title("Length of Stay") name(los_models24, replace) ///
> subtitle("from Treatment Completion to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models45.gph, replace
. tw (rarea flos_from3a_4_lci flos_from3a_4_uci timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from3b_4_lci flos_from3b_4_uci timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from3a_4 timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from3b_4 timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)1080, /// *730
> angle(h) format(%4.0f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO") pos(11) ring(0) c(1)) ///
> title("Length of Stay") name(los_models34, replace) ///
> subtitle("from Discharge w/o Completion to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models45.gph, replace
=============================================================================
Figures of Expected Lengths of Stay (RP)
=============================================================================
. tw (rarea flos_from1a_2_lci_rp flos_from1a_2_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from1b_2_lci_rp flos_from1b_2_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from1a_2_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from1b_2_rp timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)365, /// *730
> angle(h) format(%4.0f)) legend(order(1 "Outpatient CI95%" 2 "Residential CI95%" 3 "Outpatient (RP)" 4 "Residential (RP)") pos(11) ring(0) c(1)) ///
> title("Length of Stay (RP)") name(los_models12_rp, replace) ///
> subtitle("from Admission to Treatment Completion") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
variable flos_from1a_2_lci_rp not found
r(111);
. *graph export rott_los_models12_rp.gph, replace
. tw (rarea flos_from1a_3_lci_rp flos_from1a_3_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from1b_3_lci_rp flos_from1b_3_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from1a_3_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from1b_3_rp timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)400, /// *730
> angle(h) format(%4.0f)) legend(order(1 "Outpatient CI95%" 2 "Residential CI95%" 3 "Outpatient" 4 "Residential") pos(11) ring(0) c(1)) ///
> title("Length of Stay (RP)") name(los_models13_rp, replace) ///
> subtitle("from Admission to Discharge w/o Completion") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
variable flos_from1a_3_lci_rp not found
r(111);
. *graph export rott_los_models13_rp.gph, replace
. tw (rarea flos_from1a_4_lci_rp flos_from1a_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from1b_4_lci_rp flos_from1b_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from1a_4_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from1b_4_rp timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)540, /// *730
> angle(h) format(%4.0f)) legend(order(1 "Outpatient CI95%" 2 "Residential CI95%" 3 "Outpatient" 4 "Residential") pos(11) ring(0) c(1)) ///
> title("Length of Stay (RP)") name(los_models14_rp, replace) ///
> subtitle("from Admission to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
variable flos_from1a_4_lci_rp not found
r(111);
. *graph export rott_los_models14_rp.gph, replace
. tw (rarea flos_from2a_4_lci_rp flos_from2a_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from2b_4_lci_rp flos_from2b_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from2a_4_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from2b_4_rp timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)540, /// *730
> angle(h) format(%4.0f)) legend(order(1 "Outpatient CI95%" 2 "Residential CI95%" 3 "Outpatient" 4 "Residential") pos(11) ring(0) c(1)) ///
> title("Length of Stay (RP)") name(los_models24_rp, replace) ///
> subtitle("from Treatment Completion to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models34_rp.gph, replace
. tw (rarea flos_from3a_4_lci_rp flos_from3a_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) color(gs8%35)) ///
> (rarea flos_from3b_4_lci_rp flos_from3b_4_uci_rp timevar01, sort connect(stairstep stairstep) lcolor(black) color(black%35)) ///
> (line flos_from3a_4_rp timevar01, sort connect(stairstep stairstep) lcolor(gs8) lpattern("-")) ///
> (line flos_from3b_4_rp timevar01, sort connect(stairstep stairstep) lcolor(black)), ///
> xtitle(Time (years)) ytitle(LOS (years)) xlab(0(.3)5,labsize(small)) ylab(0(90)1080, /// *730
> angle(h) format(%4.0f)) legend(order(1 "GP CI95%" 2 "WO CI95%" 3 "GP" 4 "WO") pos(11) ring(0) c(1)) ///
> title("Length of Stay (RP)") name(los_models34_rp, replace) ///
> subtitle("from Treatment w/o Completion to Readmission") /// * text(0.2 4 "Outpatient", size(small)) /// * text(0.29 4 "Residential", size(small)) ///
> graphregion(col(white))
. *graph export rott_los_models45_rp.gph, replace
=============================================================================
Parametric models of differences & ratios- At 0 days
=============================================================================
. estread "${pathdata2}parmodels_m2_ago_c.sters"
-----------------------------------------------------------------------------
name | command depvar npar title
-------------+---------------------------------------------------------------
m2_1_exp | merlin no depvar 11 Fixed effects regression model
m2_1_weib | merlin no depvar 12 Fixed effects regression model
m2_1_gom | merlin no depvar 12 Fixed effects regression model
m2_1_logl | merlin no depvar 12 Fixed effects regression model
m2_1_logn | merlin no depvar 12 Fixed effects regression model
m2_1_ggam | merlin no depvar 13 Fixed effects regression model
m2_1_rp2 | merlin no depvar 13 Fixed effects regression model
m2_1_rp3 | merlin no depvar 14 Fixed effects regression model
m2_1_rp4 | merlin no depvar 15 Fixed effects regression model
m2_1_rp5 | merlin no depvar 16 Fixed effects regression model
m2_1_rp6 | merlin no depvar 17 Fixed effects regression model
m2_1_rp7 | merlin no depvar 18 Fixed effects regression model
m2_1_rp8 | merlin no depvar 19 Fixed effects regression model
m2_1_rp9 | merlin no depvar 20 Fixed effects regression model
m2_1_rp10 | merlin no depvar 21 Fixed effects regression model
m2_2_exp | merlin no depvar 11 Fixed effects regression model
m2_2_weib | merlin no depvar 12 Fixed effects regression model
m2_2_gom | merlin no depvar 12 Fixed effects regression model
m2_2_logl | merlin no depvar 12 Fixed effects regression model
m2_2_logn | merlin no depvar 12 Fixed effects regression model
m2_2_ggam | merlin no depvar 13 Fixed effects regression model
m2_2_rp2 | merlin no depvar 13 Fixed effects regression model
m2_2_rp3 | merlin no depvar 14 Fixed effects regression model
m2_2_rp4 | merlin no depvar 15 Fixed effects regression model
m2_2_rp5 | merlin no depvar 16 Fixed effects regression model
m2_2_rp6 | merlin no depvar 17 Fixed effects regression model
m2_2_rp7 | merlin no depvar 18 Fixed effects regression model
m2_2_rp8 | merlin no depvar 19 Fixed effects regression model
m2_2_rp9 | merlin no depvar 20 Fixed effects regression model
m2_2_rp10 | merlin no depvar 21 Fixed effects regression model
m2_3_exp | merlin no depvar 11 Fixed effects regression model
m2_3_weib | merlin no depvar 12 Fixed effects regression model
m2_3_gom | merlin no depvar 12 Fixed effects regression model
m2_3_logl | merlin no depvar 12 Fixed effects regression model
m2_3_logn | merlin no depvar 12 Fixed effects regression model
m2_3_ggam | merlin no depvar 13 Fixed effects regression model
m2_3_rp2 | merlin no depvar 13 Fixed effects regression model
m2_3_rp3 | merlin no depvar 14 Fixed effects regression model
m2_3_rp4 | merlin no depvar 15 Fixed effects regression model
m2_3_rp5 | merlin no depvar 16 Fixed effects regression model
m2_3_rp6 | merlin no depvar 17 Fixed effects regression model
m2_3_rp7 | merlin no depvar 18 Fixed effects regression model
m2_3_rp8 | merlin no depvar 19 Fixed effects regression model
m2_3_rp9 | merlin no depvar 20 Fixed effects regression model
m2_3_rp10 | merlin no depvar 21 Fixed effects regression model
m2_4_exp | merlin no depvar 12 Fixed effects regression model
m2_4_weib | merlin no depvar 13 Fixed effects regression model
m2_4_gom | merlin no depvar 13 Fixed effects regression model
m2_4_logl | merlin no depvar 13 Fixed effects regression model
m2_4_logn | merlin no depvar 13 Fixed effects regression model
m2_4_ggam | merlin no depvar 14 Fixed effects regression model
m2_4_rp2 | merlin no depvar 14 Fixed effects regression model
m2_4_rp3 | merlin no depvar 15 Fixed effects regression model
m2_4_rp4 | merlin no depvar 16 Fixed effects regression model
m2_4_rp5 | merlin no depvar 17 Fixed effects regression model
m2_4_rp6 | merlin no depvar 18 Fixed effects regression model
m2_4_rp7 | merlin no depvar 19 Fixed effects regression model
m2_4_rp8 | merlin no depvar 20 Fixed effects regression model
m2_4_rp9 | merlin no depvar 21 Fixed effects regression model
m2_4_rp10 | merlin no depvar 22 Fixed effects regression model
m2_5_exp | merlin no depvar 12 Fixed effects regression model
m2_5_weib | merlin no depvar 13 Fixed effects regression model
m2_5_gom | merlin no depvar 13 Fixed effects regression model
m2_5_logl | merlin no depvar 13 Fixed effects regression model
m2_5_logn | merlin no depvar 13 Fixed effects regression model
m2_5_ggam | merlin no depvar 14 Fixed effects regression model
m2_5_rp2 | merlin no depvar 14 Fixed effects regression model
m2_5_rp3 | merlin no depvar 15 Fixed effects regression model
m2_5_rp4 | merlin no depvar 16 Fixed effects regression model
m2_5_rp5 | merlin no depvar 17 Fixed effects regression model
m2_5_rp6 | merlin no depvar 18 Fixed effects regression model
m2_5_rp7 | merlin no depvar 19 Fixed effects regression model
m2_5_rp8 | merlin no depvar 20 Fixed effects regression model
m2_5_rp9 | merlin no depvar 21 Fixed effects regression model
m2_5_rp10 | merlin no depvar 22 Fixed effects regression model
-----------------------------------------------------------------------------
. predictms, transmatrix(mat_four_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logl m2_5_ggam) ///
> timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n(100000) prob diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(1) ci reset bootstrap m(500) /// * es una forma de hacer bootstrap
> at1(tipo_de_programa_2 1 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0) ///
> at2(tipo_de_programa_2 2 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0)
Calculating CIs via parametric bootstrap for at1() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
..........................probabilities must be in [0,1]
rdiscrete(): 3300 argument out of range
predictms_get_newstate(): - function returned error
predictms_sim_cr(): - function returned error
predictms_sim(): - function returned error
predictms_core(): - function returned error
predictms(): - function returned error
<istmt>: - function returned error
r(3300);
.
. drop _prob_at1_1_*_lci _prob_at1_1_*_uci _prob_at2_1_*_lci _prob_at2_1_*_uci
variable _prob_at1_1_*_lci not found
r(111);
. drop _prob_at1_1_* _prob_at2_1_*
variable _prob_at1_1_* not found
r(111);
. rename (_ratio_prob_at2_1_*_lci) (fratio_from1_*_lci_0)
variable _ratio_prob_at2_1_*_lci not found
r(111);
. rename (_ratio_prob_at2_1_*_uci) (fratio_from1_*_uci_0)
variable _ratio_prob_at2_1_*_uci not found
r(111);
. rename (_ratio_prob_at2_1_*) (fratio_from1_*_0)
variable _ratio_prob_at2_1_* not found
r(111);
. rename (_diff_prob_at2_1_*_lci) (fdiff_prob_from1_*_lci_0)
variable _diff_prob_at2_1_*_lci not found
r(111);
. rename (_diff_prob_at2_1_*_uci) (fdiff_prob_from1_*_uci_0)
variable _diff_prob_at2_1_*_uci not found
r(111);
. rename (_diff_prob_at2_1_*) (fdiff_prob_from1_*_0)
variable _diff_prob_at2_1_* not found
r(111);
.
. estread "${pathdata2}parmodels_m2_ago_c.sters"
-----------------------------------------------------------------------------
name | command depvar npar title
-------------+---------------------------------------------------------------
m2_1_exp | merlin no depvar 11 Fixed effects regression model
m2_1_weib | merlin no depvar 12 Fixed effects regression model
m2_1_gom | merlin no depvar 12 Fixed effects regression model
m2_1_logl | merlin no depvar 12 Fixed effects regression model
m2_1_logn | merlin no depvar 12 Fixed effects regression model
m2_1_ggam | merlin no depvar 13 Fixed effects regression model
m2_1_rp2 | merlin no depvar 13 Fixed effects regression model
m2_1_rp3 | merlin no depvar 14 Fixed effects regression model
m2_1_rp4 | merlin no depvar 15 Fixed effects regression model
m2_1_rp5 | merlin no depvar 16 Fixed effects regression model
m2_1_rp6 | merlin no depvar 17 Fixed effects regression model
m2_1_rp7 | merlin no depvar 18 Fixed effects regression model
m2_1_rp8 | merlin no depvar 19 Fixed effects regression model
m2_1_rp9 | merlin no depvar 20 Fixed effects regression model
m2_1_rp10 | merlin no depvar 21 Fixed effects regression model
m2_2_exp | merlin no depvar 11 Fixed effects regression model
m2_2_weib | merlin no depvar 12 Fixed effects regression model
m2_2_gom | merlin no depvar 12 Fixed effects regression model
m2_2_logl | merlin no depvar 12 Fixed effects regression model
m2_2_logn | merlin no depvar 12 Fixed effects regression model
m2_2_ggam | merlin no depvar 13 Fixed effects regression model
m2_2_rp2 | merlin no depvar 13 Fixed effects regression model
m2_2_rp3 | merlin no depvar 14 Fixed effects regression model
m2_2_rp4 | merlin no depvar 15 Fixed effects regression model
m2_2_rp5 | merlin no depvar 16 Fixed effects regression model
m2_2_rp6 | merlin no depvar 17 Fixed effects regression model
m2_2_rp7 | merlin no depvar 18 Fixed effects regression model
m2_2_rp8 | merlin no depvar 19 Fixed effects regression model
m2_2_rp9 | merlin no depvar 20 Fixed effects regression model
m2_2_rp10 | merlin no depvar 21 Fixed effects regression model
m2_3_exp | merlin no depvar 11 Fixed effects regression model
m2_3_weib | merlin no depvar 12 Fixed effects regression model
m2_3_gom | merlin no depvar 12 Fixed effects regression model
m2_3_logl | merlin no depvar 12 Fixed effects regression model
m2_3_logn | merlin no depvar 12 Fixed effects regression model
m2_3_ggam | merlin no depvar 13 Fixed effects regression model
m2_3_rp2 | merlin no depvar 13 Fixed effects regression model
m2_3_rp3 | merlin no depvar 14 Fixed effects regression model
m2_3_rp4 | merlin no depvar 15 Fixed effects regression model
m2_3_rp5 | merlin no depvar 16 Fixed effects regression model
m2_3_rp6 | merlin no depvar 17 Fixed effects regression model
m2_3_rp7 | merlin no depvar 18 Fixed effects regression model
m2_3_rp8 | merlin no depvar 19 Fixed effects regression model
m2_3_rp9 | merlin no depvar 20 Fixed effects regression model
m2_3_rp10 | merlin no depvar 21 Fixed effects regression model
m2_4_exp | merlin no depvar 12 Fixed effects regression model
m2_4_weib | merlin no depvar 13 Fixed effects regression model
m2_4_gom | merlin no depvar 13 Fixed effects regression model
m2_4_logl | merlin no depvar 13 Fixed effects regression model
m2_4_logn | merlin no depvar 13 Fixed effects regression model
m2_4_ggam | merlin no depvar 14 Fixed effects regression model
m2_4_rp2 | merlin no depvar 14 Fixed effects regression model
m2_4_rp3 | merlin no depvar 15 Fixed effects regression model
m2_4_rp4 | merlin no depvar 16 Fixed effects regression model
m2_4_rp5 | merlin no depvar 17 Fixed effects regression model
m2_4_rp6 | merlin no depvar 18 Fixed effects regression model
m2_4_rp7 | merlin no depvar 19 Fixed effects regression model
m2_4_rp8 | merlin no depvar 20 Fixed effects regression model
m2_4_rp9 | merlin no depvar 21 Fixed effects regression model
m2_4_rp10 | merlin no depvar 22 Fixed effects regression model
m2_5_exp | merlin no depvar 12 Fixed effects regression model
m2_5_weib | merlin no depvar 13 Fixed effects regression model
m2_5_gom | merlin no depvar 13 Fixed effects regression model
m2_5_logl | merlin no depvar 13 Fixed effects regression model
m2_5_logn | merlin no depvar 13 Fixed effects regression model
m2_5_ggam | merlin no depvar 14 Fixed effects regression model
m2_5_rp2 | merlin no depvar 14 Fixed effects regression model
m2_5_rp3 | merlin no depvar 15 Fixed effects regression model
m2_5_rp4 | merlin no depvar 16 Fixed effects regression model
m2_5_rp5 | merlin no depvar 17 Fixed effects regression model
m2_5_rp6 | merlin no depvar 18 Fixed effects regression model
m2_5_rp7 | merlin no depvar 19 Fixed effects regression model
m2_5_rp8 | merlin no depvar 20 Fixed effects regression model
m2_5_rp9 | merlin no depvar 21 Fixed effects regression model
m2_5_rp10 | merlin no depvar 22 Fixed effects regression model
-----------------------------------------------------------------------------
.
. predictms, transmatrix(mat_four_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logl m2_5_ggam) ///
> timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n(100000) prob diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(2) ci reset bootstrap m(500) /// * es una forma de hacer bootstrap
> at1(tipo_de_programa_2 1 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0) ///
> at2(tipo_de_programa_2 2 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0)
Calculating CIs via parametric bootstrap for at1() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
Calculating CIs via parametric bootstrap for at2() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
.
. drop _prob_at1_2_*_lci _prob_at1_2_*_uci _prob_at2_2_*_lci _prob_at2_2_*_uci
. drop _prob_at1_2_* _prob_at2_2_*
. rename (_ratio_prob_at2_2_*_lci) (fratio_from2_*_lci_0)
. rename (_ratio_prob_at2_2_*_uci) (fratio_from2_*_uci_0)
. rename (_ratio_prob_at2_2_*) (fratio_from2_*_0)
. rename (_diff_prob_at2_2_*_lci) (fdiff_prob_from2_*_lci_0)
. rename (_diff_prob_at2_2_*_uci) (fdiff_prob_from2_*_uci_0)
. rename (_diff_prob_at2_2_*) (fdiff_prob_from2_*_0)
.
. estread "${pathdata2}parmodels_m2_ago_c.sters"
-----------------------------------------------------------------------------
name | command depvar npar title
-------------+---------------------------------------------------------------
m2_1_exp | merlin no depvar 11 Fixed effects regression model
m2_1_weib | merlin no depvar 12 Fixed effects regression model
m2_1_gom | merlin no depvar 12 Fixed effects regression model
m2_1_logl | merlin no depvar 12 Fixed effects regression model
m2_1_logn | merlin no depvar 12 Fixed effects regression model
m2_1_ggam | merlin no depvar 13 Fixed effects regression model
m2_1_rp2 | merlin no depvar 13 Fixed effects regression model
m2_1_rp3 | merlin no depvar 14 Fixed effects regression model
m2_1_rp4 | merlin no depvar 15 Fixed effects regression model
m2_1_rp5 | merlin no depvar 16 Fixed effects regression model
m2_1_rp6 | merlin no depvar 17 Fixed effects regression model
m2_1_rp7 | merlin no depvar 18 Fixed effects regression model
m2_1_rp8 | merlin no depvar 19 Fixed effects regression model
m2_1_rp9 | merlin no depvar 20 Fixed effects regression model
m2_1_rp10 | merlin no depvar 21 Fixed effects regression model
m2_2_exp | merlin no depvar 11 Fixed effects regression model
m2_2_weib | merlin no depvar 12 Fixed effects regression model
m2_2_gom | merlin no depvar 12 Fixed effects regression model
m2_2_logl | merlin no depvar 12 Fixed effects regression model
m2_2_logn | merlin no depvar 12 Fixed effects regression model
m2_2_ggam | merlin no depvar 13 Fixed effects regression model
m2_2_rp2 | merlin no depvar 13 Fixed effects regression model
m2_2_rp3 | merlin no depvar 14 Fixed effects regression model
m2_2_rp4 | merlin no depvar 15 Fixed effects regression model
m2_2_rp5 | merlin no depvar 16 Fixed effects regression model
m2_2_rp6 | merlin no depvar 17 Fixed effects regression model
m2_2_rp7 | merlin no depvar 18 Fixed effects regression model
m2_2_rp8 | merlin no depvar 19 Fixed effects regression model
m2_2_rp9 | merlin no depvar 20 Fixed effects regression model
m2_2_rp10 | merlin no depvar 21 Fixed effects regression model
m2_3_exp | merlin no depvar 11 Fixed effects regression model
m2_3_weib | merlin no depvar 12 Fixed effects regression model
m2_3_gom | merlin no depvar 12 Fixed effects regression model
m2_3_logl | merlin no depvar 12 Fixed effects regression model
m2_3_logn | merlin no depvar 12 Fixed effects regression model
m2_3_ggam | merlin no depvar 13 Fixed effects regression model
m2_3_rp2 | merlin no depvar 13 Fixed effects regression model
m2_3_rp3 | merlin no depvar 14 Fixed effects regression model
m2_3_rp4 | merlin no depvar 15 Fixed effects regression model
m2_3_rp5 | merlin no depvar 16 Fixed effects regression model
m2_3_rp6 | merlin no depvar 17 Fixed effects regression model
m2_3_rp7 | merlin no depvar 18 Fixed effects regression model
m2_3_rp8 | merlin no depvar 19 Fixed effects regression model
m2_3_rp9 | merlin no depvar 20 Fixed effects regression model
m2_3_rp10 | merlin no depvar 21 Fixed effects regression model
m2_4_exp | merlin no depvar 12 Fixed effects regression model
m2_4_weib | merlin no depvar 13 Fixed effects regression model
m2_4_gom | merlin no depvar 13 Fixed effects regression model
m2_4_logl | merlin no depvar 13 Fixed effects regression model
m2_4_logn | merlin no depvar 13 Fixed effects regression model
m2_4_ggam | merlin no depvar 14 Fixed effects regression model
m2_4_rp2 | merlin no depvar 14 Fixed effects regression model
m2_4_rp3 | merlin no depvar 15 Fixed effects regression model
m2_4_rp4 | merlin no depvar 16 Fixed effects regression model
m2_4_rp5 | merlin no depvar 17 Fixed effects regression model
m2_4_rp6 | merlin no depvar 18 Fixed effects regression model
m2_4_rp7 | merlin no depvar 19 Fixed effects regression model
m2_4_rp8 | merlin no depvar 20 Fixed effects regression model
m2_4_rp9 | merlin no depvar 21 Fixed effects regression model
m2_4_rp10 | merlin no depvar 22 Fixed effects regression model
m2_5_exp | merlin no depvar 12 Fixed effects regression model
m2_5_weib | merlin no depvar 13 Fixed effects regression model
m2_5_gom | merlin no depvar 13 Fixed effects regression model
m2_5_logl | merlin no depvar 13 Fixed effects regression model
m2_5_logn | merlin no depvar 13 Fixed effects regression model
m2_5_ggam | merlin no depvar 14 Fixed effects regression model
m2_5_rp2 | merlin no depvar 14 Fixed effects regression model
m2_5_rp3 | merlin no depvar 15 Fixed effects regression model
m2_5_rp4 | merlin no depvar 16 Fixed effects regression model
m2_5_rp5 | merlin no depvar 17 Fixed effects regression model
m2_5_rp6 | merlin no depvar 18 Fixed effects regression model
m2_5_rp7 | merlin no depvar 19 Fixed effects regression model
m2_5_rp8 | merlin no depvar 20 Fixed effects regression model
m2_5_rp9 | merlin no depvar 21 Fixed effects regression model
m2_5_rp10 | merlin no depvar 22 Fixed effects regression model
-----------------------------------------------------------------------------
.
. predictms, transmatrix(mat_four_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logl m2_5_ggam) ///
> timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n(100000) prob diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(3) ci reset bootstrap m(500) /// * es una forma de hacer bootstrap
> at1(tipo_de_programa_2 1 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0) ///
> at2(tipo_de_programa_2 2 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0)
Calculating CIs via parametric bootstrap for at1() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
Calculating CIs via parametric bootstrap for at2() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
.
. drop _prob_at1_3_*_lci _prob_at1_3_*_uci _prob_at2_3_*_lci _prob_at2_3_*_uci
. drop _prob_at1_3_* _prob_at2_3_*
. rename (_ratio_prob_at2_3_*_lci) (fratio_from3_*_lci_0)
. rename (_ratio_prob_at2_3_*_uci) (fratio_from3_*_uci_0)
. rename (_ratio_prob_at2_3_*) (fratio_from3_*_0)
. rename (_diff_prob_at2_3_*_lci) (fdiff_prob_from3_*_lci_0)
. rename (_diff_prob_at2_3_*_uci) (fdiff_prob_from3_*_uci_0)
. rename (_diff_prob_at2_3_*) (fdiff_prob_from3_*_0)
.
=============================================================================
(RP) Parametric models of differences & ratios- At 0 days
=============================================================================
. estread "${pathdata2}parmodels_m2_ago_c.sters"
-----------------------------------------------------------------------------
name | command depvar npar title
-------------+---------------------------------------------------------------
m2_1_exp | merlin no depvar 11 Fixed effects regression model
m2_1_weib | merlin no depvar 12 Fixed effects regression model
m2_1_gom | merlin no depvar 12 Fixed effects regression model
m2_1_logl | merlin no depvar 12 Fixed effects regression model
m2_1_logn | merlin no depvar 12 Fixed effects regression model
m2_1_ggam | merlin no depvar 13 Fixed effects regression model
m2_1_rp2 | merlin no depvar 13 Fixed effects regression model
m2_1_rp3 | merlin no depvar 14 Fixed effects regression model
m2_1_rp4 | merlin no depvar 15 Fixed effects regression model
m2_1_rp5 | merlin no depvar 16 Fixed effects regression model
m2_1_rp6 | merlin no depvar 17 Fixed effects regression model
m2_1_rp7 | merlin no depvar 18 Fixed effects regression model
m2_1_rp8 | merlin no depvar 19 Fixed effects regression model
m2_1_rp9 | merlin no depvar 20 Fixed effects regression model
m2_1_rp10 | merlin no depvar 21 Fixed effects regression model
m2_2_exp | merlin no depvar 11 Fixed effects regression model
m2_2_weib | merlin no depvar 12 Fixed effects regression model
m2_2_gom | merlin no depvar 12 Fixed effects regression model
m2_2_logl | merlin no depvar 12 Fixed effects regression model
m2_2_logn | merlin no depvar 12 Fixed effects regression model
m2_2_ggam | merlin no depvar 13 Fixed effects regression model
m2_2_rp2 | merlin no depvar 13 Fixed effects regression model
m2_2_rp3 | merlin no depvar 14 Fixed effects regression model
m2_2_rp4 | merlin no depvar 15 Fixed effects regression model
m2_2_rp5 | merlin no depvar 16 Fixed effects regression model
m2_2_rp6 | merlin no depvar 17 Fixed effects regression model
m2_2_rp7 | merlin no depvar 18 Fixed effects regression model
m2_2_rp8 | merlin no depvar 19 Fixed effects regression model
m2_2_rp9 | merlin no depvar 20 Fixed effects regression model
m2_2_rp10 | merlin no depvar 21 Fixed effects regression model
m2_3_exp | merlin no depvar 11 Fixed effects regression model
m2_3_weib | merlin no depvar 12 Fixed effects regression model
m2_3_gom | merlin no depvar 12 Fixed effects regression model
m2_3_logl | merlin no depvar 12 Fixed effects regression model
m2_3_logn | merlin no depvar 12 Fixed effects regression model
m2_3_ggam | merlin no depvar 13 Fixed effects regression model
m2_3_rp2 | merlin no depvar 13 Fixed effects regression model
m2_3_rp3 | merlin no depvar 14 Fixed effects regression model
m2_3_rp4 | merlin no depvar 15 Fixed effects regression model
m2_3_rp5 | merlin no depvar 16 Fixed effects regression model
m2_3_rp6 | merlin no depvar 17 Fixed effects regression model
m2_3_rp7 | merlin no depvar 18 Fixed effects regression model
m2_3_rp8 | merlin no depvar 19 Fixed effects regression model
m2_3_rp9 | merlin no depvar 20 Fixed effects regression model
m2_3_rp10 | merlin no depvar 21 Fixed effects regression model
m2_4_exp | merlin no depvar 12 Fixed effects regression model
m2_4_weib | merlin no depvar 13 Fixed effects regression model
m2_4_gom | merlin no depvar 13 Fixed effects regression model
m2_4_logl | merlin no depvar 13 Fixed effects regression model
m2_4_logn | merlin no depvar 13 Fixed effects regression model
m2_4_ggam | merlin no depvar 14 Fixed effects regression model
m2_4_rp2 | merlin no depvar 14 Fixed effects regression model
m2_4_rp3 | merlin no depvar 15 Fixed effects regression model
m2_4_rp4 | merlin no depvar 16 Fixed effects regression model
m2_4_rp5 | merlin no depvar 17 Fixed effects regression model
m2_4_rp6 | merlin no depvar 18 Fixed effects regression model
m2_4_rp7 | merlin no depvar 19 Fixed effects regression model
m2_4_rp8 | merlin no depvar 20 Fixed effects regression model
m2_4_rp9 | merlin no depvar 21 Fixed effects regression model
m2_4_rp10 | merlin no depvar 22 Fixed effects regression model
m2_5_exp | merlin no depvar 12 Fixed effects regression model
m2_5_weib | merlin no depvar 13 Fixed effects regression model
m2_5_gom | merlin no depvar 13 Fixed effects regression model
m2_5_logl | merlin no depvar 13 Fixed effects regression model
m2_5_logn | merlin no depvar 13 Fixed effects regression model
m2_5_ggam | merlin no depvar 14 Fixed effects regression model
m2_5_rp2 | merlin no depvar 14 Fixed effects regression model
m2_5_rp3 | merlin no depvar 15 Fixed effects regression model
m2_5_rp4 | merlin no depvar 16 Fixed effects regression model
m2_5_rp5 | merlin no depvar 17 Fixed effects regression model
m2_5_rp6 | merlin no depvar 18 Fixed effects regression model
m2_5_rp7 | merlin no depvar 19 Fixed effects regression model
m2_5_rp8 | merlin no depvar 20 Fixed effects regression model
m2_5_rp9 | merlin no depvar 21 Fixed effects regression model
m2_5_rp10 | merlin no depvar 22 Fixed effects regression model
-----------------------------------------------------------------------------
. predictms, transmatrix(mat_four_states) models(m2_1_rp6 m2_2_rp5 m2_3_rp4 m2_4_rp7 m2_5_rp3) ///
> timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n(100000) prob diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(1) ci reset bootstrap m(500) /// * es una forma de hacer bootstrap
> at1(tipo_de_programa_2 1 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0) ///
> at2(tipo_de_programa_2 2 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0)
Calculating CIs via parametric bootstrap for at1() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.............probabilities must be in [0,1]
rdiscrete(): 3300 argument out of range
predictms_get_newstate(): - function returned error
predictms_sim_cr(): - function returned error
predictms_sim(): - function returned error
predictms_core(): - function returned error
predictms(): - function returned error
<istmt>: - function returned error
r(3300);
.
. drop _prob_at1_1_*_lci _prob_at1_1_*_uci _prob_at2_1_*_lci _prob_at2_1_*_uci
variable _prob_at1_1_*_lci not found
r(111);
. drop _prob_at1_1_* _prob_at2_1_*
variable _prob_at1_1_* not found
r(111);
. rename (_ratio_prob_at2_1_*_lci) (fratio_from1_*_lci_rp_0)
variable _ratio_prob_at2_1_*_lci not found
r(111);
. rename (_ratio_prob_at2_1_*_uci) (fratio_from1_*_uci_rp_0)
variable _ratio_prob_at2_1_*_uci not found
r(111);
. rename (_ratio_prob_at2_1_*) (fratio_from1_*_rp_0)
variable _ratio_prob_at2_1_* not found
r(111);
. rename (_diff_prob_at2_1_*_lci) (fdiff_prob_from1_*_lci_rp_0)
variable _diff_prob_at2_1_*_lci not found
r(111);
. rename (_diff_prob_at2_1_*_uci) (fdiff_prob_from1_*_uci_rp_0)
variable _diff_prob_at2_1_*_uci not found
r(111);
. rename (_diff_prob_at2_1_*) (fdiff_prob_from1_*_rp_0)
variable _diff_prob_at2_1_* not found
r(111);
.
. estread "${pathdata2}parmodels_m2_ago_c.sters"
-----------------------------------------------------------------------------
name | command depvar npar title
-------------+---------------------------------------------------------------
m2_1_exp | merlin no depvar 11 Fixed effects regression model
m2_1_weib | merlin no depvar 12 Fixed effects regression model
m2_1_gom | merlin no depvar 12 Fixed effects regression model
m2_1_logl | merlin no depvar 12 Fixed effects regression model
m2_1_logn | merlin no depvar 12 Fixed effects regression model
m2_1_ggam | merlin no depvar 13 Fixed effects regression model
m2_1_rp2 | merlin no depvar 13 Fixed effects regression model
m2_1_rp3 | merlin no depvar 14 Fixed effects regression model
m2_1_rp4 | merlin no depvar 15 Fixed effects regression model
m2_1_rp5 | merlin no depvar 16 Fixed effects regression model
m2_1_rp6 | merlin no depvar 17 Fixed effects regression model
m2_1_rp7 | merlin no depvar 18 Fixed effects regression model
m2_1_rp8 | merlin no depvar 19 Fixed effects regression model
m2_1_rp9 | merlin no depvar 20 Fixed effects regression model
m2_1_rp10 | merlin no depvar 21 Fixed effects regression model
m2_2_exp | merlin no depvar 11 Fixed effects regression model
m2_2_weib | merlin no depvar 12 Fixed effects regression model
m2_2_gom | merlin no depvar 12 Fixed effects regression model
m2_2_logl | merlin no depvar 12 Fixed effects regression model
m2_2_logn | merlin no depvar 12 Fixed effects regression model
m2_2_ggam | merlin no depvar 13 Fixed effects regression model
m2_2_rp2 | merlin no depvar 13 Fixed effects regression model
m2_2_rp3 | merlin no depvar 14 Fixed effects regression model
m2_2_rp4 | merlin no depvar 15 Fixed effects regression model
m2_2_rp5 | merlin no depvar 16 Fixed effects regression model
m2_2_rp6 | merlin no depvar 17 Fixed effects regression model
m2_2_rp7 | merlin no depvar 18 Fixed effects regression model
m2_2_rp8 | merlin no depvar 19 Fixed effects regression model
m2_2_rp9 | merlin no depvar 20 Fixed effects regression model
m2_2_rp10 | merlin no depvar 21 Fixed effects regression model
m2_3_exp | merlin no depvar 11 Fixed effects regression model
m2_3_weib | merlin no depvar 12 Fixed effects regression model
m2_3_gom | merlin no depvar 12 Fixed effects regression model
m2_3_logl | merlin no depvar 12 Fixed effects regression model
m2_3_logn | merlin no depvar 12 Fixed effects regression model
m2_3_ggam | merlin no depvar 13 Fixed effects regression model
m2_3_rp2 | merlin no depvar 13 Fixed effects regression model
m2_3_rp3 | merlin no depvar 14 Fixed effects regression model
m2_3_rp4 | merlin no depvar 15 Fixed effects regression model
m2_3_rp5 | merlin no depvar 16 Fixed effects regression model
m2_3_rp6 | merlin no depvar 17 Fixed effects regression model
m2_3_rp7 | merlin no depvar 18 Fixed effects regression model
m2_3_rp8 | merlin no depvar 19 Fixed effects regression model
m2_3_rp9 | merlin no depvar 20 Fixed effects regression model
m2_3_rp10 | merlin no depvar 21 Fixed effects regression model
m2_4_exp | merlin no depvar 12 Fixed effects regression model
m2_4_weib | merlin no depvar 13 Fixed effects regression model
m2_4_gom | merlin no depvar 13 Fixed effects regression model
m2_4_logl | merlin no depvar 13 Fixed effects regression model
m2_4_logn | merlin no depvar 13 Fixed effects regression model
m2_4_ggam | merlin no depvar 14 Fixed effects regression model
m2_4_rp2 | merlin no depvar 14 Fixed effects regression model
m2_4_rp3 | merlin no depvar 15 Fixed effects regression model
m2_4_rp4 | merlin no depvar 16 Fixed effects regression model
m2_4_rp5 | merlin no depvar 17 Fixed effects regression model
m2_4_rp6 | merlin no depvar 18 Fixed effects regression model
m2_4_rp7 | merlin no depvar 19 Fixed effects regression model
m2_4_rp8 | merlin no depvar 20 Fixed effects regression model
m2_4_rp9 | merlin no depvar 21 Fixed effects regression model
m2_4_rp10 | merlin no depvar 22 Fixed effects regression model
m2_5_exp | merlin no depvar 12 Fixed effects regression model
m2_5_weib | merlin no depvar 13 Fixed effects regression model
m2_5_gom | merlin no depvar 13 Fixed effects regression model
m2_5_logl | merlin no depvar 13 Fixed effects regression model
m2_5_logn | merlin no depvar 13 Fixed effects regression model
m2_5_ggam | merlin no depvar 14 Fixed effects regression model
m2_5_rp2 | merlin no depvar 14 Fixed effects regression model
m2_5_rp3 | merlin no depvar 15 Fixed effects regression model
m2_5_rp4 | merlin no depvar 16 Fixed effects regression model
m2_5_rp5 | merlin no depvar 17 Fixed effects regression model
m2_5_rp6 | merlin no depvar 18 Fixed effects regression model
m2_5_rp7 | merlin no depvar 19 Fixed effects regression model
m2_5_rp8 | merlin no depvar 20 Fixed effects regression model
m2_5_rp9 | merlin no depvar 21 Fixed effects regression model
m2_5_rp10 | merlin no depvar 22 Fixed effects regression model
-----------------------------------------------------------------------------
.
. predictms, transmatrix(mat_four_states) models(m2_1_rp6 m2_2_rp5 m2_3_rp4 m2_4_rp7 m2_5_rp3) ///
> timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n(100000) prob diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(2) ci reset bootstrap m(500) /// * es una forma de hacer bootstrap
> at1(tipo_de_programa_2 1 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0) ///
> at2(tipo_de_programa_2 2 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0)
Calculating CIs via parametric bootstrap for at1() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
Calculating CIs via parametric bootstrap for at2() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
.
. drop _prob_at1_2_*_lci _prob_at1_2_*_uci _prob_at2_2_*_lci _prob_at2_2_*_uci
. drop _prob_at1_2_* _prob_at2_2_*
. rename (_ratio_prob_at2_2_*_lci) (fratio_from2_*_lci_rp_0)
. rename (_ratio_prob_at2_2_*_uci) (fratio_from2_*_uci_rp_0)
. rename (_ratio_prob_at2_2_*) (fratio_from2_*_rp_0)
. rename (_diff_prob_at2_2_*_lci) (fdiff_prob_from2_*_lci_rp_0)
. rename (_diff_prob_at2_2_*_uci) (fdiff_prob_from2_*_uci_rp_0)
. rename (_diff_prob_at2_2_*) (fdiff_prob_from2_*_rp_0)
.
. estread "${pathdata2}parmodels_m2_ago_c.sters"
-----------------------------------------------------------------------------
name | command depvar npar title
-------------+---------------------------------------------------------------
m2_1_exp | merlin no depvar 11 Fixed effects regression model
m2_1_weib | merlin no depvar 12 Fixed effects regression model
m2_1_gom | merlin no depvar 12 Fixed effects regression model
m2_1_logl | merlin no depvar 12 Fixed effects regression model
m2_1_logn | merlin no depvar 12 Fixed effects regression model
m2_1_ggam | merlin no depvar 13 Fixed effects regression model
m2_1_rp2 | merlin no depvar 13 Fixed effects regression model
m2_1_rp3 | merlin no depvar 14 Fixed effects regression model
m2_1_rp4 | merlin no depvar 15 Fixed effects regression model
m2_1_rp5 | merlin no depvar 16 Fixed effects regression model
m2_1_rp6 | merlin no depvar 17 Fixed effects regression model
m2_1_rp7 | merlin no depvar 18 Fixed effects regression model
m2_1_rp8 | merlin no depvar 19 Fixed effects regression model
m2_1_rp9 | merlin no depvar 20 Fixed effects regression model
m2_1_rp10 | merlin no depvar 21 Fixed effects regression model
m2_2_exp | merlin no depvar 11 Fixed effects regression model
m2_2_weib | merlin no depvar 12 Fixed effects regression model
m2_2_gom | merlin no depvar 12 Fixed effects regression model
m2_2_logl | merlin no depvar 12 Fixed effects regression model
m2_2_logn | merlin no depvar 12 Fixed effects regression model
m2_2_ggam | merlin no depvar 13 Fixed effects regression model
m2_2_rp2 | merlin no depvar 13 Fixed effects regression model
m2_2_rp3 | merlin no depvar 14 Fixed effects regression model
m2_2_rp4 | merlin no depvar 15 Fixed effects regression model
m2_2_rp5 | merlin no depvar 16 Fixed effects regression model
m2_2_rp6 | merlin no depvar 17 Fixed effects regression model
m2_2_rp7 | merlin no depvar 18 Fixed effects regression model
m2_2_rp8 | merlin no depvar 19 Fixed effects regression model
m2_2_rp9 | merlin no depvar 20 Fixed effects regression model
m2_2_rp10 | merlin no depvar 21 Fixed effects regression model
m2_3_exp | merlin no depvar 11 Fixed effects regression model
m2_3_weib | merlin no depvar 12 Fixed effects regression model
m2_3_gom | merlin no depvar 12 Fixed effects regression model
m2_3_logl | merlin no depvar 12 Fixed effects regression model
m2_3_logn | merlin no depvar 12 Fixed effects regression model
m2_3_ggam | merlin no depvar 13 Fixed effects regression model
m2_3_rp2 | merlin no depvar 13 Fixed effects regression model
m2_3_rp3 | merlin no depvar 14 Fixed effects regression model
m2_3_rp4 | merlin no depvar 15 Fixed effects regression model
m2_3_rp5 | merlin no depvar 16 Fixed effects regression model
m2_3_rp6 | merlin no depvar 17 Fixed effects regression model
m2_3_rp7 | merlin no depvar 18 Fixed effects regression model
m2_3_rp8 | merlin no depvar 19 Fixed effects regression model
m2_3_rp9 | merlin no depvar 20 Fixed effects regression model
m2_3_rp10 | merlin no depvar 21 Fixed effects regression model
m2_4_exp | merlin no depvar 12 Fixed effects regression model
m2_4_weib | merlin no depvar 13 Fixed effects regression model
m2_4_gom | merlin no depvar 13 Fixed effects regression model
m2_4_logl | merlin no depvar 13 Fixed effects regression model
m2_4_logn | merlin no depvar 13 Fixed effects regression model
m2_4_ggam | merlin no depvar 14 Fixed effects regression model
m2_4_rp2 | merlin no depvar 14 Fixed effects regression model
m2_4_rp3 | merlin no depvar 15 Fixed effects regression model
m2_4_rp4 | merlin no depvar 16 Fixed effects regression model
m2_4_rp5 | merlin no depvar 17 Fixed effects regression model
m2_4_rp6 | merlin no depvar 18 Fixed effects regression model
m2_4_rp7 | merlin no depvar 19 Fixed effects regression model
m2_4_rp8 | merlin no depvar 20 Fixed effects regression model
m2_4_rp9 | merlin no depvar 21 Fixed effects regression model
m2_4_rp10 | merlin no depvar 22 Fixed effects regression model
m2_5_exp | merlin no depvar 12 Fixed effects regression model
m2_5_weib | merlin no depvar 13 Fixed effects regression model
m2_5_gom | merlin no depvar 13 Fixed effects regression model
m2_5_logl | merlin no depvar 13 Fixed effects regression model
m2_5_logn | merlin no depvar 13 Fixed effects regression model
m2_5_ggam | merlin no depvar 14 Fixed effects regression model
m2_5_rp2 | merlin no depvar 14 Fixed effects regression model
m2_5_rp3 | merlin no depvar 15 Fixed effects regression model
m2_5_rp4 | merlin no depvar 16 Fixed effects regression model
m2_5_rp5 | merlin no depvar 17 Fixed effects regression model
m2_5_rp6 | merlin no depvar 18 Fixed effects regression model
m2_5_rp7 | merlin no depvar 19 Fixed effects regression model
m2_5_rp8 | merlin no depvar 20 Fixed effects regression model
m2_5_rp9 | merlin no depvar 21 Fixed effects regression model
m2_5_rp10 | merlin no depvar 22 Fixed effects regression model
-----------------------------------------------------------------------------
.
. predictms, transmatrix(mat_four_states) models(m2_1_rp6 m2_2_rp5 m2_3_rp4 m2_4_rp7 m2_5_rp3) ///
> timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n(100000) prob diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(3) ci reset bootstrap m(500) /// * es una forma de hacer bootstrap
> at1(tipo_de_programa_2 1 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0) ///
> at2(tipo_de_programa_2 2 edad_al_ing_grupos 2 escolaridad_rec 2 sus_principal_mod 1 freq_cons_sus_prin 5 compromiso_biopsicosocial 2 ten
> encia_de_la_vivienda_mod 5 num_otras_sus_mod 2 numero_de_hijos_mod_rec 2 tipo_de_plan_res 2 _start 0)
Calculating CIs via parametric bootstrap for at1() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
Calculating CIs via parametric bootstrap for at2() (500)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.................................................. 50
.................................................. 100
.................................................. 150
.................................................. 200
.................................................. 250
.................................................. 300
.................................................. 350
.................................................. 400
.................................................. 450
.................................................. 500
.
. drop _prob_at1_3_*_lci _prob_at1_3_*_uci _prob_at2_3_*_lci _prob_at2_3_*_uci
. drop _prob_at1_3_* _prob_at2_3_*
. rename (_ratio_prob_at2_3_*_lci) (fratio_from3_*_lci_rp_0)
. rename (_ratio_prob_at2_3_*_uci) (fratio_from3_*_uci_rp_0)
. rename (_ratio_prob_at2_3_*) (fratio_from3_*_rp_0)
. rename (_diff_prob_at2_3_*_lci) (fdiff_prob_from3_*_lci_rp_0)
. rename (_diff_prob_at2_3_*_uci) (fdiff_prob_from3_*_uci_rp_0)
. rename (_diff_prob_at2_3_*) (fdiff_prob_from3_*_rp_0)
.
=============================================================================
Saved at= 16:18:20 26 Jul 2021