ddply summary issue

Slavek · January 30, 2019, 5:39pm

I have this data set;

data.frame(stringsAsFactors=FALSE,
             URN = c("9x013122510", "9x013680557", "9x014145627",
                     "9x015094359", "9x013122478", "9x012544645",
                     "9x013122457", "9x015094333", "9x014145611", "9x014963276",
                     "9x015104994", "9x013233865", "9x015094394", "9x013918536",
                     "9x013122546", "9x013122646", "9x013359154", "9x013122574",
                     "9x013122548", "9x015104981", "9x013918605", "9x013122618",
                     "9x013918581", "9x013122341", "9x013122547", "9x013122698",
                     "9x012654879", "9x013703410", "9x013122485", "9x014853449",
                     "9x013122486", "9x014032625", "9x013122697", "9x013122342",
                     "9x012545899", "9x014145624", "9x013233864", "9x013703438",
                     "9x013918532", "9x014635481", "9x012544643", "9x013918584",
                     "9x013122284", "9x012544539", "9x013680544", "9x013918572",
                     "9x013918545", "9x014032615", "9x014635460", "9x015094372"),
   InterviewDate = c("2018-06-04 16:42:00", "2018-05-21 08:30:00",
                     "2018-10-15 13:04:32", "2018-12-17 11:03:42",
                     "2018-06-04 21:25:00", "2017-03-20 20:54:00", "2018-06-05 10:23:00",
                     "2018-12-23 12:01:17", "2018-10-15 19:02:13",
                     "2018-12-10 09:02:05", "2018-12-22 11:01:14", "2018-05-26 15:10:00",
                     "2018-12-20 22:01:27", "2018-05-18 08:59:00", "2018-06-04 14:57:00",
                     "2018-05-31 21:02:00", "2018-05-25 16:05:00",
                     "2018-06-04 12:27:00", "2018-06-04 14:28:00", "2018-12-29 21:01:06",
                     "2018-09-25 14:01:32", "2018-06-03 22:41:00",
                     "2018-10-01 09:01:57", "2017-01-18 08:24:00", "2018-06-04 14:49:00",
                     "2018-05-28 13:46:00", "2017-01-09 11:02:00", "2018-05-20 09:19:00",
                     "2018-06-04 19:45:00", "2018-12-01 12:01:20",
                     "2018-06-04 19:14:00", "2018-10-08 08:02:09", "2018-05-28 14:14:00",
                     "2017-01-12 09:18:00", "2017-01-01 10:33:00",
                     "2018-10-15 13:04:32", "2018-05-26 16:21:00", "2018-05-18 19:16:00",
                     "2018-05-18 10:54:00", "2018-11-16 20:01:54", "2017-02-02 10:51:00",
                     "2018-09-25 12:02:14", "2017-01-15 14:32:00",
                     "2017-01-28 00:11:00", "2018-05-21 09:50:00", "2018-09-26 11:02:27",
                     "2018-05-17 19:40:00", "2018-10-02 11:02:39",
                     "2018-11-24 13:01:10", "2018-12-19 15:01:21"),
              A1 = c(100, 100, 90, 100, 80, 50, 90, 80, 70, 90, 90, 90, 100,
                     100, 100, 90, 90, 100, 100, 100, 100, 100, 100, 100, 80,
                     90, 80, 100, 100, 80, 90, 80, 90, 80, 80, 80, 100, 90, 100,
                     90, 100, 90, 90, 90, 100, 100, 100, 80, 100, 80),
              A2 = c(100, 100, 90, 100, 80, 40, 100, 80, 70, 90, 90, 90, 100,
                     100, 100, 90, 90, 100, 100, 100, 100, 100, 100, 100, 80,
                     90, 80, 100, 100, 80, 90, 80, 90, 90, 90, 80, 100, 90, 100,
                     90, 100, 90, 90, 90, 100, 100, 100, 80, 100, 80),
              B1 = c(100, 100, 90, 100, 90, 60, 90, 80, 80, 90, 100, 90, 90,
                     100, 100, 90, 90, 90, 100, 100, 100, 100, 100, 100, 80,
                     80, 80, 100, 100, 80, 80, 80, 80, 80, 90, 100, 90, 80, 100,
                     80, 80, 90, 80, 90, 100, 100, 100, 90, 80, 80),
              B2 = c(NA, NA, NA, NA, NA, 60, NA, 80, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 70, 90, NA,
                     NA, 80, 80, 80, 80, 90, NA, NA, NA, 80, NA, 80, 100, NA, 80,
                     NA, NA, NA, NA, NA, 80, 80),
              B3 = c(NA, NA, NA, NA, NA, 60, NA, 80, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 80, 80, NA,
                     NA, 80, 80, 80, 80, 90, NA, NA, NA, 70, NA, 80, 100, NA, 80,
                     NA, NA, NA, NA, NA, 90, 80),
              B4 = c(NA, NA, NA, NA, NA, 50, NA, 80, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 80, 80, NA,
                     NA, 80, 90, 80, 80, 90, NA, NA, NA, 70, NA, 80, 100, NA, 80,
                     NA, NA, NA, NA, NA, 80, 80),
              B5 = c(NA, NA, NA, NA, NA, 50, NA, NA, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 70, 80, NA,
                     NA, 80, 90, 80, 90, 80, NA, NA, NA, 70, NA, 90, 100, NA, 80,
                     NA, NA, NA, NA, NA, 90, 80),
              B6 = c(NA, NA, NA, NA, NA, NA, NA, 80, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 10, 90, NA,
                     NA, 80, 80, 80, 60, 70, NA, NA, NA, 70, NA, 100, 100, NA, 80,
                     NA, NA, NA, NA, NA, 80, 80),
            OS21 = c(NA, NA, NA, NA, NA, 50, NA, 80, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 80, 90, NA,
                     NA, 80, 90, 80, 80, 80, NA, NA, NA, 80, NA, 90, 100, NA, 80,
                     NA, NA, NA, NA, NA, 100, 80),
            OS22 = c(NA, NA, NA, NA, NA, 50, NA, 80, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 80, NA,
                     NA, 80, 90, 80, 80, 90, NA, NA, NA, 90, NA, 90, 40, NA, 80,
                     NA, NA, NA, NA, NA, 80, 80),
              C1 = c(100, 100, 100, 100, 80, 50, 90, 90, 80, 100, 90, 90, 90,
                     100, 100, 100, 90, 100, 100, 100, 100, 100, 100, 100, 80,
                     100, 90, 100, 100, 80, 90, 80, 90, 90, 80, 80, 100, 90, 100,
                     90, 100, 90, 90, 90, 100, 100, 100, 90, 90, 90),
              C2 = c(100, 100, 100, 100, 90, 50, 100, 90, 70, 100, 100, 80,
                     100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100,
                     100, 80, 100, 80, 100, 100, 80, 90, 80, 90, 90, 100, 80, 100,
                     100, 100, 100, 100, 90, 90, 60, 100, 100, 100, 100, 100, 80),
              C3 = c(100, 100, 100, 100, 90, 60, 100, 90, 80, 100, 90, 90, 90,
                     100, 100, 100, 90, 100, 100, 100, 100, 100, 100, 100, 80,
                     90, 80, 100, 100, 80, 90, 80, 90, 90, 100, 80, 100, 80, 100,
                     90, 100, 90, 90, 80, 100, 100, 100, 90, 100, 80),
              C4 = c(NA, NA, NA, NA, 90, 60, NA, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, NA, NA,
                     NA, 80, NA, 80, NA, NA, 100, 90, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              C5 = c(NA, NA, NA, NA, 90, 60, NA, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, NA, NA,
                     NA, 80, NA, 80, NA, NA, 100, 80, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              C6 = c(NA, NA, NA, NA, 90, 60, NA, NA, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, NA, NA,
                     NA, 80, NA, 80, NA, NA, 100, 90, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              C7 = c(NA, NA, NA, NA, 90, 60, NA, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, NA, NA,
                     NA, 80, NA, 80, NA, NA, 100, 100, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              C8 = c(NA, NA, NA, NA, 90, 60, NA, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, NA, NA,
                     NA, 80, NA, 80, NA, NA, 100, 90, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
            OS44 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
            OS45 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              D1 = c(100, 100, 90, 100, 80, 0, 90, NA, 100, NA, NA, 100, NA,
                     100, 100, 90, 80, 90, 100, 100, 100, 100, NA, 100, 0, 100,
                     80, 100, 100, 0, 90, 60, 90, 100, 80, 90, 100, 90, 100, NA,
                     70, NA, NA, 0, 100, 100, 100, 80, NA, 80),
              E1 = c(100, 100, 100, 100, 90, 70, 90, 90, 80, 100, 100, 80, 100,
                     100, 100, 90, 90, 100, 100, 100, 100, 100, 100, 100, 80,
                     80, 80, 100, 100, 80, 90, 80, 90, 90, 90, 80, 100, 80, 100,
                     90, 100, 100, 90, 90, 100, 100, 100, 80, 90, 80),
              E2 = c(100, 100, 100, 100, 90, 80, 50, 90, 90, 100, 100, 90, 100,
                     100, 100, 90, 90, 100, 100, 80, 100, 100, 100, 100, 80,
                     100, 80, 100, 100, 80, 90, 80, 90, 100, 100, 90, 100, 90, 100,
                     90, 100, 90, 80, 70, 100, 100, 100, 80, 100, 80),
              E3 = c(NA, NA, NA, NA, NA, 70, NA, NA, 80, NA, NA, 90, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 100, 80, NA,
                     NA, 70, NA, 80, NA, NA, NA, 90, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 80, NA, 80),
              E4 = c(NA, NA, NA, NA, NA, 60, NA, NA, 90, NA, NA, 80, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, 80, 80, NA,
                     NA, 70, NA, 80, NA, NA, NA, 90, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 70, NA, 80),
              E5 = c(100, 100, NA, 100, 80, NA, NA, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 80, NA, 100, 100, 100, NA, 80, NA, NA, NA,
                     100, NA, 90, NA, 90, NA, NA, NA, NA, NA, NA, NA, NA, 90,
                     NA, NA, NA, 100, NA, NA, NA, NA),
              E6 = c(NA, NA, NA, NA, NA, 50, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, NA, NA,
                     NA, NA, NA, 80, NA, NA, NA, NA, NA, 90, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, 80),
              E7 = c(1, 1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1,
                     1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1,
                     1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2),
            OS17 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              F1 = c(100, 100, 100, 100, 90, 60, 90, 90, 80, 100, 100, 90, 100,
                     100, 100, 90, 90, 90, 100, 100, 100, 100, 100, 100, 80,
                     100, 70, 100, 100, 80, 100, 80, 90, 80, 100, 90, 100, 90, 100,
                     90, 100, 90, 90, 90, 100, 100, 100, 80, 100, 80),
              F2 = c(NA, NA, NA, NA, NA, 80, NA, NA, 100, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 80, NA,
                     NA, 80, NA, 80, NA, 100, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 90, NA, 80),
              F3 = c(NA, NA, NA, NA, NA, 80, NA, NA, 80, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 80, NA,
                     NA, 80, NA, 80, NA, 100, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 90, NA, 80),
              F4 = c(NA, NA, NA, NA, NA, 60, NA, NA, 100, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 80, NA,
                     NA, 80, NA, 80, NA, 100, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 80, NA, 80),
              F5 = c(NA, NA, NA, NA, NA, 60, NA, NA, 60, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 80, NA,
                     NA, 80, NA, 80, NA, 90, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 80, NA, 70),
              F6 = c(NA, NA, NA, NA, NA, 70, NA, NA, 60, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 10, NA,
                     NA, 50, NA, 80, NA, 90, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 70, NA, 70),
              F7 = c(NA, NA, NA, NA, NA, 70, NA, NA, 60, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80, NA, 60, NA,
                     NA, 80, NA, 80, NA, 90, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, 60, NA, 60),
            OS46 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
            OS32 = c(1, 1, 2, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1,
                     1, 1, 2, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 1,
                     1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 2, 2),
              G1 = c(100, 100, 0, 100, 0, 60, 90, 0, 90, 100, 100, 60, 0, 100,
                     100, 90, 70, 100, 100, 100, 100, 0, 100, 100, 0, 100, 0,
                     100, 0, 80, 0, 80, 90, 80, 0, 0, 100, 90, 100, 90, 100, 0, 0,
                     90, 100, 100, 100, 0, 0, 0),
              G3 = c(2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
                     2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2,
                     2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2),
              G4 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, 1, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA),
              G5 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, 80, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
                     NA, NA, NA, NA, NA, NA, NA),
              H1 = c(100, 100, 100, 100, 90, 10, 100, 100, 80, 100, 100, 80,
                     100, 100, 100, 90, 90, 100, 100, 100, 100, 100, 100, 100,
                     80, 90, 80, 100, 100, 80, 90, 80, 90, 90, 90, 80, 100, 80,
                     100, 90, 60, 90, 60, 90, 100, 100, 100, 50, 80, 80),
              H3 = c(100, 100, 90, 100, 90, 100, 100, 100, 80, 100, 100, 80,
                     100, 100, 100, 90, 90, 100, 80, 100, 100, 100, 100, 100,
                     80, 90, 80, 100, 100, 80, 90, 100, 100, 90, 100, 90, 100, 80,
                     100, 90, 60, 90, 60, 90, 100, 100, 100, 50, 80, 80),
              H2 = c(100, 100, 90, 100, 90, 100, 100, 100, 90, 100, 100, 80,
                     100, 100, 100, 90, 90, 100, 100, 100, 100, 100, 100, 100,
                     80, 90, 90, 100, 100, 80, 90, 100, 100, 90, 100, 90, 100,
                     80, 100, 90, 80, 90, 60, 90, 100, 100, 100, 50, 90, 80),
              I1 = c(100, 100, 90, 100, 90, 90, 90, 90, 80, 100, 100, 80, 100,
                     90, 100, 80, 90, 100, 80, 90, 100, 100, 100, 100, 80, 90,
                     80, 100, 100, 80, 100, 80, 90, 80, 100, 90, 100, 90, 100,
                     80, 100, 90, 70, 90, 100, 100, 100, 80, 90, 80),
            OS40 = c(2, 1, 2, NA, 1, 1, 1, 1, NA, 1, 1, 2, 1, 1, 1, NA, 2, NA,
                     1, 2, 1, NA, NA, NA, NA, 1, 1, 2, 1, 1, 2, NA, 2, NA, 1,
                     1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 1, 1, NA, NA, NA),
       ModelLong = c("ccc", "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )", "ccc (2016 ~ )",
                     "ccc (2016 ~ )", "ccc (2016 ~ )", "aaa (2014 ~ )",
                     "aaa (2014 ~ )", "aaa (2014 ~ )", "aaa (2014 ~ )", "aaa (2014 ~ )",
                     "aaa (2014 ~ )", "aaa (2014 ~ )", "aaa (2014 ~ )",
                     "aaa (2014 ~ )"),
            year = c(2018, 2018, 2018, 2018, 2018, 2017, 2018, 2018, 2018,
                     2018, 2018, 2018, 2018, 2018, 2018, 2018, 2018, 2018,
                     2018, 2018, 2018, 2018, 2018, 2017, 2018, 2018, 2017, 2018,
                     2018, 2018, 2018, 2018, 2018, 2017, 2017, 2018, 2018, 2018,
                     2018, 2018, 2017, 2018, 2017, 2017, 2018, 2018, 2018, 2018,
                     2018, 2018),
           month = c("06", "05", "10", "12", "06", "03", "06", "12", "10",
                     "12", "12", "05", "12", "05", "06", "05", "05", "06",
                     "06", "12", "09", "06", "10", "01", "06", "05", "01", "05",
                     "06", "12", "06", "10", "05", "01", "01", "10", "05", "05",
                     "05", "11", "02", "09", "01", "01", "05", "09", "05", "10",
                     "11", "12"),
          ymonth = c(6, 5, 10, 12, 6, 3, 6, 12, 10, 12, 12, 5, 12, 5, 6, 5, 5,
                     6, 6, 12, 9, 6, 10, 1, 6, 5, 1, 5, 6, 12, 6, 10, 5, 1, 1,
                     10, 5, 5, 5, 11, 2, 9, 1, 1, 5, 9, 5, 10, 11, 12),
              yq = c("2018 Q2", "2018 Q2", "2018 Q4", "2018 Q4", "2018 Q2",
                     "2017 Q1", "2018 Q2", "2018 Q4", "2018 Q4", "2018 Q4",
                     "2018 Q4", "2018 Q2", "2018 Q4", "2018 Q2", "2018 Q2",
                     "2018 Q2", "2018 Q2", "2018 Q2", "2018 Q2", "2018 Q4", "2018 Q3",
                     "2018 Q2", "2018 Q4", "2017 Q1", "2018 Q2", "2018 Q2",
                     "2017 Q1", "2018 Q2", "2018 Q2", "2018 Q4", "2018 Q2", "2018 Q4",
                     "2018 Q2", "2017 Q1", "2017 Q1", "2018 Q4", "2018 Q2",
                     "2018 Q2", "2018 Q2", "2018 Q4", "2017 Q1", "2018 Q3", "2017 Q1",
                     "2017 Q1", "2018 Q2", "2018 Q3", "2018 Q2", "2018 Q4",
                     "2018 Q4", "2018 Q4"),
              hy = c(1, 1, 2, 2, 1, 1, 1, 2, 2, 2, 2, 1, 2, 1, 1, 1, 1, 1, 1,
                     2, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1,
                     1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 2),
        halfyear = c(20181, 20181, 20182, 20182, 20181, 20171, 20181, 20182,
                     20182, 20182, 20182, 20181, 20182, 20181, 20181, 20181,
                     20181, 20181, 20181, 20182, 20182, 20181, 20182, 20171, 20181,
                     20181, 20171, 20181, 20181, 20182, 20181, 20182, 20181,
                     20171, 20171, 20182, 20181, 20181, 20181, 20182, 20171, 20182,
                     20171, 20171, 20181, 20182, 20181, 20182, 20182, 20182),
       monthyear = c("2018_06", "2018_05", "2018_10", "2018_12", "2018_06",
                     "2017_03", "2018_06", "2018_12", "2018_10", "2018_12",
                     "2018_12", "2018_05", "2018_12", "2018_05", "2018_06",
                     "2018_05", "2018_05", "2018_06", "2018_06", "2018_12", "2018_09",
                     "2018_06", "2018_10", "2017_01", "2018_06", "2018_05",
                     "2017_01", "2018_05", "2018_06", "2018_12", "2018_06", "2018_10",
                     "2018_05", "2017_01", "2017_01", "2018_10", "2018_05",
                     "2018_05", "2018_05", "2018_11", "2017_02", "2018_09", "2017_01",
                     "2017_01", "2018_05", "2018_09", "2018_05", "2018_10",
                     "2018_11", "2018_12"),
            A2TB = c(100, 100, 100, 100, 0, 0, 100, 0, 0, 100, 100, 100, 100,
                     100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 0,
                     100, 0, 100, 100, 0, 100, 0, 100, 100, 100, 0, 100, 100,
                     100, 100, 100, 100, 100, 100, 100, 100, 100, 0, 100, 0)
)

I managed to prepare my mean scores summary using ddply which can be exported to Excel (thank you for your help again).

library(plyr)
halfyear.scores <- ddply(df, .(year,hy), summarize,  A2TB=mean(A2TB), B1=mean(B1), C1=mean(C1), D1=mean(D1,na.rm=TRUE), E1=mean(E1), F1=mean(F1), G1=mean(G1))
halfyear.scores

The result shows mean scores of all variables but I cannot add counts for each period (hy).
I've tried adding one more variable to the code (the only solution I found myself):

halfyear.scores <- ddply(df, .(year,hy), summarize,  A2TB=mean(A2TB), B1=mean(B1), C1=mean(C1), D1=mean(D1,na.rm=TRUE), E1=mean(E1), F1=mean(F1), G1=mean(G1), A2=n())

but sometimes it works but sometimes it does not.
Is any better solution for creating this table? There are some missing responses in D1 so ideally I would like to add A2TB and D1 counts in the end of the table.

Is it easy to do?

andresrcs · January 30, 2019, 6:28pm

Is not clear what you are trying to accomplish, is it similar to this?

library(dplyr)

df %>% 
    group_by(year, hy) %>% 
    summarise(A2TB = mean(A2TB),
              B1 = mean(B1),
              C1 = mean(C1),
              D1 = mean(D1, na.rm = TRUE),
              E1 = mean(E1),
              F1 = mean(F1),
              G1 = mean(G1)) %>%
    mutate(A2 = n())
#> # A tibble: 3 x 10
#> # Groups:   year [2]
#>    year    hy  A2TB    B1    C1    D1    E1    F1    G1    A2
#>   <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <int>
#> 1  2017     1  75    82.5  86.2  61.4  88.8  86.2  53.8     1
#> 2  2018     1  91.3  92.6  95.7  91.3  93.9  95.7  73.5     2
#> 3  2018     2  63.2  90    91.6  81.8  92.1  92.6  54.7     2

Slavek · January 31, 2019, 10:01am

Yes bur the result in the last column should be a count of responses which is:

table(df$hy,df$year)
   
    2017 2018
  1    8   23
  2    0   19

So in the last column I should get 8, 23 and 19.

Can I replace A2 in mutate by any variable used in summarise? Would it be a table conflict (two columns with the same variable)? Sorry for bothering you but I cannot check it myself because when I run your syntax I can see following in the console:

> library(dplyr)
> 
> df %>% 
+   group_by(year, hy) %>% 
+   summarise(A2TB = mean(A2TB),
+             B1 = mean(B1),
+             C1 = mean(C1),
+             D1 = mean(D1, na.rm = TRUE),
+             E1 = mean(E1),
+             F1 = mean(F1),
+             G1 = mean(G1)) %>%
+   mutate(A2 = n())
Error: This function should not be called directly
Call `rlang::last_error()` to see a backtrace
>

What could it be?

andresrcs · January 31, 2019, 2:41pm

library(dplyr)

df %>% 
    group_by(year, hy) %>% 
    summarise(A2TB = mean(A2TB),
              B1 = mean(B1),
              C1 = mean(C1),
              D1 = mean(D1, na.rm = TRUE),
              E1 = mean(E1),
              F1 = mean(F1),
              G1 = mean(G1),
              count = n())
#> # A tibble: 3 x 10
#> # Groups:   year [2]
#>    year    hy  A2TB    B1    C1    D1    E1    F1    G1    count
#>   <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <int>
#> 1  2017     1  75    82.5  86.2  61.4  88.8  86.2  53.8     8
#> 2  2018     1  91.3  92.6  95.7  91.3  93.9  95.7  73.5     23
#> 3  2018     2  63.2  90    91.6  81.8  92.1  92.6  54.7     19

andresrcs · January 31, 2019, 2:42pm

It seems that there is something really odd with your Rstudio setup

Slavek · January 31, 2019, 2:53pm

Brilliant!!!
Thank you very much

system · February 7, 2019, 2:54pm

This topic was automatically closed 7 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.