iterating_seqVsInsitu {cellOrigins} | R Documentation |
The same functionality as seqVsInsitu
but computationally less expensive if combinations of anatomical terms are tested.
The number of term combinations to test increases rapidly in seqVsInsitu
. For example with 350 anatomical terms there are 61425 combinations of 2 terms and 7207200 combinations of 3 terms. This makes the exhaustive search of seqVsInsitu
costly with depth>2.
iterating_seqVsInsitu
reduces the computational cost by initially testing the combinations of only a few terms. Then in each iteration the cardinality of the combinations is increased by one, but only the top anatomical terms of the previous iteration are used to reduce the number of tested combinations.
iterating_seqVsInsitu(seq_signature, upto_depth, use_topN = 50,
start_depth = 2, insitu = cellOrigins::BDGP_insitu_dmel_embryo,
insitu_discovery_function = discovery.log, saturate = 500,
prior = prior.temporal_proximity_is_good)
seq_signature |
A named vector containing FPKM RNAseq data. Each element name must correspond to the names used in the |
upto_depth |
Number of terms to combine in the final iteration. |
use_topN |
How many of the top results from the previous iteration to use to find the terms for the current iteration. |
start_depth |
Number of terms to combine in the first iteration. All combinations of all terms are tested at this step. |
insitu |
Matrix with RNA in situ hybridisation data. Rows are transcript names (queried by probes: same names as used for |
insitu_discovery_function |
A function that converts FPKM values to the probability of discovery by RNA in situ hybridisation. Values must be ]0..1[, 0 and 1 are not permitted.
Defaults to |
saturate |
Will be passed on to the |
prior |
A function that evaluates to the log2 prior probability of each anatomic term or combination of terms.
Defaults to |
Returns a named list that contains a matrix for each iteration like those produced by seqVsInsitu
.
## Not run:
fpath <- system.file("extdata", "vncMedianCoverage.tsv", package="cellOrigins")
vncExpression <- read.delim(file = fpath, header=FALSE, as.is=TRUE)
expression <- vncExpression$V2
names(expression) <- vncExpression$V1
oracleResponse <- iterating_seqVsInsitu(expression, 3)
head(oracleResponse[[1]])
head(oracleResponse[[2]])
diagnosticPlots(oracleResponse)
## End(Not run)