Re: [R] [BioC] Overlay Gene Expression on SNP (copy number) data

2012-04-23 Thread Tengfei Yin
On Mon, Apr 23, 2012 at 6:33 AM, Ekta Jain wrote:

> Hello,
> Can anyone please suggest any packages in R that can be used to overlay
> gene expression data on SNP (affymetrix) copy number ?
>
Hi Ekta,

If you mean visually, as Steve suggested, you could try packages like
ggbio, Gviz, Rcytoscape.. it depends on how you plan to visualize your
data, track-based? circular view? net work? and what format your data are?

for example, in ggbio, it depends on what data you are using, you can
arrange your data into GRanges manually or just provide data  that
rtracklayer supported like bed, then just use autoplot, it accepts
different objects, like GRanges, IRanges, bamfiles or character... allow
some transformation like coverage. For files like bed, it automatically use
bar to represent your data and use score as y(you can specify other y).
Function tracks() allow you to bind or overlay any graphics produced by
ggbio or ggplot2, so you could work from data.frame too, it will help you
align your plots after the graphics are produced. For genomic structure, if
you want to overlay with your data, try autoplot, TranscriptDb. And if you
want to show interaction between genes, you could try either arches in
linear view or links in circular view(layout_circle).

http://tengfei.github.com/ggbio/

this website is still under development, just to show some possible cases,
it will be re-built against R 2.15 and more case studies are coming.

HTH

Tengfei



>
> Thanks,
> Ekta
> Senior Research Associate
> Bioinformatics Department
> Jubilant Biosys Pvt Ltd,
> #96, Industrial Suburb, 2nd Stage
> Yeshwantpur, Bangalore 560 022
> Ph No : +91-80-66628346
>
> The information contained in this electronic message and in any
> attachments to this message is confidential, legally privileged and
> intended only for use by the person or entity to which this electronic
> message is addressed. If you are not the intended recipient, and have
> received this message in error, please notify the sender and system manager
> by return email and delete the message and its attachments and also you are
> hereby notified that any distribution, copying, review, retransmission,
> dissemination or other use of this electronic transmission or the
> information contained in it is strictly prohibited. Please note that any
> views or opinions presented in this email are solely those of the author
> and may not represent those of the Company or bind the Company. Any
> commitments made over e-mail are not financially binding on the company
> unless accompanied or followed by a valid purchase order. This message has
> been scanned for viruses and dangerous content by Mail Scanner, a!
>  nd is believed to be clean. The Company accepts no liability for any
> damage caused by any virus transmitted by this email.
> www.jubl.com
>
>[[alternative HTML version deleted]]
>
> ___
> Bioconductor mailing list
> bioconduc...@r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>



-- 
Tengfei Yin
MCDB PhD student
1620 Howe Hall, 2274,
Iowa State University
Ames, IA,50011-2274
Homepage: www.tengfei.name

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] [BioC] Overlay Gene Expression on SNP (copy number) data

2012-04-23 Thread Paul Shannon
Another possibility for visual display of several kinds of data is RCytoscape, 
for which an example can be seen here:

http://rcytoscape.systemsbiology.net/versions/current/gallery/TCGA/subnet.TCGA.02.0014.png

This portrays
  1) gene expression (green: under-expression; red: over-expression) 
  2) copy number (blue: high, black: low)
  3) gene type (node shape: hexagons are kinases, arrows are ligands, diamonds 
are receptors, circles for everything else)
  4) mutation status (not SNPs, but non-synonymous amino acid substitutions)
  5) in the context of gene relationships, from KEGG

A vignette for a larger version of this network is nearly complete.  The 
Cytoscape network map is created from data and R code.

 - Paul

On Apr 23, 2012, at 5:15 AM, Steve Lianoglou wrote:

> Hi,
> 
> On Mon, Apr 23, 2012 at 7:33 AM, Ekta Jain  
> wrote:
>> Hello,
>> Can anyone please suggest any packages in R that can be used to overlay gene 
>> expression data on SNP (affymetrix) copy number ?
> 
> I guess you mean visually? If so, I'd suggest skimming through the
> vignettes of the following packages to see which one might suit you
> best:
> 
> * Gviz
> * ggbio
> * GenomeGraphs
> 
> -steve
> 
> -- 
> Steve Lianoglou
> Graduate Student: Computational Systems Biology
>  | Memorial Sloan-Kettering Cancer Center
>  | Weill Medical College of Cornell University
> Contact Info: http://cbio.mskcc.org/~lianos/contact
> 
> ___
> Bioconductor mailing list
> bioconduc...@r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: 
> http://news.gmane.org/gmane.science.biology.informatics.conductor

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] [BioC] Overlay Gene Expression on SNP (copy number) data

2012-04-23 Thread Steve Lianoglou
Hi,

On Mon, Apr 23, 2012 at 7:33 AM, Ekta Jain  wrote:
> Hello,
> Can anyone please suggest any packages in R that can be used to overlay gene 
> expression data on SNP (affymetrix) copy number ?

I guess you mean visually? If so, I'd suggest skimming through the
vignettes of the following packages to see which one might suit you
best:

* Gviz
* ggbio
* GenomeGraphs

-steve

-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.