diff --git a/analysis/quantif_20131105/.RData b/analysis/quantif_20131105/.RData new file mode 100644 index 0000000..01c9cf4 Binary files /dev/null and b/analysis/quantif_20131105/.RData differ diff --git a/analysis/quantif_20131105/120518_07-degreeDist131109-213811.pdf b/analysis/quantif_20131105/120518_07-degreeDist131109-213811.pdf new file mode 100644 index 0000000..c864a16 Binary files /dev/null and b/analysis/quantif_20131105/120518_07-degreeDist131109-213811.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-20-01.pdf b/analysis/quantif_20131105/2013-11-05-13-20-01.pdf new file mode 100644 index 0000000..12a0609 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-20-01.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-22-06.pdf b/analysis/quantif_20131105/2013-11-05-13-22-06.pdf new file mode 100644 index 0000000..4d5479b Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-22-06.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-23-11.pdf b/analysis/quantif_20131105/2013-11-05-13-23-11.pdf new file mode 100644 index 0000000..73132d1 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-23-11.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-25-38.pdf b/analysis/quantif_20131105/2013-11-05-13-25-38.pdf new file mode 100644 index 0000000..6dca705 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-25-38.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-30-47.pdf b/analysis/quantif_20131105/2013-11-05-13-30-47.pdf new file mode 100644 index 0000000..a1c9d30 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-30-47.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-34-28.pdf b/analysis/quantif_20131105/2013-11-05-13-34-28.pdf new file mode 100644 index 0000000..9387bc9 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-34-28.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-38-54.pdf b/analysis/quantif_20131105/2013-11-05-13-38-54.pdf new file mode 100644 index 0000000..ec6df61 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-38-54.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-54-21.pdf b/analysis/quantif_20131105/2013-11-05-13-54-21.pdf new file mode 100644 index 0000000..42544f4 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-54-21.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-13-56-43.pdf b/analysis/quantif_20131105/2013-11-05-13-56-43.pdf new file mode 100644 index 0000000..41d7ec9 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-13-56-43.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-14-00-43.pdf b/analysis/quantif_20131105/2013-11-05-14-00-43.pdf new file mode 100644 index 0000000..90881c4 Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-14-00-43.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-14-29-10.pdf b/analysis/quantif_20131105/2013-11-05-14-29-10.pdf new file mode 100644 index 0000000..c15ee2b Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-14-29-10.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05-14-34-36.pdf b/analysis/quantif_20131105/2013-11-05-14-34-36.pdf new file mode 100644 index 0000000..13cdf8e Binary files /dev/null and b/analysis/quantif_20131105/2013-11-05-14-34-36.pdf differ diff --git a/analysis/quantif_20131105/2013-11-05_wholeBrain_graphTest.R b/analysis/quantif_20131105/2013-11-05_wholeBrain_graphTest.R new file mode 100644 index 0000000..2f48100 --- /dev/null +++ b/analysis/quantif_20131105/2013-11-05_wholeBrain_graphTest.R @@ -0,0 +1,97 @@ +require(igraph) + +edgelist<-read.delim('/Users/ackman/Documents/MATLAB/dCorrpairs-20131105-125023.txt') + +edgelist2<-subset(edgelist,rvalue > 0.1) + + +g <- graph.data.frame(edgelist2, directed=FALSE) +#V(g)$label <- V(g)$name #gets the actual vertice labels +plot(g,layout=layout.lgl,vertex.size=5) +title('layoutLgl, vertex size5') + +quartz(); +plot(g, layout=layout.fruchterman.reingold, vertex.size=3) +title('layoutFR, vertex size3') + +quartz(); +com <- spinglass.community(g, spins=5) #finds communities +V(g)$color <- com$membership+1 +plot(g, layout=layout.fruchterman.reingold) +title('spinglass, layoutFR, default vertex size') + + +walktrapCom<-walktrap.community(g) +V(g)$color <- walktrapCom$membership+1 + +quartz(); +palette(rainbow(max(V(g)$color),alpha=0.5)) +plot(g, layout=layout.fruchterman.reingold) +palette("default") +title('walktrap, layoutFR, default vertex size,alpha0.5') + + +quartz(); +palette(rainbow(max(V(g)$color),alpha=0.5)) +plot(g, layout=layout.lgl) +palette("default") +title('walktrap, layoutLgl, default vertex size,alpha0.5') + + +quartz(); +palette(rainbow(max(V(g)$color),alpha=0.5)) +plot(g, layout=layout.kamada.kawai) +palette("default") +title('walktrap, layout.kamada.kawai, default vertex size,alpha0.5') + +g <- graph.data.frame(edgelist2, directed=FALSE) +fastgreedyCom<-fastgreedy.community(g) +V(g)$color <- fastgreedyCom$membership+1 +quartz(); +palette(rainbow(max(V(g)$color),alpha=0.5)) +plot(g, layout=layout.kamada.kawai) +palette("default") +title('fastgreedy, layout.kamada.kawai, default vertex size,alpha0.5') + + + + +edgelist2<-subset(edgelist,rvalue > 0.1) +g <- graph.data.frame(edgelist2, directed=FALSE) +fastgreedyCom<-fastgreedy.community(g) +V(g)$color <- fastgreedyCom$membership+1 +quartz(); +palette(rainbow(max(V(g)$color),alpha=0.5)) +plot(g, layout=layout.fruchterman.reingold) +palette("default") +title('fastgreedy, layout.fruchterman.reingold, default vertex size,alpha0.5') + + + + +# Edgelist > 0.2 + +edgelist2<-subset(edgelist,rvalue > 0.2) +g <- graph.data.frame(edgelist2, directed=FALSE) +fastgreedyCom<-fastgreedy.community(g) +V(g)$color <- fastgreedyCom$membership+1 +quartz(); +palette(rainbow(max(V(g)$color),alpha=0.5)) +plot(g, layout=layout.fruchterman.reingold) +palette("default") +title('fastgreedy, layout.fruchterman.reingold, rvalue >0.2') + + +degree(g) +degree.distribution(g) +degree.distribution(g,cumulative = TRUE) + + + + +#------Histogram of degree distribution------------------------------------------------------------- +df <- data.frame(degree(g)) +colnames(df) <- c("degree") +p <- ggplot(df, aes(x=degree)) + xlab("degree") + theme_bw() +p + geom_histogram(binwidth = 2) + scale_colour_brewer(palette="Set1") + opts(aspect.ratio=1) #raw counts +ggsave(file=paste("120518_07-degreeDist", format(Sys.time(),"%y%m%d-%H%M%S"), ".pdf",sep=""))