r - 定义三向非线性交互的峰值

Question

我们有一个显着的三向二次交互（所有 3 个连续变量），现在想定义“临界点”。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 之间的哪个交叉点出现峰值/谷值。有什么建议么？