inst/pdBuilders/runPDIBSnp5.R

In benilton/pdInfoBuilder: Platform Design Information Package Builder

library("pdInfoBuilder")

targ = "../pd.genomewidesnp.5"

if (file.exists(targ)) stop(paste("cannot have folder", targ,
   "present if you want to build there"))
parms_pref = "../../parms_store/pd.genomewidesnp.5"  # should be in parms_store
anno_pref = "../../anno_store/pd.genomewidesnp.5"
pat = function(x) paste(parms_pref, x, sep="/")
pat2 = function(x) paste(anno_pref, x, sep="/")

# next 4 lines gunzip contents of anno_store, if needed
curd = getwd()
setwd(anno_pref)
try(system( paste(c("gunzip", dir(patt="gz$")), collapse=" ")))
setwd(curd)

# REMEMBER TO GUNZIP for NOW ... eventually use gzfiles?
cdfFile <- pat2("GenomeWideSNP_5.cdf")
csvAnno <- pat2("GenomeWideSNP_5.na24.annot.csv")
csvSeq <- pat2("GenomeWideSNP_5.probe_tab")
#../../anno_store/pd.genomewidesnp.6/GenomeWideSNP_5.probe_tab
csvAnnoCnv <- pat2("GenomeWideSNP_5.cn.na24.annot.csv")
csvSeqCnv <- pat2("GenomeWideSNP_5.CN_probe_tab")
spline <- pat("pd.genomewidesnp.5.spline.params.rda")
refd <- pat("pd.genomewidesnp.5Ref.rda")
crlmmInf <- pat("pd.genomewidesnp.5CrlmmInfo.rda")

pkg <- new("AffySNPCNVPDInfoPkgSeed",
           version="0.3.5",
           author="Vince Carey", email="stvjc@channing.harvard.edu",
           biocViews="AnnotationData",
           genomebuild="NCBI Build 36",
           cdfFile=cdfFile, csvAnnoFile=csvAnno, csvSeqFile=csvSeq,
           csvAnnoFileCnv=csvAnnoCnv, csvSeqFileCnv=csvSeqCnv,
           splineParamFile=spline,
           crlmmInfoFile=crlmmInf, referenceDistFile=refd)

makePdInfoPackage(pkg, destDir="..")

# next 4 lines gzip contents of anno_store
curd = getwd()
setwd(anno_pref)
try(system( paste(c("gzip", dir()), collapse=" ")))
setwd(curd)