Hello everybody,

I am using the Mfuzz package (v. 2.30.0) for analysing the effect of a time series on gene expression and I therefore manage to cluster the genes from this dataset, based on their expression profiles. I would like to know how it is possible to actually identify which genes are part of each cluster, in order to export these lists and make some gene ontology analysis. Is it possible with a function of this package or should I use another package?

I thank you in advance

Sincerely

Tom

The function acore() can create a list object with all the clusters and their genes with their gene membership values.

You should run something like that:

acore.list <- acore(CountTable,cl=cl)

temp   = sapply(acore.list,'[[',"NAME")

temp is than a list of vectors you can use for any downstream analysis