1 | getTargets(db, table, mir.fam)
|
db |
|
table |
|
mir.fam |
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 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (db, table, mir.fam)
{
require("data.table")
db = merge(db, mir.fam[, c(1, 4)], by.x = "miRNA", by.y = "miRBaseID")
db = as.data.table(db)
var1 <- unique(as.character(db[miRNA %in% table$mir1, Gene.Symbol]))
var2 <- unique(as.character(db[miRNA %in% table$mir2, Gene.Symbol]))
overlap <- length(intersect(var1, var2))
var1 <- length(var1)
var2 <- length(var2)
miRFamily.mir1 <- unique(as.character(db[miRNA %in% table$mir1,
miRFamily]))
miRFamily.mir2 <- unique(as.character(db[miRNA %in% table$mir2,
miRFamily]))
res <- data.frame(miRFamily.mir1 = miRFamily.mir1, mir1.targets = var1,
miRFamily.mir2 = miRFamily.mir2, mir2.targets = var2,
overlap = overlap, mir1.targets.perc = round(overlap/var1 *
100, 2), mir2.targets.perc = round(overlap/var2 *
100, 2))
return(res)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.