Displaying 2 results from an estimated 2 matches for "wfun".
Did you mean:
fun
2007 Nov 21
1
equivalent of Matlab robustfit?
...rithm uses iteratively
reweighted least squares with the bisquare weighting function. By
default, ROBUSTFIT adds a column of ones to X, corresponding to a
constant term in the first element of B. Do not enter a column of ones
directly into the X matrix.
B = ROBUSTFIT(X,Y,'WFUN',TUNE) uses the weighting function 'WFUN' and
tuning constant TUNE. 'WFUN' can be any of 'andrews', 'bisquare',
'cauchy', 'fair', 'huber', 'logistic', 'talwar', or 'welsch'.
Alternatively 'WFUN'...
2009 Nov 30
3
bug or bizarre feature?
....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, it is only 115x115.
> udim
[1] 116
> length(uvec)
[1] 116...