これは、SE に関する以前の投稿の一種のフォローアップです: https://stats.stackexchange.com/questions/70858/right-censored-survival-fit-with-jags
しかし、ここでは、完全な R スクリプト (最初から最後まで) が JAGS の右打ち切りデータに対して生存分析を実行することを望みます。私が見つけたすべてのサイトでは、JAGS に非常に高いレベルの習熟度が求められるため、あるコード行から別のコード行に移動する方法を理解するのは困難です。私はこれが多くの質問であることを知っています...
とにかく、生存データの例をいくつか示します。グループは t1、t2、t3 です。NA は、右打ち切りデータを参照します (検閲カットオフ = 3)。
t1 <- c(1.73, NA, NA, NA, NA,0.77, NA, NA, NA, NA, NA, NA,0.5,1.06, NA, NA, NA, NA, NA,0.46, NA)
t2 <- c(1.42, NA, NA, NA, NA, NA,0.69,1.84, NA, NA, NA,1.47,1.6,1.8, NA, NA, NA, NA, NA,0.73, NA,1.28,3,2.97)
t3 <- c(0.88, NA, NA,1.65,1.75, NA, NA,1.01,1.46,1.95, NA, NA,2.93, NA,0.78,1.05,1.52, NA)
#Specify model in BUGS language
sink("model.txt")
cat("
model
{
}
",fill = TRUE)
sink()
#Bundle data
data<- list()
#Parameters monitored
parameters<-c()
#Initial values
inits <- list(
# MCMC settings
ni <-
nt <-
nb <-
nc <-
fit <- jags(data, inits, parameters, "model.txt", n.iter=ni, n.thin=nt, n.burnin=nb, n.chains=nc, working.directory = getwd())
質問が多いことは承知していますが、何かをつなぎ合わせようと何日も費やしており、迷子になったり混乱したりしています。この種の分析を実行するパッケージがあることは知っていますが、これをゼロから自分で構築する方法を本当に学びたいです! 読者の皆さん、ありがとう！