Запланированные контрасты с использованием Wilcox, надежного ANOVA в R (пакет WRS2)

Question

Я запускаю двухсторонний ANOVA и пост-тесты на моем наборе данных с использованием WRS2. Тем не менее, в рамках контрастов я не нахожу сравнения, которые мне нужны. Можно ли планировать мои собственные контрасты?

Answer 1

К сожалению, похоже, что контрасты жестко закодированы в функции mcp2atm. Некоторые отрывки из кишки этой функции (которую вы можете увидеть, набрав mcp2atm в командной строке: вы также можете использовать dump() для внешнего файла).

... 
temp <- con2way(J, K) 
conA <- temp$conA 
conB <- temp$conB 
conAB <- temp$conAB 
if (!op) { ## note: op is hardcoded to FALSE earlier in the function 
    Factor.A <- lincon1(x, con = conA, tr = tr, alpha = alpha) 
    Factor.B <- lincon1(x, con = conB, tr = tr, alpha = alpha) 
    Factor.AB <- lincon1(x, con = conAB, tr = tr, alpha = alpha) 
} 
## ... stuff for generating contrast names ... 
contrasts <- as.data.frame(cbind(conA, conB, conAB)) 
colnames(contrasts) <- c(dnamesA, dnamesB, dnamesAB) 
rownames(contrasts) <- colnames(dataWide) 
... 

con2way и lincon1 являются внутренними WRS2 функции (вы можете увидеть их с помощью WRS2::con2way и WRS2:::lincon1). (Если вы действительно хотите вникнуть в исходный код, вероятно, имеет смысл загрузить исходный пакет и распаковать его.)

Чтобы выполнить пользовательские контрасты, вам, вероятно, придется взять его с собой сопровождающий WRS2 (maintainer("WRS2"), или возможно, попробуйте issues tracker on R-forge ...) Я не знаю, предоставит ли какой-либо другой механизм надежной статистики (см. robust stats task view) альтернативный маршрут ...