# \donttest{
# Sample dataset in long format
set.seed(123)
group1 <- rnorm(30, mean=50, sd=2)
group2 <- rnorm(30, mean=51, sd=3)
df <- data.frame(value=c(group1, group2), group=rep(c("A", "B"), each=30))
# Compare the groups using various descriptive statistics
result <- groupcompare(df, cl=0.95, alternative="two.sided",
q=c(0.25, 0.5, 0.75), qt=0, R=500, out=FALSE, verbose=TRUE)
result
# Compare the groups using Huber's M-estimator of location with bootstrap
bshubermean <- bootstrap(df, statistic=calchubermeandif,
alternative="two.sided", alpha=0.05, R=500)
bshubermean
permhubermean <- permtest(df, statistic=calchubermeandif,
alternative="two.sided", R=500)
permhubermean$pval
# }
Run the code above in your browser using DataLab