Don't like ads? PRO users don't see any ads ;-)
Public Pastes
Guest

Untitled

By: a guest on May 9th, 2012  |  syntax: None  |  size: 0.51 KB  |  hits: 15  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. # Define the function
  2. ggd.qqplot = function(pvector, main=NULL, ...) {
  3.     o = -log10(sort(pvector,decreasing=F))
  4.     e = -log10( 1:length(o)/length(o) )
  5.     plot(e,o,pch=19,cex=1, main=main, ...,
  6.         xlab=expression(Expected~~-log[10](italic(p))),
  7.         ylab=expression(Observed~~-log[10](italic(p))),
  8.         xlim=c(0,max(e)), ylim=c(0,max(o)))
  9.     lines(e,e,col="red")
  10. }
  11.  
  12. #Generate some fake data that deviates from the null
  13. set.seed(42)
  14. pvalues=runif(10000)
  15. pvalues[sample(10000,10)]=pvalues[sample(10000,10)]/5000
create a new version of this paste RAW Paste Data