Displaying 1 result from an estimated 1 matches for "db1db1".
Did you mean:
dbodb1
2013 Mar 10
0
Steepest Ascent Algorithm
...2] = hess[2,1] = sum(-log(t)*mu) #dadB0
hess[2,2] = sum(-mu) #dBodBo
hess[1,3] = hess[3,1] = sum(-d*log(t)*mu) #dadB1
hess[2,3] = hess[3,2] = sum(-d*mu) #dBodB1
hess[3,3] = sum(-mu*d^2) #dB1dB1
}
list(l=l,grad=if(gcomp) grad,hess=if(hesscomp)hess)
}
#Steepest ascent algorithm----------------------------------
ascent2 <- function(data,maxit,theta) {
for (it in 1:maxit){
a =like(theta,data,gcomp=T,hesscomp=F)
theta1 = theta + a$grad # Steepest Ascent
atmp = like(theta...