Date created: 19:14:34 10 Nov 2022.
Install commands that are unavailable or out of date.
. *<< dd_do : noout > >
. clear all
.
. cap noi which predictms
c:\ado\plus\p\predictms.ado
*! version 4.4.0 22dec2021 MJC
. if _rc==111 {
. cap noi net install multistate, from("https://www.mjcrowther.co.uk/code/multistate")
. }
. cap noi which merlin
c:\ado\plus\m\merlin.ado
*! version 2.1.5 18mar2022 MJC
. if _rc==111 {
. cap noi net install merlin, from("https://www.mjcrowther.co.uk/code/merlin/")
. }
. cap noi which sumat
c:\ado\plus\s\sumat.ado
*! Part of package matrixtools v. 0.28
*! Support: Niels Henrik Bruun, niels.henrik.bruun@gmail.com
*! 2021-01-03 toxl added
. if _rc==111 {
. cap noi scc install matrixtools
. }
. cap noi which estwrite
c:\ado\plus\e\estwrite.ado
*! version 1.2.4 04sep2009
*! version 1.0.1 15may2007 (renamed from -eststo- to -estwrite-; -append- added)
*! version 1.0.0 29apr2005 Ben Jann (ETH Zurich)
. if _rc==111 {
. cap noi ssc install estwrite
. }
. cap noi which winsor2
c:\ado\plus\w\winsor2.ado
*! Inspirit of -winsor-(NJ Cox) and -winsorizeJ-(J Caskey)
*! Lian Yujun, arlionn@163.com, 2013-12-25
*! 1.1 2014.12.16
. if _rc==111 {
. cap noi ssc install winsor2
. }
.
We need to obtain the file and the work folder.
. mata : st_numscalar("OK", direxists("/volumes/sdrive/data//"))
. if scalar(OK) == 1 {
. cap noi cd "/volumes/sdrive/data//"
. global pathdata "/volumes/sdrive/data//"
. di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }
. else display "This file does not exist"
This file does not exist
.
. mata : st_numscalar("OK", direxists("E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\"))
. if scalar(OK) == 1 {
. cap noi cd "E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags
. global pathdata "E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
. global pathdata2 "E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
. di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
Location= E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date: 10 Nov 2022, considering an OS Windows for the user: andre
. }
. else display "This file does not exist"
.
. mata : st_numscalar("OK", direxists("G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\"))
. if scalar(OK) == 1 {
. cap noi cd "G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
. global pathdata "G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
. global pathdata2 "G:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
. di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }
. else display "This file does not exist"
This file does not exist
.
. mata : st_numscalar("OK", direxists("C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\"))
. if scalar(OK) == 1 {
. cap noi cd "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)"
. global pathdata "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
. global pathdata2 "C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
. di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }
. else display "This file does not exist"
This file does not exist
.
. mata : st_numscalar("OK", direxists("C:\Users\andre\Desktop\_mult_state_ags\"))
. if scalar(OK) == 1 {
. cap noi cd "C:\Users\andre\Desktop\_mult_state_ags"
. global pathdata "C:\Users\andre\Desktop\_mult_state_ags"
. global pathdata2 "C:/Users/andre/Desktop/_mult_state_ags/"
. di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }
. else display "This file does not exist"
This file does not exist
.
. mata : st_numscalar("OK", direxists("C:\Users\CISS Fondecyt\OneDrive\Documentos\"))
. if scalar(OK) == 1 {
. cap noi cd "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
. global pathdata "C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags"
. global pathdata2 "C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/"
. di "Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'"
. }
. else display "This file does not exist"
This file does not exist
.
Path data= E:\Mi unidad\Alvacast\SISTRAT 2019 (github)_mult_state_ags;
Timestamp: 10 Nov 2022, considering that is a Windows OS for the username: andre
The file is located and named as: E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/archivo_multiestado1_oct_22_corr3_cc_3y.dta
Defined a transition matrix.
. gen timevar01 = timevar0/365.25 variable timevar01 already defined r(110);
. gen _t_scaled= _t/365.25 variable _t_scaled already defined r(110);
. *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_from1b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_1st_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs12 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 6.1 | 5.8 | 6.5 | 9.2 | 8.8 | 9.7 |
| 1_yr | 11.3 | 10.7 | 11.9 | 16.5 | 15.6 | 17.4 |
| 3_yrs | 14.7 | 13.9 | 15.6 | 19.7 | 18.4 | 21.1 |
| 5_yrs | 14.8 | 13.8 | 15.8 | 18.5 | 16.9 | 20.1 |
. 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_from1b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_2nd_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs22 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 30.7 | 29.7 | 31.7 | 22.7 | 21.7 | 23.8 |
| 1_yr | 44.7 | 42.9 | 46.5 | 31.3 | 29.4 | 33.2 |
| 3_yrs | 41.5 | 39 | 44.1 | 26.6 | 24.3 | 29.1 |
| 5_yrs | 37 | 34.3 | 39.7 | 22.7 | 20.5 | 25.2 |
. 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_from2b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_3rd_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs32 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 7.5 | 6.3 | 8.8 | 13.1 | 11.1 | 15.3 |
| 1_yr | 13.7 | 11.8 | 15.9 | 22 | 19.2 | 25.1 |
| 3_yrs | 21.1 | 18.4 | 24.1 | 30.3 | 27.1 | 33.8 |
| 5_yrs | 24 | 20.9 | 27.4 | 32.7 | 29.3 | 36.3 |
. 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_from3b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_4th_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs42 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 13.3 | 11.2 | 15.8 | 24 | 20.7 | 27.6 |
| 1_yr | 21.6 | 18.8 | 24.7 | 34.1 | 30.4 | 37.9 |
| 3_yrs | 27.6 | 24.8 | 30.7 | 39.4 | 36.1 | 42.8 |
| 5_yrs | 28.8 | 26 | 31.6 | 39 | 36.1 | 42 |
. 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_from4b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_5th_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs52 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 2.1 | 1.3 | 3.5 | 3.2 | 2.3 | 4.3 |
| 1_yr | 9.6 | 7.2 | 12.6 | 12.4 | 10.6 | 14.6 |
| 3_yrs | 20.4 | 16.7 | 24.7 | 23.8 | 21 | 26.8 |
| 5_yrs | 24.1 | 19.7 | 29.1 | 27.5 | 24.1 | 31.2 |
. 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_from5b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_6th_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs62 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 4.1 | 3.4 | 5 | 3.7 | 3 | 4.4 |
| 1_yr | 13.5 | 12 | 15.1 | 12.2 | 11 | 13.6 |
| 3_yrs | 24.7 | 22.7 | 26.8 | 21.7 | 20 | 23.5 |
| 5_yrs | 27.5 | 25.2 | 30 | 24.1 | 22.1 | 26.2 |
. 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_from6b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_7th_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs72 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 4.9 | 2 | 11.4 | 3.8 | 1.8 | 7.6 |
| 1_yr | 18.5 | 11.6 | 28.2 | 16.1 | 11 | 22.9 |
| 3_yrs | 36.2 | 27.1 | 46.4 | 33.6 | 26.9 | 41.1 |
| 5_yrs | 44.5 | 34.6 | 54.9 | 42 | 34.6 | 49.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_from7b_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/pr_8th_transition_1st_pat_22_2corr3_cc.html)
The transition probabilities are presented here:
| est_bs82 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3.1 | 2.1 | 4.6 | 5 | 3.7 | 6.7 |
| 1_yr | 13.5 | 10.9 | 16.8 | 18.5 | 15.6 | 21.8 |
| 3_yrs | 33.7 | 29.4 | 38.3 | 40.8 | 36.7 | 45.1 |
| 5_yrs | 45.2 | 40.2 | 50.2 | 52.4 | 47.7 | 57 |
. 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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_1st_actual_st_3y_1st_pat_22_2corr3_cc.html)
The lengths of stay are presented here:
| est_cs12 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 91.6 | 90.9 | 92.2 | 93.6 | 93 | 94.3 |
| 1_yr | 199.3 | 196.4 | 202.3 | 207.6 | 204.6 | 210.7 |
| 3_yrs | 365.2 | 356.3 | 374.3 | 381 | 371.9 | 390.2 |
| 5_yrs | 489 | 474.5 | 504.1 | 502.5 | 487.9 | 517.6 |
. 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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_2nd_actual_st_3y_1st_pat_22_2corr3_cc.html)
The lengths of stay are presented here:
| est_cs22 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 113.5 | 112.4 | 114.6 | 109.1 | 107.4 | 110.8 |
| 1_yr | 336.8 | 331.4 | 342.4 | 312.9 | 304.4 | 321.7 |
| 3_yrs | 912.5 | 886.9 | 938.7 | 804.1 | 768.3 | 841.6 |
| 5_yrs | 1424.9 | 1373.2 | 1478.5 | 1213.3 | 1145.3 | 1285.3 |
. 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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_3rd_actual_st_3y_1st_pat_22_2corr3_cc.html)
The lengths of stay are presented here:
| est_cs32 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 109 | 107 | 110.9 | 99.6 | 96.4 | 102.8 |
| 1_yr | 312 | 302.8 | 321.4 | 270.8 | 258 | 284.3 |
| 3_yrs | 812.5 | 776.9 | 849.8 | 663.7 | 620.2 | 710.2 |
| 5_yrs | 1248.3 | 1183.7 | 1316.4 | 986.9 | 912.9 | 1066.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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/los_4th_actual_st_3y_1st_pat_22_2corr3_cc.html)
The lengths of stay are presented here:
| est_cs42 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 118.2 | 117.7 | 118.8 | 117.8 | 117.2 | 118.3 |
| 1_yr | 354.6 | 348.8 | 360.4 | 348.9 | 344.4 | 353.4 |
| 3_yrs | 940.2 | 906.8 | 974.8 | 908.7 | 886.5 | 931.4 |
| 5_yrs | 1435.8 | 1368.2 | 1506.6 | 1375.3 | 1331.8 | 1420.1 |
. *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)|inr
> ange(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_from1_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_fdi
> ff_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t1_3y_22_corr3_cc.html)
The transition probabilities are presented here:
| est_diff11 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 3.1 | 2.6 | 3.7 | -8 | -9 | -6.9 |
| 1_yr | 5.2 | 4.2 | 6.2 | -13.4 | -14.7 | -12 |
| 3_yrs | 5 | 3.6 | 6.4 | -14.9 | -16.3 | -13.4 |
| 5_yrs | 3.7 | 2.2 | 5.2 | -14.2 | -15.7 | -12.8 |
. 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_from3_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_fdi
> ff_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t2_3y_22_corr3_cc.html)
The transition probabilities are presented here:
| est_diff21 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 5.6 | 3.9 | 7.3 | 10.6 | 8.9 | 12.4 |
| 1_yr | 8.2 | 5.9 | 10.6 | 12.5 | 10.7 | 14.3 |
| 3_yrs | 9.2 | 6.3 | 12 | 11.8 | 10 | 13.6 |
| 5_yrs | 8.7 | 5.7 | 11.8 | 10.3 | 8.4 | 12.1 |
. 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_from5_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_fdi
> ff_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t3_3y_22_corr3_cc.html)
The transition probabilities are presented here:
| est_diff31 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 1 | -.2 | 2.3 | -.5 | -1.3 | .4 |
| 1_yr | 2.9 | -.3 | 6 | -1.2 | -3.1 | .7 |
| 3_yrs | 3.4 | -1 | 7.6 | -3 | -5.5 | -.4 |
| 5_yrs | 3.4 | -1.2 | 8 | -3.5 | -6.2 | -.7 |
. 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_from7_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_fdi
> ff_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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_probs_res_out_mod_t4_3y_22_corr3_cc.html)
The transition probabilities are presented here:
| est_diff41 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | -1.1 | -5.2 | 3.1 | 1.9 | .4 | 3.5 |
| 1_yr | -2.4 | -10.1 | 5.4 | 5 | 1 | 8.9 |
| 3_yrs | -2.6 | -11.4 | 6.3 | 7.2 | 1.5 | 12.8 |
| 5_yrs | -2.5 | -11.1 | 6.1 | 7.2 | 1.4 | 12.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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t1_3y_22_corr3_cc.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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t2_3y_22_corr3_cc.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.4 | -5.7 | -3 | -9.4 | -11.1 | -7.7 |
| 1_yr | -23.9 | -30.9 | -16.9 | -41.1 | -47.6 | -34.6 |
| 3_yrs | -108.4 | -137.8 | -78.9 | -148.9 | -170 | -127.7 |
| 5_yrs | -211.6 | -267.4 | -155.8 | -261.4 | -297.2 | -225.7 |
. 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_hi"
.
. esttab matrix(est_diff32_l) using "${pathdata2}diff_los_res_out_mod_t3_3y_22_corr3_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t3_3y_22_corr3_cc.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 | -.5 | -1.1 | .1 | .3 | -.2 | .8 |
| 1_yr | -5.7 | -12.5 | 1.2 | 2 | -2.4 | 6.5 |
| 3_yrs | -31.5 | -71.4 | 8.3 | 13 | -10.8 | 36.7 |
| 5_yrs | -60.5 | -139.4 | 18.4 | 25.7 | -20.7 | 72 |
. 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_cc.html", replace
(output written to E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/diff_los_res_out_mod_t4_3y_22_corr3_cc.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 | .5 | -1.7 | 2.7 | -.9 | -1.7 | -.1 |
| 1_yr | 5.2 | -13 | 23.5 | -10.3 | -18.4 | -2.2 |
| 3_yrs | 24.2 | -56.2 | 104.6 | -56.7 | -101.4 | -12 |
| 5_yrs | 43 | -100.7 | 186.7 | -109.1 | -195.3 | -22.8 |
. *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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_3y", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(start) cell(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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_pr_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_prob_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_3y", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(start) 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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("1pat_los_act_st_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_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(timevar01 , 4.999, 5.02) ///
> using trans_prob_los_22_nov.xlsx, sheet("diff_los_corr3_cc_3y", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
.
Ended (not saved) at= 19:14:42 10 Nov 2022