my.data <- filter(mydf, CTRY == 'PL')
my.fit <- cfa(cfa.data, data = my.data, meanstructure = TRUE)
summary(my.fit, standardized= TRUE, rsquare = TRUE)