Expression containing results from heteroscedastic one-way ANOVA for trimmed means

expr_anova_robust(
  data,
  x,
  y,
  subject.id = NULL,
  paired = FALSE,
  k = 2L,
  conf.level = 0.95,
  tr = 0.1,
  nboot = 100L,
  output = "expression",
  ...
)

Arguments

data

A dataframe (or a tibble) from which variables specified are to be taken. A matrix or tables will not be accepted.

x

The grouping variable from the dataframe data.

y

The response (a.k.a. outcome or dependent) variable from the dataframe data.

subject.id

In case of repeated measures design (paired = TRUE, i.e.), this argument specifies the subject or repeated measures id. Note that if this argument is NULL (which is the default), the function assumes that the data has already been sorted by such an id by the user and creates an internal identifier. So if your data is not sorted and you leave this argument unspecified, the results can be inaccurate.

paired

Logical that decides whether the experimental design is repeated measures/within-subjects or between-subjects. The default is FALSE.

k

Number of digits after decimal point (should be an integer) (Default: k = 2L).

conf.level

Scalar between 0 and 1. If unspecified, the defaults return 95% lower and upper confidence intervals (0.95).

tr

Trim level for the mean when carrying out robust tests. If you get error stating "Standard error cannot be computed because of Winsorized variance of 0 (e.g., due to ties). Try to decrease the trimming level.", try to play around with the value of tr, which is by default set to 0.1. Lowering the value might help.

nboot

Number of bootstrap samples for computing confidence interval for the effect size (Default: 100).

output

If "expression", will return expression with statistical details, while "dataframe" will return a dataframe containing the results.

...

Additional arguments (currently ignored).

Value

For more details, see- https://indrajeetpatil.github.io/statsExpressions/articles/stats_details.html

Examples

# for reproducibility set.seed(123) library(statsExpressions) # ------------------------ between-subjects ----------------------------- expr_anova_robust( data = ggplot2::midwest, x = state, y = percbelowpoverty )
#> paste(italic("F")["trimmed-means"], "(", "4", ",", "169.73", #> ") = ", "11.43", ", ", italic("p"), " = ", "3.09e-08", ", ", #> widehat(italic(xi)), " = ", "0.35", ", CI"["95%"], " [", #> "0.27", ", ", "0.43", "]", ", ", italic("n")["obs"], " = ", #> 437L)
# ------------------------ within-subjects ----------------------------- expr_anova_robust( data = iris_long, x = condition, y = value, paired = TRUE, k = 3 )
#> paste(italic("F")["trimmed-means"], "(", "1.118", ",", "132.992", #> ") = ", "559.197", ", ", italic("p"), " = ", "0e+00", ", ", #> italic("n")["pairs"], " = ", 150L)