using FEOLS vs PLM comparison

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

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.