This function is merely a thin wrapper around stats::TukeyHSD with tidying done by broom::tidy and optional formatting via pixiedust::sprinkle. Its input is not a aov model like in the original function, but instead the aov model is fit internally based on the arguments given. This is meant to enable a consistent usage between the tadaa_pairwise-functions.

tadaa_pairwise_tukey(data, response, group1, group2 = NULL, print = "df", ...)

Arguments

data

A data.frame containing the variables.

response

The response variable, i.e. the dependent numeric vector.

group1

The grouping variables, typically a factor.

group2

(Optional) second grouping variable.

print

Print method, defaults to df for data.frame output, otherwise passed to pixiedust::sprinkle_print_method.

...

Further arguments passed to stats::TukeyHSD

Value

A data.frame or pixiedust::dust object depending on print.

See also

Examples

tadaa_pairwise_tukey(data = ngo, deutsch, jahrgang, geschl)
#> # A tibble: 19 x 6 #> term estimate conf.low conf.high adj.p.value contrast #> <chr> <dbl> <dbl> <dbl> <dbl> <chr> #> 1 jahrgang 0.0500 -0.590 0.690 0.981 12-11 #> 2 jahrgang 1.17 0.387 1.95 0.00148 13-11 #> 3 jahrgang 1.12 0.337 1.90 0.00250 13-12 #> 4 geschl 1.03 0.554 1.51 0.0000300 Weiblich-Männlich #> 5 jahrgang:geschl 0.740 -0.362 1.84 0.387 12:Männlich-11:Männl… #> 6 jahrgang:geschl 1.16 -0.190 2.51 0.138 13:Männlich-11:Männl… #> 7 jahrgang:geschl 1.58 0.478 2.68 0.000736 11:Weiblich-11:Männl… #> 8 jahrgang:geschl 0.940 -0.162 2.04 0.143 12:Weiblich-11:Männl… #> 9 jahrgang:geschl 2.76 1.41 4.11 0.000000207 13:Weiblich-11:Männl… #> 10 jahrgang:geschl 0.420 -0.930 1.77 0.948 13:Männlich-12:Männl… #> 11 jahrgang:geschl 0.840 -0.262 1.94 0.246 11:Weiblich-12:Männl… #> 12 jahrgang:geschl 0.200 -0.902 1.30 0.995 12:Weiblich-12:Männl… #> 13 jahrgang:geschl 2.02 0.670 3.37 0.000354 13:Weiblich-12:Männl… #> 14 jahrgang:geschl 0.420 -0.930 1.77 0.948 11:Weiblich-13:Männl… #> 15 jahrgang:geschl -0.220 -1.57 1.13 0.997 12:Weiblich-13:Männl… #> 16 jahrgang:geschl 1.60 0.0415 3.16 0.0404 13:Weiblich-13:Männl… #> 17 jahrgang:geschl -0.640 -1.74 0.462 0.554 12:Weiblich-11:Weibl… #> 18 jahrgang:geschl 1.18 -0.170 2.53 0.125 13:Weiblich-11:Weibl… #> 19 jahrgang:geschl 1.82 0.470 3.17 0.00190 13:Weiblich-12:Weibl…
tadaa_pairwise_tukey(data = ngo, deutsch, jahrgang, print = "console")
#> Term Diff CI_low CI_high p (adj.) Comparison #> 1 jahrgang 0.05 -0.619 0.719 .983 12-11 #> 2 jahrgang 1.17 0.351 1.989 < .01 13-11 #> 3 jahrgang 1.12 0.301 1.939 < .01 13-12 #> #>