非线性泊松回归纠正代码公式的问题

目前有一组数据为一个生物模型数据，主要牵涉到剩余物种对应这每个岛屿的面积（area），以及度过的时间 t的关系。

目前所出现的问题：

2.计算非线性泊松模型回归并使用公式。

目前状况：area对应的parameter目前代码无法有效展示

把公式转换成可以处理的方式：

我们看到公式 其中 S_t是剩余物种数量，S_infinit ，c,a,k为parameter，a是岛屿面积，t是时间。

                              

如图所示 最后一行为需要表达的公式，代码在r语言中表现如下：

第46行 文件

选择1.

fit2_1 = gnm(richness~-1+ #此处为扣除beta_0

                Mult(1  #此处为s_infinito

, Const(1)         #此处为括号里面第一个数字1

+Mult(Mult(1,Exp(log(area)))#此处有问题 。表现的是we^{z*Ina},问题是z是parameter对应着Ina而不是a，所以此处的所用代码使用的是对应的 area而不是Ina，

                  +Const(-1) #此处为-1数字

,Exp(year)))   #此处为表达e^{-kt}

,family=poisson, data = x)

选择2

fit2_1 = gnm(richness~-1+

              Mult(1, Const(1)+Mult(Mult(1,Exp(log(area)))

                                    +Const(-1),Exp(year))),family=poisson, data = x)

选择3

fit2_1 = gnm(richness~-1+

              Mult(1, Const(1)+Mult(Mult(1,Exp(Exp(area)))

                                    +Const(-1),Exp(year))),family=poisson, data = x)

选择4

fit2_1 = gnm(richness~-1+

              Mult(1, Const(1)+Mult(Mult(1,Exp(Mult(1,offset(Exp(area)))))

                                    +Const(-1),Exp(year))),family=poisson, data = x)

标黄部分为所需修正的代码，需要表现的是we^{z*Ina}，z是parameter对应着是 Ina 而不是a，所以需要修正代码。

期望大神可以帮忙改一下代码。

建议你简化一下问题，举一个简单的例子，不然很难理解你的需求

希望可以添加一列 作为 Exp(area)
目前的代码是
Inarea <- Exp(x$area)
x1 <- cbind(x,Inarea)
但是出现的error是：
Error in as.data.frame.default(x[[i]], optional = TRUE) :
  cannot coerce class ‘"expression"’ to a data.frame
该如何处理。