Question

我想从矢量（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路线，但这并没有导入数据。

Answer 1

你可以做到

model  <- function(x = IC50, y = descs) 
  lm(x ~ ., data = y)

Data.frame列中的公式

1 个答案: