13

最初に言いたいのは、非線形回帰の R^2 値を計算することは正確ではない、または有効なことではないことを理解しているということです。

しかし、私は SigmaPlot での作業のほとんどを R に移行する過渡期にあり、非線形 (濃度-反応) モデルの場合、同僚はモデルに関連付けられた R^2 値を見て良さを推定することに慣れています。 -適合。

SigmaPlot は 1-(残余 SS/総 SS) を使用して R^2 を計算しますが、RI では総 SS を抽出できないようです (残余 SS は要約で報告されます)。

これを機能させるための助けをいただければ幸いです。適合度のより良い推定器を使用するように私たちを動かそうとしています。

乾杯。

4

2 に答える 2

12

総平方和を抽出する代わりに、私はそれらを計算しました：

test.mdl <- nls(ctrl.adj~a/(1((conc.calc/x0)^b)),
                data=dataSet,
                start=list(a=100,b=10,x0=40), trace=T);

1 - (deviance(test.mdl)/sum((ctrl.adj-mean(ctrl.adj))^2))

SigmaPlotを使用した場合と同じR^2が得られるので、すべて問題ないはずです。

于 2011-04-14T21:51:36.963 に答える
5

したがって、y の合計変動は (n-1)*var(y) のようになり、比率は説明されていません。私のモデルは次のsum(residuals(fit)^2)ようになります。1-(sum(residuals(fit)^2)/((n-1)*var(y)) )

于 2011-04-14T21:04:07.493 に答える