R/HR.R

HR <-
function(dat,t,h=NULL,alpha=1/5,bound=Inf){
	if (is.null(h)){
		h<-1/((length(dat))^(alpha))
		}
	b<-dat; c<-.InvYnBis(dat,max(dat)+1)
	sum( ((1/h)*.ker( (t-b)/h )*c)[b<bound] )
}

Try the EstSimPDMP package in your browser

Any scripts or data that you put into this service are public.

EstSimPDMP documentation built on May 2, 2019, 3:40 p.m.