3

次のように MA 正規化 ( http://en.wikipedia.org/wiki/MA_plot )を使用して、2 つの変数xyRでレスを実行しようとしています。

> 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$residualslloessxy

4

1 に答える 1