我是 R 的新手。目前我正在为我拥有的一些生存数据拟合对数正态分布,但是在尝试计算中位数和平均值等统计数据时我陷入了困境。这是我到目前为止使用的代码,谁能告诉我接下来应该输入什么来找到平均值?
# rm(list=ls(all=TRUE))
library(survival)
data<-read.table("M:\\w2k\\Diss\\Hoyle And Henley True IPD with number at risk known.txt",header=T)
attach(data)
data
times_start <-c( rep(start_time_censor, n_censors), rep(start_time_event, n_events) )
times_end <-c( rep(end_time_censor, n_censors), rep(end_time_event, n_events) )
model <- survreg(Surv(times_start, times_end, type="interval2")~1, dist="lognormal")
intercept <- summary(model)$table[1]
log_scale <- summary(model)$table[2]
这是我卡住的地方,我试过:
meantime<-exp(intercept+log_scale/2)
但这似乎并没有给出一个现实的意思。