Как написать формулу линейной модели со 100 переменными в R

Есть ли простой способ в R создать линейную регрессию по модели с 100 параметрами в R? Допустим, у нас есть вектор Y с 10 значениями и фрейм данных X с 10 столбцами и 100 строками. В математической записи я бы написал Y = X[[1]] + X[[2]] + ... + X[[100]]. Как мне написать что-то подобное в синтаксисе R?

Попробуй это

df<-data.frame(y=rnorm(10),x1=rnorm(10),x2=rnorm(10))
lm(y~.,df)

Отличные ответы!

Я бы добавил, что по умолчанию вызов formulaa data.frameсоздает аддитивную формулу для регрессии первого столбца на другие.

Так что в случае ответа @ danas.zuokas вы можете даже сделать

lm(df)

который интерпретируется правильно.

$Y$$n$$X$$n \times p$$p=100$

Z = as.data.frame(cbind(Y,X))
lm(Y ~ .,data=Z)

Если есть другие столбцы, которые вы не хотите включать в качестве предикторов, вам придется удалить их, Xпрежде чем использовать этот трюк или использовать -в формуле модели, чтобы исключить их. Например, если вы хотите исключить 67-й предиктор (имеющий соответствующее имя x67), вы можете написать

lm(Y ~ .-x67,data=Z)

Кроме того, если вы хотите включить взаимодействия и т. Д., Вам нужно будет добавить их вручную как (например)

lm(Y ~ .+X[,1]*X[,2],data=Z)

или убедитесь, что они введены в виде столбцов X.

Вы также можете использовать комбинацию из formulaи pasteфункций.

Данные настройки : Давайте представим , что мы быть_наст data.frame, содержащие предикторы x1к x100и наши зависимые переменных y, но есть также неприятность переменный asdfasdf. Кроме того, переменные предиктора расположены в таком порядке, что они не являются непрерывными в data.frame.

Data <- data.frame(matrix(rnorm(102 * 200), ncol=102))
names(Data) <- c(paste("x", 1:50, sep=""), 
    "asdfasdf", "y", paste("x", 51:100, sep=""))

Представьте также, что у вас есть строка, содержащая имена переменных-предикторов. В этом случае это можно легко создать с помощью pasteфункции, но в других ситуациях grepили для получения этой строки можно использовать другой подход.

PredictorVariables <- paste("x", 1:100, sep="")

Применить подход : мы можем затем построить формулу следующим образом:

Formula <- formula(paste("y ~ ", 
     paste(PredictorVariables, collapse=" + ")))
lm(Formula, Data)

• то collapseаргумент вставка +между предикторами
• formulaпреобразует строку в объект формулы класса, подходящей для lmфункции.

В более общем смысле я использую следующую функцию довольно регулярно, когда я хочу предоставить переменные предиктора в качестве вектора имен переменных.

regression <- function(dv, ivs, data) {
    # run a linear model with text arguments for dv and ivs
    iv_string <- paste(ivs, collapse=" + ")
    regression_formula <- as.formula(paste(dv, iv_string, sep=" ~ "))
    lm(regression_formula, data)

}

Например,

regression("y", PredictorVariables, Data)