inst/eval/GSE11976/01.createDataSets.R

filenames <- sprintf("%s,b=%s.xdr", simName, 1:B)
regA <- data.frame(region = c("(1,1)", "(0,2)", "(1,2)", "(0,1)", "(0,3)", "(0,0)", "(2,2)"),
                   freq = c(0.90/4, 0.90/4, 0.90/4, 0.90/4, 0.10/3, 0.10/3, 0.10/3))
for (bb in 1:B) {
    filename <- filenames[bb]
    print(filename)
    pathname <- file.path(spath, filename)
    if (!file.exists(pathname) || simForce) {
        simTag100 <- sprintf("ROC,n=%s,K=%s,regSize=%s,minL=%s,pct=100", len, K, regSize, minL)
        simName100 <- sprintf("%s,%s", dataSet, simTag100)
        filename100 <- sprintf("%s,b=%s.xdr", simName100, 1:B)[bb]
        spath100 <- file.path(simPath, simName100)
        pathname100 <- file.path(spath100, filename100)
        if(file.exists(pathname100)){
            sim100 <- loadObject(pathname100)
            bkp100 <- sim100$bkp
            regions100 <- sim100$region
        }else{
            sim100 <- NULL
            bkp100 <- NULL
            regions100 <- NULL
        }
        sim <- getCopyNumberDataByResampling(len, K, bkp=bkp100,
                                             regions=regions100, minLength=minL,
                                             regionSize=regSize, regData=dat,
                                             connex=TRUE, regAnnot=regA)
        saveObject(sim, file=pathname)
    }
}

Try the jointseg package in your browser

Any scripts or data that you put into this service are public.

jointseg documentation built on May 2, 2019, 6:10 a.m.