search for: svysmooth

Hi all- So sorry to bother you all with something pretty basic. I am trying to add the lines method output from svysmooth to a svyplot with style="grayhex". However, the line either appears in the wrong place or if I am running in R Studio it causes the system to crash. I know this is something to do with Lattice graphics, but for the life of me I can not figure out how. Dr. Lumley in his excellent book on...

Greetings, I am using library(survey) to analyze some complex sample data. After fitting a model I tried to use termplot() with make.panel.svysmooth(), but I received an error (see below). Could someone help me interpret the error message so I can make the necessary corrections? The make.panel.svysmooth() function seems to work fine, and termplot() worked fine after I dropped the smoother. This led me to believe that the error is coming fr...

...row=c(2,2)) #Chart 1 options( survey.lonely.psu = "adjust" ) svyhist (~age_p, subset (nhis, xsmoke=='Current SMK'), breaks=MyBreaks, ylim = c(0,0.040), main= " ", col="red", xlab="Current Smoker's Age ") lines (svysmooth(~age_p, bandwidth=5,subset(nhis, xsmoke=='Current SMK')), lwd=2) #Chart 2 options( survey.lonely.psu = "adjust" ) svyhist (~age_p, subset (nhis, xsmoke=='Former SMK'), breaks=MyBreaks, ylim = c(0,0.040), main= " ", col=&quot...

...dclus2$variables$pc3 <- with(apiclus2, api99 * mycoef["api99"] + api00 * mycoef["api00"]) svycor(~ pc1 + pc2 + pc3, dclus2)[, ] # pc1 pc2 pc3 # pc1 1.0000000 0.2078789 0.2078789 # pc2 0.2078789 1.0000000 1.0000000 # pc3 0.2078789 1.0000000 1.0000000 plot(svysmooth(~ pc1, dclus2), xlim = c(-2.5, 5), ylim = 0:1) lines(svysmooth(~ pc2, dclus2), col = 2) lines(svysmooth(~ pc3, dclus2), col = 3) legend("topright", legend = c('pc$x[, "PC1"]', 'predict(pc, apiclus2)[, "PC1"]', 'Reconstituted indicator'), col = 1...

...t=TRUE) MyBreaks <- c(18,35,45,55,65,75,85,95) png("svyhist_age_at_death.png") svyhist (~dthage, subset (nhis, xspd2==2), breaks=MyBreaks, main= " ", col="grey80", xlab="Age at Death Distribution" ) lines (svysmooth(~dthage, bandwidth=5,subset(nhis, xspd2==2)), lwd=2) dev.off () Pradip K. Muhuri Statistician Substance Abuse & Mental Health Services Administration The Center for Behavioral Health Statistics and Quality Division of Population Surveys 1 Choke Cherry Road, Room 2-1071 Rockville, MD 2085...

...ain new features are - loglinear models with svyloglin() - database-backed designs now allow new variables to be created, support ODBC in addition to DBI database interfaces, and support multiple imputation - there is an approximation (Brewer,2002) for PPS sampling without replacement - svysmooth() now does quantile as well as mean smoothing - barplot() and dotchart() methods for many objects. More information at the package home page: http://faculty.washington.edu/tlumley/survey/ -thomas Thomas Lumley Assoc. Professor, Biostatistics tlumley at u.washington.edu University of Wash...

...ain new features are - loglinear models with svyloglin() - database-backed designs now allow new variables to be created, support ODBC in addition to DBI database interfaces, and support multiple imputation - there is an approximation (Brewer,2002) for PPS sampling without replacement - svysmooth() now does quantile as well as mean smoothing - barplot() and dotchart() methods for many objects. More information at the package home page: http://faculty.washington.edu/tlumley/survey/ -thomas Thomas Lumley Assoc. Professor, Biostatistics tlumley at u.washington.edu University of Wash...

I'm doing some work with a weighted data set (the CPS) and I'd like to use nonparametric techniques. Is there an R package that supports this? The np package doesn't appear to have provisions for weighted data. (Or have I missed something?) I need to perform both density estimation and regressions. thanks, bw