我目前正在学习使用R进行响应面建模。我成功地使用了带有函数ccd
的rsm
包来创建中心复合材料设计。现在我的问题是,我如何将我收集的数据添加到其中,即在进行实验时获得的响应变量的实际数据?我想用我的数据填充下例中的y
,这些数据在数据帧中作为数字。
这是我的代码:
library(rsm)
ccd <- ccd(y ~ x1+x2+x3, n0 = c(0,3), alpha = "rotatable", randomize = F, oneblock = T
ccd
run.order std.order x1.as.is x2.as.is x3.as.is y
1 1 1 -1.000000 -1.000000 -1.000000 NA
2 2 2 1.000000 -1.000000 -1.000000 NA
3 3 3 -1.000000 1.000000 -1.000000 NA
4 4 4 1.000000 1.000000 -1.000000 NA
5 5 5 -1.000000 -1.000000 1.000000 NA
6 6 6 1.000000 -1.000000 1.000000 NA
7 7 7 -1.000000 1.000000 1.000000 NA
8 8 8 1.000000 1.000000 1.000000 NA
9 1 1 -1.681793 0.000000 0.000000 NA
10 2 2 1.681793 0.000000 0.000000 NA
11 3 3 0.000000 -1.681793 0.000000 NA
12 4 4 0.000000 1.681793 0.000000 NA
13 5 5 0.000000 0.000000 -1.681793 NA
14 6 6 0.000000 0.000000 1.681793 NA
15 7 7 0.000000 0.000000 0.000000 NA
16 8 8 0.000000 0.000000 0.000000 NA
17 9 9 0.000000 0.000000 0.000000 NA
Data are stored in coded form using these coding formulas ...
x1 ~ x1.as.is
x2 ~ x2.as.is
x3 ~ x3.as.is
事实证明,最简单的解决方案往往是正确的。将数据存储在相同长度的数据帧df
中,使用简单添加
ccd$y <- df$y