Please don't post in HTML. It makes your code unreadable. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium
To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey 2015-11-02 13:13 GMT+01:00 Axel Urbiz <axel.ur...@gmail.com>: > Sorry, this is *related* to a recent post, but not the same. I'd appreciate > your help in getting the same results with the two methods below (the first > using plyr and the second using dplyr. The former works, but not the > latter.) > > > ### Sample data > set.seed(4)df <- data.frame(pred = rnorm(100), y = sample(c(0,1), 100, > replace = TRUE), models = gl(2, 50, 100, labels = > c("model1", "model2"))) > ### using plyr > fooFun <- function(x) { xcuts <- unique(x$pred) x$bin <- cut(x$pred, > breaks = xcuts, include.lowest = TRUE) x <- > dplyr::summarize(dplyr::group_by(x, bin), sumY = sum(y)) > x}head(plyr::ddply(df, plyr::.(models), fooFun)) > > ### Using dplyr > fooFun2 <- function(pred, y) { xcuts <- unique(pred) bin <- cut(pred, > breaks = xcuts, include.lowest = TRUE) dft <- data.frame(bin, pred, y) > dft <- dplyr::summarize(dplyr::group_by(dft, bin), sumY = sum(y)) dft} > res_dplyr <- dplyr::mutate(dplyr::group_by(df, models), fooFun2(pred, > y))Error: incompatible size (2), expecting 50 (the group size) or 1 > head(res_dplyr) > Thank you. > Axel. > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.