Re: [R] Power curves

From: Chuck Cleland <ccleland_at_optonline.net>
Date: Mon, 14 Apr 2008 06:35:15 -0400

On 4/14/2008 6:21 AM, Louisa Hay wrote:
>
> Hi,
> I am trying to create a power curve to show how the power of a t-test varies depending on the mean. Any ideas how I should go about this?
> Louisa

   Something like this maybe?

mydiffs <- seq(.05,1,.05)

mypower <- vector("numeric", 20)

for(i in 1:20){

   mypower[i] <- power.t.test(n=100, sd=1, delta=mydiffs[i])$power }

plot(mymeans, mypower, xlab="delta", ylab="Power", type="b")

   See ?power.t.test

> _________________________________________________________________
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.

-- 
Chuck Cleland, Ph.D.
NDRI, Inc. (www.ndri.org)
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212) 845-4495 (Tu, Th)
tel: (732) 512-0171 (M, W, F)
fax: (917) 438-0894

______________________________________________
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.
Received on Mon 14 Apr 2008 - 10:37:24 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 14 Apr 2008 - 13:30:28 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.

list of date sections of archive