hello I want to confirm two code is the same:. I am clustering standard errors by "parents".
#1 using PLM from PLM package
plm(kids ~ treatment, data = dataset_x, index = c("love", "parents", "marriage"),
model = "within", effect = "twoways",
vcov = function(x) vcovHC(x, cluster = "parents"))
#2 using FEOLS from fixest package
feols(data = dataset_x, kids ~ treatment | love + parents + marriage,
cluster = "parents")
is this the same code? am I missing something? thank you all. I have 3 fixed effects(I want all 3 suppressed)
thank you! <3 <3