Displaying 1 result from an estimated 1 matches for "m2b2".
Did you mean:
a2b2
2006 Jan 02
0
boostrap astronomy problem
...Function
glm=function(x){
e<-x[1]
k<-x[2]
c(1/e + (k+1)*(n/(e-t)) - k*sum(1/(e-s[b])),-n/(k+1) +
n*log(e-t) - sum(log(e-s[b])))
}
optim(c(480.,2.),lm,glm,method="BFGS",control=list(maxit=10000000))$par
}
#Compute Bootstrap replicates of escape velocity and kr
m2B2=boot(vg,mystat,5000)
------------------------------------------------------------------------
Does this appear to be correct for what I'd like to achieve? I have
looked at the distribution and it appears to be about Normal, but can I
say that this is true for the sampling distribution as well?...