git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/crlmm@52846 bc3139a8-67e5-0310-9ffc-ced21a209358
... | ... |
@@ -109,9 +109,6 @@ assign these identifiers to the variable \Robject{batch}. |
109 | 109 |
celFiles <- list.celfiles(pathToCels, full.names=TRUE, pattern=".CEL") |
110 | 110 |
celFiles <- celFiles[substr(basename(celFiles), 13, 13) %in% c("C", "Y")] |
111 | 111 |
batch <- as.factor(substr(basename(celFiles), 13, 13)) |
112 |
-sample.index <- c(sample(which(batch == "C"), 15), sample(which(batch == "Y"), 11)) |
|
113 |
-celFiles <- celFiles[sample.index] |
|
114 |
-batch <- batch[sample.index] |
|
115 | 112 |
@ |
116 | 113 |
|
117 | 114 |
The preprocessing steps for copy number estimation includes quantile |
... | ... |
@@ -193,10 +190,7 @@ The \Rfunction{crlmmCopynumber} performs the following steps: |
193 | 190 |
|
194 | 191 |
<<LDS_copynumber>>= |
195 | 192 |
GT.CONF.THR <- 0.8 |
196 |
-cnSet <- checkExists("cnSet", .path=outdir, .FUN=crlmmCopynumber, object=gtSet, GT.CONF.THR=GT.CONF.THR, |
|
197 |
- MIN.SAMPLES=12) |
|
198 |
-cnSet <- checkExists("cnSet", .path=outdir, .FUN=crlmmCopynumber, object=gtSet, GT.CONF.THR=GT.CONF.THR, |
|
199 |
- MIN.SAMPLES=11) |
|
193 |
+cnSet <- checkExists("cnSet", .path=outdir, .FUN=crlmmCopynumber, object=gtSet, GT.CONF.THR=GT.CONF.THR) |
|
200 | 194 |
@ |
201 | 195 |
|
202 | 196 |
In an effort to reduce I/O, the \Rpackage{crlmmCopynumber} function no |