help
Hi, the thing is, you didn't really have any questions in your email. What is the problem with your code? Is it not giving the expected results? More below. On Wed, Aug 22, 2012 at 8:31 AM, Yanis El Omari
<yanis.elomari at esma.europa.eu> wrote:
[...]
Rm(list())
What is this this supposed to do? Do you mean 'rm(list=ls())' ?
library(R.matlab)
library(igraph)
path<-("D:/adj_matrices/nonweighted_GOVSOV")
setwd(path)
vector.files<-paste("nonweighted_GOVSOV",1:265, ".mat",sep="")
dat<-list()
for(i in 1:265) {
dat[[i]] <-readMat(file.path(path,vector.files[i]))
adj<-do.call(rbind, dat[[i]])
g<-graph.adjacency(adj, mode=c("directed"), weighted=T) # the graph object of the network
a<-length(V(g)) #the number of nodes
indeg<-rep(0,a) #the vector where I want to write the degree of each node for week i
The above two lines are not needed, you can remove them.
indeg<-degree(g) #degree calculation for each node of the network write(indeg, "indeg[i]".txt, sep="\n")
This is probably wrong, the results are written to a file named
"indeg[i].txt", I think you want something like
write(indeg, sprintf("indeg[%i].txt", i), sep="\n")
} Is there some obvious mistake in my program, or does anyone have a solution to do this?
Maybe, but what is your problem with it? (Apart from the incorrect naming of the output files.) Gabor [...]
Gabor Csardi <csardi at rmki.kfki.hu> MTA KFKI RMKI