В книге Добсона и Барни «Введение в обобщенные линейные модели - третье издание» на стр. 161 (база данных находится здесь [bit.ly/2NtJFAw]);Я пытаюсь повторить некоторые результаты, но знаки оценок меняются местами и их значения одинаковы.
library(readxl) #para usar la función "read_excel"
car_pref <- read_excel('car.xls')
car_pref <- car_pref[-1,]
car_pref <- car_pref[-1,]
colnames(car_pref) <- c("sex","age","response","frequency")
car_pref$frequency <- as.numeric(car_pref$frequency)
car_pref$sex <- factor(car_pref$sex)
car_pref$age <- factor(car_pref$age)
car_pref$response <- factor(car_pref$response)
car_pref$sex <- relevel(car_pref$sex, ref="women")
car_pref$age <- relevel(car_pref$age, ref="24-40")
car_pref$age <- relevel(car_pref$age, ref="18-23")
car_pref$response <- relevel(car_pref$response, ref="very important")
library(MASS)
res.polr <- polr(factor(response)~factor(age) + factor(sex), weights=frequency,data=car_pref)
summary(res.polr)
Изображение книги
Мои результаты
Обратите внимание, что мои оценки имеют противоположный знак.Я подозреваю, что это значения повторной параметризации, но я не вижу matrix.model
для подтверждения.Любая идея?