[R] Help for Power analysis
Ass.Prof. Nikom Thanomsieng
nikom at kku.ac.th
Sun Dec 9 18:29:12 CET 2001
Dear colleague,
I not sure this R code is correctly ? I would to show
the number of Sample Size at Sample Size Axis that line
draw from Power Axis (80%) from R code.
How I show this and select the most appropriate of
this power (.79955687 - 80983575).
Thank for your help and answer.
Best Regards,
Nikom Thanomsieng,
Email: nikom at kku.ac.th
....
#Power analysis: Sample size for Chi-Square 2x2, RxC: R software
#Concept from SAS program to calculate power of ANOVA F-test
#http://www2.tltc.ttu.edu/Westfall/images/5347/power_analysis_of_anova_f_test.htm
# and Cohen,J (1977). Statistical Power Analysis for Behavioral
Sciences.
#New York: Academic Press.
#Asst. Prof. Nikom Thanomsieng. 29/09/2000
#Department of Biostatistics & Demography. Faculty of Public Health.
#Khon Kaen University. Thailand.
#Email: nikom at kku.ac.th
#Modify data value of the first two line
x1 <- c(6,9)
x2 <- c(6,6)
nc <-cbind(x1)
nr <-rbind(x2)
data1 <- rbind(x1,x2)
chi2<- chisq.test(data1,correct=F)$statistic
Ntotal<-sum(data1)
df<- ncol(nc-1)*nrow(nr-1)
ifelse(df==1,w<- sqrt(chi2/Ntotal), w<-sqrt(chi2/(chi2+Ntotal)))
Ntotal1<-900 #change this if power not enough
alpha <-0.05 #change this for One tailed =0.05
ncp<-0
chicrit<-NULL
power<-NULL
n<-NULL
samplesize<-NULL
for (i in 1:Ntotal1){
ncp[i] <- w^2 * i
chicrit<-qchisq(1-alpha,df)
power[i] <- 1-(pchisq(chicrit , df, ncp[i]))
n[i]<-i
samplesize<-cbind(n, ncp,power) }
samplesize
plot(n,power,type="l",col="red", lwd=1,
panel.first = grid(10,10),
main="Power as a function of Sample Size",
xlab="Sample Size",
ylab="Power" )
segments(785, .8, 785, 0, col ="pink")
segments(785, .8, 0, 0.8, col ="pink")
mtext("Chi-Square", side = 3, line = 0.35,
outer = FALSE, at = mean(par("usr")[1:2]), cex = 1, font = 4)
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list