AlgDesign中的quad函数是用于构建二次模型的,但是该函数不是AlgDesign自带的。因此,需要安装quadprog包来提供quad函数。
示例代码:
# 安装quadprog包
install.packages("quadprog")
# 加载AlgDesign包和quadprog包
library(AlgDesign)
library(quadprog)
# 构建二次模型
x1 <- runif(10, 0, 1)
x2 <- runif(10, 0, 1)
design <- data.frame(x1, x2)
model <- lm(y ~ . + .^2, data = design)
# 生成D-最优设计
opt <- optFederov(~ . + .^2, data = design, center = TRUE, max.iter = 1000, method = "sequential")
# 查看结果
opt