Re: [R] binom.test [Broadcast]

From: Liaw, Andy <andy_liaw_at_merck.com>
Date: Sat 21 Oct 2006 - 08:25:38 GMT


To quote one of the previous answers you've got: The formula you're using is the TV. The one binom.test() uses is the ballpark. Take your pick.

Andy

From: Ethan Johnsons
>
> Thank you for the info. It helps.
>
> After all, it would be:
>
> > 0.1304348-1.96*(sqrt((0.1304348*(1-0.1304348))/46))
> [1] 0.03310968
> > 0.1304348+1.96*(sqrt((0.1304348*(1-0.1304348))/46))
> [1] 0.2277599
>
> Does R have a function for the calculation above?
>
> ej
>
>
> On 10/20/06, Francisco J. Zagmutt <gerifalte28@hotmail.com> wrote:
> > Ethan,
> >
> > You need to explain why you think this is "not the right
> function to
> > use". R is doing exactly what you are asking it to do. Now
> is up to
> > you to choose the methodology you feel is correct.
> > For a good discussion on your particular issue I recommend you the
> > following
> > reference:
> >
> > A. Agresti and B. A. Coull, "Approximate is better than "exact" for
> > interval estimation of binomial proportions," The American
> Statistician, vol. 52, no.
> > 2, pp. 119-126, 1998.
> >
> > Once you figure out the "right function to use" see if the
> function is
> > available in R. If not readily available, and if after
> searching through
> > R's documentation and the forum archives you still can't
> find a way to
> > perform the calculation, then is time to get back to this forum.
> >
> > Regards,
> >
> > Francisco
> >
> >
> > Dr. Francisco J. Zagmutt
> > College of Veterinary Medicine and Biomedical Sciences
> Colorado State
> > University
> >
> >
> >
> >
> > >From: "Ethan Johnsons" <ethan.johnsons@gmail.com>
> > >To: r-help@stat.math.ethz.ch
> > >Subject: [R] binom.test
> > >Date: Fri, 20 Oct 2006 17:18:02 -0400
> > >
> > >A quick question, please.
> > >
> > >46 e coli lab samples are tested, 6 of them returned positive.
> > >
> > >So, the best point estimate for p is 6/46 = 0.1304348.
> > >
> > >For a 95% CI for p, I thought binom.test would give me the correct
> > >result, but it seems it is not the right function to use. What is
> > >the R function for this?
> > >
>
> ______________________________________________
> R-help@stat.math.ethz.ch 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.
>
>
>



Notice: This e-mail message, together with any attachments,...{{dropped}}

R-help@stat.math.ethz.ch 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. Received on Sat Oct 21 18:32:42 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 21 Oct 2006 - 09:30:11 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.