Displaying 1 result from an estimated 1 matches for "dbodbo".
Did you mean:
dbodbc
2013 Mar 10
0
Steepest Ascent Algorithm
...sum(d*w-d*mu) #B
}
if(hesscomp==TRUE){
#Recall mu = (t^a)*exp(Bo+d*B1)
hess = matrix(0,3,3)
hess[1,1] = sum(-(mu*log(t)^2)-w/(a^2)) #dada
hess[1,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----------------------...