Получение совокупных значений риска - PullRequest
1 голос
/ 27 сентября 2019

Рассмотрим следующий игрушечный пример:

use https://data.princeton.edu/pop509/justices2.dta, clear
stset tenure, fail(event == 1)
stcrreg age year, compete (event == 2)
stcurve, cif

Я хочу построить кумулятивную кривую заболеваемости, как описано выше, но затем я хочу сохранить эти значения с их 95% доверительными интервалами.Однако мне не ясно, как получить доступ к ним / сохранить их как переменные.

Перекрестная публикация на Statalist .

1 Ответ

2 голосов
/ 27 сентября 2019

Используйте параметр outfile() команды stcurve:

stcurve, cif outfile(stdata)

use stdata

list in 1/10

     +---------------------+
     |      ci1         _t |
     |---------------------|
  1. | .0465373   5.691992 |
  2. |        0   1.045859 |
  3. | .2600816    20.6078 |
  4. | .1169629   8.876112 |
  5. | .0465373   5.724846 |
     |---------------------|
  6. | .1249585   9.440109 |
  7. |        0   .4462697 |
  8. | .1574731   13.49213 |
  9. | .1991083   15.36756 |
 10. | .0232038   4.769336 |
     +---------------------+
...