################################################### ### chunk number 1: ################################################### require(GLAD) ################################################### ### chunk number 2: ################################################### data(snijders) profileCGH <- as.profileCGH(gm13330) res <- glad(profileCGH, mediancenter=FALSE, smoothfunc="lawsglad", bandwidth=10, round=1.5, model="Gaussian", lkern="Exponential", qlambda=0.999, base=FALSE, lambdabreak=8, lambdacluster=8, lambdaclusterGen=40, type="tricubic", param=c(d=6), alpha=0.001, msize=5, method="centroid", nmax=8, verbose=FALSE) ################################################### ### chunk number 3: ################################################### plotProfile(res, unit=3, Bkp=TRUE, labels=FALSE, plotband=FALSE, Smoothing="Smoothing") ################################################### ### chunk number 4: ################################################### data(veltman) profileCGH <- as.profileCGH(P9) profileCGH <- daglad(profileCGH, mediancenter=FALSE, normalrefcenter=FALSE, genomestep=FALSE, smoothfunc="lawsglad", lkern="Exponential", model="Gaussian", qlambda=0.999, bandwidth=10, base=FALSE, round=1.5, lambdabreak=8, lambdaclusterGen=40, param=c(d=6), alpha=0.001, msize=5, method="centroid", nmin=1, nmax=8, amplicon=1, deletion=-5, deltaN=0.2, forceGL=c(-0.3,0.3), nbsigma=3, MinBkpWeight=0.35, CheckBkpPos=TRUE) ################################################### ### chunk number 5: ################################################### plotProfile(profileCGH, Smoothing="Smoothing", Bkp=TRUE, plotband=FALSE) abline(h=c(-0.3,-0.2,0.2,0.3),col=c("green","black","black","red")) axis(2,at=c(-0.3,-0.2,0.2,0.3), labels=c("forceGL[1]","-deltaN","+deltaN","forceGL[2]"), las=2) ################################################### ### chunk number 6: ################################################### data(veltman) profileCGH <- as.profileCGH(P9) profileCGH <- daglad(profileCGH, mediancenter=FALSE, normalrefcenter=FALSE, genomestep=FALSE, smoothfunc="lawsglad", lkern="Exponential", model="Gaussian", qlambda=0.999, bandwidth=10, base=FALSE, round=1.5, lambdabreak=8, lambdaclusterGen=40, param=c(d=6), alpha=0.001, msize=5, method="centroid", nmin=1, nmax=8, amplicon=1, deletion=-5, deltaN=0.10, forceGL=c(-0.15,0.15), nbsigma=3, MinBkpWeight=0.35, CheckBkpPos=TRUE) ################################################### ### chunk number 7: ################################################### plotProfile(profileCGH, Smoothing="Smoothing", Bkp=TRUE, plotband=FALSE) abline(h=c(-0.15,-0.1,0.1,0.15),col=c("green","black","black","red")) axis(2,at=c(-0.15,-0.1,0.1,0.15), labels=c("forceGL[1]","-deltaN","+deltaN","forceGL[2]"), las=2) ################################################### ### chunk number 8: ################################################### data(arrayCGH) # object of class arrayCGH array <- list(arrayValues=array2, arrayDesign=c(4,4,21,22)) class(array) <- "arrayCGH" ################################################### ### chunk number 9: ################################################### arrayPlot(array,"Log2Rat", bar="none") ################################################### ### chunk number 10: ################################################### arrayPersp(array,"Log2Rat", box=FALSE, theta=110, phi=40, bar=FALSE) ################################################### ### chunk number 11: ################################################### data(snijders) profileCGH <- as.profileCGH(gm13330) res <- glad(profileCGH, mediancenter=FALSE, smoothfunc="lawsglad", bandwidth=10, round=2, model="Gaussian", lkern="Exponential", qlambda=0.999, base=FALSE, lambdabreak=8, lambdacluster=8, lambdaclusterGen=40, type="tricubic", param=c(d=6), alpha=0.001, msize=5, method="centroid", nmax=8, verbose=FALSE) ################################################### ### chunk number 12: ################################################### plotProfile(res, unit=3, Bkp=TRUE, labels=FALSE, Smoothing="Smoothing", plotband=FALSE) ################################################### ### chunk number 13: ################################################### plotProfile(res, unit=3, Bkp=TRUE, labels=FALSE, Smoothing="Smoothing") ################################################### ### chunk number 14: ################################################### text <- list(x=c(90000,200000),y=c(0.15,0.3),labels=c("NORMAL","GAIN"), cex=2) plotProfile(res, unit=3, Bkp=TRUE, labels=TRUE, Chromosome=1, Smoothing="Smoothing", plotband=FALSE, text=text) ################################################### ### chunk number 15: ################################################### text <- list(x=c(90000,200000),y=c(0.15,0.3),labels=c("NORMAL","GAIN"), cex=2) plotProfile(res, unit=3, Bkp=TRUE, labels=TRUE, Chromosome=1, Smoothing="Smoothing", text=text, main="Chromosome 1")