[R] binom.test [Broadcast]
Liaw, Andy
andy_liaw at merck.com
Sat Oct 21 10:25:38 CEST 2006
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 at 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 at gmail.com>
> > >To: r-help at 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 at 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}}
More information about the R-help
mailing list