估计R中的模型参数

Question

我有一组值对（x，y）我想让它们适合下面的模型，以便最小化平方误差。在R中最简单的方法是什么：

Y =（A *（N + 1-X）^ b）/ X ^ a，其中N是max（y），A，a，b是各种参数

以上是秩排序分布（见http://www.plosone.org/article/fetchObject.action?uri=info:doi/10.1371/journal.pone.0004791&representation=PDF）

任何帮助将不胜感激。

Answer 1

在我看来，你愿意做一个线性回归：log Y = b log A + b log (N+1-X) - a log X

在这种情况下，我会创建说YY = log Y，X1 = log(N+1-X)和X2 = log X并与lm(YY ~ X1 + X2)

一起使用