/usr/lib/R/site-library/matrixStats/benchmarking/meanOver.md.rsp is in r-cran-matrixstats 0.14.2-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | <%@include file="includes/setup.md.rsp"%>
<%@string fcnname="meanOver"%>
<% fcnname <- "<%@string name="fcnname"%>" %>
<%@meta title="${fcnname}() benchmarks"%>
<%@meta author="Henrik Bengtsson"%>
<%@meta date="2014-11-02"%>
<%@include file="${header}"%>
# <%@meta name="title"%>
This report benchmark the performance of <%=fcnname%>() against alternative methods.
## Alternative methods
* mean() + [()
* mean.default() + [() - avoids method dispatching
as below
```r
<%=withCapture({
meanOver_R_v1 <- function(x, na.rm=FALSE, idxs) {
mean(x[idxs], na.rm=na.rm)
}
})%>
```
and
```r
<%=withCapture({
meanOver_R_v2 <- function(x, na.rm=FALSE, idxs) {
mean.default(x[idxs], na.rm=na.rm)
}
})%>
```
<% for (mode in c("integer", "double")) { %>
## Data type "<%=mode%>"
### Data
```r
<%=withCapture({
<%@include file="R/random-vectors.R"%>
data <- rvectors(mode=mode)
##data <- data[1:3]
})%>
```
### Results
<% for (ii in seq_along(data)) { %>
<%
dataLabel <- names(data)[ii]
mprintf("%s: %s\n", mode, dataLabel)
x <- data[[dataLabel]]
gc()
%>
### <%=dataLabel%> vector
#### All elements
```r
<%=withCapture({
x <- data[[.dataLabel.]]
gc()
stats <- microbenchmark(
"meanOver" = meanOver(x, refine=TRUE),
"meanOver_no_refine" = meanOver(x, refine=FALSE),
"mean" = mean(x),
"mean.default" = mean.default(x),
unit = "ms"
)
})%>
```
<% benchmarkResults(stats, tags=c(dataLabel, "all")) %>
<% for (subset in c(0.2, 0.4, 0.8)) { %>
#### A <%=sprintf("%g", 100*subset)%>% subset
```r
<%=withCapture({
x <- data[[.dataLabel.]]
subset
idxs <- sort(sample(length(x), size=subset*length(x), replace=FALSE))
gc()
stats <- microbenchmark(
"meanOver" = meanOver(x, idxs=idxs, refine=TRUE),
"meanOver_no_refine" = meanOver(x, idxs=idxs, refine=FALSE),
"mean+[()" = meanOver_R_v1(x, idxs=idxs),
"mean.default+[()" = meanOver_R_v2(x, idxs=idxs),
unit = "ms"
)
})%>
```
<% benchmarkResults(stats, tags=c(mode, dataLabel, subset)) %>
<% } # for (subset in ...) %>
<% } # for (ii ...) %>
<% } # for (mode ...) %>
<%@include file="${footer}"%>
<%---------------------------------------------------------------------------
HISTORY:
2014-11-02
o Created.
---------------------------------------------------------------------------%>
|