Learn to search the Archives (or at least one of the forms of the
Archives)
http://search.r-project.org/cgi-bin/namazu.cgi?query=confidence+interval+logit+proportion&max=100&result=normal&sort=score&idxname=functions&idxname=Rhelp08&idxname=Rhelp10&idxname=Rhelp02
(That is constructed with an URL I use for the purpose that brings up
the RSiteSearch page with different parameters than are currently in
the default implementation. This is an R function that does the same:
> rhelpSearch
function(string,
restrict = c("Rhelp10", "Rhelp08", "Rhelp02",
"functions" ),
matchesPerPage = 100, ...)
RSiteSearch(string=string, restrict = restrict,
matchesPerPage = matchesPerPage, ...)
--
David
On Jun 20, 2011, at 8:50 AM, Muhuri, Pradip (SAMHSA/CBHSQ) wrote:
Dear Patrick,
I do agree with you that it is a very simple problem. Actually, I
do have the following SAS program written to compute the
asymmetrical confidence interval.
As a new user of R, I just wanted to see the corresponding codes in
R if they already exist.
Thanks,
Pradip
*SAS program begins here;
/********************************************************
MEAN = prevalence rate
PLOWER = lower 95% confidence limit for the rate
PPER = upper 95% confidence limit for the rate
TLOWER = lower 95% confidence limit for the total
TUPPER = upper 95% confidence limit for the total
Calculate the 95% CI FOR PREVALENCE RATES AND TOTALS
********************************************************/
IF MEAN=0 OR MEAN=1 THEN DO;
L=.;
NUMBER=.;
A=.; B=.;
PLOWER=.; PUPPER=.; TLOWER=.; TUPPER=.;
END;
ELSE DO;
L=LOG(MEAN/(1-MEAN));
NUMBER=SEMEAN/(MEAN*(1-MEAN));
A=L-1.96*NUMBER;
B=L+1.96*NUMBER;
PLOWER=1/(1+EXP(-A)); PUPPER=1/(1+EXP(-B));
TLOWER=WSUM*PLOWER; TUPPER=WSUM*PUPPER;
END;
RUN;
*SAS program ends here:
Pradip K. Muhuri, PhD
Statistician
Substance Abuse & Mental Health Services Administration
The Center for Behavioral Health Statistics and Quality
Division of Population Surveys
1 Choke Cherry Road, Room 7-1023
Rockville, MD 20857
Tel: 240-276-1070
Fax: 240-276-1260
e-mail: pradip.muh...@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
-----Original Message-----
From: Patrick Connolly [mailto:p_conno...@slingshot.co.nz]
Sent: Sunday, June 19, 2011 1:57 AM
To: Muhuri, Pradip (SAMHSA/CBHSQ)
Cc: r-help@r-project.org; 'tlum...@u.washington.edu'
Subject: Re: [R] Asymetrical Confidence Interval
On Thu, 16-Jun-2011 at 04:43PM -0400, Muhuri, Pradip (SAMHSA/CBHSQ)
wrote:
|>
|> Dear List,
|>
|> I wanted to calculate the asymmetrical confidence interval based on
|> the sample statistic and standard error that available from the
|> published report (complex survey-based).
|> The calculation details can be seen from pages 17-18 of the
|> document at the following link:
|> http://www.oas.samhsa.gov/nsduh/2k5MRB/2k5statInference.pdf.
|>
|> Could someone tell me whether R has any function included in it
|> "survey" or other contributed package of R.
There might be one in a package somewhere, but it's so trivial to make
your own function by using the information you already have.
This is sounding suspiciously like a homework question.
--
~
.~
.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.
___ Patrick Connolly
{~._.~} Great minds discuss ideas
_( Y )_ Average minds discuss events
(:_~*~_:) Small minds discuss people
(_)-(_) ..... Eleanor Roosevelt
~
.~
.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.
______________________________________________
R-help@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.
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@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.
