inst/doc/dks.R

In dks: The double Kolmogorov-Smirnov package for evaluating multiple testing procedures.

### R code from vignette source 'dks.Rnw'

###################################################
### code chunk number 1: dks.Rnw:73-77
###################################################
library(dks)
library(cubature)
data(dksdata)
dim(P)


###################################################
### code chunk number 2: dks1
###################################################
dks1 <- dks(P,plot=TRUE)


###################################################
### code chunk number 3: pprob
###################################################
 ## Calculate the posterior distribution
 delta <- 0.1
  dist1 <- pprob.dist(P[,1])

  ## Plot the posterior distribution
  alpha <- seq(0.1,10,by=delta)
  beta <- seq(0.1,10,by=delta)
  pprobImage = image(log10(alpha),log10(beta),dist1,xaxt="n",yaxt="n",xlab="Alpha",ylab="Beta")
  axis(1,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
  axis(2,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
  points(0,0,col="blue",cex=1,pch=19)	


###################################################
### code chunk number 4: credSet
###################################################
 ## Calculate a 80% credible set
  cred1 <- cred.set(dist1,delta=0.1, level=0.80)

  ## Plot the posterior and the credible set
  
  alpha <- seq(0.1,10,by=delta)
  beta <- seq(0.1,10,by=delta)

  credImage=image(log10(alpha),log10(beta),cred1$cred,xaxt="n",yaxt="n",xlab="Alpha",ylab="Beta")
  axis(1,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
  axis(2,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
  points(0,0,col="blue",cex=1,pch=19)