【发布时间】：2016-03-04 07:42:32
【问题描述】：
我想从一个向量 (IC50) 针对许多不同的分子描述符（A、B、C、D 等）创建一个回归模型。
我想用，
model <- lm (IC50 ~ A + B + C + D)
分子描述符位于 data.frame 的列中。我想使用一个函数，它将 IC50 向量和适当的子集 data.frame 作为输入。
我的问题是我无法将列转换为模型的公式。
谁能帮忙。
样本数据和微不足道的尝试，
IC50 <- c(0.1,0.2,0.55,0.63,0.005)
descs <- data.frame(A=c(0.002,0.2,0.654,0.851,0.654),
B=c(56,25,89,55,60),
C=c(0.005,0.006,0.004,0.009,0.007),
D=c(189,202,199,175,220))
model <- function(x=IC50,y=descs) {
a <- lm(x ~ y)
return(a)
}
我沿着substitute/deparse 路线走，但这并没有导入数据。
【问题讨论】：
-
a
-
除了您使用
x和y令人困惑，因为它与正常用法相反，为什么不简单地将单个向量作为列添加到具有已知列名，然后你的公式就是known_column_name ~ . -
@VenYao 你能做到这六十三栏吗？
-
@joran 我不知道点公式。谢谢。
标签： r