2つの変数があり、どちらもクラス「数値」からのものです。
> head(y)
[1] 0.4651804 0.6185849 0.3766175 0.5489810 0.3695258 0.4002567
> head(x)
[1] 59.32820 68.46436 80.76974 132.90824 216.75995 153.25551
プロットしたので、指数モデルをデータに適合させたい(そしてそれをデータに追加したい)プロット)しかし、Rの多変量データにモデルを適合させるための情報が見つかりません!データを単変量にするためだけに、誰かが助けることができますか?どこから始めればいいのかわからない…ありがとう!
コメント
回答
あなたの用語がオフになっているため、あなたが何を求めているのか完全にはわかりません。しかし、変数が互いに独立していないと仮定します (もしそうなら、彼らは「見つけることとは関係ありません)私はそれを試してみます。 x
が独立(または予測)変数であり、y
が従属(または応答)変数である場合、これは機能するはずです。
# generate data beta <- 0.05 n <- 100 temp <- data.frame(y = exp(beta * seq(n)) + rnorm(n), x = seq(n)) # plot data plot(temp$x, temp$y) # fit non-linear model mod <- nls(y ~ exp(a + b * x), data = temp, start = list(a = 0, b = 0)) # add fitted curve lines(temp$x, predict(mod, list(x = temp$x)))
コメント
- ご回答ありがとうございます。 independent "ご指摘のとおり、'意味がありませんでした。データにコードを使用すると、モデルを適合させることができますが、結果として、グラフには1本ではなく数十本の線が表示されます。理由はわかりますか?
- @ sbg-いいえ、申し訳ありませんが、'理由がわかりません。
nls()
はモデルに適合しますか? - そう思います:非線形回帰モデルモデル:y〜exp(a + b * x)データ:DF ab -0.535834 -0.002024残差平方和:18.62収束までの反復回数:6達成された収束許容誤差:8.08e-06
- @sbg
x
変数:lines(sort(temp$x),predict(mod, list(x=sort(temp$x)))
独立した"変数(私は"予測子"ですが、'は重要ではありません)。"依存" / " response "変数?これらが両方とも応答変数である場合、パラメトリックな2変量確率分布(ありまたはなし)を当てはめることを想像できます。分布'のパラメーターが依存する予測変数)-または2Dカーネル密度推定。おそらくコンテキストをもう少し説明できます。(PSwhoever '質問に賛成する人は、それが何を意味するのかを知っている必要があります…誰かがチャイムを鳴らしますか?)