Exercise 1: Solutions


load("~/public/BIOL0021/dream4-multifact4.Rdata")
#ex1
m<-cor(dataset,method="spearman")^2
diag(m)<-0
library(minet)
#ex2
conf<-validate(m,true.net)
curve<-show.pr(conf,col="green",type="l")
#ex3
cm<-clr(m)
confclr<-validate(cm,true.net)
show.pr(confclr,col="red",device=curve,type="l")
#ex4
cm[cm<0.9]<-0
cm[cm>=0.9]<-1
ig<-graph_from_adjacency_matrix(cm,mode="undirected")
#ex5
ecount(ig)
edge_density(ig)
diameter(ig)
#ex6
fastgreedy.community(ig)->comm
membership(comm)