inst/testData/testData.R

library(bumphunter)
  
## Generate dummy data
set.seed(123)
dat <- dummyData()

## loessByCluster      
loessSmoothed <- loessByCluster(cluster=dat$cluster, y=dat$mat[,1],
                              bpSpan = 1000, minNum=7, minInSpan=5, maxSpan=1)
save(loessSmoothed, file="../unitTests/loessSmoothed.rda")

## runmedByCluster
runmedSmoothed <- runmedByCluster(cluster=dat$cluster, y=dat$mat[,1],  k=5, endrule="constant")
save(runmedSmoothed, file="../unitTests/runmedSmoothed.rda")

## bumphunter test 1
set.seed(123)
bumpsTest1 <- bumphunter (dat$mat, design=dat$design, chr=dat$chr, pos=dat$pos, cluster=dat$cluster,
                          coef=2, cutoff=0.28, smooth=TRUE, B=500,
                          verbose=TRUE, smoothFunction=loessByCluster)
save(bumpsTest1, file="../unitTests/bumpsTest1.rda")

## bumphunter test 2: chr unspecified
set.seed(123)
bumpsTest2 <- bumphunter (dat$mat, design=dat$design, pos=dat$pos, cluster=dat$cluster,
                          coef=2, cutoff=0.28, maxGap=500, smooth=TRUE, B=500,
                          verbose=TRUE, smoothFunction=loessByCluster)
save(bumpsTest2, file="../unitTests/bumpsTest2.rda")
rafalab/bumphunter documentation built on March 20, 2024, 6:22 a.m.