### J Appl. Stat 2026 paper: A multi-model procedure for testing versatile alternatives 
### in multi-sample medical studies using a maximum-type test
### Eco-tox example Table 10
library(Rfit)
data(quail)
mod1<-lm(ldl~treat, data=quail)
library(coin)
quail$AB<- ansari_trafo(quail$ldl, ties.method ="mid-ranks")
quail$SA<-savage_trafo(quail$ldl, ties.method ="mid-ranks")
quail$Rldl<- rank(quail$ldl) 
mod2<-lm(Rldl~treat, data=quail)
mod3<-lm(AB~treat, data=quail)
mod4<-lm(SA~treat, data=quail)
library(multcomp)
Joint2 <- glht(mmm(location = mod2, scale= mod3, shape=mod4), mlf(mcp(treat ="GrandMean")),df=Glob$df[2])
summary(Joint2)