如何在R中的nlmer中实现自己的非线性函数？

问题描述

我正在尝试实现一个新的非线性函数，以在lme4软件包的nlmer函数中使用。但是我不确定是什么问题。这是我第一次尝试使用nlmer，但我会遵循在互联网上找到的所有说明。第一个错误是关于我的数据框。

  data <- read.csv(paste("C:/Users/oguz/Desktop/Runs4SiteModels/db/","DB4NLSiteModel",Periods[i],".txt",sep=""),sep = "",header = TRUE)
  
  psa_rock <- data$PSAr
  
  nparams <- c("c")
  nonl_fn <- deriv(~ log(( psa_rock + c)/c),namevec = c("c"),function.arg=c("c",psa_rock))
  
  fm <- nlmer(log(data$PSAm) ~ nonl_fn(c,psa_rock) ~ 1 + data$M1 + data$M3 + data$M85 + data$Nflag + data$Rflag + data$FDepth + 
              data$Dist1 + data$Dist3 + data$VN + (exp(-1*exp(2*log(data$Vs)- 11)) *  log((data$PSAr + c) / c) ) +
              (1|data$EQID) + (1|data$STID),data=data,start=c(c=0.1))

运行此代码时，出现以下错误：

Error in model.frame.default(data = data,drop.unused.levels = TRUE,formula = log(data$PSAm) ~  : 
  invalid type (list) for variable 'data'

在使用lmer函数（当然没有非线性函数）时我没有得到它。这就是为什么我认为我的问题与我的数据帧无关。

我不禁要思考的其他问题，即固定效应中的一部分：

(exp(-1*exp(2*log(data$Vs)- 11)) *  log((data$PSAr + c) / c) )

您可以看到我的非线性函数也参与了我的固定效应公式，但我不确定如何实现。我希望我的方法是正确的，但是由于我的第一个问题，我找不到机会进行测试。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）