请大佬解释下生存分析代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

n <- 1000<br>
set.seed(731)<br>
age <- 50 + 12*rnorm(n)<br>
label(age) <- "Age"<br>
sex <- factor(sample(c('Male','Female'), n,<br>
rep=TRUE, prob=c(.6, .4)))<br>
cens <- 15*runif(n)<br>
h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))<br>
dt <- -log(runif(n))/h<br>
label(dt) <- 'Follow-up Time'<br>
e <- ifelse(dt <= cens,1,0)<br>
dt <- pmin(dt, cens)<br>
units(dt) <- "Year"
主要是Age和h看不懂

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：生存分析 Female Sample factor follow

R初探者 发表于 2019-1-5 11:22
n <- 1000<br>
set.seed(731)<br>
age <- 50 + 12*rnorm(n)<br>

n &lt;- 1000<br>
set.seed(731)<br>
age &lt;- 50 + 12*rnorm(n)<br>
label(age) &lt;- \"Age\"<br>
sex &lt;- factor(sample(c(\'Male\',\'Female\'), n,<br>
rep=TRUE, prob=c(.6, .4)))<br>
cens &lt;- 15*runif(n)<br>
h &lt;- .02*exp(.04*(age-50)+.8*(sex==\'Female\'))<br>
dt &lt;- -log(runif(n))/h<br>
label(dt) &lt;- \'Follow-up Time\'<br>
e &lt;- ifelse(dt &lt;= cens,1,0)<br>
dt &lt;- pmin(dt, cens)<br>
units(dt) &lt;- \"Year\"

R初探者 发表于 2019-1-5 11:22
n &amp;amp;lt;- 1000&amp;lt;br&amp;gt;
set.seed(731)&amp;lt;br&amp;gt;
age &amp;amp;lt;- 50 + 12*rnorm(n)&amp;lt;br&amp;gt;

age和h都在随机生成观测值。
age服从均值50标准差12的正态分布。
h这里应该是hazard吧，是0.2*exp(0.04*(age-50) + 0.08*female)。
Female是二分类变量，女性时候为1，男性为0，是随机生成的，男性女性概率均为0.5。<br>