0

我正在尝试实现一个新的非线性函数以在 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) )

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

4

0 回答 0