[R] vcov error in svyby and svytable functions
Hi, I'm trying to compute survey svytable statistic on subsets by using the svyby function. Here is the code: b-svyby(~V024+V751, by=~V025, design=strat2, svytable, round=TRUE) The vars, V024, V751 and V025 are factors. The by var has 2 levels, and hence there will be two subsets. strat2 is created by the svydesign function. It's giving me the following error: b-svyby(~V024+V751, by=~V025, design=strat2, svytable, round=TRUE) Error in vcov(object, ...) : no applicable method for vcov I can't understand what vcov is and why it is giving this error because both svyby and svytable functions do not have vcov as one of their arguements. I tried svytable without the svyby and it works fine. Any help on this would be appreciated. Thanks, Debarchana. -- Debarchana Ghosh Research Assistant, Department of Geography University of Minnesota. PH: 8143607580 www.tc.umn.edu/~ghos0033/ [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list 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.
[R] Lattice package par.settings/trellis.par.settings questions
Hi All, I'm trying to modify some of the default graphic parameters in a conditional histogram. While I was able to change the default grey background to white, I couldn't change the axis.font or the xlab font. I used the following code: /histogram(~V751|V013+V025, finalbase, xlab=Heard of HIV/AIDS (No/Yes), col=c(cyan,magenta), par.settings=list(background=white)) /The arguments for example like /axis.font=2/, or /cex=2/ are not working in the /par.settings(). /I also tried to read the manual of /trellis.par.settings()/ but didn't understand how to use it and where exactly to put it. Any help with this will be appreciated. Thanks, Debarchana. -- Debarchana Ghosh Research Assistant Department of Geography University of Minnesota PH: 8143607580 email to: [EMAIL PROTECTED] www.tc.umn.edu/~ghos0033 __ R-help@stat.math.ethz.ch mailing list 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.
[R] Survey-weighted ordered logistic regression
Hi, I am trying to fit a model with an ordered response variable (3 levels) and 13 predictor variables. The sample has complex survey design and I've used 'svydesign' command from the survey package to specify the sampling design. After reading the manual of 'svyglm' command, I've found that you can fit a logistic regression (binary response variable) by specifying the family=binomial in svyglm function. However I'm unable to fit an ordered logistic model in 'svyglm function'. Any help with this will be appreciated. Thanks, Debarchana. -- Debarchana Ghosh Research Assistant, Department of Geography University of Minnesota. PH: 8143607580 www.tc.umn.edu/~ghos0033/ [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list 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.
[R] Problem with ordered logistic regression using polr function.
Hi, I'm trying to fit a ordered logistic regression. The response variable (y) has three levels (0,1,2). The command I've used is: /ordlog-polr(y~x1+x2+x3+x4, data=finalbase, subset=heard, weight=wt, na.action=na.omit) / (There are no NA's in y but there are NA's in X's) The error I'm getting is: /Warning messages: 1: non-integer #successes in a binomial glm! in: eval(expr, envir, enclos) 2: design appears to be rank-deficient, so dropping some coefs in: polr(right.ans ~ S107 + children + work + rel + media + V013 + / Also, if I write /summary(ordlog)/ I'm getting the following error: /Re-fitting to get Hessian Error in if (all(pr 0)) -sum(wt * log(pr)) else Inf : missing value where TRUE/FALSE needed /Could anybody point out the problem? Thanks, Debarchana. -- Debarchana Ghosh Research Assistant Department of Geography University of Minnesota PH: 8143607580 email to: [EMAIL PROTECTED] www.tc.umn.edu/~ghos0033 __ R-help@stat.math.ethz.ch mailing list 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.
[R] Problem with ordered logistic regression using polr function
Hi, I'm trying to fit a ordered logistic regression. The response variable (y) has three levels (0,1,2). The command I've used is: ordlog-polr(y~x1+x2+x3+x4, data=finalbase, subset=heard, weight=wt, na.action=na.omit) (There are no NA's in y but there are NA's in X's) The error I'm getting is: Warning messages: 1: non-integer #successes in a binomial glm! in: eval(expr, envir, enclos) 2: design appears to be rank-deficient, so dropping some coefs in: polr(right.ans ~ S107 + children + work + rel + media + V013 + Also, if I write summary(ordlog) I'm getting the following error: Re-fitting to get Hessian Error in if (all(pr 0)) -sum(wt * log(pr)) else Inf : missing value where TRUE/FALSE needed Could anybody point out the problem? Thanks, Debarchana. -- Debarchana Ghosh Research Assistant Department of Geography University of Minnesota PH: 8143607580 email to: [EMAIL PROTECTED] www.tc.umn.edu/~ghos0033 __ R-help@stat.math.ethz.ch mailing list 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.
[R] Ordered Logistic Regression in survey command
Hi, How can I do ordered logistic regression in svyglm? Thanks, D. -- Debarchana Ghosh Research Assistant Department of Geography University of Minnesota PH: 8143607580 email to: [EMAIL PROTECTED] www.tc.umn.edu/~ghos0033 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
