r - 设置指数衰减截距

Question

我正在尝试将我的数据拟合到指数模型

y=exp(a-b*x)

x=0 中的y 的值必须正好是 13.75，所以我想在拟合数据之前设置截距。

a=log(13.75)

我正在使用函数 gnm。我读过我应该使用 constrain 或 constrainTo 来设置系数；但是，我不明白该怎么做。

x=seq(0,3120,120)

y= c(13.75, 13.74, 13.73, 13.71, 13.70, 13.69, 13.68,
     13.66, 13.65,13.64, 13.64, 13.63, 13.63, 13.62,
     13.62, 13.61, 13.61, 13.60,13.60, 13.60, 13.60,
     13.59, 13.59, 13.59, 13.59, 13.59, 13.59)

xl<-gnm(y ~ Exp(1 - x), verbose = FALSE,constrain=????)

#

大家好。非常感谢你。我一直在按照你的建议工作，我终于找到了一个非常合适的人选。现在我将使用重量来抛光它。我这样做了:

xl<-gnm(y ~ Exp(1 + x) -1, verbose = FALSE,
constrain="Exp(.+x).Intercept",constrainTo=a,start=c(a,-0.05),
prior.weights=W,weights=W)   

还有一个大问题。我需要从 y 中减去一个数字，以便将我的数据的水平渐近线置于零。

Answer 1

在 Exp 中使用 -1 告诉模式不要在指数中设置截距项。尝试使用这个模型:Exp(1 + x) -1，删除实际的 y 截距。

要设置指数截距，你可以尝试:

xl<-gnm(y ~ Exp(1 + x) -1, verbose = FALSE,constrain="Exp(.+x).Intercept",constrainTo=a,start=c(a,-0.05))
xl

必须为要估计的模型输入起始参数要交互式地选择约束参数，您还可以这样做:

xl<-gnm(y ~ Exp(1 + x) -1, verbose = FALSE,constrain="[?]",constrainTo=a,start=c(a,-0.05))