Database (step 2b)

<head>
  <link rel="stylesheet" type="text/css" href="stmarkdown.css">
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/javascript" async
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_CHTML">
</script>
</head>
<head>
  <link rel="stylesheet" type="text/css" href="stmarkdown.css">
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/javascript" async
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_CHTML">
</script>
</head>
<head>
  <link rel="stylesheet" type="text/css" href="stmarkdown.css">
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/javascript" async
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_CHTML">
</script>
</head>
<head>
  <link rel="stylesheet" type="text/css" href="stmarkdown.css">
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/javascript" async
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_CHTML">
</script>
</head>
<head>
  <link rel="stylesheet" type="text/css" href="stmarkdown.css">
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/javascript" async
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_CHTML">
</script>
</head>
<h1><a href="#database-step-2b" id="database-step-2b">Database (step 2b)</a></h1>
<p>Date created:  09:35:54 23 Oct 2022.</p>
<p>Install commands that are unavailable or out of date.</p>
<pre><code>. *&lt;&lt; dd_do : noout &gt; &gt;
. 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(&quot;https://www.mjcrowther.co.uk/code/multistate&quot;) 
.         }

. 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(&quot;https://www.mjcrowther.co.uk/code/merlin/&quot;) 
.         }

. 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
.         }       

.         
</code></pre>
<p>We need to obtain the file and the work folder.</p>
<pre><code>. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;/volumes/sdrive/data//&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;/volumes/sdrive/data//&quot;
.         global pathdata &quot;/volumes/sdrive/data//&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'&quot;
. }

. else display &quot;This file does not exist&quot;
This file does not exist

. 
. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
.         global pathdata &quot;E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
.         global pathdata2 &quot;E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'&quot;
. }

. else display &quot;This file does not exist&quot;
This file does not exist

. 
. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
.         global pathdata &quot;G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
.         global pathdata2 &quot;G:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'&quot;
. }

. else display &quot;This file does not exist&quot;
This file does not exist

.                 
. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)&quot;
C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)
.         global pathdata &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
.         global pathdata2 &quot;C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'&quot;
Location= C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date: 23 Oct 2022, considering an OS Windows for
&gt;  the user: CISS Fondecyt
. }

. else display &quot;This file does not exist&quot;

. 
. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;C:\Users\andre\Desktop\_mult_state_ags\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;C:\Users\andre\Desktop\_mult_state_ags&quot;
.         global pathdata &quot;C:\Users\andre\Desktop\_mult_state_ags&quot;
.         global pathdata2 &quot;C:/Users/andre/Desktop/_mult_state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'&quot;
. }

. else display &quot;This file does not exist&quot;
This file does not exist

. 
. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;C:\Users\CISS Fondecyt\OneDrive\Documentos\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags
.         global pathdata &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
.         global pathdata2 &quot;C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the user: `c(username)'&quot;
Location= C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date: 23 Oct 2022, considering an OS Windows for
&gt;  the user: CISS Fondecyt
. }

. else display &quot;This file does not exist&quot;

. 
</code></pre>
<head>
  <link rel="stylesheet" type="text/css" href="stmarkdown.css">
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/javascript" async
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_CHTML">
</script>
</head>
<p>Path data= C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)_mult_state_ags;</p>
<p>Timestamp: 23 Oct 2022, considering that is a Windows OS for the username: CISS Fondecyt</p>
<p>The file is located and named as: C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/archivo_multiestado0_oct_22_corr3_3y.dta</p>
<p>timevar0 not found</p>
<p>Defined a transition matrix.</p>
<p>=============================================================================</p>
<h2><a href="#transition-probabilities-length-of-stay-parametric-models-rp" id="transition-probabilities-length-of-stay-parametric-models-rp">Transition Probabilities &amp; Length of Stay, Parametric Models (RP)</a></h2>
<p>=============================================================================</p>
<p>(70,463 missing values generated)</p>
<pre><code>. *We corrected the first transition due to errors in the probabilities. Given that it showed an error, we decided to calculate the estimat
&gt; es with the parameter `latent'.
. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         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
&gt; trap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_1_*lci) (fprob_from1a_*lci_3y_rp)

. rename (_prob_at1_1_*uci) (fprob_from1a_*uci_3y_rp)

. rename (_prob_at1_1_*) (fprob_from1a_*_3y_rp)

. 
. rename (_los_at1_1_*lci) (flos_from1a_*lci_3y_rp)

. rename (_los_at1_1_*uci) (flos_from1a_*uci_3y_rp)

. rename (_los_at1_1_*) (flos_from1a_*_3y_rp)

. 
. rename (_prob_at2_1_*lci) (fprob_from1b_*lci_3y_rp)

. rename (_prob_at2_1_*uci) (fprob_from1b_*uci_3y_rp)

. rename (_prob_at2_1_*) (fprob_from1b_*_3y_rp)

. 
. rename (_los_at2_1_*lci) (flos_from1b_*lci_3y_rp)

. rename (_los_at2_1_*uci) (flos_from1b_*uci_3y_rp)

. rename (_los_at2_1_*) (flos_from1b_*_3y_rp)

. 
. rename (_ratio_prob_at2_1_*_lci) (fratio_prob_from1_*_lci_3y_rp)

. rename (_ratio_prob_at2_1_*_uci) (fratio_prob_from1_*_uci_3y_rp)

. rename (_ratio_prob_at2_1_*) (fratio_prob_from1_*_3y_rp)

. rename (_diff_prob_at2_1_*_lci) (fdiff_prob_from1_*_lci_3y_rp)

. rename (_diff_prob_at2_1_*_uci) (fdiff_prob_from1_*_uci_3y_rp)

. rename (_diff_prob_at2_1_*) (fdiff_prob_from1_*_3y_rp) 

. 
. rename (_ratio_los_at2_1_*_lci) (fratio_los_from1_*_lci_3y_rp)

. rename (_ratio_los_at2_1_*_uci) (fratio_los_from1_*_uci_3y_rp)

. rename (_ratio_los_at2_1_*) (fratio_los_from1_*_3y_rp)

. rename (_diff_los_at2_1_*_lci) (fdiff_los_from1_*_lci_3y_rp)

. rename (_diff_los_at2_1_*_uci) (fdiff_los_from1_*_uci_3y_rp)

. rename (_diff_los_at2_1_*) (fdiff_los_from1_*_3y_rp)

. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         from(2) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una for
&gt; ma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_2_*lci) (fprob_from2a_*lci_3y_rp)

. rename (_prob_at1_2_*uci) (fprob_from2a_*uci_3y_rp)

. rename (_prob_at1_2_*) (fprob_from2a_*_3y_rp)

. 
. rename (_los_at1_2_*lci) (flos_from2a_*lci_3y_rp)

. rename (_los_at1_2_*uci) (flos_from2a_*uci_3y_rp)

. rename (_los_at1_2_*) (flos_from2a_*_3y_rp)

. 
. rename (_prob_at2_2_*lci) (fprob_from2b_*lci_3y_rp)

. rename (_prob_at2_2_*uci) (fprob_from2b_*uci_3y_rp)

. rename (_prob_at2_2_*) (fprob_from2b_*_3y_rp)

. 
. rename (_los_at2_2_*lci) (flos_from2b_*lci_3y_rp)

. rename (_los_at2_2_*uci) (flos_from2b_*uci_3y_rp)

. rename (_los_at2_2_*) (flos_from2b_*_3y_rp)

. 
. rename (_ratio_prob_at2_2_*_lci) (fratio_prob_from2_*_lci_3y_rp)

. rename (_ratio_prob_at2_2_*_uci) (fratio_prob_from2_*_uci_3y_rp)

. rename (_ratio_prob_at2_2_*) (fratio_prob_from2_*_3y_rp)

. rename (_diff_prob_at2_2_*_lci) (fdiff_prob_from2_*_lci_3y_rp)

. rename (_diff_prob_at2_2_*_uci) (fdiff_prob_from2_*_uci_3y_rp)

. rename (_diff_prob_at2_2_*) (fdiff_prob_from2_*_3y_rp) 

. 
. rename (_ratio_los_at2_2_*_lci) (fratio_los_from2_*_lci_3y_rp)

. rename (_ratio_los_at2_2_*_uci) (fratio_los_from2_*_uci_3y_rp)

. rename (_ratio_los_at2_2_*) (fratio_los_from2_*_3y_rp)

. rename (_diff_los_at2_2_*_lci) (fdiff_los_from2_*_lci_3y_rp)

. rename (_diff_los_at2_2_*_uci) (fdiff_los_from2_*_uci_3y_rp)

. rename (_diff_los_at2_2_*) (fdiff_los_from2_*_3y_rp)

. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         from(3) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una for
&gt; ma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_3_*lci) (fprob_from3a_*lci_3y_rp)

. rename (_prob_at1_3_*uci) (fprob_from3a_*uci_3y_rp)

. rename (_prob_at1_3_*) (fprob_from3a_*_3y_rp)

. 
. rename (_los_at1_3_*lci) (flos_from3a_*lci_3y_rp)

. rename (_los_at1_3_*uci) (flos_from3a_*uci_3y_rp)

. rename (_los_at1_3_*) (flos_from3a_*_3y_rp)

. 
. rename (_prob_at2_3_*lci) (fprob_from3b_*lci_3y_rp)

. rename (_prob_at2_3_*uci) (fprob_from3b_*uci_3y_rp)

. rename (_prob_at2_3_*) (fprob_from3b_*_3y_rp)

. 
. rename (_los_at2_3_*lci) (flos_from3b_*lci_3y_rp)

. rename (_los_at2_3_*uci) (flos_from3b_*uci_3y_rp)

. rename (_los_at2_3_*) (flos_from3b_*_3y_rp)

. 
. rename (_ratio_prob_at2_3_*_lci) (fratio_prob_from3_*_lci_3y_rp)

. rename (_ratio_prob_at2_3_*_uci) (fratio_prob_from3_*_uci_3y_rp)

. rename (_ratio_prob_at2_3_*) (fratio_prob_from3_*_3y_rp)

. rename (_diff_prob_at2_3_*_lci) (fdiff_prob_from3_*_lci_3y_rp)

. rename (_diff_prob_at2_3_*_uci) (fdiff_prob_from3_*_uci_3y_rp)

. rename (_diff_prob_at2_3_*) (fdiff_prob_from3_*_3y_rp) 

. 
. rename (_ratio_los_at2_3_*_lci) (fratio_los_from3_*_lci_3y_rp)

. rename (_ratio_los_at2_3_*_uci) (fratio_los_from3_*_uci_3y_rp)

. rename (_ratio_los_at2_3_*) (fratio_los_from3_*_3y_rp)

. rename (_diff_los_at2_3_*_lci) (fdiff_los_from3_*_lci_3y_rp)

. rename (_diff_los_at2_3_*_uci) (fdiff_los_from3_*_uci_3y_rp)

. rename (_diff_los_at2_3_*) (fdiff_los_from3_*_3y_rp)

. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         from(4) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una for
&gt; ma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_4_*lci) (fprob_from4a_*lci_3y_rp)

. rename (_prob_at1_4_*uci) (fprob_from4a_*uci_3y_rp)

. rename (_prob_at1_4_*) (fprob_from4a_*_3y_rp)

. 
. rename (_los_at1_4_*lci) (flos_from4a_*lci_3y_rp)

. rename (_los_at1_4_*uci) (flos_from4a_*uci_3y_rp)

. rename (_los_at1_4_*) (flos_from4a_*_3y_rp)

. 
. rename (_prob_at2_4_*lci) (fprob_from4b_*lci_3y_rp)

. rename (_prob_at2_4_*uci) (fprob_from4b_*uci_3y_rp)

. rename (_prob_at2_4_*) (fprob_from4b_*_3y_rp)

. 
. rename (_los_at2_4_*lci) (flos_from4b_*lci_3y_rp)

. rename (_los_at2_4_*uci) (flos_from4b_*uci_3y_rp)

. rename (_los_at2_4_*) (flos_from4b_*_3y_rp)

. 
. rename (_ratio_prob_at2_4_*_lci) (fratio_prob_from4_*_lci_3y_rp)

. rename (_ratio_prob_at2_4_*_uci) (fratio_prob_from4_*_uci_3y_rp)

. rename (_ratio_prob_at2_4_*) (fratio_prob_from4_*_3y_rp)

. rename (_diff_prob_at2_4_*_lci) (fdiff_prob_from4_*_lci_3y_rp)

. rename (_diff_prob_at2_4_*_uci) (fdiff_prob_from4_*_uci_3y_rp)

. rename (_diff_prob_at2_4_*) (fdiff_prob_from4_*_3y_rp) 

. 
. rename (_ratio_los_at2_4_*_lci) (fratio_los_from4_*_lci_3y_rp)

. rename (_ratio_los_at2_4_*_uci) (fratio_los_from4_*_uci_3y_rp)

. rename (_ratio_los_at2_4_*) (fratio_los_from4_*_3y_rp)

. rename (_diff_los_at2_4_*_lci) (fdiff_los_from4_*_lci_3y_rp)

. rename (_diff_los_at2_4_*_uci) (fdiff_los_from4_*_uci_3y_rp)

. rename (_diff_los_at2_4_*) (fdiff_los_from4_*_3y_rp)

. 
</code></pre>
<pre><code>. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         from(5) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una for
&gt; ma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_5_*lci) (fprob_from5a_*lci_3y_rp)

. rename (_prob_at1_5_*uci) (fprob_from5a_*uci_3y_rp)

. rename (_prob_at1_5_*) (fprob_from5a_*_3y_rp)

. 
. rename (_los_at1_5_*lci) (flos_from5a_*lci_3y_rp)

. rename (_los_at1_5_*uci) (flos_from5a_*uci_3y_rp)

. rename (_los_at1_5_*) (flos_from5a_*_3y_rp)

. 
. rename (_prob_at2_5_*lci) (fprob_from5b_*lci_3y_rp)

. rename (_prob_at2_5_*uci) (fprob_from5b_*uci_3y_rp)

. rename (_prob_at2_5_*) (fprob_from5b_*_3y_rp)

. 
. rename (_los_at2_5_*lci) (flos_from5b_*lci_3y_rp)

. rename (_los_at2_5_*uci) (flos_from5b_*uci_3y_rp)

. rename (_los_at2_5_*) (flos_from5b_*_3y_rp)

. 
. rename (_ratio_prob_at2_5_*_lci) (fratio_prob_from5_*_lci_3y_rp)

. rename (_ratio_prob_at2_5_*_uci) (fratio_prob_from5_*_uci_3y_rp)

. rename (_ratio_prob_at2_5_*) (fratio_prob_from5_*_3y_rp)

. rename (_diff_prob_at2_5_*_lci) (fdiff_prob_from5_*_lci_3y_rp)

. rename (_diff_prob_at2_5_*_uci) (fdiff_prob_from5_*_uci_3y_rp)

. rename (_diff_prob_at2_5_*) (fdiff_prob_from5_*_3y_rp) 

. 
. rename (_ratio_los_at2_5_*_lci) (fratio_los_from5_*_lci_3y_rp)

. rename (_ratio_los_at2_5_*_uci) (fratio_los_from5_*_uci_3y_rp)

. rename (_ratio_los_at2_5_*) (fratio_los_from5_*_3y_rp)

. rename (_diff_los_at2_5_*_lci) (fdiff_los_from5_*_lci_3y_rp)

. rename (_diff_los_at2_5_*_uci) (fdiff_los_from5_*_uci_3y_rp)

. rename (_diff_los_at2_5_*) (fdiff_los_from5_*_3y_rp)

. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         from(6) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una for
&gt; ma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_6_*lci) (fprob_from6a_*lci_3y_rp)

. rename (_prob_at1_6_*uci) (fprob_from6a_*uci_3y_rp)

. rename (_prob_at1_6_*) (fprob_from6a_*_3y_rp)

. 
. rename (_los_at1_6_*lci) (flos_from6a_*lci_3y_rp)

. rename (_los_at1_6_*uci) (flos_from6a_*uci_3y_rp)

. rename (_los_at1_6_*) (flos_from6a_*_3y_rp)

. 
. rename (_prob_at2_6_*lci) (fprob_from6b_*lci_3y_rp)

. rename (_prob_at2_6_*uci) (fprob_from6b_*uci_3y_rp)

. rename (_prob_at2_6_*) (fprob_from6b_*_3y_rp)

. 
. rename (_los_at2_6_*lci) (flos_from6b_*lci_3y_rp)

. rename (_los_at2_6_*uci) (flos_from6b_*uci_3y_rp)

. rename (_los_at2_6_*) (flos_from6b_*_3y_rp)

. 
. rename (_ratio_prob_at2_6_*_lci) (fratio_prob_from6_*_lci_3y_rp)

. rename (_ratio_prob_at2_6_*_uci) (fratio_prob_from6_*_uci_3y_rp)

. rename (_ratio_prob_at2_6_*) (fratio_prob_from6_*_3y_rp)

. rename (_diff_prob_at2_6_*_lci) (fdiff_prob_from6_*_lci_3y_rp)

. rename (_diff_prob_at2_6_*_uci) (fdiff_prob_from6_*_uci_3y_rp)

. rename (_diff_prob_at2_6_*) (fdiff_prob_from6_*_3y_rp) 

. 
. rename (_ratio_los_at2_6_*_lci) (fratio_los_from6_*_lci_3y_rp)

. rename (_ratio_los_at2_6_*_uci) (fratio_los_from6_*_uci_3y_rp)

. rename (_ratio_los_at2_6_*) (fratio_los_from6_*_3y_rp)

. rename (_diff_los_at2_6_*_lci) (fdiff_los_from6_*_lci_3y_rp)

. rename (_diff_los_at2_6_*_uci) (fdiff_los_from6_*_uci_3y_rp)

. rename (_diff_los_at2_6_*) (fdiff_los_from6_*_3y_rp)

. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22_corr3.sters&quot;

. 
. predictms, transmatrix(mat_nine_states) models(m2_1_rp9 m2_2_rp10 m2_3_rp10 m2_4_rp7 m2_5_rp3 m2_6_rp4 m2_7_rp4 m2_8_rp2) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt()/obs()
&gt;         seed(2125) n($sim) prob los diff ratio /// * si agrego bootstrap, at#() limit reached, or unrecognised option
&gt;         from(7) ci at1(tipo_de_plan_res_1 0  _start 731) at2(tipo_de_plan_res_1 1  _start 731) reset bootstrap m($boots)   //* es una for
&gt; ma de hacer bootstrap 

Calculating CIs via parametric bootstrap for at1() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

Calculating CIs via parametric bootstrap for at2() (1000)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000

. 
</code></pre>
<pre><code>. rename (_prob_at1_7_*lci) (fprob_from7a_*lci_3y_rp)

. rename (_prob_at1_7_*uci) (fprob_from7a_*uci_3y_rp)

. rename (_prob_at1_7_*) (fprob_from7a_*_3y_rp)

. 
. rename (_los_at1_7_*lci) (flos_from7a_*lci_3y_rp)

. rename (_los_at1_7_*uci) (flos_from7a_*uci_3y_rp)

. rename (_los_at1_7_*) (flos_from7a_*_3y_rp)

. 
. rename (_prob_at2_7_*lci) (fprob_from7b_*lci_3y_rp)

. rename (_prob_at2_7_*uci) (fprob_from7b_*uci_3y_rp)

. rename (_prob_at2_7_*) (fprob_from7b_*_3y_rp)

. 
. rename (_los_at2_7_*lci) (flos_from7b_*lci_3y_rp)

. rename (_los_at2_7_*uci) (flos_from7b_*uci_3y_rp)

. rename (_los_at2_7_*) (flos_from7b_*_3y_rp)

. 
. rename (_ratio_prob_at2_7_*_lci) (fratio_prob_from7_*_lci_3y_rp)

. rename (_ratio_prob_at2_7_*_uci) (fratio_prob_from7_*_uci_3y_rp)

. rename (_ratio_prob_at2_7_*) (fratio_prob_from7_*_3y_rp)

. rename (_diff_prob_at2_7_*_lci) (fdiff_prob_from7_*_lci_3y_rp)

. rename (_diff_prob_at2_7_*_uci) (fdiff_prob_from7_*_uci_3y_rp)

. rename (_diff_prob_at2_7_*) (fdiff_prob_from7_*_3y_rp) 

. 
. rename (_ratio_los_at2_7_*_lci) (fratio_los_from7_*_lci_3y_rp)

. rename (_ratio_los_at2_7_*_uci) (fratio_los_from7_*_uci_3y_rp)

. rename (_ratio_los_at2_7_*) (fratio_los_from7_*_3y_rp)

. rename (_diff_los_at2_7_*_lci) (fdiff_los_from7_*_lci_3y_rp)

. rename (_diff_los_at2_7_*_uci) (fdiff_los_from7_*_uci_3y_rp)

. rename (_diff_los_at2_7_*) (fdiff_los_from7_*_3y_rp)

</code></pre>
<p>Saved at= 17:09:35 24 Oct 2022</p>
<pre><code>
</code></pre>