Displaying 1 result from an estimated 1 matches for "ckern".
Did you mean:
kern
2009 Nov 30
3
bug or bizarre feature?
...5)
> xdim <- 100
> mobility <- 0.1
> slope <- 1.16
> urange <- slope*xrange
> udim <- max(slope*xdim,5)
> du <- (urange[2]-urange[1])/udim
> uvec <- urange[1]+(1:udim-0.5)*du
> # type dependent weight function
> ckern <- array(0,dim=c(udim,udim))
> diag(ckern) = wfun(uvec,slope,mobility)
Error in `diag<-`(`*tmp*`, value = c(0.992300064325398, 0.990746129315703, :
replacement diagonal has wrong length
It turns out that the array ckern has the wrong size for some reason.
Instead of 116x116, i...