Für die Lösung eines Bradley – Terry Modells (logistisches Modell) liegt eine Syntax in R vor sowie der gesamte Lösungsweg und die zugehörigen Daten. Problem: Wie würde die äquivalente SPSS-Befehlssyntax dafür lauten? Ich bin leider daran gescheitert, die Ergebnisse in SPSS zu reproduzieren. Wie würde denn die Syntax dazu lauten, bzw. zumindest der Lösungsweg. Hatte versucht über GENLIN das Modell zu berechnen, bin allerdings in Log. Modellen zu wenig eingearbeitet.
Der Quellcode in R:
Code: Alles auswählen
ybb <- c(7,9,7,7,9,11,7,5,11,9,9,7,7,8,12,6,7,10,7,12,6)/13
xbb <- matrix(c(1,-1,0,0,0,0,0,
+ 1,0,-1,0,0,0,0,
+ 1,0,0,-1,0,0,0,
+ 1,0,0,0,-1,0,0,
+ 1,0,0,0,0,-1,0,
+ 1,0,0,0,0,0,-1,
+ 0,1,-1,0,0,0,0,
+ 0,1,0,-1,0,0,0,
+ 0,1,0,0,-1,0,0,
+ 0,1,0,0,0,-1,0,
+ 0,1,0,0,0,0,-1,
+ 0,0,1,-1,0,0,0,
+ 0,0,1,0,-1,0,0,
+ 0,0,1,0,0,-1,0,
+ 0,0,1,0,0,0,-1,
+ 0,0,0,1,-1,0,0,
+ 0,0,0,1,0,-1,0,
+ 0,0,0,1,0,0,-1,
+ 0,0,0,0,1,-1,0,
+ 0,0,0,0,1,0,-1,
+ 0,0,0,0,0,1,-1), nrow=21, ncol=7, byrow=TRUE)
btbb <- glm(ybb ~ xbb - 1, family = binomial, weights = c(rep(13, 21)))
Call: glm(formula = ybb ~ xbb - 1, family = binomial, weights = c(rep(13, 21)))Code: Alles auswählen
Coefficients:
xbb1 xbb2 xbb3 xbb4 xbb5 xbb6 xbb7
1.5814 1.4364 1.2945 1.2476 1.1077 0.6839 NA
SPSS Versuch (Variablennamen xbb1-xbb7 in Städtnamen geändert):
Code: Alles auswählen
GENLIN ybb BY Milwaukee Detroit Toronto NewYork Boston Cleveland Baltimore Vergleich zeile
(ORDER=ASCENDING)
/MODEL Milwaukee Detroit Toronto NewYork Boston Cleveland Baltimore Vergleich zeile INTERCEPT=YES
DISTRIBUTION=POISSON LINK=LOG
/CRITERIA METHOD=FISHER(1) SCALE=1 COVB=MODEL MAXITERATIONS=100 MAXSTEPHALVING=5
PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 CITYPE=WALD
LIKELIHOOD=FULL
/MISSING CLASSMISSING=EXCLUDE
/PRINT SOLUTION (EXPONENTIATED).Sämtliche Daten (SPSS / R ) und Darstellung des Modells (pdf)
http://www.frag-einen-statistiker.de/Br ... modell.zip
