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

1012

收藏 2019-01-05

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

扫码加我拉你入群

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

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

全部回复

R初探者

2019-1-5 11:23:37

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

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

扫码加我拉你入群

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

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

caimiao0714

2019-1-5 20:39:24

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

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

扫码加我拉你入群

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

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

栏目导航

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群