Tidy summarizes information about the components of a model. A model component might be a single term in a regression, a single hypothesis, a cluster, or a class. Exactly what tidy considers to be a model component varies cross models but is usually self-evident. If a model has several distinct types of components, you will need to specify which components to return.
# S3 method for manova
tidy(x, test = "Pillai", ...)
A manova
object return from stats::manova()
.
One of "Pillai" (Pillai's trace), "Wilks" (Wilk's lambda), "Hotelling-Lawley" (Hotelling-Lawley trace) or "Roy" (Roy's greatest root) indicating which test statistic should be used. Defaults to "Pillai".
Arguments passed on to stats::summary.manova
An object of class "manova"
or an aov
object with multiple responses.
The name of the test statistic to be used. Partial matching is used so the name can be abbreviated.
logical. If TRUE
, the intercept term is
included in the table.
tolerance to be used in deciding if the residuals are
rank-deficient: see qr
.
A tibble::tibble with columns:
\item{term}{Term in design} \item{statistic}{Approximate F statistic} \item{num.df}{Degrees of freedom} \item{p.value}{P-value}
Depending on which test statistic is specified, one of the following columns is also included:
\item{pillai}{Pillai's trace} \item{wilks}{Wilk's lambda} \item{hl}{Hotelling-Lawley trace} \item{roy}{Roy's greatest root}
tidy()
, stats::summary.manova()
Other anova tidiers: tidy.TukeyHSD
,
tidy.anova
, tidy.aovlist
,
tidy.aov
# NOT RUN {
npk2 <- within(npk, foo <- rnorm(24))
m <- manova(cbind(yield, foo) ~ block + N * P * K, npk2)
tidy(m)
# }
Run the code above in your browser using DataLab