在指數衰減中設置截距

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))