library(cisPath)
# Generate the identifier mapping file
input <- system.file("extdata", "uniprot_sprot_human10.dat", package="cisPath")
mappingFile <- file.path(tempdir(), "mappingFile.txt")
getMappingFile(input, output=mappingFile, taxonId="9606")
# Format the file downloaded from STRING database
output <- file.path(tempdir(), "PINA2PPI.txt")
fileFromPINA2 <- system.file("extdata", "Homo_sapiens_PINA2.sif", package="cisPath")
formatSIFfile(fileFromPINA2, mappingFile, output)
source("http://bioconductor.org/biocLite.R")
biocLite("R.utils")
library(R.utils)
outputDir <- file.path(getwd(), "cisPath_test")
dir.create(outputDir, showWarnings=FALSE, recursive=TRUE)
# Generate the identifier mapping file
fileFromUniProt <- file.path(outputDir, "uniprot_sprot_human.dat")
mappingFile <- file.path(outputDir, "mappingFile.txt")
getMappingFile(fileFromUniProt, output=mappingFile)
# Download PINA2 PPI (SIF format) for humans only (~2.8M)
destfile <- file.path(outputDir, "Homo_sapiens.sif")
cat("Downloading...\n")
download.file("http://cbg.garvan.unsw.edu.au/pina/download/Homo%20sapiens-20140521.sif", destfile)
# Format PINA2 PPI
fileFromPINA2 <- file.path(outputDir, "Homo_sapiens.sif")
PINA2PPI <- file.path(outputDir, "PINA2PPI.txt")
formatSIFfile(fileFromPINA2, mappingFile, output=PINA2PPI)Run the code above in your browser using DataLab