次のように MA 正規化 ( http://en.wikipedia.org/wiki/MA_plot )を使用して、2 つの変数x
とy
Rでレスを実行しようとしています。
> x = rnorm(100) + 5
> y = x + 0.6 + rnorm(100)*0.8
> m = log2(x/y)
> a = 0.5*log(x*y)
標準の MA 正規化のように、平均m
が 0 になるように x と y を正規化し、正しい x と y の値を逆算します。MA で最初にレスを実行する:
> l = loess(m ~ a)
補正m
値を取得する方法は何ですか?これは正しいです?
> mc <- predict(l, a)
# original MA plot
> plot(a,m)
# corrected MA plot
> plot(a,m-mc)
predict
オブジェクトの場合に実際に何をするのか、返されたオブジェクトloess
で使用するのとどう違うのか、私にはわかりません-誰かが説明できますか? 最後に、この修正に基づいて新しい値と値を逆算するにはどうすればよいですか?l$residuals
l
loess
x
y