Nothing
VarBandHazEst<-
#Adaptive variable bandwidth hazard rate estimator
function(xin, xout, kfun, h1, h2, ci)
{
n<- length(xin)
nn<-length(xout)
n1<-1:n
xin.use<-sort(xin)
pilot.est<-HazardRateEst(xin.use, xin.use, kfun, h1, ci) #Pilot estimate (W-L)
sqrt.pilot<-sqrt(pilot.est)
vhre<-sapply(1:nn, function(i, xin.use, xout, h2, kfun, n, n1,ci, sqrt.pilot)
{
sum( (sqrt.pilot * kfun( ((xout[i]-xin.use)/h2) * sqrt.pilot) )*ci/(n-n1+1))
}, xin.use, xout, h2, kfun, n, n1,ci, sqrt.pilot)
vhre/h2
}
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.