Displaying 1 result from an estimated 1 matches for "out33".
Did you mean:
out3
2009 Oct 30
2
Efficient way to code using optim()
...for (i in 2:length(realized_prob)){
realized_prob[i] = exp(theta_multiple[cum_count[(i-1)]+1]) / sum(exp(theta_multiple[((cum_count[(i-1)]+1):cum_count[i])]))
}
-sum(log(realized_prob))
}
initial = rep(0,32)
out33 = optim(initial, simple_logit, method="BFGS", hessian=TRUE)
Many thanks in advance!!!
_________________________________________________________________
[[alternative HTML version deleted]]