我正在尝试使用nlme R软件包拟合具有多个随机效应的非线性模型。该模型具有3个参数k1,k2和k3,对于级别,它具有YEAR
和COUNTRY
,其中YEAR
嵌套在COUNTRY
中。
y
是因变量,人口和月份是常数。
这是我尝试过的
mymodel <- nlme(y ~ population*k1*(1 - exp(-(months/k2)^k3))
, fixed = list(k1~1, k2~1, k3 ~ 1)
, random = list(k1 ~1 | COUNTRY/YEAR, k2 ~ 1 | COUNTRY, k3 ~ 1 | COUNTRY)
, weights = varPower(fixed=0.5)
, data=mydata, start = start.vals)
运行此命令会给我这个错误。
解析时出错(文本= paste(“〜”,paste(nVal,折叠=“ /”))): :2:0:输入意外结束 1:〜 ^
我在这里怎么可能做错了?