[R] how to extract p values in svyglm

[David Winsemius]

On Feb 12, 2012, at 6:39 PM, Tanu Soni wrote:

> summary(result)
>
> Call:
> svyglm(Injury ~ seat, sD, family = quasibinomial(link = "logit"))
>
> Survey design:
> svydesign(~1, prob = NULL, strata = Data[, 1], weights = Data[,
>    4], data = Data, fpc = ~fPc)
>
> Coefficients:
>             Estimate Std. Error t value Pr(>|t|)
> (Intercept) -4.256875   0.001421 -2996.7   <2e-16 ***
> seatbad      0.681504   0.001689   403.4   <2e-16 ***
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
> (Dispersion parameter for quasibinomial family taken to be 1.000095)
>
> Number of Fisher Scoring iterations: 6
>

Have you examined the output of str(summary(result))? Most versions of  
summary will retrun as one of theor list elements a matrix and you can  
generally access those matrices by name with "[[" and name, followed  
by extracting that result with "[" using a column name of "Pr(>|t|)".

-- 
David Winsemius, MD
West Hartford, CT