What is the equivalent command in SPSS for he following R-command?
btbb <- glm(ybb ~ xbb - 1, family = binomial, weights = c(rep(13, 21)))
(no intercept)
Here is the data (SPSS)
http://www.datest.de/frag_einen_statistiker/Bradley_Terry/Beispiel_Pietsch.savand here is the R syntax:
> 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, 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)))
> btbb
Call: glm(formula = ybb ~ xbb - 1, family = binomial, weights = c(rep(13, 21)))
Coefficients:
xbb1 xbb2 xbb3 xbb4 xbb5 xbb6 xbb7
1.5814 1.4364 1.2945 1.2476 1.1077 0.6839 NA
Degrees of Freedom: 21 Total (i.e. Null); 15 Residual
Null Deviance: 49.7
Residual Deviance: 15.74 AIC: 87.32
Dr. Frank Gaeth