re_fe_with_hidden_collider

library(jtools)
library(plm)

i <- rep(1:100, each = 100)
v <- rep(rnorm(100), each = 100)
a <- rnorm(10000)
t <- a + v + rnorm(10000)
y <- a + v + t + rnorm(10000)
u <- t + y + rnorm(10000)

df <- data.frame(i, y, t, a, v, u)

ols <- lm(y ~ t + a, data = df)
re <- plm(y ~ t + a, data = df, model = "random", index = "i")
fe <- plm(y ~ t + a, data = df, model = "within", index = "i")

export_summs(ols, re, fe)