> # basic model > model_basic <- nlme(wt ~ a(1+b*exp(-k*t))^m, + data=df_phen_group, + fixed=a+b+k+m ~ 1, + random=a+k ~ 1, + start=c(20, 10, 0.03, 0.5)) Error in a(1 + b * exp(-k * t)) : could not find function "a"

在这段代码中，a(1+b*exp(-k*t))^m 并不是一个合法的R表达式，因此R会报错。如果你想使用幂函数的形式来拟合生长曲线，可以使用以下代码：

model_basic <- nlme(wt ~ a * (1 + b * exp(-k * t))^m, 
                     data=df_phen_group, 
                     fixed=a+b+k+m ~ 1, 
                     random=a+k ~ 1, 
                     start=c(20, 10, 0.03, 0.5))

在这里，我们将幂函数的形式改为了a*(1+b*exp(-k*t))^m的形式，并且在拟合时需要提供一个初始值m。这样就可以使用nlme函数来拟合基本模型了。