Displaying 1 result from an estimated 1 matches for "g_noise".
Did you mean:
_noise
2011 Oct 07
1
loess question
...o use loess to smooth a 2D image, and also obtain the
standard error for every pixel. I see that the standard error does not
make sense. For example, running the following:
library(stats)
x <- array(c(1:100), dim=c(100,100))
y <- t(x)
v <- exp(-((x-50)^2+(y-50)^2)/30^2)
s <- v*0.02
g_noise <- rnorm(10000, mean = 0, sd = s)
f <- v + g_noise
f.loess <- loess(f ~ x + y, span=0.1, data.frame(x=c(x),y=c(y),f=c(f)))
f.predict <- predict(f.loess, data = data.frame(x = c(x), y = c(y), f =
c(f)), span = 0.1,se=TRUE)
image(1:100,1:100,matrix(f.predict$se,nrow=100))
I get an image...