使用R中的“cuminc”函数定义累积发生率曲线的特定时间点

gmxoilav  于 2023-04-18  发布在  其他
关注(0)|答案(1)|浏览(425)

我是一个R语言的博士生初学者,我在使用R语言中的cuminc()函数获得的累积发生率函数中定义特定时间点时遇到了一些问题。
下面是与我的数据集对应的模拟数据集:

# Set the sample size
n<- 8076
# CReate a variable for follow-up time
time<- c(rep(180,6533),sample(1:179,n-6533,replace = TRUE))
# Create a variable for status
status<-ifelse(time==180,0,1) # O= alive/censored
                          # 1 = death
# Create age
age<-sample(18:90,n,replace = TRUE)
# Create gender
sex<-sample(c("male","female"),n,replace = TRUE)
# Combine
df<-data.frame(time,status,age,sex) `

下面是计算累积关联函数和绘制累积关联函数曲线的代码:

`# Cumulative incidence function
fit<-cuminc(df$time,df$status,df$sex)
print(fit)
# Plot the cumulative incidence curves
plot(fit,xlab = "Days",ylab = "Cumulative incidence",ylim=c(0,0.4))`

获得的时间点是0、50、100、150,但我希望时间点是0、30、90和180天(在图表上和print(fit)提供的输出上)。据我所知,“Cuminc”中的参数“timepoint”需要用于定义时间点,但我不知道如何正确使用此参数。
非常感谢您的帮助提前:)
我尝试定义一个向量,如:times〈-c(0,30,90,180),我尝试在“cuminc”代码中添加这个向量。我还尝试通过使用xlim=c(0,30,90,180)直接修改图中的x轴,但我收到一条错误消息,而对于y-lim,它可以正常工作。我还尝试在cuminc函数中使用timepoints参数:时间点(c(0,30,90,180))但它不起作用..

oxiaedzo

oxiaedzo1#

你几乎已经完成了。你只需要使用同一个包中的timepoints()函数,并指定你想要的时间点:

timepoints(fit, c(0,30,90,180))
$est
         0         30         90       180
female 1 0 0.02985075 0.09676617 0.1878109
male 1   0 0.03402367 0.10108481 0.1942801

$var
         0           30           90          180
female 1 0 7.205719e-06 2.174761e-05 3.795530e-05
male 1   0 8.105104e-06 2.240886e-05 3.860418e-05

相关问题