## if(CHR.X) {
## sample.index <- which(object$gender==2)
## batches <- split(sample.index, as.character(batch(object))[sample.index])
## } else {
## batches <- split(seq_along(batch(object)), as.character(batch(object)))
## }
git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/crlmm@52814 bc3139a8-67e5-0310-9ffc-ced21a209358
Rob Scharp authored on 16/02/2011 15:57:01 | ... | ... |
@@ -1430,12 +1430,13 @@ summarizeSnps <- function(strata, |
| 1430 | 1430 |
} |
| 1431 | 1431 |
if(verbose) message("Probe stratum ", strata, " of ", length(index.list))
|
| 1432 | 1432 |
index <- index.list[[strata]] |
| 1433 |
- if(CHR.X) {
|
|
| 1434 |
- sample.index <- which(object$gender==2) |
|
| 1435 |
- batches <- split(sample.index, as.character(batch(object))[sample.index]) |
|
| 1436 |
- } else {
|
|
| 1437 |
- batches <- split(seq_along(batch(object)), as.character(batch(object))) |
|
| 1438 |
- } |
|
| 1433 |
+## if(CHR.X) {
|
|
| 1434 |
+## sample.index <- which(object$gender==2) |
|
| 1435 |
+## batches <- split(sample.index, as.character(batch(object))[sample.index]) |
|
| 1436 |
+## } else {
|
|
| 1437 |
+## batches <- split(seq_along(batch(object)), as.character(batch(object))) |
|
| 1438 |
+## } |
|
| 1439 |
+ batches <- split(seq_along(batch(object)), as.character(batch(object))) |
|
| 1439 | 1440 |
batchnames <- batchNames(object) |
| 1440 | 1441 |
nr <- length(index) |
| 1441 | 1442 |
nc <- length(batchnames) |