Nothing
#This file belongs to
#netmes: NETMES, <http://www.bioconductor.org/packages/release/Software.html>
#This R package allows local network based assesments for
##inferring regulatory networks from expression data.
## Copyright (C) August 2009 Frank Emmert-Streib and Gokmen Altay
##<f.emmert-streib@qub.ac.uk>
## This program is free software; you can redistribute it and/or
## modify it under the terms of the GNU LESSER GENERAL PUBLIC LICENSE
## as published by the Free Software Foundation; either version 3
## of the License, or any later version.
##
## This program is distributed WITHOUT ANY WARRANTY;
## You can get a copy of the GNU LESSER GENERAL PUBLIC LICENSE
## from
## http://www.r-project.org/Licenses/LGPL-3
readSifNet <- function(infilenet, infilepath){
file.adr <- paste(infilepath, "1/mim.exp", sep = "");
cc <-read.table(file.adr)
G <- dim(cc)[1] # number of genes is assigned inside the local function
dat <- matrix(scan(file = infilenet, what = "character", skip = 0), ncol = 3, byrow = TRUE)
N <- dim(dat)[1]
names <- union(dat[,1], dat[,3])
net <- matrix(0, nrow = G, ncol = G)
for(i in 1:N){
aux <- dat[i,]
ind1 <- which(names == aux[1])
ind2 <- which(names == aux[3])
if(aux[2] == "ac"){
net[ind1, ind2] <- 1
}
else{
net[ind1, ind2] <- -1
}
}
diag(net) <- c(1:G)*0 # delete self-loops if existing
res <- new.env()
assign("net", net, envir=res)
assign("names", names, envir=res)
res
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.