Ich führe eine zweiseitige ANOVA- und eine posthoc-Prüfung an meinem Datensatz mit WRS2 durch. Aber innerhalb der Kontraste finde ich nicht die Vergleiche, die ich brauche. Ist es möglich, meine eigenen Kontraste zu planen?Geplante Kontraste unter Verwendung der robusten ANOVA von Wilcox in R (WRS2-Paket)
Leider sind die Kontraste in der mcp2atm Funktion fest codiert. Einige Auszüge aus dem Inneren dieser Funktion (die Sie sehen können, indem Sie mcp2atm an der Eingabeaufforderung eingeben: Sie können auch dump() die Funktion zu einer externen Datei).
...
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 und lincon1 sind interne WRS2 Funktionen (Sie sie mit WRS2::con2way und WRS2:::lincon1 sehen). (Wenn Sie wirklich in den Quellcode graben wollen macht es wahrscheinlich sinnvoll, das Quellcode-Paket herunterladen und entpacken.)
Um benutzerdefinierte zu tun Kontraste würden Sie wahrscheinlich mit dem Betreuer des WRS2 nehmen müssen bis (maintainer("WRS2") oder möglicherweise versuchen Sie die issues tracker on R-forge ...) Ich weiß nicht, ob irgendeine der anderen robusten Statistikmaschinerie (sehen Sie die robust stats task view) eine alternative Route zur Verfügung stellen würde ...