Hi
I have a gene presence-absence matrix of accessory genes, for multiple strains of bacteria. I would like to produce a good figure with this data, and maybe include a tree as well . I was wondering what the best software programmes are available to do this? As I have >4000 genes, I don't each gene to be labelled!
Thanks
•
written
by
natasha
Maybe try a PCA on your data. 4000 observations is quiet big to plot in a figure..
or you could do both (a. PCA ; b. heatmap)
Perhaps a heatmap would do ?
For 4000 observations? I doubt whether that would send a clear message. Although a heatmap would be a nice combination with a hierarchical clustering tree.
You want to produce a good figure, but you should think first which message you want that figure to contain.
That's still a lot. What is the message your figure should send? You could try as suggested a PCA with an additional hierarchical clustering (tree) and a PCA and see how that looks. There are many R packages that can do that for you.
How about a gene presence/absence network, or project the presence/absence onto a known phylogeny, like Figures 2 and 3 in the CoPAP paper? http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3692100/
Log in to answer this question.
how many strains do you have ?
50 strains. I make a mistake, I have 1800 genes to plot not 4000.