[R] Significance of Svyrepdesign Object Warning

Courtney Benjamin cbenjami at BTBOCES.ORG
Sun Oct 23 17:56:33 CEST 2016


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 matched the published standard errors of the data set.

When doing a logistic regression through the svyglm call, I receive the following warning:

In object$survey.design$pweights * presid^2 :
  longer object length is not a multiple of shorter object length?
I have search around quite a bit online and have not been able to find any good interpretation of its meaning.  I want to be sure that I am not making some type of mistake that is causing this warning to be produced.  Any advisement is greatly appreciated.
The following is an MRE that can be pasted into the R console:
library(RCurl)
library(survey)
data <- getURL("https://raw.githubusercontent.com/cbenjamin1821/careertech-ed/master/elsq1adj.csv")
elsq1ch <- read.csv(text = data)
#Specifying the svyrepdesign object which applies the BRR weights
elsq1ch_brr<-svrepdesign(variables = elsq1ch[,1:16], repweights = elsq1ch[,18:217], weights = elsq1ch[,17], combined.weights = TRUE, type = "BRR")
elsq1ch_brr
#Logistic regression call which yields a warning regarding svyrepdesign object
svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1,na.action=na.exclude)
allCC <- summary(svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1,na.action=na.exclude))
allCC

#Session Info
#R version 3.3.1 (2016-06-21)
#Platform: x86_64-w64-mingw32/x64 (64-bit)
#Running under: Windows >= 8 x64 (build 9200)

#locale:
#  [1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252
#[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
#[5] LC_TIME=English_United States.1252
#attached base packages:
#  [1] grid      stats     graphics  grDevices utils     datasets  methods   base
#other attached packages:
#[1] survey_3.31-2   survival_2.39-4 Matrix_1.2-6    RCurl_1.95-4.8  bitops_1.0-6
#loaded via a namespace (and not attached):
#[1] tools_3.3.1     splines_3.3.1   knitr_1.14      lattice_0.20-33


Courtney Benjamin

Broome-Tioga BOCES

Automotive Technology II Teacher

Located at Gault Toyota

Doctoral Candidate-Educational Theory & Practice

State University of New York at Binghamton

cbenjami at btboces.org<mailto:cbenjami at btboces.org>

607-763-8633

	[[alternative HTML version deleted]]



More information about the R-help mailing list