Допустим, у меня есть такой набор данных:
library(dplyr)
library(MatchIt)
set.seed(1)
df <- data_frame(outcome = sample(c(rep("TRUE", 50), rep("FALSE", 50))),
age = rnorm(100, mean = 35, sd = 15),
gender = sample(c(rep("MALE", 50), rep("FEMALE", 50))),
var1 = rnorm(100, mean = 1, sd = 0.3),
var2 = rnorm(100, mean = 1000, sd = 125),
var3 = rnorm(100, mean = 0, sd = 300))
Я хочу контролировать возраст и пол и определить влияние вариаций 1, 2 и 3 на результат. Я полагаю, что могу использовать функцию MatchIt для контроля возраста и пола следующим образом:
match_it <- matchit(formula = outcome ~ age + gender,
data = df,
method = "nearest")
Но как только это будет сделано, как я могу определить влияние var1 на результат? Спасибо за любую помощь.