Date created: 18:32:55 4 Nov 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: 4 Nov 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: 4 Nov 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: 4 Nov 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_multiestado1_oct_22_corr3_3y.dta
variable _t_scaled already defined
Defined a transition matrix.
. cap gen timevar01 = timevar0/365.25
. cap gen _t_scaled= _t/365.25
. *a= Ambulatory . *b= residential . . *.24640657 =90 days . *.48537963 = 6 months . *.98987835 = 1 year . *3.0078732 = 3 years . *4.9993157 = 5 years
. foreach var of varlist fprob_from1a_2_3y fprob_from1a_2_lci_3y fprob_from1a_2_uci_3y fprob_from1b_2_3y fprob_from1b_2_lci_3y fprob_from1b
> _2_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs12 = (ests_fprob_from1a_2_3y, ests_fprob_from1a_2_lci_3y, ests_fprob_from1a_2_uci_3y, ests_fprob_from1b_2_3y, ests_fprob_fro
> m1b_2_lci_3y, ests_fprob_from1b_2_uci_3y)
. matrix colnames est_bs12 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs12) using "${pathdata2}pr_1st_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_1st_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs12 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 6.3 | 6 | 6.6 | 9.7 | 9.3 | 10.1 |
| 1_yr | 11.7 | 11.2 | 12.3 | 17.6 | 16.9 | 18.4 |
| 3_yrs | 15.8 | 15 | 16.5 | 22.1 | 21 | 23.3 |
| 5_yrs | 16.4 | 15.6 | 17.2 | 22 | 20.8 | 23.3 |
. foreach var of varlist fprob_from1a_3_3y fprob_from1a_3_lci_3y fprob_from1a_3_uci_3y fprob_from1b_3_3y fprob_from1b_3_lci_3y fprob_from1b
> _3_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs22 = (ests_fprob_from1a_3_3y, ests_fprob_from1a_3_lci_3y, ests_fprob_from1a_3_uci_3y, ests_fprob_from1b_3_3y, ests_fprob_fro
> m1b_3_lci_3y, ests_fprob_from1b_3_uci_3y)
. matrix colnames est_bs22 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs22) using "${pathdata2}pr_2nd_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_2nd_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs22 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 31.4 | 30.7 | 32.2 | 24.1 | 23.3 | 25 |
| 1_yr | 47.6 | 46.3 | 48.9 | 34.8 | 33.3 | 36.4 |
| 3_yrs | 47.2 | 45.4 | 49 | 32.4 | 30.5 | 34.4 |
| 5_yrs | 43.6 | 41.6 | 45.6 | 29 | 27.1 | 31.1 |
. foreach var of varlist fprob_from2a_4_3y fprob_from2a_4_lci_3y fprob_from2a_4_uci_3y fprob_from2b_4_3y fprob_from2b_4_lci_3y fprob_from2b
> _4_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs32 = (ests_fprob_from2a_4_3y, ests_fprob_from2a_4_lci_3y, ests_fprob_from2a_4_uci_3y, ests_fprob_from2b_4_3y, ests_fprob_fro
> m2b_4_lci_3y, ests_fprob_from2b_4_uci_3y)
. matrix colnames est_bs32 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs32) using "${pathdata2}pr_3rd_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_3rd_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs32 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 5.4 | 4.5 | 6.5 | 11.1 | 9.5 | 12.9 |
| 1_yr | 9.9 | 8.5 | 11.6 | 18.2 | 16 | 20.5 |
| 3_yrs | 15.1 | 13.2 | 17.3 | 24 | 21.7 | 26.6 |
| 5_yrs | 17 | 14.9 | 19.3 | 25.1 | 22.8 | 27.6 |
. foreach var of varlist fprob_from3a_5_3y fprob_from3a_5_lci_3y fprob_from3a_5_uci_3y fprob_from3b_5_3y fprob_from3b_5_lci_3y fprob_from3b
> _5_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs42 = (ests_fprob_from3a_5_3y, ests_fprob_from3a_5_lci_3y, ests_fprob_from3a_5_uci_3y, ests_fprob_from3b_5_3y, ests_fprob_fro
> m3b_5_lci_3y, ests_fprob_from3b_5_uci_3y)
. matrix colnames est_bs42 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs42) using "${pathdata2}pr_4th_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_4th_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs42 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 9.8 | 8.3 | 11.4 | 17.8 | 15.5 | 20.3 |
| 1_yr | 16.1 | 14.2 | 18.3 | 25.9 | 23.2 | 28.8 |
| 3_yrs | 21.6 | 19.4 | 23.9 | 31 | 28.4 | 33.7 |
| 5_yrs | 22.8 | 20.7 | 25 | 31.5 | 29 | 34 |
. foreach var of varlist fprob_from4a_6_3y fprob_from4a_6_lci_3y fprob_from4a_6_uci_3y fprob_from4b_6_3y fprob_from4b_6_lci_3y fprob_from4b
> _6_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs52 = (ests_fprob_from4a_6_3y, ests_fprob_from4a_6_lci_3y, ests_fprob_from4a_6_uci_3y, ests_fprob_from4b_6_3y, ests_fprob_fro
> m4b_6_lci_3y, ests_fprob_from4b_6_uci_3y)
. matrix colnames est_bs52 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs52) using "${pathdata2}pr_5th_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_5th_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs52 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3 | 2.1 | 4.2 | 3.6 | 2.8 | 4.5 |
| 1_yr | 9.7 | 7.6 | 12.3 | 11.3 | 9.9 | 12.9 |
| 3_yrs | 20.8 | 17.2 | 25 | 22.6 | 20.3 | 25.2 |
| 5_yrs | 25.6 | 21.3 | 30.4 | 27.5 | 24.6 | 30.7 |
. foreach var of varlist fprob_from5a_7_3y fprob_from5a_7_lci_3y fprob_from5a_7_uci_3y fprob_from5b_7_3y fprob_from5b_7_lci_3y fprob_from5b
> _7_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs62 = (ests_fprob_from5a_7_3y, ests_fprob_from5a_7_lci_3y, ests_fprob_from5a_7_uci_3y, ests_fprob_from5b_7_3y, ests_fprob_fro
> m5b_7_lci_3y, ests_fprob_from5b_7_uci_3y)
. matrix colnames est_bs62 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs62) using "${pathdata2}pr_6th_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_6th_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs62 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3.7 | 3.1 | 4.4 | 3.7 | 3.1 | 4.3 |
| 1_yr | 12.4 | 11.1 | 13.7 | 12 | 10.9 | 13.1 |
| 3_yrs | 23.1 | 21.4 | 24.9 | 21.4 | 20 | 22.9 |
| 5_yrs | 25.9 | 24 | 28 | 23.8 | 22.1 | 25.5 |
. foreach var of varlist fprob_from6a_8_3y fprob_from6a_8_lci_3y fprob_from6a_8_uci_3y fprob_from6b_8_3y fprob_from6b_8_lci_3y fprob_from6b
> _8_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs72 = (ests_fprob_from6a_8_3y, ests_fprob_from6a_8_lci_3y, ests_fprob_from6a_8_uci_3y, ests_fprob_from6b_8_3y, ests_fprob_fro
> m6b_8_lci_3y, ests_fprob_from6b_8_uci_3y)
. matrix colnames est_bs72 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs72) using "${pathdata2}pr_7th_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_7th_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs72 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3.3 | 1.3 | 8 | 3.3 | 1.8 | 6.1 |
| 1_yr | 15.1 | 9.6 | 22.9 | 15 | 10.8 | 20.5 |
| 3_yrs | 32.2 | 24.6 | 40.7 | 32.1 | 26.4 | 38.4 |
| 5_yrs | 40.2 | 32 | 49 | 40.1 | 33.8 | 46.8 |
. foreach var of varlist fprob_from7a_9_3y fprob_from7a_9_lci_3y fprob_from7a_9_uci_3y fprob_from7b_9_3y fprob_from7b_9_lci_3y fprob_from7b
> _9_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_bs82 = (ests_fprob_from7a_9_3y, ests_fprob_from7a_9_lci_3y, ests_fprob_from7a_9_uci_3y, ests_fprob_from7b_9_3y, ests_fprob_fro
> m7b_9_lci_3y, ests_fprob_from7b_9_uci_3y)
. matrix colnames est_bs82 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs82) using "${pathdata2}pr_8th_transition_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_8th_transition_1st_pat_22_2corr3.html
> )
The transition probabilities are presented here:
| est_bs82 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3.4 | 2.4 | 4.7 | 4.8 | 3.7 | 6.2 |
| 1_yr | 13.9 | 11.4 | 16.9 | 17.7 | 15.2 | 20.6 |
| 3_yrs | 33.7 | 29.8 | 37.8 | 39 | 35.4 | 42.7 |
| 5_yrs | 44.7 | 40.3 | 49.2 | 50.4 | 46.3 | 54.4 |
. foreach var of varlist flos_from1a_1_3y flos_from1a_1_lci_3y flos_from1a_1_uci_3y flos_from1b_1_3y flos_from1b_1_lci_3y flos_from1b_1_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs12 = (ests_flos_from1a_1_3y, ests_flos_from1a_1_lci_3y, ests_flos_from1a_1_uci_3y, ests_flos_from1b_1_3y, ests_flos_from1b_1
> _lci_3y, ests_flos_from1b_1_uci_3y)
. matrix colnames est_cs12 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs12) using "${pathdata2}los_1st_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_1st_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs12 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 91.6 | 91 | 92.2 | 93.2 | 92.6 | 93.7 |
| 1_yr | 198.9 | 196.4 | 201.4 | 205 | 202.4 | 207.6 |
| 3_yrs | 361.8 | 354.3 | 369.4 | 368.4 | 360.8 | 376.3 |
| 5_yrs | 482.2 | 469.9 | 494.8 | 480.4 | 468.2 | 493 |
. foreach var of varlist flos_from2a_2_3y flos_from2a_2_lci_3y flos_from2a_2_uci_3y flos_from2b_2_3y flos_from2b_2_lci_3y flos_from2b_2_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs22 = (ests_flos_from2a_2_3y, ests_flos_from2a_2_lci_3y, ests_flos_from2a_2_uci_3y, ests_flos_from2b_2_3y, ests_flos_from2b_2
> _lci_3y, ests_flos_from2b_2_uci_3y)
. matrix colnames est_cs22 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs22) using "${pathdata2}los_2nd_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_2nd_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs22 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 115.2 | 114.4 | 116 | 110.6 | 109.1 | 112.1 |
| 1_yr | 346.3 | 342.2 | 350.6 | 322.1 | 315.3 | 329.1 |
| 3_yrs | 964.4 | 945.5 | 983.8 | 860.4 | 833 | 888.6 |
| 5_yrs | 1539.3 | 1501.9 | 1577.5 | 1342.2 | 1291.6 | 1394.8 |
. foreach var of varlist flos_from3a_3_3y flos_from3a_3_lci_3y flos_from3a_3_uci_3y flos_from3b_3_3y flos_from3b_3_lci_3y flos_from3b_3_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs32 = (ests_flos_from3a_3_3y, ests_flos_from3a_3_lci_3y, ests_flos_from3a_3_uci_3y, ests_flos_from3b_3_3y, ests_flos_from3b_3
> _lci_3y, ests_flos_from3b_3_uci_3y)
. matrix colnames est_cs32 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs32) using "${pathdata2}los_3rd_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_3rd_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs32 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 111.7 | 110.3 | 113 | 104.7 | 102.4 | 106.9 |
| 1_yr | 327.8 | 321.4 | 334.3 | 296.2 | 286.8 | 305.9 |
| 3_yrs | 885.1 | 859.8 | 911.2 | 765.8 | 732.8 | 800.4 |
| 5_yrs | 1390.6 | 1343.9 | 1438.9 | 1177 | 1119 | 1237.9 |
. foreach var of varlist flos_from4a_4_3y flos_from4a_4_lci_3y flos_from4a_4_uci_3y flos_from4b_4_3y flos_from4b_4_lci_3y flos_from4b_4_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs42 = (ests_flos_from4a_4_3y, ests_flos_from4a_4_lci_3y, ests_flos_from4a_4_uci_3y, ests_flos_from4b_4_3y, ests_flos_from4b_4
> _lci_3y, ests_flos_from4b_4_uci_3y)
. matrix colnames est_cs42 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs42) using "${pathdata2}los_4th_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_4th_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs42 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 117.6 | 117 | 118.2 | 117.3 | 116.7 | 117.8 |
| 1_yr | 352.6 | 347.5 | 357.9 | 349.2 | 345.4 | 353 |
| 3_yrs | 942.5 | 912.2 | 973.9 | 922 | 903.3 | 941 |
| 5_yrs | 1441.7 | 1378 | 1508.4 | 1398.7 | 1360.6 | 1437.9 |
. foreach var of varlist flos_from5a_5_3y flos_from5a_5_lci_3y flos_from5a_5_uci_3y flos_from5b_5_3y flos_from5b_5_lci_3y flos_from5b_5_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs52 = (ests_flos_from5a_5_3y, ests_flos_from5a_5_lci_3y, ests_flos_from5a_5_uci_3y, ests_flos_from5b_5_3y, ests_flos_from5b_5
> _lci_3y, ests_flos_from5b_5_uci_3y)
. matrix colnames est_cs52 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs52) using "${pathdata2}los_5th_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_5th_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs52 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 117.3 | 116.9 | 117.7 | 117.4 | 117 | 117.7 |
| 1_yr | 347.5 | 344.4 | 350.5 | 347.6 | 344.9 | 350.4 |
| 3_yrs | 912.6 | 897.3 | 928.2 | 915.3 | 901.8 | 929.1 |
| 5_yrs | 1392.1 | 1362.1 | 1422.7 | 1396.7 | 1370.2 | 1423.7 |
. foreach var of varlist flos_from6a_6_3y flos_from6a_6_lci_3y flos_from6a_6_uci_3y flos_from6b_6_3y flos_from6b_6_lci_3y flos_from6b_6_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs62 = (ests_flos_from6a_6_3y, ests_flos_from6a_6_lci_3y, ests_flos_from6a_6_uci_3y, ests_flos_from6b_6_3y, ests_flos_from6b_6
> _lci_3y, ests_flos_from6b_6_uci_3y)
. matrix colnames est_cs62 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs62) using "${pathdata2}los_6th_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_6th_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs62 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 117.9 | 116.3 | 119.4 | 117.9 | 116.8 | 119 |
| 1_yr | 346 | 331.5 | 361 | 345.9 | 335.1 | 357 |
| 3_yrs | 891.4 | 825 | 963.1 | 892.4 | 842.7 | 945 |
| 5_yrs | 1353.2 | 1231.3 | 1487.2 | 1354.6 | 1263.7 | 1452.1 |
. foreach var of varlist flos_from7a_7_3y flos_from7a_7_lci_3y flos_from7a_7_uci_3y flos_from7b_7_3y flos_from7b_7_lci_3y flos_from7b_7_uci
> _3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001),.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001),.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001),.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001),.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_cs72 = (ests_flos_from7a_7_3y, ests_flos_from7a_7_lci_3y, ests_flos_from7a_7_uci_3y, ests_flos_from7b_7_3y, ests_flos_from7b_7
> _lci_3y, ests_flos_from7b_7_uci_3y)
. matrix colnames est_cs72 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_cs72) using "${pathdata2}los_7th_actual_st_3y_1st_pat_22_2corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_7th_actual_st_3y_1st_pat_22_2corr3.h
> tml)
The lengths of stay are presented here:
| est_cs72 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 117.6 | 117 | 118.2 | 116.9 | 116.1 | 117.6 |
| 1_yr | 347.6 | 341.9 | 353.3 | 339.8 | 333.9 | 345.9 |
| 3_yrs | 893.1 | 862.7 | 924.5 | 850.8 | 822 | 880.8 |
| 5_yrs | 1332.2 | 1273 | 1394.2 | 1249.5 | 1194.9 | 1306.5 |
. *list fdiff_prob_from1_2_3y fdiff_prob_from1_2_lci_3y fdiff_prob_from1_2_uci_3y if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999
> , 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(timevar01 , 4.999, 5.02)
.
. foreach var of varlist fdiff_prob_from1_2_3y fdiff_prob_from1_2_lci_3y fdiff_prob_from1_2_uci_3y fdiff_prob_from1_3_3y fdiff_prob_from
> 1_3_lci_3y fdiff_prob_from1_3_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_diff11 = (ests_fdiff_prob_from1_2_3y, ests_fdiff_prob_from1_2_lci_3y, ests_fdiff_prob_from1_2_uci_3y, ests_fdiff_prob_from1_3_
> 3y, ests_fdiff_prob_from1_3_lci_3y, ests_fdiff_prob_from1_3_uci_3y)
. matrix colnames est_diff11 = "Completion" "Completion_lo" "Completion_hi" "Non-completion" "Non-completion_lo" "Non-completion_hi"
.
. esttab matrix(est_diff11) using "${pathdata2}diff_probs_res_out_mod_t1_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t1_3y_22_corr3.ht
> ml)
The transition probabilities are presented here:
| est_diff11 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 3.4 | 2.9 | 3.9 | -7.3 | -8.3 | -6.4 |
| 1_yr | 5.9 | 5 | 6.7 | -12.7 | -14 | -11.5 |
| 3_yrs | 6.4 | 5.2 | 7.6 | -14.8 | -16 | -13.5 |
| 5_yrs | 5.6 | 4.3 | 6.9 | -14.6 | -15.8 | -13.3 |
. foreach var of varlist fdiff_prob_from2_4_3y fdiff_prob_from2_4_lci_3y fdiff_prob_from2_4_uci_3y fdiff_prob_from3_5_3y fdiff_prob_from
> 3_5_lci_3y fdiff_prob_from3_5_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_diff21 = (ests_fdiff_prob_from2_4_3y, ests_fdiff_prob_from2_4_lci_3y, ests_fdiff_prob_from2_4_uci_3y, ests_fdiff_prob_from3_5_
> 3y, ests_fdiff_prob_from3_5_lci_3y, ests_fdiff_prob_from3_5_uci_3y)
. matrix colnames est_diff21 = "Completion" "Completion_lo" "Completion_hi" "Non-completion" "Non-completion_lo" "Non-completion_hi"
.
. esttab matrix(est_diff21) using "${pathdata2}diff_probs_res_out_mod_t2_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t2_3y_22_corr3.ht
> ml)
The transition probabilities are presented here:
| est_diff21 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 5.7 | 4.4 | 7 | 8 | 6.7 | 9.3 |
| 1_yr | 8.2 | 6.4 | 10 | 9.8 | 8.3 | 11.2 |
| 3_yrs | 8.9 | 6.8 | 11 | 9.4 | 7.9 | 10.9 |
| 5_yrs | 8.1 | 5.8 | 10.4 | 8.7 | 7.2 | 10.2 |
. foreach var of varlist fdiff_prob_from4_6_3y fdiff_prob_from4_6_lci_3y fdiff_prob_from4_6_uci_3y fdiff_prob_from5_7_3y fdiff_prob_from
> 5_7_lci_3y fdiff_prob_from5_7_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_diff31 = (ests_fdiff_prob_from4_6_3y, ests_fdiff_prob_from4_6_lci_3y, ests_fdiff_prob_from4_6_uci_3y, ests_fdiff_prob_from5_7_
> 3y, ests_fdiff_prob_from5_7_lci_3y, ests_fdiff_prob_from5_7_uci_3y)
. matrix colnames est_diff31 = "Completion" "Completion_lo" "Completion_hi" "Non-completion" "Non-completion_lo" "Non-completion_hi"
.
. esttab matrix(est_diff31) using "${pathdata2}diff_probs_res_out_mod_t3_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t3_3y_22_corr3.ht
> ml)
The transition probabilities are presented here:
| est_diff31 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | .6 | -.5 | 1.7 | 0 | -.8 | .7 |
| 1_yr | 1.6 | -1.1 | 4.3 | -.4 | -2 | 1.2 |
| 3_yrs | 1.8 | -2.5 | 6.1 | -1.6 | -3.8 | .5 |
| 5_yrs | 1.9 | -2.8 | 6.6 | -2.2 | -4.5 | .2 |
. foreach var of varlist fdiff_prob_from6_8_3y fdiff_prob_from6_8_lci_3y fdiff_prob_from6_8_uci_3y fdiff_prob_from7_9_3y fdiff_prob_from
> 7_9_lci_3y fdiff_prob_from7_9_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*100,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*100,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*100,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*100,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
.
. matrix est_diff41 = (ests_fdiff_prob_from6_8_3y, ests_fdiff_prob_from6_8_lci_3y, ests_fdiff_prob_from6_8_uci_3y, ests_fdiff_prob_from7_9_
> 3y, ests_fdiff_prob_from7_9_lci_3y, ests_fdiff_prob_from7_9_uci_3y)
. matrix colnames est_diff41 = "Completion" "Completion_lo" "Completion_hi" "Non-completion" "Non-completion_lo" "Non-completion_hi"
.
. esttab matrix(est_diff41) using "${pathdata2}diff_probs_res_out_mod_t4_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t4_3y_22_corr3.ht
> ml)
The transition probabilities are presented here:
| est_diff41 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 0 | -3.2 | 3.2 | 1.4 | 0 | 2.9 |
| 1_yr | 0 | -6.7 | 6.6 | 3.7 | .2 | 7.3 |
| 3_yrs | -.1 | -7.8 | 7.6 | 5.3 | .3 | 10.4 |
| 5_yrs | -.1 | -7.5 | 7.4 | 5.7 | .5 | 10.9 |
. foreach var of varlist fdiff_los_from1_2_3y fdiff_los_from1_2_lci_3y fdiff_los_from1_2_uci_3y fdiff_los_from1_3_3y fdiff_los_from1_3_lci_
> 3y fdiff_los_from1_3_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*1,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*1,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*1,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*1,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
. matrix est_diff12_l = (ests_fdiff_los_from1_2_3y, ests_fdiff_los_from1_2_lci_3y, ests_fdiff_los_from1_2_uci_3y, ests_fdiff_los_from1_3_3y
> , ests_fdiff_los_from1_3_lci_3y, ests_fdiff_los_from1_3_uci_3y)
. matrix colnames est_diff12_l = "to_Completion" "to_Completion_lo" "to_Completion_hi" "to_Non_Completion" "to_Non_Completion_lo" "to_Non_
> Completion_hi"
.
. esttab matrix(est_diff12_l) using "${pathdata2}diff_los_res_out_mod_t1_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t1_3y_22_corr3.html
> )
The lengths of stay are presented here:
| est_diff12_l | ||||||
| Not_Completed | Not_Completed_lo | Not_Completed_hi | Completed | Completed_lo | Completed_hi | |
| 3_mths | .3 | .2 | .4 | .3 | .2 | .4 |
| 1_yr | 5.1 | 4.2 | 6.1 | 5.7 | 4.7 | 6.6 |
| 3_yrs | 42 | 35.3 | 48.6 | 43.6 | 36.9 | 50.3 |
| 5_yrs | 87.5 | 73.1 | 101.8 | 88.8 | 74.4 | 103.1 |
. foreach var of varlist fdiff_los_from2_2_3y fdiff_los_from2_2_lci_3y fdiff_los_from2_2_uci_3y fdiff_los_from3_3_3y fdiff_los_from3_3_lci_
> 3y fdiff_los_from3_3_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*1,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*1,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*1,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*1,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
. matrix est_diff22_l = (ests_fdiff_los_from2_2_3y, ests_fdiff_los_from2_2_lci_3y, ests_fdiff_los_from2_2_uci_3y, ests_fdiff_los_from3_3_3y
> , ests_fdiff_los_from3_3_lci_3y, ests_fdiff_los_from3_3_uci_3y)
. matrix colnames est_diff22_l = "Completion" "Completion_lo" "Completion_hi" "Non-Completion" "Non-Completion_lo" "Non-Completion_hi"
.
. esttab matrix(est_diff22_l) using "${pathdata2}diff_los_res_out_mod_t2_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t2_3y_22_corr3.html
> )
The lengths of stay are presented here:
| est_diff22_l | ||||||
| Completion | Completion_lo | Completion_hi | Non-Completion | Non-Completion_lo | Non-Completion_hi | |
| 3_mths | -4.6 | -5.7 | -3.5 | -7 | -8.3 | -5.8 |
| 1_yr | -24.2 | -29.6 | -18.8 | -31.6 | -36.5 | -26.7 |
| 3_yrs | -104.1 | -125.7 | -82.4 | -119.3 | -135.6 | -103 |
| 5_yrs | -197.1 | -237.3 | -156.8 | -213.6 | -241.6 | -185.6 |
. foreach var of varlist fdiff_los_from4_4_3y fdiff_los_from4_4_lci_3y fdiff_los_from4_4_uci_3y fdiff_los_from5_5_3y fdiff_los_from5_5_lci_
> 3y fdiff_los_from5_5_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*1,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*1,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*1,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*1,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
. matrix est_diff32_l = (ests_fdiff_los_from4_4_3y, ests_fdiff_los_from4_4_lci_3y, ests_fdiff_los_from4_4_uci_3y, ests_fdiff_los_from5_5_3y
> , ests_fdiff_los_from5_5_lci_3y, ests_fdiff_los_from5_5_uci_3y)
. matrix colnames est_diff32_l = "Readmission" "Readmission_lo" "Readmission_hi" "Readmission_TNC" "Readmission_TNC_lo" "Readmission_TNC_h
> i"
.
. esttab matrix(est_diff32_l) using "${pathdata2}diff_los_res_out_mod_t3_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t3_3y_22_corr3.html
> )
The lengths of stay are presented here:
| est_diff32_l | ||||||
| Readmission | Readmission_lo | Readmission_hi | Readmission_TNC | Readmission_TNC_lo | Readmission_TNC_hi | |
| 3_mths | -.3 | -.9 | .3 | 0 | -.4 | .5 |
| 1_yr | -3.4 | -9.3 | 2.4 | .1 | -3.7 | 4 |
| 3_yrs | -20.6 | -56.8 | 15.7 | 2.7 | -17.6 | 23.1 |
| 5_yrs | -43 | -119.2 | 33.1 | 4.7 | -34.9 | 44.2 |
. foreach var of varlist fdiff_los_from6_6_3y fdiff_los_from6_6_lci_3y fdiff_los_from6_6_uci_3y fdiff_los_from7_7_3y fdiff_los_from7_7_lci_
> 3y fdiff_los_from7_7_uci_3y {
2. scalar variable = "`var'"
3. qui summarize `var' if inrange(timevar01, .32, .34)
4. scalar e3m_`var' = round(round(r(mean),.001)*1,.1)
5. qui summarize `var' if inrange(timevar01, .999, 1.02)
6. scalar e1y_`var' = round(round(r(mean),.001)*1,.1)
7. qui summarize `var' if inrange(timevar01, 2.99, 3.02)
8. scalar e3y_`var' = round(round(r(mean),.001)*1,.1)
9. qui summarize `var' if inrange(timevar01, 4.999, 5.02)
10. scalar e5y_`var' = round(round(r(mean),.001)*1,.1)
11. cap noi matrix ests_`var' = (`=scalar(e3m_`var')'\ `=scalar(e1y_`var')'\ `=scalar(scalar(e3y_`var'))'\ `=scalar(scalar(e5y_`
> var'))')
12. matrix colnames ests_`var' = `var'
13. matrix rownames ests_`var' = 3_mths 1_yr 3_yrs 5_yrs
14. }
. matrix est_diff42_l = (ests_fdiff_los_from6_6_3y, ests_fdiff_los_from6_6_lci_3y, ests_fdiff_los_from6_6_uci_3y, ests_fdiff_los_from7_7_3y
> , ests_fdiff_los_from7_7_lci_3y, ests_fdiff_los_from7_7_uci_3y)
. matrix colnames est_diff42_l = "Readmission2" "Readmission2_lo" "Readmission2_hi" "Readmission2_TNC" "Readmission2_TNC_lo" "Readmission2
> _TNC_hi"
.
. esttab matrix(est_diff42_l) using "${pathdata2}diff_los_res_out_mod_t4_3y_22_corr3.html", replace
(output written to C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t4_3y_22_corr3.html
> )
The lengths of stay are presented here:
| est_diff42_l | ||||||
| Readmission2 | Readmission2_lo | Readmission2_hi | Readmission2_TNC | Readmission2_TNC_lo | Readmission2_TNC_hi | |
| 3_mths | 0 | -1.5 | 1.5 | -.7 | -1.5 | 0 |
| 1_yr | -.1 | -15.1 | 15 | -7.7 | -15.2 | -.3 |
| 3_yrs | 1 | -68.2 | 70.2 | -42.3 | -82.6 | -1.9 |
| 5_yrs | 1.4 | -122.9 | 125.7 | -82.8 | -160.5 | -5 |
. *https://www.statalist.org/forums/forum/general-stata-discussion/general/1594605-twoway-graph-serset
.
.
. *####################################################################
. *# Probabilities (In the next state) (1st pat)
. *###################################################################
.
. local x `" "1a_2_3y" "1a_2_lci_3y" "1a_2_uci_3y" "1b_2_3y" "1b_2_lci_3y" "1b_2_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(start) ce
> ll(A1)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "1a_3_3y" "1a_3_lci_3y" "1a_3_uci_3y" "1b_3_3y" "1b_3_lci_3y" "1b_3_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2a_4_3y" "2a_4_lci_3y" "2a_4_uci_3y" "2b_4_3y" "2b_4_lci_3y" "2b_4_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3a_5_3y" "3a_5_lci_3y" "3a_5_uci_3y" "3b_5_3y" "3b_5_lci_3y" "3b_5_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4a_6_3y" "4a_6_lci_3y" "4a_6_uci_3y" "4b_6_3y" "4b_6_lci_3y" "4b_6_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5a_7_3y" "5a_7_lci_3y" "5a_7_uci_3y" "5b_7_3y" "5b_7_lci_3y" "5b_7_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6a_8_3y" "6a_8_lci_3y" "6a_8_uci_3y" "6b_8_3y" "6b_8_lci_3y" "6b_8_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7a_9_3y" "7a_9_lci_3y" "7a_9_uci_3y" "7b_9_3y" "7b_9_lci_3y" "7b_9_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fprob_from`l'"
3. }
. *di "`wx'"
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_3y", modify) cell(B30) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
.
.
. *####################################################################
. *# Differences in Probabilities (In the next state)
. *###################################################################
.
. local x `""1_2_3y" "1_2_lci_3y" "1_2_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(start)
> cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "1_3_3y" "1_3_lci_3y" "1_3_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2_4_3y" "2_4_lci_3y" "2_4_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3_5_3y" "3_5_lci_3y" "3_5_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4_6_3y" "4_6_lci_3y" "4_6_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5_7_3y" "5_7_lci_3y" "5_7_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6_8_3y" "6_8_lci_3y" "6_8_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7_9_3y" "7_9_lci_3y" "7_9_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_prob_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_3y", modify) cell(B30) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
.
.
. *####################################################################
. *# Length of Stay (In the actual state)
. *###################################################################
.
. local x `" "1a_1_3y" "1a_1_lci_3y" "1a_1_uci_3y" "1b_1_3y" "1b_1_lci_3y" "1b_1_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(s
> tart) cell(A1)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2a_2_3y" "2a_2_lci_3y" "2a_2_uci_3y" "2b_2_3y" "2b_2_lci_3y" "2b_2_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3a_3_3y" "3a_3_lci_3y" "3a_3_uci_3y" "3b_3_3y" "3b_3_lci_3y" "3b_3_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4a_4_3y" "4a_4_lci_3y" "4a_4_uci_3y" "4b_4_3y" "4b_4_lci_3y" "4b_4_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5a_5_3y" "5a_5_lci_3y" "5a_5_uci_3y" "5b_5_3y" "5b_5_lci_3y" "5b_5_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6a_6_3y" "6a_6_lci_3y" "6a_6_uci_3y" "6b_6_3y" "6b_6_lci_3y" "6b_6_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7a_7_3y" "7a_7_lci_3y" "7a_7_uci_3y" "7b_7_3y" "7b_7_lci_3y" "7b_7_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' flos_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_3y", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
.
. *####################################################################
. *# Differences in Length of Stay (In the actual state)
. *###################################################################
.
. local x `""1_1_3y" "1_1_lci_3y" "1_1_uci_3y""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(start)
> cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2_2_3y" "2_2_lci_3y" "2_2_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3_3_3y" "3_3_lci_3y" "3_3_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4_4_3y" "4_4_lci_3y" "4_4_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5_5_3y" "5_5_lci_3y" "5_5_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6_6_3y" "6_6_lci_3y" "6_6_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7_7_3y" "7_7_lci_3y" "7_7_uci_3y" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fdiff_los_from`l'"
3. }
. export excel timevar01 `wx' if inrange(timevar01 , .32, .34)|inrange(timevar01 , .999, 1.02)|inrange(timevar01 , 2.99, 3.02)|inrange(tim
> evar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_3y", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
.
Ended (not saved) at= 18:33:07 4 Nov 2022