inst/eval/GSE29172/02e.PSCBS.R

filenames <- sprintf("%s,b=%s.xdr", simName, 1:B)
for (bb in 1:B) {
  filename <- filenames[bb]
  print(filename)
  pathname <- file.path(spath, filename)
  sim <- loadObject(pathname)
  if (!is.na(normFrac)) {
    dat <- setNormalContamination(sim$profile, normFrac)
  } else {
    dat <- sim$profile
  }
  ## drop outliers
  CNA.object <- CNA(dat$c,rep(1,len),1:len)
  smoothed.CNA.obj <- smooth.CNA(CNA.object)
  dat$c <- smoothed.CNA.obj$Sample.1
  
  stat <- c("c", "b", "genotype")
  dpStat <- c("c", "d")
  methTag <- "PSCBS"
  filename <- sprintf("%s,b=%s,%s.xdr", simNameNF, bb, methTag)
  pathname <- file.path(bpath, filename)
  if (!file.exists(pathname) || segForce) {
    res <- PSSeg(dat, method="PSCBS", stat=stat, profile=FALSE, verbose=TRUE)
    print(res$prof[, "time"])
    saveObject(res, 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, 5:20 p.m.