In attempting to use the svyglm call in the R Survey Package, I am receiving the error: Error in pwt[i] : invalid subscript type 'list' I have not been able to find a lot of information on how to resolve the error; one source advised it was related to how the subsetting command was executed. This was my initial attempt: mc1 <- svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset(elsq1ch_brr,BYSCTRL==1 & G10COHRT==1),na.action) summary(mc1) This was my second approach trying to change up how I had subsetted the data: summary(mc1) samp1 <- subset(elsq1ch_brr,BYSCTRL==1 & G10COHRT==1) dim(samp1) mc1 <- svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset=samp1,na.action) summary(mc1)? Both attempts resulted in the same error stated above. Any advisement in how to resolve this error would be greatly appreciated. Sincerely, Courtney Benjamin ? 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]]
hi could you make this a minimal reproducible example? On Sep 24, 2016 12:03 PM, "Courtney Benjamin" <cbenjami at btboces.org> wrote:> In attempting to use the svyglm call in the R Survey Package, I am > receiving the error: Error in pwt[i] : invalid subscript type 'list' > > I have not been able to find a lot of information on how to resolve the > error; one source advised it was related to how the subsetting command was > executed. > > This was my initial attempt: > > mc1 <- svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset(elsq1ch_brr,BYSCTRL==1 > & G10COHRT==1),na.action) > summary(mc1) > This was my second approach trying to change up how I had subsetted the > data: > summary(mc1) > samp1 <- subset(elsq1ch_brr,BYSCTRL==1 & G10COHRT==1) > dim(samp1) > mc1 <- svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset> samp1,na.action) > summary(mc1)? > > Both attempts resulted in the same error stated above. Any advisement in > how to resolve this error would be greatly appreciated. > Sincerely, > Courtney Benjamin > > ? > > > > 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]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/ > posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]
Dear Courtney,
I think that you're confused about how to use the subset argument to
svyglm() and about what the subset() function returns. The subset argument
should be a logical expression, evaluating to TRUE or FALSE for each case;
subset() returns a data set (e.g., a "survey.design" object).
Here's an example adapted from ?svyglm:
------------------- snip-------------------
data(api)
dstrat<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)
summary(svyglm(api00~ell+meals+mobility, design=dstrat))
summary(svyglm(api00~ell+meals+mobility, design=dstrat, subset = pcttest >
90))
summary(svyglm(api00~ell+meals+mobility, design=subset(dstrat, subset = pcttest
> 90)))
------------------- snip-------------------
The last two commands produce the same results.
I hope this helps,
John
-----------------------------
John Fox, Professor
McMaster University
Hamilton, Ontario
Canada L8S 4M4
Web: socserv.mcmaster.ca/jfox
> -----Original Message-----
> From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Courtney
> Benjamin
> Sent: September 23, 2016 11:01 PM
> To: r-help at r-project.org
> Subject: [R] Svyglm Error in Survey Package
>
> In attempting to use the svyglm call in the R Survey Package, I am
receiving the
> error: Error in pwt[i] : invalid subscript type 'list'
>
> I have not been able to find a lot of information on how to resolve the
error;
> one source advised it was related to how the subsetting command was
> executed.
>
> This was my initial attempt:
>
> mc1 <-
> svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset(elsq1ch_brr,
> BYSCTRL==1 & G10COHRT==1),na.action)
> summary(mc1)
> This was my second approach trying to change up how I had subsetted the
> data:
> summary(mc1)
> samp1 <- subset(elsq1ch_brr,BYSCTRL==1 & G10COHRT==1)
> dim(samp1)
> mc1 <-
> svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset=samp1,na.ac
> tion)
> summary(mc1)?
>
> Both attempts resulted in the same error stated above. Any advisement in
how
> to resolve this error would be greatly appreciated.
> Sincerely,
> Courtney Benjamin
>
> ?
>
>
>
> 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]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.
Dear Courtney, You're confusing a function call, na.action(na.omit), with an argument specification, na.action=na.omit (and, in any event, there is no na.action() function). But you don't have to specify na.action=na.omit, because na.omit (which produces a complete-case analysis) is the default na.action. See ?glm for more information. If you haven't already done so, you might want to read something about how statistical modeling functions in R -- and possibly R more generally -- work. In the long run, that likely will save you some time. I'm cc'ing this to the r-help email list, where you posed your original question. It's generally a good idea to keep responses on the list. I hope this helps, John> -----Original Message----- > From: Courtney Benjamin [mailto:cbenjami at BTBOCES.ORG] > Sent: September 24, 2016 10:55 PM > To: Fox, John <jfox at mcmaster.ca> > Subject: Re: Svyglm Error in Survey Package > > Hello Dr. Fox, > Thank you very much for your explanation; I am in better shape now with the > subset argument and I am no longer getting the error. Now I am not having > success with specifying the na.action argument. I would like to exclude NAs > like I have done in other commands with na.rm=TRUE. > The following are my attempts: > > > summary(svyglm(F3ATTAINMENT~F1PARED+F1SES2QU,design=elsq1ch_brr,su > bset > > =BYSCTRL==1 & G10COHRT==1,na.action)) > Error in glm(formula = F3ATTAINMENT ~ F1PARED + F1SES2QU, subset > BYSCTRL == : > argument "na.action" is missing, with no default > > > > summary(svyglm(F3ATTAINMENT~F1PARED+F1SES2QU,design=elsq1ch_brr,su > bset > > =BYSCTRL==1 & G10COHRT==1,na.action(na.omit)) > + > Any guidance you would be willing to provide is greatly appreciated; I have > only been using R for about 6 months. > Sincerely > Courtney > > 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 > 607-763-8633 > > ________________________________________ > From: Fox, John <jfox at mcmaster.ca> > Sent: Saturday, September 24, 2016 9:38 AM > To: Courtney Benjamin > Cc: r-help at r-project.org > Subject: RE: Svyglm Error in Survey Package > > Dear Courtney, > > I think that you're confused about how to use the subset argument to svyglm() > and about what the subset() function returns. The subset argument should be a > logical expression, evaluating to TRUE or FALSE for each case; subset() returns > a data set (e.g., a "survey.design" object). > > Here's an example adapted from ?svyglm: > > ------------------- snip------------------- > data(api) > dstrat<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc) > summary(svyglm(api00~ell+meals+mobility, design=dstrat)) > summary(svyglm(api00~ell+meals+mobility, design=dstrat, subset = pcttest > > 90)) summary(svyglm(api00~ell+meals+mobility, design=subset(dstrat, subset > = pcttest > 90))) > ------------------- snip------------------- > > The last two commands produce the same results. > > I hope this helps, > John > > ----------------------------- > John Fox, Professor > McMaster University > Hamilton, Ontario > Canada L8S 4M4 > Web: socserv.mcmaster.ca/jfox > > > > > -----Original Message----- > > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of > > Courtney Benjamin > > Sent: September 23, 2016 11:01 PM > > To: r-help at r-project.org > > Subject: [R] Svyglm Error in Survey Package > > > > In attempting to use the svyglm call in the R Survey Package, I am > > receiving the > > error: Error in pwt[i] : invalid subscript type 'list' > > > > I have not been able to find a lot of information on how to resolve > > the error; one source advised it was related to how the subsetting > > command was executed. > > > > This was my initial attempt: > > > > mc1 <- > > > svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset(elsq1ch_brr, > > BYSCTRL==1 & G10COHRT==1),na.action) > > summary(mc1) > > This was my second approach trying to change up how I had subsetted > > the > > data: > > summary(mc1) > > samp1 <- subset(elsq1ch_brr,BYSCTRL==1 & G10COHRT==1) > > dim(samp1) > > mc1 <- > > > svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset=samp1,na.ac > > tion) > > summary(mc1)? > > > > Both attempts resulted in the same error stated above. Any advisement > > in how to resolve this error would be greatly appreciated. > > Sincerely, > > Courtney Benjamin > > > > ? > > > > > > > > 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]] > > > > ______________________________________________ > > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide http://www.R-project.org/posting- > > guide.html and provide commented, minimal, self-contained, > > reproducible code.
Hello Dr. Fox, Thank you very much for your recommendations and I will take action on following them. I appreciate your willingness to provide guidance on my novice questions. Sincerely, Courtney ________________________________________ From: Fox, John <jfox at mcmaster.ca> Sent: Sunday, September 25, 2016 9:21 AM To: Courtney Benjamin Cc: r-help at r-project.org Subject: RE: Svyglm Error in Survey Package Dear Courtney, You're confusing a function call, na.action(na.omit), with an argument specification, na.action=na.omit (and, in any event, there is no na.action() function). But you don't have to specify na.action=na.omit, because na.omit (which produces a complete-case analysis) is the default na.action. See ?glm for more information. If you haven't already done so, you might want to read something about how statistical modeling functions in R -- and possibly R more generally -- work. In the long run, that likely will save you some time. I'm cc'ing this to the r-help email list, where you posed your original question. It's generally a good idea to keep responses on the list. I hope this helps, John> -----Original Message----- > From: Courtney Benjamin [mailto:cbenjami at BTBOCES.ORG] > Sent: September 24, 2016 10:55 PM > To: Fox, John <jfox at mcmaster.ca> > Subject: Re: Svyglm Error in Survey Package > > Hello Dr. Fox, > Thank you very much for your explanation; I am in better shape now with the > subset argument and I am no longer getting the error. Now I am not having > success with specifying the na.action argument. I would like to exclude NAs > like I have done in other commands with na.rm=TRUE. > The following are my attempts: > > > summary(svyglm(F3ATTAINMENT~F1PARED+F1SES2QU,design=elsq1ch_brr,su > bset > > =BYSCTRL==1 & G10COHRT==1,na.action)) > Error in glm(formula = F3ATTAINMENT ~ F1PARED + F1SES2QU, subset > BYSCTRL == : > argument "na.action" is missing, with no default > > > > summary(svyglm(F3ATTAINMENT~F1PARED+F1SES2QU,design=elsq1ch_brr,su > bset > > =BYSCTRL==1 & G10COHRT==1,na.action(na.omit)) > + > Any guidance you would be willing to provide is greatly appreciated; I have > only been using R for about 6 months. > Sincerely > Courtney > > 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 > 607-763-8633 > > ________________________________________ > From: Fox, John <jfox at mcmaster.ca> > Sent: Saturday, September 24, 2016 9:38 AM > To: Courtney Benjamin > Cc: r-help at r-project.org > Subject: RE: Svyglm Error in Survey Package > > Dear Courtney, > > I think that you're confused about how to use the subset argument to svyglm() > and about what the subset() function returns. The subset argument should be a > logical expression, evaluating to TRUE or FALSE for each case; subset() returns > a data set (e.g., a "survey.design" object). > > Here's an example adapted from ?svyglm: > > ------------------- snip------------------- > data(api) > dstrat<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc) > summary(svyglm(api00~ell+meals+mobility, design=dstrat)) > summary(svyglm(api00~ell+meals+mobility, design=dstrat, subset = pcttest > > 90)) summary(svyglm(api00~ell+meals+mobility, design=subset(dstrat, subset > = pcttest > 90))) > ------------------- snip------------------- > > The last two commands produce the same results. > > I hope this helps, > John > > ----------------------------- > John Fox, Professor > McMaster University > Hamilton, Ontario > Canada L8S 4M4 > Web: socserv.mcmaster.ca/jfox > > > > > -----Original Message----- > > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of > > Courtney Benjamin > > Sent: September 23, 2016 11:01 PM > > To: r-help at r-project.org > > Subject: [R] Svyglm Error in Survey Package > > > > In attempting to use the svyglm call in the R Survey Package, I am > > receiving the > > error: Error in pwt[i] : invalid subscript type 'list' > > > > I have not been able to find a lot of information on how to resolve > > the error; one source advised it was related to how the subsetting > > command was executed. > > > > This was my initial attempt: > > > > mc1 <- > > > svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset(elsq1ch_brr, > > BYSCTRL==1 & G10COHRT==1),na.action) > > summary(mc1) > > This was my second approach trying to change up how I had subsetted > > the > > data: > > summary(mc1) > > samp1 <- subset(elsq1ch_brr,BYSCTRL==1 & G10COHRT==1) > > dim(samp1) > > mc1 <- > > > svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset=samp1,na.ac > > tion) > > summary(mc1)? > > > > Both attempts resulted in the same error stated above. Any advisement > > in how to resolve this error would be greatly appreciated. > > Sincerely, > > Courtney Benjamin > > > > ? > > > > > > > > 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]] > > > > ______________________________________________ > > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide http://www.R-project.org/posting- > > guide.html and provide commented, minimal, self-contained, > > reproducible code.