Datensatz einlesen und Variablen spezifizieren
# Datensatz 'data': pro Versuchsperson eine Zeile
# Häufigkeitstabellen zuerst umwandeln, vgl. Umwandlung von Datenstrukturen
Vollpräparat <- rep(c("gering", "stark"), c(14, 24))
Leerpräparat <- rep(c("gering", "stark", "gering", "stark"), c(10, 4, 15, 9))
data <- data.frame(Vollpräparat, Leerpräparat)
# Variablen spezifizieren
Faktor1 <- "Vollpräparat" # Name des ersten Faktors eingeben
Faktor2 <- "Leerpräparat" # Name des zweiten Faktors eingeben
Häufigkeitstabelle und Signifikanztest
# Häufigkeitstabelle
data2 <- na.omit(data[, c(Faktor1, Faktor2), drop=FALSE])
tab <- xtabs(paste("~", Faktor1, "+", Faktor2), data2)
# Signifikantest
test <- mcnemar.test(tab)
list(Häufigkeitstabelle=tab, Signifikanztest=test)
## $Häufigkeitstabelle
## Leerpräparat
## Vollpräparat gering stark
## gering 10 4
## stark 15 9
##
## $Signifikanztest
##
## McNemar's Chi-squared test with continuity correction
##
## data: tab
## McNemar's chi-squared = 5.2632, df = 1, p-value = 0.02178