2013-04-29 5 views
4

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)를 사용하는 경우에는

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

. glm 모델에서 stepwise와 같은 모델을 비교하는 데 사용할 수있는 것은 무엇입니까?

감사 help(glm) 가입일 나탈리

답변

2

는 :

비 NULL 가중치 다른 관찰 (가중치 분산액에 반비례의 값으로) 다른 분산액을 나타 내기 위해 사용될 수있다; 또는 동등하게, 가중치의 요소가 양의 정수 w_i 일 때, 각 응답 y_i는 w_i 단위 중량 관측의 평균입니다. 이항 GLM의 경우, 이전의 가중치는 응답이 성공의 비율 일 때 시행 횟수를 나타 내기 위해 사용됩니다. Poisson GLM에는 거의 사용되지 않습니다.

나는 당신이 그 무게를 찾고 있다고 생각하지 않습니다. 귀하의 사례를 보면 성층화 설문 조사를 다루고있는 것으로 보입니다. 당신은 확실히 surveyglm을 사용해야합니다.

+0

감사합니다. Florian !,하지만 견적서 사용자인 svyglm을 알고 계십니까? 또는 내가 책에서 발달을 찾을 수있는 곳? – Natuk

+1

'svyglm'은 디자인 (무게, 지층 등)과 관련이 있지만 결국에는'glm'을 호출합니다. 그래서'? glm'을 살펴 봐야합니다. 근본적으로 당신의 예제에서 당신이 사용하고자하는 링크 함수 ("logit"또는 "probit"등에 해당)를 설정하기 위해'family'를 지정함으로써. 간결한 설명은 Kleiber & Zeileis 2008 (p122)에서 제공합니다. 또한 설문 조사 패키지의 [웹 사이트] (http://staff.washington.edu/tlumley/survey/)를 확인하십시오. 도움이된다고 생각한다면 제 대답을 받아주십시오. –

+1

예, 이미 다른 책과 조사 패키지를 체크인했지만 svyglm이 사용하는 견적 도구를 찾을 수 없습니다. 나는 svyglm이 최대 가능성에 의해 평가되지 않는다고 생각합니다. Horvitz Thomson 견적서가 될 수 있습니까? 자세히 설명하는 텍스트를 찾고 싶습니다. 감사! – Natuk