[R] getting the smoother matrix from smooth.spline
Simon Wood
sw283 at maths.bath.ac.uk
Fri Jun 30 18:07:40 CEST 2006
smooth.matrix = function(x, df){
n = length(x);
A = matrix(0, n, n);
for(i in 1:n){
y = rep(0, n); y[i]=1;
yi = predict(smooth.spline(x, y, df=df),x)$y;
A[,i]= yi;
}
(A+t(A))/2;
}
>- Simon Wood, Mathematical Sciences, University of Bath, Bath BA2 7AY
>- +44 (0)1225 386603 www.maths.bath.ac.uk/~sw283/
On Sat, 24 Jun 2006, Gregory Gentlemen wrote:
> Can anyone tell me the trick for obtaining the smoother matrix from smooth.spline when there are non-unique values for x. I have the following code but, of course, it only works when all values of x are unique.
>
> ## get the smoother matrix (x having unique values
> smooth.matrix = function(x, df){
> n = length(x);
> A = matrix(0, n, n);
> for(i in 1:n){
> y = rep(0, n); y[i]=1;
> yi = smooth.spline(x, y, df=df)$y;
> A[,i]= yi;
> }
> (A+t(A))/2;
> }
>
>
> Thanks for any assistance,
> Gregory
>
>
> ---------------------------------
>
> ---------------------------------
> Get a sneak peak at messages with a handy reading pane.
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
More information about the R-help
mailing list