similar to: problem with svyglm

Until a weeks ago I used stata for everything. Now I'm learning R and trying to move. But, in this stage I'm testing R trying to do the same things than I used to do in stata whit the same outputs. I have a problem with the logit, applying weights. in stata I have this output . svy: logit bach job2 mujer i.egp4 programa delay mdeo i.str evprivate (running logit on estimation sample)

Dear R-help, I'd like ask for your opinion on choosing the "right" strategy for a particular dataset. We conducted 24-hour electric field measurements on 90 subjects. They are grouped by job (2 categories) and location (3 categories). There are four exposure metrics assigned to each subject. An excerpt from the data: n job location M OA UE all 0 job1 dist_200 0.297 0.072 0.171

Greetings all, I see the following issue during snap creation under IO: Transaction commit calls btrfs_start_delalloc_inodes() that locks the delalloc_inodes list, fetches the first inode, unlocks the list, triggers btrfs_alloc_delalloc_work/btrfs_queue_worker for this inode and then locks the list again. Then it checks the head of the list again. In my case, this is always exactly the same

I have been doing some experiments with rsync on btrfs, a copy-on-write file system that is approaching or having just achieved production-ready status depending on your requirements. For my purposes the reliability appears by almost all accounts to be there, and the compression alone makes it very compelling. However the following two experiments show rsync behaviors that are disappointing to

Sorry for the limited detail here. To make matters worse, I don't have access to the machines in question any more. But perhaps someone would be able to give me some pointers to useful avenues to pursue. I'm trying to use Wine on FreeBSD (versions 4.2 and 4.4), to act as a Lightwave 5.6 render node, using Lightwave's ScreamerNet software. The software is for NT 4.0. This software

foreign package does not support dataset saved at stata 13. anyone knows any wayt to make it works? [[alternative HTML version deleted]]

I'm performing item response theory with eRm packages I am excluding the persons that doesn't fit in the infit/outfit persons. for that I created a condition. then I have to create a new subset or matrix but with the condition. So: ORIGINAL matrix<-cbind(item1, item2, item3, item4) IF I PERFORM A head(matrix) item1 item2 item3 item4 3 2 3 1 3 1

Hi, I have a problem with the predict function of svyglm: When I try to predict with an object of svyglm I get the error message that R can't allocate a vector about 3 GB big. Has anybody a hint how I can reduce the memory used for this? I have about 25000 observation and use the weigths parameter. Where is so much data created? Thanks, Werner

summary(result) Call: svyglm(Injury ~ seat, sD, family = quasibinomial(link = "logit")) Survey design: svydesign(~1, prob = NULL, strata = Data[, 1], weights = Data[, 4], data = Data, fpc = ~fPc) Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -4.256875 0.001421 -2996.7 <2e-16 *** seatbad 0.681504 0.001689 403.4 <2e-16 *** ---

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

Dear list, I am using svyglm in the survey library to fit a binomial logistic regression accounting for sample design. The documentation says the models are not fit by maximum likelihood, so my question is what is the fitting method? Pseudo likelihood? Generalized least squares? Thanks for any insight. Best, Corey Corey Sparks Assistant Professor Department of Demography and

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)

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.

Hello I want to know if the summary of the logistic model with survey Pr (> | t |) to test if the coefficient of the model is significant, ie is the p_valor wald test for the model coefficients, for I am interested to know if the three levels of the variable educational level are significant to the model (significance of handling 0.2), I present below the results of my model

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]]

Hello, I want to know how to perform stepwise elimination of variables to svyglm thanks [[alternative HTML version deleted]]

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

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

Dear Professor Lumley I am relatively new to using R and also to logistic regression. We have analysed our Dudley Health Survey using the survey package. I am now trying to look at associations using svyglm but I am unsure of how to interpret the output and present the resulting model or whether there are any other things I should do to check the validity of the model. Below is an example of

Hi, I am trying to run analyzes incorporating sample weight, strata and cluster (three-stage sample) with PNS data (national health survey) and is giving error. I describe below the commands used. I could not make the code reproducible properly. Thanks, ################################################# library(survey)####change to 0 and 1 variable