Displaying 1 result from an estimated 1 matches for "xeval".
Did you mean:
eval
2001 Nov 05
1
Why doesn't outer work?
...<- p.vec[11];
fac1<-sqrt(2*pi)*sqrt(kidsize.var);
fac2<-((y-kidsize.mean)^2)/(2*kidsize.var);
f<-fx(x)*nkids*exp(-fac2)/fac1;
return(f);
}
# Part (III)
################################################################
# THE KERNEL K(y,x)
# note the chopoff
Kyx<-function(y,x) {
xeval<-max(x,minsize); xeval<-min(xeval,maxsize);
yeval<-max(y,minsize); yeval<-min(yeval,maxsize);
#return(pxy(xeval,yeval)+fxy(xeval,yeval))
return(pxy(x,y)+fxy(x,y))
};
############## The 'big matrix' M of size n x n this uses outer and doesn't
work, don't...