Ich möchte einen Cochran-Mantel-Haenszel Chi-Squared Test durchführen und verzweifle gerade daran, meine Daten in das dazu nötige Table-Format zu bringen. Ich hab mir den Datensatz UCBAdmissions angeschaut, wie der aussieht und versucht, das bei mir nachzubauen:
Code: Alles auswählen
mydata <- structure(list(H3K27ac_Common_I_clade_1 = structure(c(90L, 177L,
1125L, 586L), .Dim = c(2L, 2L), .Dimnames = list(c("observed_CAGE_enhancers",
"observed_other_enhancers"), c("in clade", "in cluster"))), H3K27ac_Common_I_clade_2 = structure(c(83L,
290L, 1132L, 473L), .Dim = c(2L, 2L), .Dimnames = list(c("observed_CAGE_enhancers",
"observed_other_enhancers"), c("in clade", "in cluster"))), H3K27ac_Common_I_clade_3 = structure(c(109L,
20L, 1106L, 743L), .Dim = c(2L, 2L), .Dimnames = list(c("observed_CAGE_enhancers",
"observed_other_enhancers"), c("in clade", "in cluster"))), H3K27ac_Common_I_clade_4 = structure(c(292L,
52L, 923L, 711L), .Dim = c(2L, 2L), .Dimnames = list(c("observed_CAGE_enhancers",
"observed_other_enhancers"), c("in clade", "in cluster"))), H3K27ac_Common_I_clade_5 = structure(c(151L,
12L, 1064L, 751L), .Dim = c(2L, 2L), .Dimnames = list(c("observed_CAGE_enhancers",
"observed_other_enhancers"), c("in clade", "in cluster"))), H3K27ac_Common_I_clade_6 = structure(c(132L,
115L, 1083L, 648L), .Dim = c(2L, 2L), .Dimnames = list(c("observed_CAGE_enhancers",
"observed_other_enhancers"), c("in clade", "in cluster")))), .Names = c("H3K27ac_Common_I_clade_1",
"H3K27ac_Common_I_clade_2", "H3K27ac_Common_I_clade_3", "H3K27ac_Common_I_clade_4",
"H3K27ac_Common_I_clade_5", "H3K27ac_Common_I_clade_6"))
Code: Alles auswählen
as.table(mydata)
Fehler in as.table.default(mydata) : kann nicht in eine Tabelle umwandeln
Kann man das doch irgendwie einfach konvertieren oder soll ich lieber die Liste gar nicht erst so bauen?
Viele Grüße und Danke
Matthias