R中的三变量线性模型
提问于 2022-01-21 07:45:33
我想得到一个与不同化学物质之间的协同/拮抗作用相关的线性模型的系数。
化学品X，Y，Z.系数b0...b7.
0= b0 + b1x + b2y + b3z + b4xy + b5xz + b6yz + b7xyz
某些组合的X，Y，Z将杀死50%的细胞(b0代表这种总有效性)，而高阶项的符号/大小代表着化学物质之间的相互作用。
考虑到实际的数据点，我想要适应这个模型。
编辑:我已经摆脱了琐碎的解决方案，在开始时添加了一个强制值。测试数据：
x1 <- c(0,1,2,3,4)
y1 <- c(0,2,1,5,4)
z1 <- c(0,1,-0.66667,-6,-7.25)
q <- c(-1,0,0,0,0)
model <- lm(q ~ x1*y1*z1)该集合具有系数：-30，12，6，4，1，-1，0，0.5
编辑:从以前取得的进展，将尝试更多的数据点。前四个系数看起来很好(乘以30)：
Coefficients:
(Intercept) x1 y1 z1 x1:y1 x1:z1 y1:z1 x1:y1:z1
-1.00000 0.47826 0.24943 0.13730 -0.05721 NA NA NA 编辑:添加更多的数据点到目前为止还没有成功，不确定是否需要一定的最低数量才能准确。
我把事情安排好了吗？一旦我有了系数，我想要解z，这样我就可以绘制一个三维曲面。谢谢!
回答 1
Stack Overflow用户
回答已采纳
发布于 2022-01-21 11:44:12
我只使用了16个任意数据点就可以得到系数，并附加了一个点来排除琐碎的答案：
x1 <- c(0,1,2,3,4,1,2,3,4,1,2,3,4,5,6,7,8)
y1 <- c(0,2,1,5,4,3,7,5,8,6,2,1,5,5,3,5,7)
z1 <- c(0,1,-0.66667,-6,-7.25,-0.66667,-5.55556,-6,-6.125,-4,-2.5,-6,-6.8,-7.3913,-11.1429,-8.2069,-6.83333)
q <- c(-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)
model <- lm(q ~ x1*y1*z1)页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/70803964
复制相关文章
腾讯云开发者
