Nothing
`pomdev.ts` <-
function(object1,object2,eps=10^-30,nrange=1000)
{
if(!is.numeric(object1)|!is.numeric(object2))
stop("objects must be numeric ")
if(any(!is.finite(object1))|any(!is.finite(object2)))
stop("objects must have finite values ")
if(!(nrow(object2)==length(object1) & ncol(object2)>2))
stop("simulation data must have the form of a table (matrix) of nrow=length(field data) and ncol>2")
#fix a band width of the kernel according to the temporal distribution of the field data
r1<-range(object1)
r<-extendrange(r1,f=0.01)
bw1<-bw.nrd0(object1)
x<-0
for(rowi in 1:nrow(object2))
{
d2<-density(object2[rowi,],bw=bw1,from=r[1],to=r[2],n=nrange)
o2<-d2$y
w <- o2 < eps
if (any(w)) o2[w] <- eps
pdf2 <- approxfun( d2$x, o2, yleft=eps, yright=eps)
x <- x - 2 * log(pdf2(object1[rowi]))
}
x
}
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.