0

R / S / 非線形回帰に関連する問題があり、私は R プログラマーではないので、助けが必要です。

tt と td の 2 つの配列があります。

非線形関数の最小二乗和が最小になるように、パラメーター a、b、c を見つける必要があります。

td / tt - a * exp( b * tt ) + c

これを行う方法がわかりません。nls()関数を試しnls2()ましたが、運がありませんでした...

前もって感謝します。

編集：

私のデータ:

td <-as.array(0.2, 0.4, 0.8, 1.5, 3);

tt <-as.array(0.016, 0.036, 0.0777, 0.171, 0.294);

以下の回答の方法を使用すると、ランダム データの値は問題ありませんが、使用しているデータはモデル メッセージの評価時に生成された欠損値または無限大を返します。

すぐにデータを提供できなくて申し訳ありません。

4

1 に答える 1

4

あなたのデータ：

n <- 100
td <- runif(n)
tt <- runif(n)
data <- data.frame(td = td, tt = tt)

関数の構成結果

a <- 0.5
b <- 2
c <- 5
y <- jitter(td / tt - a * exp( b * tt ) + c)

(実際には、a、b、c が何なのかは後でわかります。ここでは、それらを使用して答えと比較します。)

フィッティング:

nls(
  y ~ td / tt - a * exp( b * tt ) + c, 
  data = data, 
  start = list(a = 1, b = 1, c = 1)
)

答え：

Nonlinear regression model
  model:  y ~ td/tt - a * exp(b * tt) + c 
   data:  data 
     a      b      c 
0.4996 2.0008 4.9994 
 residual sum-of-squares: 0.0001375

Number of iterations to convergence: 7 
Achieved convergence tolerance: 1.604e-06
于 2011-07-27T13:24:59.840 に答える