R语言做cox回归

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在用R语言做cox回归，用什么函数来生成服从某一概率分布的随机数（生存时间）

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Cox R语言 概率分布 随机数 统计

这个给顶起来。。。。。。

谢谢帮顶的兄弟哈  

n = 10000
beta1 = 2; beta2 = -1
lambdaT = .002 # baseline hazard
lambdaC = .004  # hazard of censoring

x1 = rnorm(n,0)
x2 = rnorm(n,0)
# true event time
T = rweibull(n, shape=1, scale=lambdaT*exp(-beta1*x1-beta2*x2))
C = rweibull(n, shape=1, scale=lambdaC)  #censoring time
time = pmin(T,C)  #observed time is minof censored and true
event = time==T   # set to 1 if event isobserve

DM小菜鸟 发表于 2014-12-29 17:58
谢谢帮顶的兄弟哈  

n = 10000

lambdaT = .002 # baseline hazard
lambdaC = .004  # hazard of censoring
lamda数值是自己随便设置吗？
T = rweibull(n, shape=1,scale=lambdaT*exp(-beta1*x1-beta2*x2))里面不用加随机误差项吗？
求解惑，，谢谢

客户流失快？Cox回归帮你寻找真相https://bbs.pinggu.org/thread-10847167-1-1.html