search for: tricubic

Hi I am trying to get the tricube weights from the loess outputs as I need to calculate an error function which requires the weight. So I have used the following example from the R: cars.lo <- loess(dist ~ speed, cars, span=0.5, degree=1, family="symmetric") Then i try to get the weights: cars.lo$weights [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

...near polynomial fit is + \file{src/library/stats/src/lowess.doc}. Normally a local linear polynomial fit is used, but under some circumstances (see the file) a local constant fit can be used. \sQuote{Local} is defined by the distance to the \code{floor(f*n)}th nearest neighbour, and tricubic weighting is used

Hello, I need help with the details of loess prediction algorithm. I would like to get it implemented as a part of a measurement system programmed in LabView. My job is provide a detailed description of the algorithm. This is a simple one-dimensional problem - smoothing an (x, y) data set. I found quite a detailed description of the fitting procedure in the "white book". It is also

...stance from x (with differences in `parametric' variables being ignored when computing the distance). The size of the neighbourhood is controlled by &alpha; (set by span or enp.target). For &alpha; < 1, the neighbourhood includes proportion &alpha; of the points, and these have tricubic weighting (proportional to (1 - (dist/maxdist)^3)^3. For &alpha; > 1, all points are used, with the `maximum distance' assumed to be &alpha; times the actual maximum distance. </SNIP> I understand the situation when alpha is less than one, presuming it uses the same conven...

...stance from x (with differences in `parametric' variables being ignored when computing the distance). The size of the neighbourhood is controlled by &alpha; (set by span or enp.target). For &alpha; < 1, the neighbourhood includes proportion &alpha; of the points, and these have tricubic weighting (proportional to (1 - (dist/maxdist)^3)^3. For &alpha; > 1, all points are used, with the `maximum distance' assumed to be &alpha; times the actual maximum distance. </SNIP> I understand the situation when alpha is less than one, presuming it uses the same conven...

...;m looking at is more complex so I need a small span value to allow for a close fit. I foresee that changing the "weigthing" is the way to go but I do not really understand how the "weight" option is used (I tried to change it and nothing happened), and also the embedded "tricubic weighting" does not seem changeable. So any idea would be very welcome. Emmanuel

For each element in w I want to find a good match (subscript number) of an element in x. x and w can be long. Instead of just finding the closest match I want to use weighted multinomial sampling (which I've already figured out once I have the probabilities) where the probabilities come from the tricube function of absolute differences between donor and target values, but normalized to sum

SAS proc rank has ties options of high and low that would allow producing ranks of the type found in the sports pages, e.g., rank (c(1,1,2,2,2,2,3)) == 1 1 3 3 3 3 7 Could R support these ties.methods?

...nfidence intervals for a loess smoother (function: loess()), but have been thus far unsuccessful. The CI for the loess predicted values, yhat, are apparently yhat +- t*s * sqrt(w^2), where s is the residual sum of squares and w is the weight function Correct me of I'm wrong, but R uses the tricubic function (1-abs(z)^3)^3, where z = (x-xi)/h, where h is the the number of neighboring values within the bandwidth. Assuming that's correct, here is my code (for the first observation in x1: x1 <- 1:156 # x set.seed(123) y <- arima.sim(list(order=c(2,0,0), ar=c(1,-.1)), n=...

I received some great ideas (see below) from a number of people to whom I am grateful. Here is the code I put together from many of their suggestions: lx <- length(x) lw <- length(w) z <- matrix(abs( rep( x , lw ) - rep( w, each = lx ) ), nrow=lw, ncol=lx, byrow=TRUE) s <- pmax( abs( w - min(x) ), abs( w - max(x) ) ) z <- (1 - (z/rep(s,length=lx*lw))^3)^3 sums <-