Hi guys,
I'm playing around with cummeRbund, and I have some questions I would like to ask you.
First, I would like to cluster a subset of genes and to plot them in a heatmap. I found on this forum how to create a list of genes, and how to draw a heatmap.
> myIDs<-read.table("text.txt")
> myIDs<-as.vector(myIDs$V1)
> myGenes<-getGenes(cuff,myIDs)
> h<-csHeatmap(myGenes,cluster='both')
> h
Now that gives a heatmap with the genes I want, but I somehow think that the clustering is not adequate. I tried to do the partitioning thing like so :
> ic<-csCluster(myGenes,k=4)
And then, I figured I could only draw a heatmap with the output of csCluster
> hic<-csHeatmap(ic)
And of course it doesn't work.
Do you guys know how I could plot a heatmap from the csCluster?
Thanks!
I'm playing around with cummeRbund, and I have some questions I would like to ask you.
First, I would like to cluster a subset of genes and to plot them in a heatmap. I found on this forum how to create a list of genes, and how to draw a heatmap.
> myIDs<-read.table("text.txt")
> myIDs<-as.vector(myIDs$V1)
> myGenes<-getGenes(cuff,myIDs)
> h<-csHeatmap(myGenes,cluster='both')
> h
Now that gives a heatmap with the genes I want, but I somehow think that the clustering is not adequate. I tried to do the partitioning thing like so :
> ic<-csCluster(myGenes,k=4)
And then, I figured I could only draw a heatmap with the output of csCluster
> hic<-csHeatmap(ic)
And of course it doesn't work.
Do you guys know how I could plot a heatmap from the csCluster?
Thanks!
Comment