Add function for computing cluster standard errors for linear FE panel data model