[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