The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
n and
trim (when n < 2 * floor(k) + 3, with
k = trim * n, and for n = 2 or 3
with k < 1) the sum of the interior observations used a
range a:b with a > b, which counts
backwards in R instead of being empty. Examples fixed:
sgmean(1:9, 0.45) returned 27.22 (now 5),
sgmean(1:10, 0.45) returned 16.5 (now 5.5),
sgmean(1:2, 0.3) returned 5.25 (now 1.5). Results for usual
settings (small trim, moderate n) are
unchanged.[k, n - k], which removes the case analysis
and the : ranges.na.rm = FALSE with missing values now returns
NA, as mean() does (before,
sort() silently dropped them). Invalid trim
(length other than one, NA) now gives a clear error.testthat suite, including an independent
reference implementation and regression tests for the cases above.Critical fix: Corrected the proportional trimmed
mean formula for cases where k = trim * n >= 1. The
previous implementation (v0.1.0) incorrectly applied a
(1 - k) discount to x[1] and x[n]
for all cases, which is only correct when k < 1 (Type A
distortion). For k >= 1, this produced incorrect
results.
For exact integer k >= 2, sgmean()
now correctly returns the same result as mean(..., trim)
(as documented and expected).
For non-integer k > 1 (Type B distortion),
sgmean() now correctly applies a fractional discount
(1 - delta) to x[floor(k) + 1] and
x[n - floor(k)], consistent with expression (4) of the
companion article and verified against Statgraphics output.
The fix was discovered and verified during peer review of the companion article submitted to The R Journal. The corrected implementation produces results identical to Statgraphics for all tested cases.
mean(..., trim)mean(..., trim).These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.