1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
library(TxDb.Hsapiens.UCSC.hg19.knownGene)
txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene
test_rename_seqlevels <- function()
{
txdb <- restoreSeqlevels(txdb)
new_seqlevels <- as.character(seq_along(seqlevels(txdb)))
seqlevels(txdb) <- new_seqlevels
checkIdentical(new_seqlevels, seqlevels(txdb))
}
test_restrict_seqlevels <- function()
{
## This should work
txdb <- restoreSeqlevels(txdb)
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5")
checkEquals(length(seqinfo(txdb)), 1)
## This should work
txdb <- restoreSeqlevels(txdb)
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5", chr6="6", chr4="4")
checkTrue(length(seqinfo(txdb)) == 3)
checkIdentical(c("5", "6", "4"), seqlevels(txdb))
checkTrue(seqlengths(txdb)[2] == min(seqlengths(txdb)))
checkTrue(seqlengths(txdb)[3] == max(seqlengths(txdb)))
## And this should NOT work
txdb <- restoreSeqlevels(txdb)
checkException(seqlevels(txdb, pruning.mode="coarse") <- c(foo="2"))
}
test_seqinfo_setter <- function()
{
txdb <- restoreSeqlevels(txdb)
new_seqinfo <- seqinfo(txdb)
seqnames(new_seqinfo) <- paste0("NEW_", seqnames(new_seqinfo))
seqinfo(txdb, new2old=seq_along(new_seqinfo)) <- new_seqinfo
checkIdentical(new_seqinfo, seqinfo(txdb))
txdb <- restoreSeqlevels(txdb)
new_seqinfo <- seqinfo(txdb)
seqlengths(new_seqinfo) <- 5 * seqlengths(new_seqinfo)
checkException(seqinfo(txdb) <- new_seqinfo)
txdb <- restoreSeqlevels(txdb)
new_seqinfo <- seqinfo(txdb)
isCircular(new_seqinfo) <- rep(TRUE, length(new_seqinfo))
checkException(seqinfo(txdb) <- new_seqinfo)
txdb <- restoreSeqlevels(txdb)
new_seqinfo <- seqinfo(txdb)
genome(new_seqinfo) <- "foo"
seqinfo(txdb) <- new_seqinfo
checkIdentical(new_seqinfo, seqinfo(txdb))
}
test_transcripts_accessor <- function()
{
txdb <- restoreSeqlevels(txdb)
txs1 <- transcripts(txdb)
seqlevels(txs1, pruning.mode="coarse") <- c(chr5="5")
## Then change seqlevels for txdb
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5")
txs2 <- transcripts(txdb)
checkIdentical(txs1, txs2)
}
test_exons_accessor <- function()
{
txdb <- restoreSeqlevels(txdb)
exs1 <- exons(txdb)
seqlevels(exs1, pruning.mode="coarse") <- c(chr5="5")
## Then change seqlevels for txdb
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5")
exs2 <- exons(txdb)
checkIdentical(exs1, exs2)
}
test_cds_accessor <- function()
{
txdb <- restoreSeqlevels(txdb)
cds1 <- cds(txdb)
seqlevels(cds1, pruning.mode="coarse") <- c(chr5="5")
## Then change seqlevels for txdb
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5")
cds2 <- cds(txdb)
checkIdentical(cds1, cds2)
}
test_promoters_accessor <- function()
{
txdb <- restoreSeqlevels(txdb)
prm1 <- promoters(txdb)
seqlevels(prm1, pruning.mode="coarse") <- c(chr5="5")
## Then change seqlevels for txdb
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5")
prm2 <- promoters(txdb)
checkIdentical(prm1, prm2)
}
test_transcriptsBy_accessors <- function()
{
## This one is a "fun" one.
## There are issues because some genes are annotated as being on
## TWO different chromosomes. Such genes are filtered for txs3,
## but NOT for txs4... Hmmmm.
txdb <- restoreSeqlevels(txdb)
txs3 <- transcriptsBy(txdb, by="gene")
seqlevels(txs3, pruning.mode="coarse") <- c(chr5="5")
## Then change seqlevels for txdb
seqlevels(txdb, pruning.mode="coarse") <- c(chr5="5")
txs4 <- transcriptsBy(txdb, by="gene")
## checkIdentical(txs3, txs4) ## TROUBLE!!
}
## What to do about this? The reason for the difference is because of order of operations. txs3 gets all the ranges and then removes any that are not kosher (this is correct), txs4 OTOH gets only ranges from chr5 (efficient!), but then fails to filter out things that have hybrid seqnames (as they were pre-filtered). I think I have to make the query less efficient to fix this, but I want to discuss it with Herve 1st to get a 2nd opinion.
|
