Functions for calculating M- and MM-estimators for location given values and associated standard errors or standard uncertainties.

1 2 3 4 5 6 7 8 9 | ```
MM.estimate(x, ...)
## Default S3 method:
MM.estimate(x, u, c = 4.685, ...)
huber.estimate(x, ...)
## Default S3 method:
huber.estimate(x, u, k= 1.345, ...)
``` |

`x` |
numeric vector of mean values for groups |

`u` |
numeric vector of standard deviations or standard uncertainties associated with the values |

`c, k` |
Tuning parameters passed to other functions (see |

`...` |
Parameters passed to other functions. |

These functions are wrappers for robust estimation using `rlm`

. All simply
call `rlm`

with the formula `x~1`

and weights `1/u^2`

.

An object of class ‘loc.est’.

S. L. R. Ellison s.ellison@lgc.co.uk

None, yet.

1 2 3 4 5 6 7 8 | ```
## Cd heat of vapourisation example
## from Paule, R. C. and Mandel, J. (1982) - see ?mpaule
x2<-c(27.044, 26.022, 26.340, 26.787, 26.796)
v<-c(3, 76, 464, 3, 14)*1e-3
MM.estimate(x2, sqrt(v))
huber.estimate(x2, sqrt(v))
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.