Code: Alles auswählen
A<-c("A1","A1","A5","A1","A1","A1","A1","A5","A1","A5","A5","A5","A5","A6","A5","A5","A5","A1","A5","A1","A6","A4","A5","A4","A5","A5","A4","A4","A4","A2","A5","A5","A2","A1","A1","A1","A1","A5","A6","A6","A6","A5","A5","A5","A5","A5","A6","A6","A6","A1","A4","A4","A5","A5","A4","A4","A4","A4","A4","A4","A5","A5","A5","A5","A5","A5","A1","A5","A5","A1","A5","A1","A5","A1","A1","A5","A1","A5","A5","A5","A5","A5","A5","A1","A1","A5","A5","A5","A5","A6","A6","A6","A6","A6","A6","A6","A4","A4","A4","A4","A5","A5","A4","A5","A4","A4","A4","A4","A4","A4","A5","A4","A4","A4","A4","A4","A4","A4","A4","A4","A2","A4","A1","A1","A5","A1","A5","A1","A1","A1","A3","A1","A1","A5","A1","A1","A6","A4","A4","A4","A4","A4","A4","A4","A5","A4","A4","A4","A5","A4","A4","A5","A1","A1","A5","A2","A1","A5","A1","A2","A1","A1","A1","A1","A1","A1","A1","A6","A6","A6","A6","A6","A5","A2","A4","A4","A4","A5","A2","A2","A5","A5","A4","A5","A4","A5","A5","A4","A5","A5","A4","A4","A5","A1","A3","A1","A5","A1","A1","A4","A1","A1","A2","A5","A1","A1","A1","A1","A5","A5","A6","A6","A6","A6","A4","A4","A4","A4","A4","A4","A5","A4","A4","A4","A4","A4","A4","A4","A4","A4","A4","A4","A5","A2","A4","A4","A4","A1","A5","A1","A5","A4","A2","A1","A5","A5","A4","A5","A5","A5","A1","A1","A5","A5","A5","A1")
B<-c(„B2","B2","B2","B2","B2","B2","B2","B2","B2","B2","B2","B2","B2","B4","B2","B2","B2","B2","B2","B2","B4","B3","B3","B3","B3","B3","B3","B3","B3","B1","B1","B1","B1","B1","B1","B1","B2","B2","B4","B4","B4","B3","B3","B3","B1","B1","B4","B4","B4","B4","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B2","B2","B2","B2","B2","B2","B2","B2","B2","B2","B2","B4","B4","B4","B4","B4","B4","B4","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B2","B4","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B2","B4","B4","B4","B4","B4","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B2","B2","B4","B4","B4","B4","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B3","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1","B1")
C<-c(0,0,0,0,0,1,0,1,0,0,0,0,1,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,4,0,1,1,0,0,0,0,0,0,26,3,58,30,7,0,0,1,0,2,3,14,0,5,10,2,3,8,11,7,2,9,22,31,27,10,4,9,1,4,1,10,1,2,35,1,11,5,3,0,2,6,6,4,1,5,5,2,6,0,0,1,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,7,1,0,0,101,29,15,1,37,25,35,27,21,2,14,40,4,22,27,3,0,2,7,1,4,0,0,0,0,3,9,2,3,3,0,1,13,5,8,52,13,45,10,19,5,10,20,0,4,13,1,0,0,0,0,0,0,0,5,1,0,0,1,4,0,1,1,0,2,1,1,0,1,0,0,14,14,4,7,5,14,1,6,16,2,4,30,5,6,3,4,6,0,8,0,0,0,0,0,0,0,0,0,1,1,0,0,0,0,0,1,0,0,0,0,1,1,0,0,0,0,1,0,1,2,0,1,0,0,0,1,1,0,1,0,0,1,0,0,0)
mydata<-data.frame(A,B,C)
Im negativ-binomialen Modell nutze ich nur die Zähldaten, die nicht 0 sind. Damit ist n=139. Im logit Modell ist n=256.
Am Beispiel des negbin Modells ist A2 (p=0.001), A5 (p=0.006) und A6 (p=0.007) gegenüber A1 signifikant (Das NA in Variable B4 kommt daher, weil alle A6 nur in B4 beobachtet wurden und umgekehrt bis auf eine Ausnahme auch. Somit sind beide Level kollinear. Das musste ich leider hier dulden und diskutieren)
Wenn ich jetzt mithilfe von lsmeans alle 15 Level-Kombinationen der 6 Level von A Holm-korrigiert haben möchte
dann ist nur noch A2 (p=0.022) signifikant, während A5 (p=0.0819) und A6 (NA) rausfallen. Das NA bei A6 rührt wahrscheinlich daher, dass A1 nur ein einziges mal mit A6 im selben B-Level (B4) vorkamen. Die anderen p-Werte sind jetzt um den Faktor 11 bis 15 erhöht. Katastrophe!
Daher meine Frage, ob ich lsmeans richtig verwende bzw. ob man die Modelle nicht mehrmals mit verschiedenen Referenzleveln rechnen kann und dann den p-Wert separat korrigiert?