Question: what Markov clustering and random walk with restart need as input ?
0
gravatar for Chaimaa
2.5 years ago by
Chaimaa150
Chaimaa150 wrote:

HI friends, I have a dense unweighted biological network and I want to cluster this network to subnetworks by using Markov cluster or random walk with restart .

Igot the adjacency matrix from this network but when i tried markov cluster it produce for me NaN as a result .

In Random walk with restart i don't have p0 and pt. plz anyone can help me to !!!.

networks matlab rwr mcl • 866 views
ADD COMMENTlink modified 2.5 years ago by Jean-Karim Heriche21k • written 2.5 years ago by Chaimaa150
0
gravatar for Jean-Karim Heriche
2.5 years ago by
EMBL Heidelberg, Germany
Jean-Karim Heriche21k wrote:

Check that the format of the input is suitable for the implementation of the algorithm you're using. If you're using the original MCL implementation then read the manual, in particular this page about input formats.
For your problem with random walk with restart, you'll have to give more info: What are you trying to do ? What implementation are you using ?

ADD COMMENTlink written 2.5 years ago by Jean-Karim Heriche21k

@Jean-Karim Heriche friend, the only input i have is adjacency matrix in which i get it from the network; i already build this adjacency matrix which has only 1 and 0. MCL algorithm need as input this adjacency matrix, and so i get NaN as output, then i turned to use a Random walk with restart which required W, p0, and pt as inputs but sadly i only have the W(adjacency matrix). My goal is to cluster my dense network into subnetworks. I hope anyone already tried these methods can explain me how to deal with them. I'm using Matlab.

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by Chaimaa150

I know you have the adjacency matrix but did you make sure it is in the right format for input to whatever program you're using. Without more details, I don't know what p0 and pt are. Matlab is not very much used in this community so you may have better luck getting answers on Matlab mailing lists or by contacting the authors of the code you're using.

ADD REPLYlink written 2.5 years ago by Jean-Karim Heriche21k

ok i will trying to read more about thess methods

ADD REPLYlink written 2.5 years ago by Chaimaa150

@Jean-Karim Heriche friend, do you know any methods can cluster or partitioned a dense unweighted biological network into small subnetworks?

ADD REPLYlink written 2.5 years ago by Chaimaa150
1

It depends on the particular type of data and what criterion you want to work with. MCL can work well in some instances but I find the clusters too unbalanced. I often get good results with normalized cut spectral clustering. However, you don't have to restrict yourself to graph partitioning methods. In general, you can derive a similarity/distance matrix from the adjacency matrix (for example using a kernel on graph nodes) and apply a clustering algorithm to this new matrix (this is actually one way of viewing spectral clustering).

ADD REPLYlink written 2.5 years ago by Jean-Karim Heriche21k

@Jean-Karim Heriche; I really appreciate your help bro, everytime your comments make so many idea more clear for me, Thank you too too much.

ADD REPLYlink written 2.5 years ago by Chaimaa150

@Jean-Karim Heriche, MCL needed a transition matrix. So, i have updated the adjacency matrix which i got from the network into transition matrix (dividing every value in the adjacency matrix by the column sum) and then i assign a value of 1 to the diagonal of this matrix to avoid NaN [since so many columns in my matrix have all 0's and no 1's, and when Matlab start running produces NaN (0/0)]. Then i run the MCL code on Matlab and got a new matrix after many expansion and inflation iterations. But, i see the new matrix have some self-loops, is it right?? I need your suggestions! Also, i checked Spectral clustering; it takes as input a square adjacency matrix and a number of clusters to look for. I don't know how many clusters i should take ?

ADD REPLYlink written 2.5 years ago by Chaimaa150

If you add 1 to the diagonal of a transition matrix then you don't have a transition matrix anymore because the resulting matrix is not stochastic. You have to add the loops to your graph before you derive the transition matrix. As for choosing the number of clusters for spectral clustering, I usually get good results with the eigengap heuristic as long as there's some structure in the graph. You may want to have a look at this tutorial on spectral clustering.

ADD REPLYlink written 2.5 years ago by Jean-Karim Heriche21k

So, if I didn't add 1 to the diagonal, is the transition matrix (dividing every value in the adjacency matrix by the column sum) I got right? Plz confirm me, and so sorry if I bother you. I'm doing research alone, and getting confused !

ADD REPLYlink written 2.5 years ago by Chaimaa150
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 729 users visited in the last hour