(0、10)上の100の等しいポイントでP(以下を参照)を使用して、XとY(データmから)の関係を推定したいと思います。Rを使用してこれを行うにはどうすればよいですか?
> m
X Y
1 0.5 6.0
2 1.5 5.0
3 2.5 6.7
4 3.5 7.1
5 4.5 6.1
6 5.5 8.1
7 6.5 8.0
8 7.5 7.9
9 8.5 8.0
10 9.5 6.0
(1)
P=function(x,X,Y,sigma){
# x = point to evaluate our estimate.
# X = vector of observation X values
# Y = vector of observation Y values
# sigma = standard deviation.
weights = rep(0,length(X))
sumweights = 0
smooth = 0
for(i in 1:length(X)){
weights[i] = dnorm(x,mean=X[i],sd=sigma)
sumweights = sumweights + weights[i]
smooth = smooth + weights[i]*Y[i]
}
return( smooth/sumweights )
}