Re: [R] Significance of Svyrepdesign Object Warning

2016-10-27 Thread Courtney Benjamin
7-763-8633 From: William Dunlap <wdun...@tibco.com> Sent: Thursday, October 27, 2016 9:59 PM To: Courtney Benjamin Cc: Anthony Damico; r-help@r-project.org; Thomas Lumley Subject: Re: [R] Significance of Svyrepdesign Object Warning For, now I would just use na.action

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-27 Thread William Dunlap via R-help
mp; Practice > > State University of New York at Binghamton > > cbenj...@btboces.org > > 607-763-8633 > ------ > *From:* William Dunlap <wdun...@tibco.com> > *Sent:* Sunday, October 23, 2016 2:24 PM > *To:* Anthony Damico > *Cc:* Cou

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-27 Thread Courtney Benjamin
William Dunlap <wdun...@tibco.com> Sent: Sunday, October 23, 2016 2:24 PM To: Anthony Damico Cc: Courtney Benjamin; r-help@r-project.org; Thomas Lumley Subject: Re: [R] Significance of Svyrepdesign Object Warning The immediate problem could be solved by changing the following lines in surv

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-23 Thread Courtney Benjamin
607-763-8633 From: William Dunlap <wdun...@tibco.com> Sent: Sunday, October 23, 2016 2:24 PM To: Anthony Damico Cc: Courtney Benjamin; r-help@r-project.org; Thomas Lumley Subject: Re: [R] Significance of Svyrepdesign Object Warning The immediate pro

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-23 Thread William Dunlap via R-help
The immediate problem could be solved by changing the following lines in survey:::summary.svrepglm from presid <- resid(object, "pearson") dispersion <- sum(object$survey.design$pweights * presid^2, na.rm = TRUE)/sum(object$survey.design$pweights) to presid <- resid(object,

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-23 Thread Anthony Damico
hi, great example. i am ccing survey package author/maintainer dr. lumley. why do you have `na.action=na.exclude`? if you remove it, things work as expected-- library(RCurl) library(survey) data <- getURL("

[R] Significance of Svyrepdesign Object Warning

2016-10-23 Thread Courtney Benjamin
Hello R Users, I am using Lumley's Survey Package in R to analyze complex survey data that involves 200 balanced repeated replicate (BRR) weight variables. I have ensured that my svyrepdesign object that specifies the application of the BRR weights to the data set is accurate and I have