similar to: Svy function doesn't work nested in user-defined function

Hello, I?m hoping you have a few minutes to help out someone very new to R. I?ve done some searching, but cannot find this particular issue. I have survey data from two different time periods (years). Both years are stratified samples and have the same variables (and variable names), but are different people in the community answering in the different years. Everything loads into the survey

Dear Forum I a running svychisq from the survey package and get errors with the number of dimensions, errors that I do not understand and do not know how to fix. I ask you kindly to help me out. The eror message follows with some information below. I hope there are enough information to help me to fix the problem if not please let me know what is needed Best Stefan Hrafn Jonsson > QHISQ2

Dear All I am trying to use svychisq with a two-dimensional table 4 x 5. The command I am using is summary(svytable(~dietperception+dietstatus,dudleyls1rake,na.rm=TRUE),"C hisq") It is throwing up an error message as follows: Error in NCOL(y) : only 0's may be mixed with negative subscripts In addition: Warning messages: 1: is.na() applied to non-(list or vector) in:

Hello forum, I want to do a test of independence with svychisq, but I get an error, then this my code: am18 <- read.spss("C:/Users/diana/Dropbox/Semestre 10/Tesis 10/Tesis Diana/AMcomuna18-29MAR2012.sav", use.value.labels=TRUE, max.value.labels=Inf, to.data.frame=TRUE) b<-matrix(c(am18$N6_MANZANA),ncol=1) c<-matrix(c(am18$PM1_1_PONDEMUESTRA),ncol=1)

Dear Md Kamruzzaman, To answer your second question first, you could just use the svychisq() function. The difference-of-proportion test is equivalent to a chisquare test for the 2-by-2 table. You don't say how you computed the confidence intervals for the two separate proportions, but if you have their standard errors (and if not, you should be able to infer them from the confidence

[url]http://svy.mk/g7RhjY Please take 10 minutes to complete this wine survey. Your feedback will inform a Columbia University graduate student's thesis project for boxed wine. Your answers will be completely anonymous.If you would like to be entered into a prize draw for a $50 Amazon gift card, please provide your email (optional) when you complete the survey. Winner will be notified

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

Hello, excuse me if this is trivial. I have some survey data with the following design MyDesign<-(id=ident, weights=~fwgt, fcp=~nval1+nval2, data=hh_data1) I would like to run simple frequency tables such as those provided by the freq() function of rgrs to check the impact of the sampling design on my Ns. Anyway I can do this? thanks, V [[alternative HTML version deleted]]

Dear all, I am looking for an equivalent of the pairwise.t.test() function using a complex survey design. Any suggestions on this or how to account for weights using the pairwise.t.test itself? Thanks, V [[alternative HTML version deleted]]

Hello forum, I want to find in survey how to test V-cramer and coefficient Goodman-Kruskal gamma, to test the collinearity between varibles to enter a model. For the V-Cramer function is used with the package vcd assocstats but without considering the survey package. I appreciate the cooperation. Att. Diana Martinez [[alternative HTML version deleted]]

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

Hi all, I am running a set of logistic regressions, where we want to use some weights, and I am not sure whether what I am doing is reasonable or not. The dependent variable is turnout in an election - i.e. survey respondents were asked whether or not they voted. The percentage of those who say they voted is much higher than the actual turnout, probably due both to non-response bias and social

Hello Everyone, I was analysing big survey data using survey packages on RStudio. Survey package allows survey data analysis with the design effect.The survey package included functions for all other statistical analysis except two-proportion z tests. I was trying to calculate the difference in prevalence of Diabetes and Prediabetes between the year 2011 and 2017 (with 95%CI). I was able to

Dear Md Kamruzzaman, I've copied this response to the r-help list, where you originally asked your question. That way, other people can follow the conversation, if they're interested and there will be a record of the solution. Please keep r-help in the loop See below: On 2024-01-17 9:47 p.m., Md. Kamruzzaman wrote: > > Caution: External email. > > > Dear John >

Dear All, I have a problem understanding the difference between the outcome of a fisher exact test and a chi-square test (with simulated p.value). For some sample data (see below), fisher reports p=.02337. The normal chi-square test complains about "approximation may be incorrect", because there is a column with cells with very small values. I therefore tried the chi-square with

I am trying to use the reshape package for the first time. I have two waves of a survey, so the id variables include a subject identification number and a variable denoting the wave of the survey. I used the following arguments: library(reshape) svy.melt <- melt(svy, id=c("id", "WAVE")) svy.wide <- cast(svy.melt, id ~ WAVE + ...) and got the following error:

I have two waves of a survey given to students at various middle schools and high schools, with student id numbers for each student. I am having difficulty reshaping the file from long to wide. My code is below: library(foreign) svy <- read.spss("studsur4.SAV") svy.wide <- reshape(svy, timevar="WAVE", idvar="id", direction="wide")

¡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

I'm studying the calibration function in the survey package in preparation for raking some survey data. Results from the rake function below agree with other sources. When I run calibrate, I get a warning message and the M and F weights seem to be reversed. Even allowing for that, the deviation between calibrated and raked weights is much more than I expected. I see that in the calibrate

Dear R users, Can anyone please tell me how to generate a large number of samples in R, given certain distribution and size. For example, if I want to generate 1000 samples of size n=100, with a N(0,1) distribution, how should I proceed? (Since I dont want to do "rnorm(100,0,1)" in R for 1000 times) Thanks for help Debbie