-6

(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 )
}
4

1 に答える 1