Let\'s consider data following :
library(pglm) data(\'UnionWage\', package = \'pglm\') df1 <- data.frame(UnionWage[1:2],\'wage\' = UnionWage$wage, \'exper\