Nothing
elltheta <-
function(theta, p0, r, tol=1e-7, maxit=500){
n<-length(r)
del<-1
it<-0
p<-p0
while(del>tol && it<maxit){
res<-H.Binv(theta,p, r)
p.nu<-p-res$Binv%*%res$H
del<-sum(abs(res$H))+sum(abs(p-p.nu))
p<-p.nu
#cat(it, del, "theta=", theta, "sum of p=", sum(p), "\n")
it<-it+1
}
ell<-res$ell
Fisher<-res$Fisher
list(ell=ell, p=p)
}
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.