1<%@include file="includes/setup.md.rsp"%> 2 3<%@string fcnname="weightedMean"%> 4<% fcnname <- "<%@string name="fcnname"%>" %> 5<%@meta title="${fcnname}() benchmarks"%> 6<%@meta author="Henrik Bengtsson"%> 7<%@meta date="2014-12-09"%> 8 9<%@include file="${header}"%> 10<% 11weighted.mean <- stats:::weighted.mean 12weighted.mean.default <- stats:::weighted.mean.default 13%> 14 15# <%@meta name="title"%> 16 17This report benchmark the performance of <%=fcnname%>() against alternative methods. 18 19## Alternative methods 20 21* stats::weighted.mean() 22* stats:::weighted.mean.default() 23 24 25<% for (mode in c("integer", "double")) { %> 26 27## Data type "<%=mode%>" 28### Data 29```r 30<%=withCapture({ 31<%@include file="R/random-vectors.R"%> 32data <- rvectors(mode = mode) 33data <- data[1:4] 34})%> 35``` 36 37## Results 38 39<% for (ii in seq_along(data)) { %> 40<% 41dataLabel <- names(data)[ii] 42mprintf("%s: %s\n", mode, dataLabel) 43x <- data[[dataLabel]] 44gc() 45%> 46### <%=dataLabel%> vector 47 48```r 49<%=withCapture({ 50x <- data[[.dataLabel.]] 51w <- runif(length(x)) 52gc() 53 54stats <- microbenchmark( 55 "weightedMean" = weightedMean(x, w = w, na.rm = FALSE), 56 "stats::weighted.mean" = weighted.mean(x, w = w, na.rm = FALSE), 57 "stats:::weighted.mean.default" = weighted.mean.default(x, w = w, na.rm = FALSE), 58 unit = "ms" 59) 60})%> 61``` 62 63<% benchmarkResults(stats, tags=c(mode, dataLabel)) %> 64 65<% } # for (ii ...) %> 66 67<% } # for (mode ...) %> 68 69<%@include file="${footer}"%> 70 71 72<%--------------------------------------------------------------------------- 73HISTORY: 742014-12-09 75o Created. 76---------------------------------------------------------------------------%> 77