svyboxplot - library (survey)
Hello,
I understand that svyhist () provides density histograms with density values on the y-axis (R code shown below). Is there a way one can have relative relative frequency histograms with relative freqencies on the y-axis?
Any advice/help would be appreciated.
Thanks,
Pradip Muhuri
###### svyhist - Density Histogram
options( survey.lonely.psu = "adjust" )
svyhist (~age_p,
subset (nhis, xspd2=='SPD'), breaks=MyBreaks,
ylim = c(0,0.040),
main= " ",
col="red",
xlab="Age at Interview (SPD Category)"
)
lines (svysmooth(~age_p, bandwidth=5,subset(nhis, xspd2=='SPD')), lwd=2)
From: Anthony Damico [ajdamico at gmail.com]
Sent: Monday, October 01, 2012 10:07 AM
To: Muhuri, Pradip (SAMHSA/CBHSQ)
Cc: R help
Subject: Re: [R] svyboxplot - library (survey)
Sent: Monday, October 01, 2012 10:07 AM
To: Muhuri, Pradip (SAMHSA/CBHSQ)
Cc: R help
Subject: Re: [R] svyboxplot - library (survey)
using a slight modification of the example shown in ?svyboxplot
# load survey library
library(survey)
# load example data
data(api)
# create an example svydesign
dstrat <- svydesign(id = ~1, strata = ~stype, weights = ~pw, data = apistrat,
fpc = ~fpc)
# set the plot window to display 1 plot x 2 plots
par(mfrow=c(1,2))
# generate two example boxplots
svyboxplot(enroll~stype,dstrat,all.outliers=TRUE)
svyboxplot(enroll~1,dstrat)
# done
# alternative: not as nice
# set the plot window to display 2 plots x 1 plot
par(mfrow=c(2,1))
# generate two example boxplots
svyboxplot(enroll~stype,dstrat,all.outliers=TRUE)
svyboxplot(enroll~1,dstrat)
# done
On Mon, Oct 1, 2012 at 9:50 AM, Muhuri, Pradip (SAMHSA/CBHSQ) <Pradip.Muhuri at samhsa.hhs.gov<mailto:Pradip.Muhuri at samhsa.hhs.gov>> wrote:
Hello,
I have used the library (survey) package for boxplots using the following code.
Could anyone please tell me why I am getting only 1 boxplot instead of 2 boxplots (1-SPD, 2-No SPD).
What changes in the following code would be required to get 2 boxplots in the same plot frame?
Thanks,
Pradip
###################################################
nhis <- svydesign (id=~psu, strat=~stratum, weights=~wt8,
data=tor, nest=TRUE)
svyboxplot (dthage~xspd2, subset (nhis, mortstat==1), col="gray80",
varwidth=TRUE, ylab="Age at Death", xlab="SPD Status: 1-SPD, 2=No SPD")
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 20857
Tel: 240-276-1070
Fax: 240-276-1260
e-mail: Pradip.Muhuri at samhsa.hhs.gov<mailto:Pradip.Muhuri at samhsa.hhs.gov>
The Center for Behavioral Health Statistics and Quality your feedback. Please click on the following link to complete a brief customer survey: http://cbhsqsurvey.samhsa.gov
vide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help at r-project.org<mailto:R-help at r-project.org> mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.