similar to: svytable and na's

hi all, after fitting a multiple linear regression model <- lm(y ~ a + b+ c+d) i wanted to plot diagnostics plot(model) but get the error message Error in object$coefficients : $ operator is invalid for atomic vectors. which does not make a lot of sense, since there is no "$" - i am working with an attached dataset. can anyone help, please?? thanks a lot, kat

hi, i am running R with JGR and Deducer under Ubuntu Mint Lucid. since i updated to R 2.12 i can no longer make graphs - R just crashes. the full error report was posted to launchpad (https://bugs.launchpad.net/ubuntu/+source/openjdk-6/+bug/675905) without any success so far. anyone with a similar problem or helpful suggestions? thanks! kat

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

Hi all, I want to know how to perform the test Breslow-Day test for homogeneity of odds ratios (OR) stratified for svytable. This test is obtained with the following code: epi.2by2 (dat = daty, method = "case.control" conf.level = 0.95, units = 100, homogeneity = "breslow.day", verbose = TRUE) where "daty" is the object type table svytable consider it, but

Dear all I'm creating two way tables based on weighted survey data using svytable and want to include NAs as a level in the table. For unweighted data I use table(..., useNA="ifany"). Is there any equivalent for svytable? If not, is there a solution using a different function? Thanks Tim -- Tim Elwell-Sutton The University of Hong Kong Office: 3906 2053 Mob: 6084 5654

Hi, I am trying to write a function to compute many cross-tabulations with the -svytable- command. Here is a simplified example of the structure of my code (adapted from the -svytable- help file): data(api) func.example<-function(variable){ dclus1<-svydesign(id=~1, weights=~pw,data=apiclus1, fpc=~fpc) svytable(~ variable, dclus1) } When I call this function with:

hi, is there a way to make r open the help file in a separate window? i'm quite new to linux(kubuntu) and used to work with r under windows... thanks, katharina

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:

hi, for some unknown reason i cannot install the package ca in R running in a ubuntu mint system. i keep getting the following error message: configure: error: missing required header GL/gl.h ERROR: configuration failed for package ?rgl? * removing ?/home/kat/R/i486-pc-linux-gnu-library/2.12/rgl? ERROR: dependency ?rgl? is not available for package ?ca? * removing

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)

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

I have survey data that I am working on. I need to make some multi-way tables and regression analyses on the data. After attaching the data, this is the code I use for tables for four variables (sweight is the weight variable): > a <- xtabs(sweight~research.area + gender + a2n2 + age) > tmp <- ftable(a) Is this correct? I don't think I need to use the strata and cluster

Hi, is there a way to analyze subpopulations (e.g. women over 50, those who answered "yes" to a particular question) in a survey using Survey package? Other packages (e.g. Stata, SUDAAN) do this with a subpopulation option to identify the subpopulation for which the analysis shoud be done. I did not see this option in the Survey package. Is there another way to do this?

Hy Together I have compied Samba3.0beta2 on Solaris 9 . When i want to start the smbd daemon the folloing Error occured [2003/07/14 08:52:29, 3] smbd/server.c:main(777) loaded services [2003/07/14 08:52:29, 3] smbd/server.c:main(792) Becoming a daemon. [2003/07/14 08:52:29, 0] lib/messages.c:message_init(106) ERROR: Failed to initialise messages database Configure Option's env

Hy Together I have compied Samba3.0beta2 on Solaris 9 . When i want to start the smbd daemon the folloing Error occured [2003/07/14 08:52:29, 3] smbd/server.c:main(777) loaded services [2003/07/14 08:52:29, 3] smbd/server.c:main(792) Becoming a daemon. [2003/07/14 08:52:29, 0] lib/messages.c:message_init(106) ERROR: Failed to initialise messages database Configure Option's env

Hello. First, I am sorry to bother you. I found your email through doing a search on sites related to Switzerland and their related sites. I am a musician (guitarist/singer) as well as other things :-) and I have a website: www.rezamusic.com which I'd like to share with you -- it's a personal site, not a commercial one. It offers many topics, some of which you might find interesting. My

Hello I'm trying to get the proportion "true" for dichotomous variable for various subgroups in a survey. This works fine, but obviously doesn't give proportions directly: svytable(~SurvYear+problem.vandal, seh.dsn, round=TRUE) problem.vandal SurvYear FALSE TRUE 1995 8906 786 1997 17164 2494 1998 17890 1921 1999 18322 1669 2001 17623 2122 ...

Dear List, I have a dataset that provides sampling weights (National Survey of Family Growth 2002). I want to produce a cross-tabulation and use the provided sampling weights to obtain representative population estimates. (I believe they are simply frequency weights but codebook is uninformative). I can reproduce results (using this data) that were reported in a recent journal article, if I

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

Hi, I'm trying to get standard errors for some of the variables in my data frame. One of the questions on my survey is whether faculty coordinate across curriculum to include Arts Education as subject matter. All the responses are coded in zeros and ones obviously. For some of the other variables I have a 2 for those that responded with "Don't Know". I'm getting NA for