Database (step 4)

<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-4" id="database-step-4">Database (step 4)</a></h1>
<p>Date created: 12:12:48  2 Aug 2022.</p>
<p>Install commands that are unavailable or out of date.</p>
<pre><code>. clear all

. 
. cap noi which predictms
c:\ado\plus\p\predictms.ado
*! version 4.3.0 14mar2021 MJC

. if _rc==111 {
.         cap noi net install multistate, from(&quot;https://www.mjcrowther.co.uk/code/multistate&quot;) 
.         }

. cap noi which merlin
c:\ado\plus\m\merlin.ado
*! version 2.0.2 19mar2021 MJC

. if _rc==111 {
.         cap noi net install merlin, from(&quot;https://www.mjcrowther.co.uk/code/merlin/&quot;) 
.         }

. cap noi which sumat
c:\ado\plus\s\sumat.ado
*! Part of package matrixtools v. 0.28
*! Support: Niels Henrik Bruun, niels.henrik.bruun@gmail.com
*! 2021-01-03 toxl added

. if _rc==111 {
.         cap noi scc install matrixtools
.         }

. cap noi which estwrite
c:\ado\plus\e\estwrite.ado
*! version 1.2.4 04sep2009
*! version 1.0.1 15may2007 (renamed from -eststo- to -estwrite-; -append- added)
*! version 1.0.0 29apr2005 Ben Jann (ETH Zurich)

. if _rc==111 {
.         cap noi ssc install estwrite
.         }       

</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 u
&gt; ser: `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_ag
&gt; s\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags&quot;
E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags
.         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 u
&gt; ser: `c(username)'&quot;
Location= E:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ags; Date:  2 Aug 2022, conside
&gt; ring an OS Windows for the user: andre
. }

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

. 
. mata : st_numscalar(&quot;OK&quot;, direxists(&quot;G:\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_state_ag
&gt; s\&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 u
&gt; ser: `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 (gi
&gt; thub)\_mult_state_ags\&quot;))

. if scalar(OK) == 1 {
.         cap noi cd &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)&quot;
.         global pathdata &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_
&gt; state_ags&quot;
.         global pathdata2 &quot;C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult
&gt; _state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the u
&gt; ser: `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\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 u
&gt; ser: `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
&gt; _ags&quot;
.         global pathdata &quot;C:\Users\CISS Fondecyt\Mi unidad\Alvacast\SISTRAT 2019 (github)\_mult_
&gt; state_ags&quot;
.         global pathdata2 &quot;C:/Users/CISS Fondecyt/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult
&gt; _state_ags/&quot;
.         di &quot;Location= ${pathdata}; Date: `c(current_date)', considering an OS `c(os)' for the u
&gt; ser: `c(username)'&quot;
. }

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

. 
</code></pre>
<p>Path data= E:\Mi unidad\Alvacast\SISTRAT 2019 (github)_mult_state_ags;</p>
<p>Timestamp:  2 Aug 2022, considering that is a Windows OS for the username: andre</p>
<p>The file is located and named as: E:/Mi unidad/Alvacast/SISTRAT 2019 (github)/_mult_state_ags/archivo_multiestado_jun_22.dta</p>
<p>Defined transition matrix.</p>
<p>=============================================================================</p>
<h2><a href="#transition-probabilities-lengths-of-stay-parametric-models-3-years-treatment-completion-at-baseline" id="transition-probabilities-lengths-of-stay-parametric-models-3-years-treatment-completion-at-baseline">Transition probabilities &amp; Lengths of Stay, Parametric Models (3 years) &amp; Treatment Completion at Baseline</a></h2>
<p>=============================================================================</p>
<pre><code>. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22.sters&quot;

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(1) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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_1_*lci) (fprob2_from1a_*lci_3y)

. rename (_prob_at1_1_*uci) (fprob2_from1a_*uci_3y)

. rename (_prob_at1_1_*) (fprob2_from1a_*_3y)

. 
. rename (_los_at1_1_*lci) (flos2_from1a_*lci_3y)

. rename (_los_at1_1_*uci) (flos2_from1a_*uci_3y)

. rename (_los_at1_1_*) (flos2_from1a_*_3y)

. 
. rename (_prob_at2_1_*lci) (fprob2_from1b_*lci_3y)

. rename (_prob_at2_1_*uci) (fprob2_from1b_*uci_3y)

. rename (_prob_at2_1_*) (fprob2_from1b_*_3y)

. 
. rename (_los_at2_1_*lci) (flos2_from1b_*lci_3y)

. rename (_los_at2_1_*uci) (flos2_from1b_*uci_3y)

. rename (_los_at2_1_*) (flos2_from1b_*_3y)

. 
. rename (_ratio_prob_at2_1_*_lci) (fratio2_prob_from1_*_lci_3y)

. rename (_ratio_prob_at2_1_*_uci) (fratio2_prob_from1_*_uci_3y)

. rename (_ratio_prob_at2_1_*) (fratio2_prob_from1_*_3y)

. rename (_diff_prob_at2_1_*_lci) (fdiff2_prob_from1_*_lci_3y)

. rename (_diff_prob_at2_1_*_uci) (fdiff2_prob_from1_*_uci_3y)

. rename (_diff_prob_at2_1_*) (fdiff2_prob_from1_*_3y)  

. 
. rename (_ratio_los_at2_1_*_lci) (fratio2_los_from1_*_lci_3y)

. rename (_ratio_los_at2_1_*_uci) (fratio2_los_from1_*_uci_3y)

. rename (_ratio_los_at2_1_*) (fratio2_los_from1_*_3y)

. rename (_diff_los_at2_1_*_lci) (fdiff2_los_from1_*_lci_3y)

. rename (_diff_los_at2_1_*_uci) (fdiff2_los_from1_*_uci_3y)

. rename (_diff_los_at2_1_*) (fdiff2_los_from1_*_3y)

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(2) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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) (fprob2_from2a_*lci_3y)

. rename (_prob_at1_2_*uci) (fprob2_from2a_*uci_3y)

. rename (_prob_at1_2_*) (fprob2_from2a_*_3y)

. 
. rename (_los_at1_2_*lci) (flos2_from2a_*lci_3y)

. rename (_los_at1_2_*uci) (flos2_from2a_*uci_3y)

. rename (_los_at1_2_*) (flos2_from2a_*_3y)

. 
. rename (_prob_at2_2_*lci) (fprob2_from2b_*lci_3y)

. rename (_prob_at2_2_*uci) (fprob2_from2b_*uci_3y)

. rename (_prob_at2_2_*) (fprob2_from2b_*_3y)

. 
. rename (_los_at2_2_*lci) (flos2_from2b_*lci_3y)

. rename (_los_at2_2_*uci) (flos2_from2b_*uci_3y)

. rename (_los_at2_2_*) (flos2_from2b_*_3y)

. 
. rename (_ratio_prob_at2_2_*_lci) (fratio2_prob_from2_*_lci_3y)

. rename (_ratio_prob_at2_2_*_uci) (fratio2_prob_from2_*_uci_3y)

. rename (_ratio_prob_at2_2_*) (fratio2_prob_from2_*_3y)

. rename (_diff_prob_at2_2_*_lci) (fdiff2_prob_from2_*_lci_3y)

. rename (_diff_prob_at2_2_*_uci) (fdiff2_prob_from2_*_uci_3y)

. rename (_diff_prob_at2_2_*) (fdiff2_prob_from2_*_3y) 

. 
. rename (_ratio_los_at2_2_*_lci) (fratio2_los_from2_*_lci_3y)

. rename (_ratio_los_at2_2_*_uci) (fratio2_los_from2_*_uci_3y)

. rename (_ratio_los_at2_2_*) (fratio2_los_from2_*_3y)

. rename (_diff_los_at2_2_*_lci) (fdiff2_los_from2_*_lci_3y)

. rename (_diff_los_at2_2_*_uci) (fdiff2_los_from2_*_uci_3y)

. rename (_diff_los_at2_2_*) (fdiff2_los_from2_*_3y) 

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(3) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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) (fprob2_from3a_*lci_3y)

. rename (_prob_at1_3_*uci) (fprob2_from3a_*uci_3y)

. rename (_prob_at1_3_*) (fprob2_from3a_*_3y)

. 
. rename (_los_at1_3_*lci) (flos2_from3a_*lci_3y)

. rename (_los_at1_3_*uci) (flos2_from3a_*uci_3y)

. rename (_los_at1_3_*) (flos2_from3a_*_3y)

. 
. rename (_prob_at2_3_*lci) (fprob2_from3b_*lci_3y)

. rename (_prob_at2_3_*uci) (fprob2_from3b_*uci_3y)

. rename (_prob_at2_3_*) (fprob2_from3b_*_3y)

. 
. rename (_los_at2_3_*lci) (flos2_from3b_*lci_3y)

. rename (_los_at2_3_*uci) (flos2_from3b_*uci_3y)

. rename (_los_at2_3_*) (flos2_from3b_*_3y)

. 
. rename (_ratio_prob_at2_3_*_lci) (fratio2_prob_from3_*_lci_3y)

. rename (_ratio_prob_at2_3_*_uci) (fratio2_prob_from3_*_uci_3y)

. rename (_ratio_prob_at2_3_*) (fratio2_prob_from3_*_3y)

. rename (_diff_prob_at2_3_*_lci) (fdiff2_prob_from3_*_lci_3y)

. rename (_diff_prob_at2_3_*_uci) (fdiff2_prob_from3_*_uci_3y)

. rename (_diff_prob_at2_3_*) (fdiff2_prob_from3_*_3y)  

. 
. rename (_ratio_los_at2_3_*_lci) (fratio2_los_from3_*_lci_3y)

. rename (_ratio_los_at2_3_*_uci) (fratio2_los_from3_*_uci_3y)

. rename (_ratio_los_at2_3_*) (fratio2_los_from3_*_3y)

. rename (_diff_los_at2_3_*_lci) (fdiff2_los_from3_*_lci_3y)

. rename (_diff_los_at2_3_*_uci) (fdiff2_los_from3_*_uci_3y)

. rename (_diff_los_at2_3_*) (fdiff2_los_from3_*_3y) 

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(4) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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) (fprob2_from4a_*lci_3y)

. rename (_prob_at1_4_*uci) (fprob2_from4a_*uci_3y)

. rename (_prob_at1_4_*) (fprob2_from4a_*_3y)

. 
. rename (_los_at1_4_*lci) (flos2_from4a_*lci_3y)

. rename (_los_at1_4_*uci) (flos2_from4a_*uci_3y)

. rename (_los_at1_4_*) (flos2_from4a_*_3y)

. 
. rename (_prob_at2_4_*lci) (fprob2_from4b_*lci_3y)

. rename (_prob_at2_4_*uci) (fprob2_from4b_*uci_3y)

. rename (_prob_at2_4_*) (fprob2_from4b_*_3y)

. 
. rename (_los_at2_4_*lci) (flos2_from4b_*lci_3y)

. rename (_los_at2_4_*uci) (flos2_from4b_*uci_3y)

. rename (_los_at2_4_*) (flos2_from4b_*_3y)

. 
. rename (_ratio_prob_at2_4_*_lci) (fratio2_prob_from4_*_lci_3y)

. rename (_ratio_prob_at2_4_*_uci) (fratio2_prob_from4_*_uci_3y)

. rename (_ratio_prob_at2_4_*) (fratio2_prob_from4_*_3y)

. rename (_diff_prob_at2_4_*_lci) (fdiff2_prob_from4_*_lci_3y)

. rename (_diff_prob_at2_4_*_uci) (fdiff2_prob_from4_*_uci_3y)

. rename (_diff_prob_at2_4_*) (fdiff2_prob_from4_*_3y)  

. 
. rename (_ratio_los_at2_4_*_lci) (fratio2_los_from4_*_lci_3y)

. rename (_ratio_los_at2_4_*_uci) (fratio2_los_from4_*_uci_3y)

. rename (_ratio_los_at2_4_*) (fratio2_los_from4_*_3y)

. rename (_diff_los_at2_4_*_lci) (fdiff2_los_from4_*_lci_3y)

. rename (_diff_los_at2_4_*_uci) (fdiff2_los_from4_*_uci_3y)

. rename (_diff_los_at2_4_*) (fdiff2_los_from4_*_3y) 

</code></pre>
<p>=============================================================================</p>
<h2><a href="#transition-probabilities-length-of-stay-parametric-models-rp-3-years-treatment-completion-at-baseline" id="transition-probabilities-length-of-stay-parametric-models-rp-3-years-treatment-completion-at-baseline">Transition Probabilities &amp; Length of Stay, Parametric Models (RP) (3 years) &amp; Treatment Completion at Baseline</a></h2>
<p>=============================================================================</p>
<pre><code>. 
. qui estread &quot;${pathdata2}parmodels_m2_jun_22.sters&quot;

. 
. predictms, transmatrix(mat_five_states) models(m2_1_rp5 m2_2_rp4 m2_3_rp2 m2_4_rp3) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(1) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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_1_*lci) (fprob2_from1a_*lci_3y_rp)

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

. rename (_prob_at1_1_*) (fprob2_from1a_*_3y_rp)

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

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

. rename (_los_at1_1_*) (flos2_from1a_*_3y_rp)

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

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

. rename (_prob_at2_1_*) (fprob2_from1b_*_3y_rp)

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

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

. rename (_los_at2_1_*) (flos2_from1b_*_3y_rp)

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

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

. rename (_ratio_prob_at2_1_*) (fratio2_prob_from1_*_3y_rp)

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

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

. rename (_diff_prob_at2_1_*) (fdiff2_prob_from1_*_3y_rp) 

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

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

. rename (_ratio_los_at2_1_*) (fratio2_los_from1_*_3y_rp)

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

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

. rename (_diff_los_at2_1_*) (fdiff2_los_from1_*_3y_rp)

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_rp5 m2_2_rp4 m2_3_rp2 m2_4_rp3) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(2) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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) (fprob2_from2a_*lci_3y_rp)

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

. rename (_prob_at1_2_*) (fprob2_from2a_*_3y_rp)

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

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

. rename (_los_at1_2_*) (flos2_from2a_*_3y_rp)

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

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

. rename (_prob_at2_2_*) (fprob2_from2b_*_3y_rp)

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

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

. rename (_los_at2_2_*) (flos2_from2b_*_3y_rp)

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

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

. rename (_ratio_prob_at2_2_*) (fratio2_prob_from2_*_3y_rp)

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

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

. rename (_diff_prob_at2_2_*) (fdiff2_prob_from2_*_3y_rp) 

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

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

. rename (_ratio_los_at2_2_*) (fratio2_los_from2_*_3y_rp)

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

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

. rename (_diff_los_at2_2_*) (fdiff2_los_from2_*_3y_rp)

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_rp5 m2_2_rp4 m2_3_rp2 m2_4_rp3) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(3) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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) (fprob2_from3a_*lci_3y_rp)

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

. rename (_prob_at1_3_*) (fprob2_from3a_*_3y_rp)

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

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

. rename (_los_at1_3_*) (flos2_from3a_*_3y_rp)

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

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

. rename (_prob_at2_3_*) (fprob2_from3b_*_3y_rp)

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

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

. rename (_los_at2_3_*) (flos2_from3b_*_3y_rp)

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

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

. rename (_ratio_prob_at2_3_*) (fratio2_prob_from3_*_3y_rp)

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

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

. rename (_diff_prob_at2_3_*) (fdiff2_prob_from3_*_3y_rp) 

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

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

. rename (_ratio_los_at2_3_*) (fratio2_los_from3_*_3y_rp)

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

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

. rename (_diff_los_at2_3_*) (fdiff2_los_from3_*_3y_rp)

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_rp5 m2_2_rp4 m2_3_rp2 m2_4_rp3) ///
&gt;         timevar(timevar0) /// mint(0) maxt(1826) timevar() cannot be specified with mint()/maxt
&gt; ()/obs()
&gt;         seed(2125) n(100000) prob los diff ratio /// * si agrego bootstrap, at#() limit reached
&gt; , or unrecognised option
&gt;         from(4) ci at1(tipo_de_plan_res_1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_d
&gt; e_plan_res_1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) reset bootstrap m(1000)   //* es una fo
&gt; rma 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) (fprob2_from4a_*lci_3y_rp)

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

. rename (_prob_at1_4_*) (fprob2_from4a_*_3y_rp)

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

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

. rename (_los_at1_4_*) (flos2_from4a_*_3y_rp)

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

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

. rename (_prob_at2_4_*) (fprob2_from4b_*_3y_rp)

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

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

. rename (_los_at2_4_*) (flos2_from4b_*_3y_rp)

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

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

. rename (_ratio_prob_at2_4_*) (fratio2_prob_from4_*_3y_rp)

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

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

. rename (_diff_prob_at2_4_*) (fdiff2_prob_from4_*_3y_rp) 

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

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

. rename (_ratio_los_at2_4_*) (fratio2_los_from4_*_3y_rp)

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

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

. rename (_diff_los_at2_4_*) (fdiff2_los_from4_*_3y_rp)

. 
</code></pre>
<p>=============================================================================</p>
<h2><a href="#ambulatory-parametric-models-of-differences-ratios-between-completion-and-non-completion-of-baseline-treatment-at-3-years" id="ambulatory-parametric-models-of-differences-ratios-between-completion-and-non-completion-of-baseline-treatment-at-3-years">Ambulatory- Parametric models of differences &amp; ratios between Completion and Non-completion of baseline treatment- At 3 years</a></h2>
<p>=============================================================================</p>
<pre><code>. qui estread &quot;${pathdata2}parmodels_m2_jun_22.sters&quot;

. *m2_1_ggam m2_2_ggam m2_3_logn m2_4_logn
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(1) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 0 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_1_*) (fprob2a_from1a_*_3y)

. rename (_prob_at2_1_*) (fprob2a_from1b_*_3y)

. 
. rename (_diff_prob_at2_1_*_lci) (fdiff2a_prob_from1_*_lci_3y)

. rename (_diff_prob_at2_1_*_uci) (fdiff2a_prob_from1_*_uci_3y)

. rename (_diff_prob_at2_1_*) (fdiff2a_prob_from1_*_3y)

. 
. rename (_diff_los_at2_1_*_lci) (fdiff2a_los_from1_*_lci_3y)

. rename (_diff_los_at2_1_*_uci) (fdiff2a_los_from1_*_uci_3y)

. rename (_diff_los_at2_1_*) (fdiff2a_los_from1_*_3y  
parentheses unbalanced
r(198);

. 
. rename (_diff_los_at2_1_*_lci) (fdiff2a_los_from1_*_lci_3y)
variable _diff_los_at2_1_*_lci not found
r(111);

. rename (_diff_los_at2_1_*_uci) (fdiff2a_los_from1_*_uci_3y)
variable _diff_los_at2_1_*_uci not found
r(111);

. rename (_diff_los_at2_1_*) (fdiff2a_los_from1_*_3y)  

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(2) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 0 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_2_*) (fprob2a_from2a_*_3y)

. rename (_prob_at2_2_*) (fprob2a_from2b_*_3y)

. 
. rename (_diff_prob_at2_2_*_lci) (fdiff2a_prob_from2_*_lci_3y)

. rename (_diff_prob_at2_2_*_uci) (fdiff2a_prob_from2_*_uci_3y)

. rename (_diff_prob_at2_2_*) (fdiff2a_prob_from2_*_3y)  

. 
. rename (_diff_los_at2_2_*_lci) (fdiff2a_los_from2_*_lci_3y)

. rename (_diff_los_at2_2_*_uci) (fdiff2a_los_from2_*_uci_3y)

. rename (_diff_los_at2_2_*) (fdiff2a_los_from2_*_3y)  

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(3) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 0 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_3_*) (fprob2a_from3a_*_3y)

. rename (_prob_at2_3_*) (fprob2a_from3b_*_3y)

. 
. rename (_diff_prob_at2_3_*_lci) (fdiff2a_prob_from3_*_lci_3y)

. rename (_diff_prob_at2_3_*_uci) (fdiff2a_prob_from3_*_uci_3y)

. rename (_diff_prob_at2_3_*) (fdiff2a_prob_from3_*_3y)  

. 
. rename (_diff_los_at2_3_*_lci) (fdiff2a_los_from3_*_lci_3y)

. rename (_diff_los_at2_3_*_uci) (fdiff2a_los_from3_*_uci_3y)

. rename (_diff_los_at2_3_*) (fdiff2a_los_from3_*_3y)  

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(4) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 0 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 0 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_4_*) (fprob2a_from4a_*_3y)

. rename (_prob_at2_4_*) (fprob2a_from4b_*_3y)

.         
. rename (_diff_prob_at2_4_*_lci) (fdiff2a_prob_from4_*_lci_3y)

. rename (_diff_prob_at2_4_*_uci) (fdiff2a_prob_from4_*_uci_3y)

. rename (_diff_prob_at2_4_*) (fdiff2a_prob_from4_*_3y)  

.         
. rename (_diff_los_at2_4_*_lci) (fdiff2a_los_from4_*_lci_3y)

. rename (_diff_los_at2_4_*_uci) (fdiff2a_los_from4_*_uci_3y)

. rename (_diff_los_at2_4_*) (fdiff2a_los_from4_*_3y)  

. 
</code></pre>
<p>=============================================================================</p>
<h2><a href="#residential-parametric-models-of-differences-ratios-between-completion-and-non-completion-of-baseline-treatment-at-3-years" id="residential-parametric-models-of-differences-ratios-between-completion-and-non-completion-of-baseline-treatment-at-3-years">Residential- Parametric models of differences &amp; ratios between Completion and Non-completion of baseline treatment- At 3 years</a></h2>
<p>=============================================================================</p>
<pre><code>. qui estread &quot;${pathdata2}parmodels_m2_jun_22.sters&quot;

. *m2_1_ggam m2_2_ggam m2_3_logn m2_4_logn
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(1) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 1 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_1_*) (fprob2b_from1a_*_3y)

. rename (_prob_at2_1_*) (fprob2b_from1b_*_3y)

. 
. rename (_diff_prob_at2_1_*_lci) (fdiff2b_prob_from1_*_lci_3y)

. rename (_diff_prob_at2_1_*_uci) (fdiff2b_prob_from1_*_uci_3y)

. rename (_diff_prob_at2_1_*) (fdiff2b_prob_from1_*_3y)  

. 
. rename (_diff_los_at2_1_*_lci) (fdiff2b_los_from1_*_lci_3y)

. rename (_diff_los_at2_1_*_uci) (fdiff2b_los_from1_*_uci_3y)

. rename (_diff_los_at2_1_*) (fdiff2b_los_from1_*_3y)  

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(2) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 1 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_2_*) (fprob2b_from2a_*_3y)

. rename (_prob_at2_2_*) (fprob2b_from2b_*_3y)

. 
. rename (_diff_prob_at2_2_*_lci) (fdiff2b_prob_from2_*_lci_3y)

. rename (_diff_prob_at2_2_*_uci) (fdiff2b_prob_from2_*_uci_3y)

. rename (_diff_prob_at2_2_*) (fdiff2b_prob_from2_*_3y)  

. 
. rename (_diff_los_at2_2_*_lci) (fdiff2b_los_from2_*_lci_3y)

. rename (_diff_los_at2_2_*_uci) (fdiff2b_los_from2_*_uci_3y)

. rename (_diff_los_at2_2_*) (fdiff2b_los_from2_*_3y)  

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(3) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 1 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_3_*) (fprob2b_from3a_*_3y)

. rename (_prob_at2_3_*) (fprob2b_from3b_*_3y)

. 
. rename (_diff_prob_at2_3_*_lci) (fdiff2b_prob_from3_*_lci_3y)

. rename (_diff_prob_at2_3_*_uci) (fdiff2b_prob_from3_*_uci_3y)

. rename (_diff_prob_at2_3_*) (fdiff2b_prob_from3_*_3y)  

. 
. rename (_diff_los_at2_3_*_lci) (fdiff2b_los_from3_*_lci_3y)

. rename (_diff_los_at2_3_*_uci) (fdiff2b_los_from3_*_uci_3y)

. rename (_diff_los_at2_3_*) (fdiff2b_los_from3_*_3y)  

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

. 
. predictms, transmatrix(mat_five_states) models(m2_1_ggam m2_2_ggam m2_3_ggam m2_4_logn) ///
&gt;          timevar(timevar0) /// mint(0) maxt(1096) timevar() cannot be specified with mint()/max
&gt; t()/obs()
&gt;          seed(2125) n(100000) prob los diff /// * si agrego bootstrap, at#() limit reached, or 
&gt; unrecognised option
&gt;          from(4) ci reset bootstrap m(1000)  /// * es una forma de hacer bootstrap 
&gt;          at1(tipo_de_plan_res_1 1 TD_1 0 TD_2 0 TD_3 0 TD_4 0 _start 1096) at2(tipo_de_plan_res
&gt; _1 1 TD_1 1 TD_2 0 TD_3 0 TD_4 0 _start 1096) 

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

.                                  
. rename (_prob_at1_4_*) (fprob2b_from4a_*_3y)

. rename (_prob_at2_4_*) (fprob2b_from4b_*_3y)

. 
. rename (_diff_prob_at2_4_*_lci) (fdiff2b_prob_from4_*_lci_3y)

. rename (_diff_prob_at2_4_*_uci) (fdiff2b_prob_from4_*_uci_3y)

. rename (_diff_prob_at2_4_*) (fdiff2b_prob_from4_*_3y)  

. 
. rename (_diff_los_at2_4_*_lci) (fdiff2b_los_from4_*_lci_3y)

. rename (_diff_los_at2_4_*_uci) (fdiff2b_los_from4_*_uci_3y)

. rename (_diff_los_at2_4_*) (fdiff2b_los_from4_*_3y)  

. 
</code></pre>
<h1><a href="#">=============================================================================</a></h1>
<p>Saved at= 20:03:00  4 Aug 2022</p>
<pre><code>
</code></pre>