0

2 つの情報セット (残差) の違いをプロットしようとしています。両方のデータ セットは類似した (ただし異なる) 特性を持ち、両方のデータ セットは 0 から同じ X 値になります。唯一の矛盾は、それらのインデックスが異なることです。つまり、最初のグラフは A ステップで X に到達しますが、2 番目のグラフは B ステップで X に到達します。したがって、あるデータ フレームの従属変数の値を別のデータ フレームから単純に減算することはできません。私は非常に一般的な言葉で話しているので、簡単な例を示しました。次のような 2 つのデータ セット間の残差をプロットします。

data1 <- data.frame(x1=c(1,2,3,4,5,6), y1=c(10,5,7,3,2,4))
data2 <- data.frame(x2=c(1,3,6), y2=c(1,3,2))

plot(data1, y1 ~ x1, type = 'l', lty = 1, col = 'blue', xlim = c(1,6), ylim = c(0,10))
points(data2$y2 ~ data2$x2, type = 'l', lty = 1, col = 'red')

だから私は私の質問は次のとおりだと思います:

  1. インデックスが異なる 2 つの関数 (上記のような) の残差をプロットするにはどうすればよいですか。2 つのデータセット間の残差を解決する関数はありますか?

EDIT1: 例に問題がありました。Spacedman がこれを修正するのを手伝ってくれました。

4

1 に答える 1

1

線形補間で十分な場合はapprox、X 座標の束で補間するために使用できます。例えば:

> xout = sort(unique(c(seq(1,6,len=100),data1$x1,data2$x2))) # include data coords (untested)
> d1 = approx(data1$x1,data1$y1,xout)
> d2 = approx(data2$x2,data2$y2,xout)
> plot(xout,d1$y-d2$y,type="l")
于 2012-08-16T14:04:47.260 に答える