[R] Unbalanced Design Power Analysis
Peter Alspach
Peter.Alspach at plantandfood.co.nz
Wed Jun 13 05:38:43 CEST 2012
Tena koe Jacob
Possibly your best option is to simulate.
HTH ....
Peter Alspach
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Jacob Warren (RIT Student)
Sent: Wednesday, 13 June 2012 4:40 a.m.
To: r-help at r-project.org
Subject: [R] Unbalanced Design Power Analysis
I have an unbalanced design I would like to run a power analysis on.
What I have been able to find has pointed me to using the pwr.f2.test
function as described below. My problem is that I don't know how to
appropriately define the numerator and demoninator df.
If someone can help here is some more info about my design.
It is an unbalanced 2^3 x 3 design where the factor with 3 levels is a
random effect.
So say there are factors A, B, C with 2 levels and D with 3 levels
A will have 77 observations at level -1 and 19 at 1
B will have 67 at -1 and 29 at 1
C will have 57 at -1 and 39 at 1
D will have 32 at -1, 0, and 1
General Linear Model
Function: pwr.f2.test
Arguments: u: Numerator degrees of freedom
v: Denominator degrees of freedom
f2: Effect Size
sig.level: Significance level
power: Power of test
______________________________________________
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.
The contents of this e-mail are confidential and may be subject to legal privilege.
If you are not the intended recipient you must not use, disseminate, distribute or
reproduce all or any part of this e-mail or attachments. If you have received this
e-mail in error, please notify the sender and delete all material pertaining to this
e-mail. Any opinion or views expressed in this e-mail are those of the individual
sender and may not represent those of The New Zealand Institute for Plant and
Food Research Limited.
More information about the R-help
mailing list