在三向非线性相互作用中定义峰值

Question

我们有一个重要的三向二次交互（所有3个连续变量），现在想要定义＆＃34;引爆点＆＃34;。 Johnson-Neyman技术仅允许双向交互。有什么建议我们如何计算峰值？

这是一些随机数据：

set.seed(1)
mydata <- data.frame( SID=sample(1:150,400,replace=TRUE),age=sample(50:70,400,replace=TRUE), sex=sample(c("Male","Female"),200, replace=TRUE),time= seq(0.7, 6.2, length.out=400), AM =rnorm(mean=-0.8, sd=0.08, n=400),DMD=rnorm(mean=1.4, sd=0.2, n=400), MEC =rnorm(mean=0, sd=2, n=400))  
mydata$time <- as.numeric(mydata$time)

以下是二次三向交互的lme：

model <- lme(MEC ~ age+sex+AM*DMD*time+I((AM*DMD*time)^2), random=~time|SID, data=mydata); summary(model)

我现在想知道AM和DMD之间的哪个交叉点出现峰值/谷值。有什么建议？