我是一个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))但它不起作用..
1条答案
按热度按时间oxiaedzo1#
你几乎已经完成了。你只需要使用同一个包中的
timepoints()
函数,并指定你想要的时间点: