On 05-Dec-04 Ted Harding wrote: > [...] > For example, adopting the ritual "sigificant == P<0.05, > power = 80%", you can see a histogram of the p-values > over the conventional "significance breaks" with > > hist(pvals,breaks=c(0,0.01,0.03,0.1,0.5,0.9,0.95,0.99,1),freq=TRUE)
Sorry for the typo! That should of course (if you were looking closely) have been hist(pvals,breaks=c(0,0.01,0.05,0.1,0.5,0.9,0.95,0.99,1),freq=TRUE) Apologies to anyone who was misled by pasting in the bad version. > and you can see your probability of getting a "significant" result > as e.g. sum(pvals < 0.05)/1000 > [...] Best wishes to all, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <[EMAIL PROTECTED]> Fax-to-email: +44 (0)870 094 0861 [NB: New number!] Date: 05-Dec-04 Time: 13:54:11 ------------------------------ XFMail ------------------------------ ______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
