PISA_wynikiPISA

wynik_matematyka <- badanie %>%
          do(data.frame(PV1 = svymean(~PV1MATH,design = .$schemat),
                        PV2 = svymean(~PV2MATH,design = .$schemat),
                        PV3 = svymean(~PV3MATH,design = .$schemat),
                        PV4 = svymean(~PV4MATH,design = .$schemat),
                        PV5 = svymean(~PV5MATH,design = .$schemat))) %>%
          ungroup() %>%
          mutate(CNT = badanie$CNT) %>%
          gather(Stat,Value,-CNT) %>%
          separate(Stat,c('PV','STAT'),sep='\\.') %>%
          group_by(CNT,STAT) %>%
          summarise(PV = mean(Value),
                    SE_test = sd(Value)) %>%
          spread(STAT,PV) %>%
          mutate(SE_all = SE + 1.2*SE_test) %>%
          select(-SE_test,-SE) %>%
          group_by(CNT) %>%
          summarise_each(funs(m=mean(.,na.rm=T)),mean,SE_all) %>%
          rename(PV=mean,SE=SE_all)
wynik_matematyka