Разница между svyglm и взвешенным glm

Я хотел бы знать, в чем разница между использованием svyglm или взвешенного glm.

Например:

M1 = glm(formula = yy ~ age + gender + country , 
         family = binomial(link = "probit"), 
         data = P2013, 
         subset = (P2013$E27>=14 & P2013$E27<=17), 
         weights = P2013$PESOANO)

или определить дизайн образца как:

diseño = svydesign(id =~ NUMERO, 
                   strata =~ ESTRATOGEO, 
                   data = p2013, 
                   weights = P2013$PESOANO)

diseño_per_1417 = subset(diseño, (P2013$E27>=14 & P2013$E27<=17))

а затем используйте svyglm:

M2 = svyglm(formula = yy ~ age + gender + country, 
            family = quasibinomial(link = "probit"),
            data = P2013, 
            subset = (stratum=!0), 
            design = diseño_per_1417)

В случае, если я использую M2 (svyglm). Что я могу использовать для сравнения моделей, как пошагово для модели glm?

Спасибо, Наталья


r survey glm
person Natuk    schedule 29.04.2013    source источник

Ответы (1)


arrow_upward
3
arrow_downward

От 1_:

Веса, отличные от NULL, могут использоваться для указания того, что разные наблюдения имеют разные дисперсии (при этом значения весов обратно пропорциональны дисперсиям); или, что то же самое, когда элементы весов являются положительными целыми числами w_i, каждый ответ y_i является средним значением наблюдений единичного веса w_i. Для биномиальной GLM априорные веса используются для определения количества испытаний, когда ответ представляет собой долю успехов: они редко используются для GLM Пуассона.

Я не думаю, что вы ищете эти веса. Из вашего примера видно, что вы имеете дело со стратифицированным опросом. вы обязательно должны использовать surveyglm.

person Florian Oswald    schedule 01.05.2013
comment
Спасибо, Флориан!, но знаете ли вы, какой оценщик использует svyglm? или где я могу найти развитие в книге? - person Natuk; 02.05.2013
comment
svyglm выполняет некоторую работу, связанную с вашим дизайном (веса, слои и т. д.), но в конечном итоге вызывает glm. Так что вам следует заглянуть в ?glm. По сути, указав family, как вы это делаете в своем примере, вы указываете, какую функцию ссылки вы хотите использовать (что соответствует логиту или пробиту и т. д.). Краткое объяснение дано Kleiber & Zeileis 2008 (стр. 122). Также посетите веб-сайт пакета опроса. Пожалуйста, примите мой ответ, если вы думаете, что это помогает. - person Florian Oswald; 02.05.2013
comment
Да, я уже проверяю разные книги и пакеты опросов, но я не могу найти оценщик, который использует svyglm. Я думаю, что svyglm не оценивается по максимальному правдоподобию, может ли это быть оценщик Horvitz Thomson? Я хотел бы найти текст, который подробно объяснит. Спасибо! - person Natuk; 02.05.2013