source("./import.R") run("..") ls() matrix.id id.pieces n n n n new.row n n n full.name c c c matrix.id c c c Q source("./import.R") source("./import.R") run("..") library(MotifDb) length(MotifDb) humans <- query(query(MotifDb, 'hsapiens'),'jaspar2016') head(humans) tbl <- read.table("./md-raw.tsv", sep = "\t", header = TRUE) head(tbl) ncbiTaxonimicCodeToBiocLinnaean(9606) ncbiTaxonimicCodeToBiocLinnaean(3702) jaspar <- query(MotifDb, 'jaspar2016') head(jaspar) head(jaspar, 40) head(tbl$ncbi.tax.code) which(tbl$matrix.id == 'MA0002.1') names(tbl) which(tbl$id == 'MA0002') which(tbl$ma.name == 'MA0002') which(tbl$ma.name == 'MA0002' & tbl$unknown == '1') which(tbl$ma.name == 'MA0002' & tbl$unknown == '2') which(tbl$ma.name == 'MA0003') run("..") blah <- load("./jaspar2016.RData") load("./jaspar2016.RData") head(tbl.md) source("./import.R") run("..") c tbl.md md id.pieces \Q Q run("..") id.pieces matrix.id head(matrix.ids) length(matrices) head(name(matrices)) head(names(matrices)) head(matrices) Q ls() head(tbl) which(tbl$ma.name == "MA0001") subset(tbl$ma.name == "MA0001") subset(tbl,ma.name == "MA0001") subset(tbl,ma.name == "MA0003") subset(tbl,ma.name == "MA0002") run("..") Q source("./import.R") run("..") traceback() nchar(NA) is.na(NA) md <- as.list(subset(tbl, ma.name == "MA0121" & unknown == "1")) md md <- as.list(subset(tbl, ma.name == "MA0102" & unknown == "1")) md ls() guessProteinIdentifierType(NA) source("./import.R") guessProteinIdentifierType(NA) moleculeName <- NA guessProteinIdentifierType(moleculeName) is.na(moleculeName) source("./import.R") is.na(moleculeName) guessProteinIdentifierType(NA) run("..") detach("package:MotifDb", unload = TRUE) library(MotifDb) length(MotifDb) q() n