Hi Tina

From the baySeq vignette:

"The calculated priors are stored in the @priors slot of the countData object
produced as before... we build an empirical prior distribution
which we record in the list object $priors of the slot @priors. Each member
of this list object corresponds to one of the models defined by the group slot
of the countData object and contains the estimated parameters for each of the
individual counts selected under the models. The vector $sampled contained
in the slot @priors describes which rows were sampled to create these sets of
parameters."

So if you're following along using the code in the vignette you should be able to get the priors with:

write.table(CDP.NBML@priors$priors, 'test.txt', sep='\t', quote=F)

Best

Iain



________________________________
 From: Tina Asante Boahene <ma08tta@brunel.ac.uk>
To: Bioconductor mailing list <bioconductor@r-project.org> 
Sent: Monday, 6 February 2012, 13:37
Subject: [BioC] baySeq

Dear All,

I just hoping if you would not mind help me answer these questions

With regard to the plotMA I have realised that this plot was tailor to the data used my the author.
However, is there a general way to approach this to incorporate with different dataset.

plotMA.CD(CD, samplesA = c(1,3,6,8,10), samplesB = c(2,4,5,7,9), col = c(rep("red",
100), rep("black", 900)))

Also is there a way to write the output of the priors to a spreadsheet eg Excel. 

Thank you.


Kind Regards

Tina
_______________________________________________
Bioconductor mailing list
Bioconductor@r-project.org
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
	[[alternative HTML version deleted]]

