About Volcano Plot
1
0
Entering edit mode
@priscila-grynberg-3196
Last seen 10.2 years ago
Dear BioCs, I'd like to make a volcano plot, but instead of B value, I want to plot the adjusted p-values. Can I do that? Fold-changes X adjusted p-values? I have oligo slides. Do you need any further information for that? Thanks, Priscila -- Priscila Grynberg, B.Sc., M.Sc. Doutoranda em Bioinformática (Bioinformatics D.Sc student) Laboratório de Genética Bioquímica Universidade Federal de Minas Gerais Tel: +55 31 3409-2628 CV: http://lattes.cnpq.br/8808643075395963 [[alternative HTML version deleted]]
oligo oligo • 1.8k views
ADD COMMENT
0
Entering edit mode
@thomas-hampton-2820
Last seen 10.2 years ago
It isn't too hard to make your own volcano, and I do it all the time. You may find it easier to just use the R plot command like this rather than trying to wangle something in your bioconductor package: plot(M, -log10(P), pch = '.') assuming M is your vector of log2 fold differences and P is a vector of whatever you prefer as a significance measure. A very nice thing about doing it yourself is that you can then highlight a set of genes that interest you: points(M[interesting], -log10(P[interesting]), col = "red", pch = 16) Of course you can add all sorts of labels in the plot command. Tom On May 9, 2009, at 9:47 AM, Priscila Grynberg wrote: > Dear BioCs, > I'd like to make a volcano plot, but instead of B value, I want to > plot the > adjusted p-values. Can I do that? Fold-changes X adjusted p-values? > > I have oligo slides. Do you need any further information for that? > > Thanks, > > Priscila > > -- > Priscila Grynberg, B.Sc., M.Sc. > Doutoranda em Bioinform?tica (Bioinformatics D.Sc student) > Laborat?rio de Gen?tica Bioqu?mica > Universidade Federal de Minas Gerais > Tel: +55 31 3409-2628 > CV: http://lattes.cnpq.br/8808643075395963 > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/ > gmane.science.biology.informatics.conductor
ADD COMMENT
0
Entering edit mode
Dear Thomas, Thank you! It worked! Have a nice weekend! Priscila On Sat, May 9, 2009 at 11:03 AM, Thomas Hampton < Thomas.H.Hampton@dartmouth.edu> wrote: > It isn't too hard to make your own volcano, and I do it > all the time. You may find it easier to just use the R > plot command like this rather than trying to wangle something in > your bioconductor package: > > plot(M, -log10(P), pch = '.') > > > assuming M is your vector of log2 fold differences and P > is a vector of whatever you prefer as a significance measure. > > A very nice thing about doing it yourself is that you can then > highlight a set of genes that interest you: > > > points(M[interesting], -log10(P[interesting]), col = "red", pch = 16) > > > Of course you can add all sorts of labels in the plot command. > > Tom > > > > > On May 9, 2009, at 9:47 AM, Priscila Grynberg wrote: > > Dear BioCs, >> I'd like to make a volcano plot, but instead of B value, I want to plot >> the >> adjusted p-values. Can I do that? Fold-changes X adjusted p-values? >> >> I have oligo slides. Do you need any further information for that? >> >> Thanks, >> >> Priscila >> >> -- >> Priscila Grynberg, B.Sc., M.Sc. >> Doutoranda em Bioinformática (Bioinformatics D.Sc student) >> Laboratório de Genética Bioquímica >> Universidade Federal de Minas Gerais >> Tel: +55 31 3409-2628 >> CV: http://lattes.cnpq.br/8808643075395963 >> >> [[alternative HTML version deleted]] >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor@stat.math.ethz.ch >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: http://news.gmane.org/ >> gmane.science.biology.informatics.conductor >> > > -- Priscila Grynberg, B.Sc., M.Sc. Doutoranda em Bioinformática (Bioinformatics D.Sc student) Laboratório de Genética Bioquímica Universidade Federal de Minas Gerais Tel: +55 31 3409-2628 CV: http://lattes.cnpq.br/8808643075395963 [[alternative HTML version deleted]]
ADD REPLY

Login before adding your answer.

Traffic: 707 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6