RWPNPBT {bootLRTpairwise} | R Documentation |
Rademacher Wild Pairwise Non-Parametric Bootstrap Test (RWPNPBT)
Description
Performs a non-parametric bootstrap test using Rademacher wild bootstrap re-sampling. This test compares all pairwise group means using a standardized distance metric, making it robust to violations of normality and heteroscedasticity.
Usage
RWPNPBT(group_list, B = 1e+05, alpha = 0.05)
Arguments
group_list |
A list where each element is a numeric vector of raw observations for a group. |
B |
Number of bootstrap re-samples. Default is 100000. |
alpha |
Significance level for the hypothesis test. Default is 0.05. |
Details
The test statistic sums the absolute standardized differences between all pairs of group means. Rademacher weights are applied to centered observations for wild bootstrapping.
Value
A list of class "RWPNPBT"
containing:
test_statistic |
Observed value of the test statistic. |
critical_value |
Bootstrap-based critical value at the given alpha level. |
decision |
Conclusion of the hypothesis test based on the critical value. |
Examples
set.seed(123)
group1 <- rexp(18, rate = 1/10) # Exponential distribution
group2 <- runif(22, min = 10, max = 18) # Uniform distribution
group3 <- rchisq(20, df = 5) + 7 # Right-skewed Chi-square + shift
RWPNPBT(list(group1, group2, group3), B = 1000)