Nothing
`trecase.b1.A` =
function(yi, ind.lst, X, twosex, sex, ni, ni0, xs, start, iphi=1, theta=1, maxiter=100, eps=1E-3, tech.ctrl){
triali = 0
par0 = start
if(twosex){
zeroes = rep(0, 2)
seq = 1:4
np = 1:11
bsi1 = c(3, 5)
bsi2 = c(4, 6)
}else{
zeroes = rep(0, 1)
seq = 1:2
np = 1:6
bsi1 = c(2, 3)
bsi2 = c(2, 3)
}
repeat{
triali = triali + 1
tag = tryCatch({
iphi_i = iphi
theta_i = theta
log.lik0 = ll.jRCI.b1.A(par0, yi=yi, ind.lst=ind.lst, X=X, twosex=twosex, sex=sex,
ni=ni, ni0=ni0, xs=xs, iphi=iphi_i, theta=theta_i)
for(i in 1:maxiter){
out = optim(par0, ll.jRCI.b1.A, yi=yi, ind.lst=ind.lst, X=X, twosex=twosex, sex=sex,
ni=ni, ni0=ni0, xs=xs, iphi=iphi_i, theta=theta_i)
par0 = out$par
log.lik1 = ll.jRCI.b1.A(par0, yi=yi, ind.lst=ind.lst, X=X, twosex=twosex, sex=sex,
ni=ni, ni0=ni0, xs=xs, iphi=iphi_i, theta=theta_i)
older = c(par0[seq], zeroes, par0[-seq], log.lik1, iphi_i, theta_i)
iphi_i = optimize(f=ll.tRCI.iphi.A, c(tech.ctrl$iphi_l, tech.ctrl$iphi_u),
yi=yi, ind.lst=ind.lst, X=X, twosex=twosex, betas=older[np])$minimum
theta_i = optimize(f=ll.aRC.theta.A, c(tech.ctrl$theta_l, tech.ctrl$theta_u),
bs1=older[bsi1], bs2=older[bsi2], ni=ni, ni0=ni0, twosex=twosex, sex=sex, xs=xs)$minimum
log.lik1 = ll.jRCI.b1.A(par0, yi=yi, ind.lst=ind.lst, X=X, twosex=twosex, sex=sex,
ni=ni, ni0=ni0, xs=xs, iphi=iphi_i, theta=theta_i)
newer = c(par0[seq], zeroes, par0[-seq], log.lik1, iphi_i, theta_i)
if((log.lik0 - log.lik1)<eps)break;
log.lik0 = log.lik1
}
if(log.lik0 < log.lik1){
ret = older
}else{
ret = newer
}
0
}, error=function(e) {
1
})
if((tag == 0) | (triali >= tech.ctrl$maxtrial))break;
par0 = rnorm(length(par0), start, 1)
}
if(tag == 1){
ret = NULL
}
return(ret)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.