Displaying 20 results from an estimated 800 matches similar to: "(svy)glm and weights question"
2009 Oct 09
1
svy / weighted regression
Dear list,
I am trying to set up a propensity-weighted regression using the
survey package. Most of my population is sampled with a sampling
probability of one (that is, I have the full population). However, for
a subset of the data I have only a 50% sample of the full population.
In previous work on the data, I analyzed these data using SAS and
STATA. In those packages I used a propensity weight
2012 Jun 28
1
SVY: variance inflation factor VIF with complex survey
Hello,
Seeking a way to get the variance inflation factor VIF for a model of regression in complex survey, I have understood that without this package (SURVEY) RGui VIF obtained as follows:
fit <- lm(mpg~disp+hp+wt+drat, data=mtcars)
vif(fit)
But I want to know if survey, Vif is obtained so
vif( svyglm(api00~ell+meals+mobility, design=dstrat))
Thank you, happy day
2009 Jan 19
1
reference category for factor in regression
Hi all,
I am struggling with a strange issue in R that I have not encountered
before and I am not sure how to resolve this.
The model looks like this, with all irrelevant variables left out:
LABOUR - a dummy variable
NONLABOUR = 1 - LABOUR
AGE - a categorical variable / factor
VOTE - a dummy variable
glm(VOTE ~ 0 + LABOUR + NONLABOUR + LABOUR : AGE + NONLABOUR : AGE,
2011 Jul 04
1
Contrastes con el paquete survey (svycontrast)
Estimados usuarios:
Estoy intentando reproducir el ejemplo 6.4 de Thomas Lumley. Complex
Survey. Editorial Wiley. 2010 (ver la página en google:
2009 Mar 15
3
read.xls question
I'm an R newbie and had a question about the read.xls function. I've heard that this is often not a reliable function to use for importing data. However, I have created numerous xls files which contain information about voter turnout and macroeconomic indicators in India. I'm writing a paper on the relationship between economic growth and voter turnout.
This is the command I use:
2009 Apr 02
1
problem with svyglm()
Hello,
I'm trying to use the function svyglm in the library survey.
I create a data survey object:
data_svy<- svydesign(id=~PSU, strata=~sample_domain,
weights=~sample_weight, data=data, nest=TRUE)
and I try to use svyglm() with little success:
R<-svyglm(data_svy[,4]~(data_svy[,iCol]==listModality[[iVar]]
[iMod]),design=data_svy, family=binomial(link="logit")
Error in
2007 Oct 19
1
plot.Design
Dear R-users:
I am trying to use the following code to reproduce the figures on page 140
of Prof. Frank Harrell's book 'Regression Modeling Strategies':
rm(list=ls())
options(width=128)
library(Hmisc)
library(Design)
getHdata(counties)
counties$older <- counties$age6574 + counties$age75
label(counties$older) <- '% age >= 65, 1990'
counties$pdensity <-
2012 Dec 10
3
Warning message: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!
Hi there
I'm trying to fit a logistic regression model to data that looks very similar to the data in the sample below. I don't understand why I'm getting this error; none of the data are proportional and the weights are numeric values. Should I be concerned about the warning about non-integer successes in my binomial glm? If I should be, how do I go about addressing it?
I'm
2012 Aug 15
2
sensitivity and specificity in svyglm??
Hello,
As obtained from a table svyglm clasificaion, sensitivity and specificity. The funtion ConfusionMatrix () of the library (caret)
gives these results but not how to apply it to svyglm.
thanks
[[alternative HTML version deleted]]
2010 Apr 28
2
Generating a model fitness when score using svyglm?
Does anyone know how to calculated a BIC score (or an equivalent model
fitness score) when using svyglm for logistic regressions?
Thanks
Brad
--
View this message in context: http://r.789695.n4.nabble.com/Generating-a-model-fitness-when-score-using-svyglm-tp2069280p2069280.html
Sent from the R help mailing list archive at Nabble.com.
2013 May 02
1
Package survey: singularities in linear regression models
Hello,
I want to specify a linear regression model in which the metric outcome
is predicted by two factors and their interaction. glm() computes
effects for each factor level and the levels of the interaction. In the
case of singularities glm() displays "NA" for the corresponding
coefficients. However, svyglm() aborts with an error message. Is there a
possibility that svyglm()
2006 Jul 18
1
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
2008 Aug 06
1
Warning when using survey:::svyglm
Howdy,
Referencing the below exchange:
https://stat.ethz.ch/pipermail/r-help/2006-April/103862.html
I am still getting the same warning ("non-integer #successes in a
binomial glm!") when using svyglm:::survey. Using the API data:
library(survey)
data(api)
#stratified sample
dstrat<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)
2011 May 17
2
can not use plot.Predict {rms} reproduce figure 7.8 from Regression Modeling Strategies (http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf)
Dear R-users,
I am using R 2.13.0 and rms 3.3-0 , but can not reproduce figure 7.8 of the
handouts *Regression Modeling Strategies* (
http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf) by the
following code. Could any one help me figure out how to solve this?
setwd('C:/Rharrell')
require(rms)
load('data/counties.sav')
older <- counties$age6574 + counties$age75
2009 Oct 20
2
Weighted Logistic Regressions using svyglm
I?m running some logistic regressions and I?ve been trying to include weights
in the equation. However, when I run the model, I get this warning message:
Here?s what it says: Warning message: In eval(expr, envir, enclos) :
non-integer #successes in a binomial glm!
I think it is because the weights are non-integer values.
What is a good way to run logistic regressions in R when using
2008 Dec 19
1
svyglm and sandwich estimator of variance
Hi,
I would like to estimate coefficients using poisson regression and then get
standard errors that are adjusted for heteroskedasticity, using a complex
sample survey data. Then I will calculate prevalence ratio and confidence
intervals.
Can sandwich estimator of variance be used when observations aren?t
independent? In my case, observations are independent across groups
(clusters), but
2012 Nov 23
1
problem with svyglm
I have this problem.
test <- svydesign(id=~1,weights=~peso)
logit <- svyglm(bach ~ job2 + mujer + egp4 + programa + delay + mdeo + str
+ evprivate, family=binomial,design=test)
then appear:
Error in svyglm.survey.design(bach ~ job2 + mujer + egp4 + programa + :
all variables must be in design= argument
I don't know what this mean...
Please help.
Pablo.
[[alternative HTML
2015 Apr 16
4
Weighted Likelihood
¡Muchas gracias Olivier!
Un saludo.
El 16 de abril de 2015, 10:44, Olivier Nuñez <onunez en unex.es> escribió:
> Mira el paquete survey.
> Un saludo. Olivier
>
> ----- Mensaje original -----
> De: "Víctor Nalda Castellet" <victor.nalda.castellet en gmail.com>
> Para: "r-help-es" <r-help-es en r-project.org>
> Enviados: Miércoles, 15 de
2012 Jul 24
3
stepwise in svyglm???
Hello,
I want to know how to perform stepwise elimination of variables to svyglm
thanks
[[alternative HTML version deleted]]
2017 Jul 05
2
Svyglm Error
Greetings,
I am revisiting code from several different files I have saved from the past and all used to run flawlessly; now when I run any of the svyglm related functions, I am coming up with an error:
Error in model.frame.default(formula = F3ATTAINB ~ F1PARED, data = data, :
the ... list does not contain 4 elements
The following is a minimal reproducible example:
library(RCurl)