data(telomerase)
model1 = cdtamodel(copula = 'fgm')
model2 = cdtamodel(copula = 'fgm',
modelargs=list(param=2,
prior.lse='normal',
par.lse1=0,
par.lse2=5,
prior.lsp='normal',
par.lsp1=0,
par.lsp2=5))
model3 = cdtamodel(copula = 'fgm',
modelargs = list(formula.se = StudyID ~ Test - 1))
if (FALSE) {
fit1 <- fit(model1,
SID='ID',
data=telomerase,
iter=2000,
warmup=1000,
thin=1,
seed=3)
plot(fit1)
}
Run the code above in your browser using DataLab