Date created: 18:34:45 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_multiestado1b_oct_22_corr3_3y.dta
Defined a transition matrix.
We corrected the first transition due to errors in the probabilities. Given that it showed an error, we decided to calculate the estimates with the parameter `latent’.
.
.
. global sim 1e5 //5e1 1e5
. global boots 1e3 //5e1 2e3
.
. /*
> *Only if the calculation did not work
> qui estread "${pathdata2}parmodels_m2_jun_22_corr3.sters"
>
> predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp5 m2_7_rp3 m2_8_rp3) ///
> timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
> seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
> from(1) ci at1(tipo_de_plan_res_1 0) at2(tipo_de_plan_res_1 1) reset latent bootstrap m($boots) //* es una forma de hacer boots
> trap
> */
> */
.
. cap rename (_prob_at1_1_*lci) (fprob_from1a_*lci_3y_rp)
. cap rename (_prob_at1_1_*uci) (fprob_from1a_*uci_3y_rp)
. cap rename (_prob_at1_1_*) (fprob_from1a_*_3y_rp)
.
. cap rename (_los_at1_1_*lci) (flos_from1a_*lci_3y_rp)
. cap rename (_los_at1_1_*uci) (flos_from1a_*uci_3y_rp)
. cap rename (_los_at1_1_*) (flos_from1a_*_3y_rp)
.
. cap rename (_prob_at2_1_*lci) (fprob_from1b_*lci_3y_rp)
. cap rename (_prob_at2_1_*uci) (fprob_from1b_*uci_3y_rp)
. cap rename (_prob_at2_1_*) (fprob_from1b_*_3y_rp)
.
. cap rename (_los_at2_1_*lci) (flos_from1b_*lci_3y_rp)
. cap rename (_los_at2_1_*uci) (flos_from1b_*uci_3y_rp)
. cap rename (_los_at2_1_*) (flos_from1b_*_3y_rp)
.
. cap rename (_ratio_prob_at2_1_*_lci) (fratio_prob_from1_*_lci_3y_rp)
. cap rename (_ratio_prob_at2_1_*_uci) (fratio_prob_from1_*_uci_3y_rp)
. cap rename (_ratio_prob_at2_1_*) (fratio_prob_from1_*_3y_rp)
. cap rename (_diff_prob_at2_1_*_lci) (fdiff_prob_from1_*_lci_3y_rp)
. cap rename (_diff_prob_at2_1_*_uci) (fdiff_prob_from1_*_uci_3y_rp)
. cap rename (_diff_prob_at2_1_*) (fdiff_prob_from1_*_3y_rp)
.
. cap rename (_ratio_los_at2_1_*_lci) (fratio_los_from1_*_lci_3y_rp)
. cap rename (_ratio_los_at2_1_*_uci) (fratio_los_from1_*_uci_3y_rp)
. cap rename (_ratio_los_at2_1_*) (fratio_los_from1_*_3y_rp)
. cap rename (_diff_los_at2_1_*_lci) (fdiff_los_from1_*_lci_3y_rp)
. cap rename (_diff_los_at2_1_*_uci) (fdiff_los_from1_*_uci_3y_rp)
. cap rename (_diff_los_at2_1_*) (fdiff_los_from1_*_3y_rp)
.
. 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_rp fprob_from1a_2_lci_3y_rp fprob_from1a_2_uci_3y_rp fprob_from1b_2_3y_rp fprob_from1b_2_lci_3y_
> rp fprob_from1b_2_uci_3y_rp {
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_rp, ests_fprob_from1a_2_lci_3y_rp, ests_fprob_from1a_2_uci_3y_rp, ests_fprob_from1b_2_3y_rp, es
> ts_fprob_from1b_2_lci_3y_rp, ests_fprob_from1b_2_uci_3y_rp)
. matrix colnames est_bs12 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs12) using "${pathdata2}pr_1st_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs12 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 1.5 | 1.3 | 1.6 | 2.2 | 2 | 2.4 |
| 1_yr | 11.3 | 10.8 | 11.9 | 18 | 17.3 | 18.7 |
| 3_yrs | 17.2 | 16.3 | 18 | 24 | 23 | 25.1 |
| 5_yrs | 16.2 | 15.3 | 17.1 | 21.8 | 20.7 | 22.9 |
. foreach var of varlist fprob_from1a_3_3y_rp fprob_from1a_3_lci_3y_rp fprob_from1a_3_uci_3y_rp fprob_from1b_3_3y_rp fprob_from1b_3_lci_3y_
> rp fprob_from1b_3_uci_3y_rp {
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_rp, ests_fprob_from1a_3_lci_3y_rp, ests_fprob_from1a_3_uci_3y_rp, ests_fprob_from1b_3_3y_rp, es
> ts_fprob_from1b_3_lci_3y_rp, ests_fprob_from1b_3_uci_3y_rp)
. matrix colnames est_bs22 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs22) using "${pathdata2}pr_2nd_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs22 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 31.8 | 31 | 32.5 | 25.9 | 25.2 | 26.6 |
| 1_yr | 53.9 | 53 | 54.8 | 43.3 | 42.4 | 44.2 |
| 3_yrs | 51.4 | 50.4 | 52.4 | 38.6 | 37.6 | 39.6 |
| 5_yrs | 48.6 | 47.5 | 49.6 | 35.1 | 34.1 | 36.2 |
. foreach var of varlist fprob_from2a_4_3y_rp fprob_from2a_4_lci_3y_rp fprob_from2a_4_uci_3y_rp fprob_from2b_4_3y_rp fprob_from2b_4_lci_3y_
> rp fprob_from2b_4_uci_3y_rp {
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_rp, ests_fprob_from2a_4_lci_3y_rp, ests_fprob_from2a_4_uci_3y_rp, ests_fprob_from2b_4_3y_rp, es
> ts_fprob_from2b_4_lci_3y_rp, ests_fprob_from2b_4_uci_3y_rp)
. matrix colnames est_bs32 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs32) using "${pathdata2}pr_3rd_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs32 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 5.7 | 4.9 | 6.6 | 9.3 | 8.1 | 10.7 |
| 1_yr | 10.7 | 9.3 | 12.3 | 17.2 | 15.2 | 19.5 |
| 3_yrs | 16 | 14 | 18.1 | 23.7 | 21.2 | 26.5 |
| 5_yrs | 18.1 | 15.8 | 20.6 | 26 | 23.3 | 28.9 |
. foreach var of varlist fprob_from3a_5_3y_rp fprob_from3a_5_lci_3y_rp fprob_from3a_5_uci_3y_rp fprob_from3b_5_3y_rp fprob_from3b_5_lci_3y_
> rp fprob_from3b_5_uci_3y_rp {
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_rp, ests_fprob_from3a_5_lci_3y_rp, ests_fprob_from3a_5_uci_3y_rp, ests_fprob_from3b_5_3y_rp, es
> ts_fprob_from3b_5_lci_3y_rp, ests_fprob_from3b_5_uci_3y_rp)
. matrix colnames est_bs42 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs42) using "${pathdata2}pr_4th_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs42 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 9.2 | 8.1 | 10.4 | 14.2 | 12.5 | 16 |
| 1_yr | 17 | 15.1 | 19 | 25.1 | 22.4 | 27.9 |
| 3_yrs | 22.4 | 20.1 | 24.8 | 31.5 | 28.5 | 34.7 |
| 5_yrs | 23.5 | 21.2 | 25.9 | 32.6 | 29.6 | 35.7 |
. foreach var of varlist fprob_from4a_6_3y_rp fprob_from4a_6_lci_3y_rp fprob_from4a_6_uci_3y_rp fprob_from4b_6_3y_rp fprob_from4b_6_lci_3y_
> rp fprob_from4b_6_uci_3y_rp {
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_rp, ests_fprob_from4a_6_lci_3y_rp, ests_fprob_from4a_6_uci_3y_rp, ests_fprob_from4b_6_3y_rp, es
> ts_fprob_from4b_6_lci_3y_rp, ests_fprob_from4b_6_uci_3y_rp)
. matrix colnames est_bs52 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs52) using "${pathdata2}pr_5th_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs52 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 2.6 | 1.9 | 3.5 | 2.9 | 2.3 | 3.8 |
| 1_yr | 10 | 8 | 12.4 | 11.5 | 9.9 | 13.3 |
| 3_yrs | 22.9 | 18.7 | 27.7 | 25.7 | 23 | 28.6 |
| 5_yrs | 24.7 | 19.6 | 30.5 | 28.2 | 24.9 | 31.7 |
. foreach var of varlist fprob_from5a_7_3y_rp fprob_from5a_7_lci_3y_rp fprob_from5a_7_uci_3y_rp fprob_from5b_7_3y_rp fprob_from5b_7_lci_3y_
> rp fprob_from5b_7_uci_3y_rp {
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_rp, ests_fprob_from5a_7_lci_3y_rp, ests_fprob_from5a_7_uci_3y_rp, ests_fprob_from5b_7_3y_rp, es
> ts_fprob_from5b_7_lci_3y_rp, ests_fprob_from5b_7_uci_3y_rp)
. matrix colnames est_bs62 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs62) using "${pathdata2}pr_6th_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs62 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3 | 2.5 | 3.5 | 3.1 | 2.6 | 3.6 |
| 1_yr | 12.7 | 11.6 | 13.9 | 12.6 | 11.5 | 13.7 |
| 3_yrs | 23.5 | 21.7 | 25.4 | 22.6 | 20.8 | 24.4 |
| 5_yrs | 24.4 | 22.3 | 26.7 | 23.2 | 21.2 | 25.3 |
. foreach var of varlist fprob_from6a_8_3y_rp fprob_from6a_8_lci_3y_rp fprob_from6a_8_uci_3y_rp fprob_from6b_8_3y_rp fprob_from6b_8_lci_3y_
> rp fprob_from6b_8_uci_3y_rp {
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_rp, ests_fprob_from6a_8_lci_3y_rp, ests_fprob_from6a_8_uci_3y_rp, ests_fprob_from6b_8_3y_rp, es
> ts_fprob_from6b_8_lci_3y_rp, ests_fprob_from6b_8_uci_3y_rp)
. matrix colnames est_bs72 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs72) using "${pathdata2}pr_7th_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs72 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 2.7 | 1.5 | 4.9 | 2.4 | 1.4 | 4 |
| 1_yr | 12.4 | 7.8 | 19.1 | 11.3 | 8.1 | 15.6 |
| 3_yrs | 32.6 | 22.4 | 44.7 | 30.1 | 23.5 | 37.7 |
| 5_yrs | 38.5 | 26.9 | 51.5 | 35.8 | 28.4 | 43.9 |
. foreach var of varlist fprob_from7a_9_3y_rp fprob_from7a_9_lci_3y_rp fprob_from7a_9_uci_3y_rp fprob_from7b_9_3y_rp fprob_from7b_9_lci_3y_
> rp fprob_from7b_9_uci_3y_rp {
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_rp, ests_fprob_from7a_9_lci_3y_rp, ests_fprob_from7a_9_uci_3y_rp, ests_fprob_from7b_9_3y_rp, es
> ts_fprob_from7b_9_lci_3y_rp, ests_fprob_from7b_9_uci_3y_rp)
. matrix colnames est_bs82 = Out Out_lci Out_uci Res Res_lci Res_uci
.
. esttab matrix(est_bs82) using "${pathdata2}pr_8th_transition_rp_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_rp_1st_pat_22_2corr3.h
> tml)
The transition probabilities are presented here:
| est_bs82 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 3.1 | 2.3 | 4.2 | 3.5 | 2.6 | 4.7 |
| 1_yr | 15.1 | 12.7 | 17.8 | 16.9 | 14.5 | 19.5 |
| 3_yrs | 34.9 | 30.6 | 39.5 | 38.7 | 34.8 | 42.7 |
| 5_yrs | 44.4 | 39.3 | 49.6 | 49 | 44.6 | 53.4 |
. foreach var of varlist flos_from1a_1_3y_rp flos_from1a_1_lci_3y_rp flos_from1a_1_uci_3y_rp flos_from1b_1_3y_rp flos_from1b_1_lci_3y_rp fl
> os_from1b_1_uci_3y_rp {
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_rp, ests_flos_from1a_1_lci_3y_rp, ests_flos_from1a_1_uci_3y_rp, ests_flos_from1b_1_3y_rp, ests_f
> los_from1b_1_lci_3y_rp, ests_flos_from1b_1_uci_3y_rp)
. 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_rp_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_rp_1st_pat_22_2corr
> 3.html)
The lengths of stay are presented here:
| est_cs12 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 99 | 98.4 | 99.5 | 102 | 101.5 | 102.5 |
| 1_yr | 207.3 | 204.9 | 209.7 | 219.1 | 216.7 | 221.6 |
| 3_yrs | 339.5 | 332.7 | 346.5 | 340.7 | 333.8 | 347.8 |
| 5_yrs | 449.8 | 438.4 | 461.5 | 436.4 | 424.9 | 448.3 |
. foreach var of varlist flos_from2a_2_3y_rp flos_from2a_2_lci_3y_rp flos_from2a_2_uci_3y_rp flos_from2b_2_3y_rp flos_from2b_2_lci_3y_rp fl
> os_from2b_2_uci_3y_rp {
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_rp, ests_flos_from2a_2_lci_3y_rp, ests_flos_from2a_2_uci_3y_rp, ests_flos_from2b_2_3y_rp, ests_f
> los_from2b_2_lci_3y_rp, ests_flos_from2b_2_uci_3y_rp)
. 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_rp_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_rp_1st_pat_22_2corr
> 3.html)
The lengths of stay are presented here:
| est_cs22 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 114.6 | 113.9 | 115.4 | 111.6 | 110.5 | 112.8 |
| 1_yr | 344.6 | 340.8 | 348.4 | 327.3 | 321.7 | 333 |
| 3_yrs | 953.9 | 935.5 | 972.6 | 869.1 | 843.1 | 895.8 |
| 5_yrs | 1516.8 | 1479.2 | 1555.4 | 1344.4 | 1293 | 1397.8 |
. foreach var of varlist flos_from3a_3_3y_rp flos_from3a_3_lci_3y_rp flos_from3a_3_uci_3y_rp flos_from3b_3_3y_rp flos_from3b_3_lci_3y_rp fl
> os_from3b_3_uci_3y_rp {
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_rp, ests_flos_from3a_3_lci_3y_rp, ests_flos_from3a_3_uci_3y_rp, ests_flos_from3b_3_3y_rp, ests_f
> los_from3b_3_lci_3y_rp, ests_flos_from3b_3_uci_3y_rp)
. 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_rp_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_rp_1st_pat_22_2corr
> 3.html)
The lengths of stay are presented here:
| est_cs32 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 111.9 | 110.9 | 112.8 | 107.7 | 106.3 | 109.2 |
| 1_yr | 327.7 | 322.6 | 332.9 | 305.3 | 297.9 | 313 |
| 3_yrs | 875.7 | 851.9 | 900.3 | 773.9 | 740.8 | 808.5 |
| 5_yrs | 1371.3 | 1324.5 | 1419.7 | 1173.6 | 1111.1 | 1239.7 |
. foreach var of varlist flos_from4a_4_3y_rp flos_from4a_4_lci_3y_rp flos_from4a_4_uci_3y_rp flos_from4b_4_3y_rp flos_from4b_4_lci_3y_rp fl
> os_from4b_4_uci_3y_rp {
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_rp, ests_flos_from4a_4_lci_3y_rp, ests_flos_from4a_4_uci_3y_rp, ests_flos_from4b_4_3y_rp, ests_f
> los_from4b_4_lci_3y_rp, ests_flos_from4b_4_uci_3y_rp)
. 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_rp_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_rp_1st_pat_22_2corr
> 3.html)
The lengths of stay are presented here:
| est_cs42 | ||||||
| Out | Out_lci | Out_uci | Res | Res_lci | Res_uci | |
| 3_mths | 117.8 | 117.3 | 118.3 | 117.6 | 117.1 | 118.1 |
| 1_yr | 353.7 | 349.4 | 358 | 351.1 | 347.6 | 354.8 |
| 3_yrs | 930.1 | 899.8 | 961.4 | 910 | 889.5 | 931 |
| 5_yrs | 1422.8 | 1356 | 1492.8 | 1376.5 | 1334 | 1420.4 |
We had to rename some columns due to the extension of the variable name.
. **rename columns
.
. local trans `" "1_2" "1_3" "2_4" "3_5" "4_6" "5_7" "6_8" "7_9" "'
. local names `" "_3y_rp" "_lci_3y_rp" "_uci_3y_rp" "'
. local fdiff `" "fdiff" "'
. foreach opts in `trans' {
2. foreach c in `names' {
3. local fd `" "fd" "'
4. foreach fdiffmod in `fdiff' {
5. gettoken fdmod fd:fd
6. local dbname `fdiffmod'_los_from`opts'`c'
7. local dbnameout `fdmod'_l_from`opts'`c'
8. local new = "From `dbname' to `dbnameout'"
9. di "`new'"
10. cap noi gen `dbnameout'=`dbname'
11. }
12. }
13. }
From fdiff_los_from1_2_3y_rp to fd_l_from1_2_3y_rp
variable fd_l_from1_2_3y_rp already defined
From fdiff_los_from1_2_lci_3y_rp to fd_l_from1_2_lci_3y_rp
variable fd_l_from1_2_lci_3y_rp already defined
From fdiff_los_from1_2_uci_3y_rp to fd_l_from1_2_uci_3y_rp
variable fd_l_from1_2_uci_3y_rp already defined
From fdiff_los_from1_3_3y_rp to fd_l_from1_3_3y_rp
variable fd_l_from1_3_3y_rp already defined
From fdiff_los_from1_3_lci_3y_rp to fd_l_from1_3_lci_3y_rp
variable fd_l_from1_3_lci_3y_rp already defined
From fdiff_los_from1_3_uci_3y_rp to fd_l_from1_3_uci_3y_rp
variable fd_l_from1_3_uci_3y_rp already defined
From fdiff_los_from2_4_3y_rp to fd_l_from2_4_3y_rp
variable fd_l_from2_4_3y_rp already defined
From fdiff_los_from2_4_lci_3y_rp to fd_l_from2_4_lci_3y_rp
variable fd_l_from2_4_lci_3y_rp already defined
From fdiff_los_from2_4_uci_3y_rp to fd_l_from2_4_uci_3y_rp
variable fd_l_from2_4_uci_3y_rp already defined
From fdiff_los_from3_5_3y_rp to fd_l_from3_5_3y_rp
variable fd_l_from3_5_3y_rp already defined
From fdiff_los_from3_5_lci_3y_rp to fd_l_from3_5_lci_3y_rp
variable fd_l_from3_5_lci_3y_rp already defined
From fdiff_los_from3_5_uci_3y_rp to fd_l_from3_5_uci_3y_rp
variable fd_l_from3_5_uci_3y_rp already defined
From fdiff_los_from4_6_3y_rp to fd_l_from4_6_3y_rp
variable fd_l_from4_6_3y_rp already defined
From fdiff_los_from4_6_lci_3y_rp to fd_l_from4_6_lci_3y_rp
variable fd_l_from4_6_lci_3y_rp already defined
From fdiff_los_from4_6_uci_3y_rp to fd_l_from4_6_uci_3y_rp
variable fd_l_from4_6_uci_3y_rp already defined
From fdiff_los_from5_7_3y_rp to fd_l_from5_7_3y_rp
variable fd_l_from5_7_3y_rp already defined
From fdiff_los_from5_7_lci_3y_rp to fd_l_from5_7_lci_3y_rp
variable fd_l_from5_7_lci_3y_rp already defined
From fdiff_los_from5_7_uci_3y_rp to fd_l_from5_7_uci_3y_rp
variable fd_l_from5_7_uci_3y_rp already defined
From fdiff_los_from6_8_3y_rp to fd_l_from6_8_3y_rp
variable fd_l_from6_8_3y_rp already defined
From fdiff_los_from6_8_lci_3y_rp to fd_l_from6_8_lci_3y_rp
variable fd_l_from6_8_lci_3y_rp already defined
From fdiff_los_from6_8_uci_3y_rp to fd_l_from6_8_uci_3y_rp
variable fd_l_from6_8_uci_3y_rp already defined
From fdiff_los_from7_9_3y_rp to fd_l_from7_9_3y_rp
variable fd_l_from7_9_3y_rp already defined
From fdiff_los_from7_9_lci_3y_rp to fd_l_from7_9_lci_3y_rp
variable fd_l_from7_9_lci_3y_rp already defined
From fdiff_los_from7_9_uci_3y_rp to fd_l_from7_9_uci_3y_rp
variable fd_l_from7_9_uci_3y_rp already defined
.
. local trans `" "1_1" "2_2" "3_3" "4_4" "5_5" "6_6" "7_7" "'
. local names `" "_3y_rp" "_lci_3y_rp" "_uci_3y_rp" "'
. local fdiff `" "fdiff" "'
. foreach opts in `trans' {
2. foreach c in `names' {
3. local fd `" "fd" "'
4. foreach fdiffmod in `fdiff' {
5. gettoken fdmod fd:fd
6. local dbname `fdiffmod'_los_from`opts'`c'
7. local dbnameout `fdmod'_l_from`opts'`c'
8. local new = "From `dbname' to `dbnameout'"
9. di "`new'"
10. cap noi gen `dbnameout'=`dbname'
11. }
12. }
13. }
From fdiff_los_from1_1_3y_rp to fd_l_from1_1_3y_rp
variable fd_l_from1_1_3y_rp already defined
From fdiff_los_from1_1_lci_3y_rp to fd_l_from1_1_lci_3y_rp
variable fd_l_from1_1_lci_3y_rp already defined
From fdiff_los_from1_1_uci_3y_rp to fd_l_from1_1_uci_3y_rp
variable fd_l_from1_1_uci_3y_rp already defined
From fdiff_los_from2_2_3y_rp to fd_l_from2_2_3y_rp
variable fd_l_from2_2_3y_rp already defined
From fdiff_los_from2_2_lci_3y_rp to fd_l_from2_2_lci_3y_rp
variable fd_l_from2_2_lci_3y_rp already defined
From fdiff_los_from2_2_uci_3y_rp to fd_l_from2_2_uci_3y_rp
variable fd_l_from2_2_uci_3y_rp already defined
From fdiff_los_from3_3_3y_rp to fd_l_from3_3_3y_rp
variable fd_l_from3_3_3y_rp already defined
From fdiff_los_from3_3_lci_3y_rp to fd_l_from3_3_lci_3y_rp
variable fd_l_from3_3_lci_3y_rp already defined
From fdiff_los_from3_3_uci_3y_rp to fd_l_from3_3_uci_3y_rp
variable fd_l_from3_3_uci_3y_rp already defined
From fdiff_los_from4_4_3y_rp to fd_l_from4_4_3y_rp
variable fd_l_from4_4_3y_rp already defined
From fdiff_los_from4_4_lci_3y_rp to fd_l_from4_4_lci_3y_rp
variable fd_l_from4_4_lci_3y_rp already defined
From fdiff_los_from4_4_uci_3y_rp to fd_l_from4_4_uci_3y_rp
variable fd_l_from4_4_uci_3y_rp already defined
From fdiff_los_from5_5_3y_rp to fd_l_from5_5_3y_rp
variable fd_l_from5_5_3y_rp already defined
From fdiff_los_from5_5_lci_3y_rp to fd_l_from5_5_lci_3y_rp
variable fd_l_from5_5_lci_3y_rp already defined
From fdiff_los_from5_5_uci_3y_rp to fd_l_from5_5_uci_3y_rp
variable fd_l_from5_5_uci_3y_rp already defined
From fdiff_los_from6_6_3y_rp to fd_l_from6_6_3y_rp
variable fd_l_from6_6_3y_rp already defined
From fdiff_los_from6_6_lci_3y_rp to fd_l_from6_6_lci_3y_rp
variable fd_l_from6_6_lci_3y_rp already defined
From fdiff_los_from6_6_uci_3y_rp to fd_l_from6_6_uci_3y_rp
variable fd_l_from6_6_uci_3y_rp already defined
From fdiff_los_from7_7_3y_rp to fd_l_from7_7_3y_rp
variable fd_l_from7_7_3y_rp already defined
From fdiff_los_from7_7_lci_3y_rp to fd_l_from7_7_lci_3y_rp
variable fd_l_from7_7_lci_3y_rp already defined
From fdiff_los_from7_7_uci_3y_rp to fd_l_from7_7_uci_3y_rp
variable fd_l_from7_7_uci_3y_rp already defined
.
. local trans `" "1_2" "1_3" "2_4" "3_5" "4_6" "5_7" "6_8" "7_9" "'
. local names `" "_3y_rp" "_lci_3y_rp" "_uci_3y_rp" "'
. local fdiff `" "fdiff" "'
. foreach opts in `trans' {
2. foreach c in `names' {
3. local fd `" "fd" "'
4. foreach fdiffmod in `fdiff' {
5. gettoken fdmod fd:fd
6. local dbname `fdiffmod'_prob_from`opts'`c'
7. local dbnameout `fdmod'_pr_from`opts'`c'
8. local new = "From `dbname' to `dbnameout'"
9. di "`new'"
10. cap noi gen `dbnameout'=`dbname'
11. }
12. }
13. }
From fdiff_prob_from1_2_3y_rp to fd_pr_from1_2_3y_rp
variable fd_pr_from1_2_3y_rp already defined
From fdiff_prob_from1_2_lci_3y_rp to fd_pr_from1_2_lci_3y_rp
variable fd_pr_from1_2_lci_3y_rp already defined
From fdiff_prob_from1_2_uci_3y_rp to fd_pr_from1_2_uci_3y_rp
variable fd_pr_from1_2_uci_3y_rp already defined
From fdiff_prob_from1_3_3y_rp to fd_pr_from1_3_3y_rp
variable fd_pr_from1_3_3y_rp already defined
From fdiff_prob_from1_3_lci_3y_rp to fd_pr_from1_3_lci_3y_rp
variable fd_pr_from1_3_lci_3y_rp already defined
From fdiff_prob_from1_3_uci_3y_rp to fd_pr_from1_3_uci_3y_rp
variable fd_pr_from1_3_uci_3y_rp already defined
From fdiff_prob_from2_4_3y_rp to fd_pr_from2_4_3y_rp
variable fd_pr_from2_4_3y_rp already defined
From fdiff_prob_from2_4_lci_3y_rp to fd_pr_from2_4_lci_3y_rp
variable fd_pr_from2_4_lci_3y_rp already defined
From fdiff_prob_from2_4_uci_3y_rp to fd_pr_from2_4_uci_3y_rp
variable fd_pr_from2_4_uci_3y_rp already defined
From fdiff_prob_from3_5_3y_rp to fd_pr_from3_5_3y_rp
variable fd_pr_from3_5_3y_rp already defined
From fdiff_prob_from3_5_lci_3y_rp to fd_pr_from3_5_lci_3y_rp
variable fd_pr_from3_5_lci_3y_rp already defined
From fdiff_prob_from3_5_uci_3y_rp to fd_pr_from3_5_uci_3y_rp
variable fd_pr_from3_5_uci_3y_rp already defined
From fdiff_prob_from4_6_3y_rp to fd_pr_from4_6_3y_rp
variable fd_pr_from4_6_3y_rp already defined
From fdiff_prob_from4_6_lci_3y_rp to fd_pr_from4_6_lci_3y_rp
variable fd_pr_from4_6_lci_3y_rp already defined
From fdiff_prob_from4_6_uci_3y_rp to fd_pr_from4_6_uci_3y_rp
variable fd_pr_from4_6_uci_3y_rp already defined
From fdiff_prob_from5_7_3y_rp to fd_pr_from5_7_3y_rp
variable fd_pr_from5_7_3y_rp already defined
From fdiff_prob_from5_7_lci_3y_rp to fd_pr_from5_7_lci_3y_rp
variable fd_pr_from5_7_lci_3y_rp already defined
From fdiff_prob_from5_7_uci_3y_rp to fd_pr_from5_7_uci_3y_rp
variable fd_pr_from5_7_uci_3y_rp already defined
From fdiff_prob_from6_8_3y_rp to fd_pr_from6_8_3y_rp
variable fd_pr_from6_8_3y_rp already defined
From fdiff_prob_from6_8_lci_3y_rp to fd_pr_from6_8_lci_3y_rp
variable fd_pr_from6_8_lci_3y_rp already defined
From fdiff_prob_from6_8_uci_3y_rp to fd_pr_from6_8_uci_3y_rp
variable fd_pr_from6_8_uci_3y_rp already defined
From fdiff_prob_from7_9_3y_rp to fd_pr_from7_9_3y_rp
variable fd_pr_from7_9_3y_rp already defined
From fdiff_prob_from7_9_lci_3y_rp to fd_pr_from7_9_lci_3y_rp
variable fd_pr_from7_9_lci_3y_rp already defined
From fdiff_prob_from7_9_uci_3y_rp to fd_pr_from7_9_uci_3y_rp
variable fd_pr_from7_9_uci_3y_rp already defined
. *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 fd_pr_from1_2_3y_rp fd_pr_from1_2_lci_3y_rp fd_pr_from1_2_uci_3y_rp fd_pr_from1_3_3y_rp fd_pr_from1_3_lci_3y_rp
> fd_pr_from1_3_uci_3y_rp {
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_fd_pr_from1_2_3y_rp, ests_fd_pr_from1_2_lci_3y_rp, ests_fd_pr_from1_2_uci_3y_rp, ests_fd_pr_from1_3_3y_rp, ests
> _fd_pr_from1_3_lci_3y_rp, ests_fd_pr_from1_3_uci_3y_rp)
. 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_rp_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_rp_22_corr3
> .html)
The transition probabilities are presented here:
| est_diff11 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | .8 | .6 | .9 | -5.8 | -6.7 | -4.9 |
| 1_yr | 6.7 | 5.8 | 7.5 | -10.6 | -11.8 | -9.4 |
| 3_yrs | 6.9 | 5.7 | 8.1 | -12.9 | -14.1 | -11.6 |
| 5_yrs | 5.6 | 4.4 | 6.8 | -13.4 | -14.7 | -12.1 |
. foreach var of varlist fd_pr_from2_4_3y_rp fd_pr_from2_4_lci_3y_rp fd_pr_from2_4_uci_3y_rp fd_pr_from3_5_3y_rp fd_pr_from3_5_lci_3y_rp
> fd_pr_from3_5_uci_3y_rp {
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_fd_pr_from2_4_3y_rp, ests_fd_pr_from2_4_lci_3y_rp, ests_fd_pr_from2_4_uci_3y_rp, ests_fd_pr_from3_5_3y_rp, ests
> _fd_pr_from3_5_lci_3y_rp, ests_fd_pr_from3_5_uci_3y_rp)
. 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_rp_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_rp_22_corr3
> .html)
The transition probabilities are presented here:
| est_diff21 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | 3.7 | 2.7 | 4.6 | 5 | 4.1 | 6 |
| 1_yr | 6.5 | 4.8 | 8.2 | 8.1 | 6.7 | 9.6 |
| 3_yrs | 7.8 | 5.4 | 10.1 | 9.1 | 7.5 | 10.8 |
| 5_yrs | 7.9 | 5.3 | 10.6 | 9.1 | 7.3 | 10.9 |
. foreach var of varlist fd_pr_from4_6_3y_rp fd_pr_from4_6_lci_3y_rp fd_pr_from4_6_uci_3y_rp fd_pr_from5_7_3y_rp fd_pr_from5_7_lci_3y_rp
> fd_pr_from5_7_uci_3y_rp {
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_fd_pr_from4_6_3y_rp, ests_fd_pr_from4_6_lci_3y_rp, ests_fd_pr_from4_6_uci_3y_rp, ests_fd_pr_from5_7_3y_rp, ests
> _fd_pr_from5_7_lci_3y_rp, ests_fd_pr_from5_7_uci_3y_rp)
. 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_rp_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_rp_22_corr3
> .html)
The transition probabilities are presented here:
| est_diff31 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | .4 | -.3 | 1 | .1 | -.2 | .5 |
| 1_yr | 1.5 | -.8 | 3.8 | -.1 | -1.3 | 1.2 |
| 3_yrs | 2.8 | -2 | 7.6 | -1 | -3.3 | 1.3 |
| 5_yrs | 3.5 | -2.3 | 9.2 | -1.3 | -4 | 1.5 |
. foreach var of varlist fd_pr_from6_8_3y_rp fd_pr_from6_8_lci_3y_rp fd_pr_from6_8_uci_3y_rp fd_pr_from7_9_3y_rp fd_pr_from7_9_lci_3y_rp
> fd_pr_from7_9_uci_3y_rp {
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_fd_pr_from6_8_3y_rp, ests_fd_pr_from6_8_lci_3y_rp, ests_fd_pr_from6_8_uci_3y_rp, ests_fd_pr_from7_9_3y_rp, ests
> _fd_pr_from7_9_lci_3y_rp, ests_fd_pr_from7_9_uci_3y_rp)
. 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_rp_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_rp_22_corr3
> .html)
The transition probabilities are presented here:
| est_diff41 | ||||||
| Completion | Completion_lo | Completion_hi | Non-completion | Non-completion_lo | Non-completion_hi | |
| 3_mths | -.3 | -1.5 | .9 | .4 | -.3 | 1.1 |
| 1_yr | -1.1 | -5.9 | 3.7 | 1.8 | -1 | 4.6 |
| 3_yrs | -2.5 | -13.1 | 8.3 | 3.8 | -1.6 | 9.2 |
| 5_yrs | -2.7 | -14.7 | 9.3 | 4.6 | -1.7 | 10.9 |
. foreach var of varlist fd_l_from1_2_3y_rp fd_l_from1_2_lci_3y_rp fd_l_from1_2_uci_3y_rp fd_l_from1_3_3y_rp fd_l_from1_3_lci_3y_rp fd_l_fr
> om1_3_uci_3y_rp {
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_fd_l_from1_2_3y_rp, ests_fd_l_from1_2_lci_3y_rp, ests_fd_l_from1_2_uci_3y_rp, ests_fd_l_from1_3_3y_rp, ests_f
> d_l_from1_3_lci_3y_rp, ests_fd_l_from1_3_uci_3y_rp)
. 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_rp_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_rp_22_corr3.h
> tml)
The lengths of stay are presented here:
| est_diff12_l | ||||||
| to_Completion | to_Completion_lo | to_Completion_hi | to_Non_Completion | to_Non_Completion_lo | to_Non_Completion_hi | |
| 3_mths | .3 | .2 | .4 | -3.6 | -4.1 | -3 |
| 1_yr | 9.2 | 8 | 10.5 | -25.6 | -28.8 | -22.3 |
| 3_yrs | 64.5 | 55.4 | 73.6 | -113.3 | -125.2 | -101.5 |
| 5_yrs | 109.6 | 92.1 | 127.2 | -209.5 | -230.4 | -188.5 |
. foreach var of varlist fd_l_from2_2_3y_rp fd_l_from2_2_lci_3y_rp fd_l_from2_2_uci_3y_rp fd_l_from3_3_3y_rp fd_l_from3_3_lci_3y_rp fd_l_fr
> om3_3_uci_3y_rp {
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_fd_l_from2_2_3y_rp, ests_fd_l_from2_2_lci_3y_rp, ests_fd_l_from2_2_uci_3y_rp, ests_fd_l_from3_3_3y_rp, ests_f
> d_l_from3_3_lci_3y_rp, ests_fd_l_from3_3_uci_3y_rp)
. 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_rp_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_rp_22_corr3.h
> tml)
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 | -3 | -3.8 | -2.2 | -4.1 | -4.9 | -3.3 |
| 1_yr | -17.3 | -21.7 | -12.9 | -22.3 | -26.4 | -18.3 |
| 3_yrs | -84.8 | -106.3 | -63.3 | -101.8 | -119.2 | -84.4 |
| 5_yrs | -172.4 | -215.7 | -129.1 | -197.7 | -230.3 | -165.1 |
. foreach var of varlist fd_l_from4_4_3y_rp fd_l_from4_4_lci_3y_rp fd_l_from4_4_uci_3y_rp fd_l_from5_5_3y_rp fd_l_from5_5_lci_3y_rp fd_l_fr
> om5_5_uci_3y_rp {
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_fd_l_from4_4_3y_rp, ests_fd_l_from4_4_lci_3y_rp, ests_fd_l_from4_4_uci_3y_rp, ests_fd_l_from5_5_3y_rp, ests_f
> d_l_from5_5_lci_3y_rp, ests_fd_l_from5_5_uci_3y_rp)
. 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_rp_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_rp_22_corr3.h
> tml)
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 | -.2 | -.5 | .1 | 0 | -.2 | .2 |
| 1_yr | -2.6 | -6.6 | 1.5 | -.5 | -2.7 | 1.7 |
| 3_yrs | -20.1 | -54.9 | 14.7 | .7 | -17.2 | 18.6 |
| 5_yrs | -46.3 | -125.1 | 32.6 | 2.1 | -36.8 | 41.1 |
. foreach var of varlist fd_l_from6_6_3y_rp fd_l_from6_6_lci_3y_rp fd_l_from6_6_uci_3y_rp fd_l_from7_7_3y_rp fd_l_from7_7_lci_3y_rp fd_l_fr
> om7_7_uci_3y_rp {
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_fd_l_from6_6_3y_rp, ests_fd_l_from6_6_lci_3y_rp, ests_fd_l_from6_6_uci_3y_rp, ests_fd_l_from7_7_3y_rp, ests_f
> d_l_from7_7_lci_3y_rp, ests_fd_l_from7_7_uci_3y_rp)
. 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_rp_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_rp_22_corr3.h
> tml)
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 | .1 | -.3 | .6 | -.1 | -.4 | .2 |
| 1_yr | 1.9 | -6.7 | 10.5 | -3 | -7.8 | 1.7 |
| 3_yrs | 16.6 | -52.6 | 85.8 | -24.2 | -61.1 | 12.6 |
| 5_yrs | 35.4 | -117.3 | 188.2 | -55.1 | -135 | 24.7 |
. *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_rp" "1a_2_lci_3y_rp" "1a_2_uci_3y_rp" "1b_2_3y_rp" "1b_2_lci_3y_rp" "1b_2_uci_3y_rp""'
. 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_rp", 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_rp" "1a_3_lci_3y_rp" "1a_3_uci_3y_rp" "1b_3_3y_rp" "1b_3_lci_3y_rp" "1b_3_uci_3y_rp""'
. 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_rp", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2a_4_3y_rp" "2a_4_lci_3y_rp" "2a_4_uci_3y_rp" "2b_4_3y_rp" "2b_4_lci_3y_rp" "2b_4_uci_3y_rp""'
. 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_rp", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3a_5_3y_rp" "3a_5_lci_3y_rp" "3a_5_uci_3y_rp" "3b_5_3y_rp" "3b_5_lci_3y_rp" "3b_5_uci_3y_rp""'
. 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_rp", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4a_6_3y_rp" "4a_6_lci_3y_rp" "4a_6_uci_3y_rp" "4b_6_3y_rp" "4b_6_lci_3y_rp" "4b_6_uci_3y_rp""'
. 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_rp", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5a_7_3y_rp" "5a_7_lci_3y_rp" "5a_7_uci_3y_rp" "5b_7_3y_rp" "5b_7_lci_3y_rp" "5b_7_uci_3y_rp""'
. 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_rp", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6a_8_3y_rp" "6a_8_lci_3y_rp" "6a_8_uci_3y_rp" "6b_8_3y_rp" "6b_8_lci_3y_rp" "6b_8_uci_3y_rp""'
. 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_rp", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7a_9_3y_rp" "7a_9_lci_3y_rp" "7a_9_uci_3y_rp" "7b_9_3y_rp" "7b_9_lci_3y_rp" "7b_9_uci_3y_rp""'
. 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_rp", 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_rp" "1_2_lci_3y_rp" "1_2_uci_3y_rp""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcell(star
> t) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "1_3_3y_rp" "1_3_lci_3y_rp" "1_3_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2_4_3y_rp" "2_4_lci_3y_rp" "2_4_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3_5_3y_rp" "3_5_lci_3y_rp" "3_5_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4_6_3y_rp" "4_6_lci_3y_rp" "4_6_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5_7_3y_rp" "5_7_lci_3y_rp" "5_7_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6_8_3y_rp" "6_8_lci_3y_rp" "6_8_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7_9_3y_rp" "7_9_lci_3y_rp" "7_9_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_pr_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_rp", 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_rp" "1a_1_lci_3y_rp" "1a_1_uci_3y_rp" "1b_1_3y_rp" "1b_1_lci_3y_rp" "1b_1_uci_3y_rp""'
. 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_rp", modify) cell(B1) firstrow(variables) //*save to Excel worksheetcel
> l(start) cell(A1)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "2a_2_3y_rp" "2a_2_lci_3y_rp" "2a_2_uci_3y_rp" "2b_2_3y_rp" "2b_2_lci_3y_rp" "2b_2_uci_3y_rp""'
. 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_rp", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3a_3_3y_rp" "3a_3_lci_3y_rp" "3a_3_uci_3y_rp" "3b_3_3y_rp" "3b_3_lci_3y_rp" "3b_3_uci_3y_rp""'
. 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_rp", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4a_4_3y_rp" "4a_4_lci_3y_rp" "4a_4_uci_3y_rp" "4b_4_3y_rp" "4b_4_lci_3y_rp" "4b_4_uci_3y_rp""'
. 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_rp", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5a_5_3y_rp" "5a_5_lci_3y_rp" "5a_5_uci_3y_rp" "5b_5_3y_rp" "5b_5_lci_3y_rp" "5b_5_uci_3y_rp""'
. 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_rp", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6a_6_3y_rp" "6a_6_lci_3y_rp" "6a_6_uci_3y_rp" "6b_6_3y_rp" "6b_6_lci_3y_rp" "6b_6_uci_3y_rp""'
. 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_rp", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7a_7_3y_rp" "7a_7_lci_3y_rp" "7a_7_uci_3y_rp" "7b_7_3y_rp" "7b_7_lci_3y_rp" "7b_7_uci_3y_rp""'
. 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_rp", 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_rp" "1_1_lci_3y_rp" "1_1_uci_3y_rp""'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", 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_rp" "2_2_lci_3y_rp" "2_2_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", modify) cell(B6) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "3_3_3y_rp" "3_3_lci_3y_rp" "3_3_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", modify) cell(B10) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "4_4_3y_rp" "4_4_lci_3y_rp" "4_4_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", modify) cell(B14) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "5_5_3y_rp" "5_5_lci_3y_rp" "5_5_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", modify) cell(B18) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "6_6_3y_rp" "6_6_lci_3y_rp" "6_6_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", modify) cell(B22) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
. local x `" "7_7_3y_rp" "7_7_lci_3y_rp" "7_7_uci_3y_rp" "'
. local wx ""
. foreach l of local x {
2. local wx "`wx' fd_l_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_rp", modify) cell(B26) //*save to Excel worksheetcell(start) cell(A6)
file trans_prob_los_22_nov.xlsx saved
.
.
Saved at= 18:34:58 4 Nov 2022