[R] Question re: interpolation

David Winsemius dwinsemius at comcast.net
Sun May 2 23:44:03 CEST 2010


On May 2, 2010, at 5:27 PM, Jessica Schedlbauer wrote:

> Hello,
>
> I have a matrix in which two variables, x and y, are used together  
> to determine z.  The variables x and y are sorted into classes.   
> Specifically, values for variable x range from 0 to 2.7 and are  
> sorted into class increments of 0.15 and variable y ranges from  
> 0-2100 with class increments of 100.  For each x-y class  
> combination, I have calculated a mean value of z from existing  
> data.  However there are gaps in the existing data (i.e., x-y  
> combinations for which there are no z data).  I would like to  
> interpolate values to fill in these missing datapoints, but have so  
> far been unable to find a straightforward way of doing so in R.   
> There are many gaps in the data and I need to repeat this process  
> for several individual datasets.  That said, I am looking for a way  
> to avoid interpolating each missing point separately.

Why not do a loess fit _before_ you do any binning and then  
predict.loess at the intervals where you want estimates?

x=rnorm(200, 1.5)
  y= (10*x)^2 +rnorm(200, 0, 100)
  plot(x,y, cex=0.5, xlim=c(-2, 5))
  points(seq(-2,5, by=0.15), predict(loess(y~x),
                      newdata=data.frame(x=seq(-2,5, by=0.15))),  
col="red", cex=2)
>
> Any suggestions would be greatly appreciated.
>
> Regards,
> Jessica Schedlbauer
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius, MD
West Hartford, CT



More information about the R-help mailing list