
medpolish(x, eps = 0.01, maxiter = 10, trace.iter = TRUE, na.rm = FALSE)
medpolish
with the following named components:
eps
or until there have been maxiter
iterations.
The sum of absolute residuals is printed at
each iteration of the fitting process, if trace.iter
is TRUE
.
If na.rm
is FALSE
the presence of any NA
value in
x
will cause an error, otherwise NA
values are ignored. medpolish
returns an object of class medpolish
(see below).
There are printing and plotting methods for this
class, which are invoked via by the generics
print
and plot
.
median
; aov
for a mean
instead of median decomposition.require(graphics)
## Deaths from sport parachuting; from ABC of EDA, p.224:
deaths <-
rbind(c(14,15,14),
c( 7, 4, 7),
c( 8, 2,10),
c(15, 9,10),
c( 0, 2, 0))
dimnames(deaths) <- list(c("1-24", "25-74", "75-199", "200++", "NA"),
paste(1973:1975))
deaths
(med.d <- medpolish(deaths))
plot(med.d)
## Check decomposition:
all(deaths ==
med.d$overall + outer(med.d$row,med.d$col, "+") + med.d$residuals)
Run the code above in your browser using DataLab