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