R

From CSBLwiki

Contents 1 Tutorials 1.1 Tips 2 Packages 2.1 ape 2.2 DBI 2.3 seqinr 2.4 Bioconductor 3 Examples 3.1 UMLS 4 Links

Tutorials

• R Bioconductor

Bioconductor course materials

2010 introductory course

• R Graphical Manual

• A Quick Guide to Teaching R Programming to Computational Biology Students - PLoS CompBio

Tips

The script editor with R.gui is a convenient tool (c.f. ESS with Emacs)

Command + Enter (OS X) or Ctrl + R (windows)

Packages

• Installed @panflam, @panpyro

ape

• an object of class "phylo" with the following components:
  • edge a two-column matrix of mode numeric where each row represents an edge of the tree; the nodes and the tips are symbolized with numbers; the tips are numbered 1, 2, ..., and the nodes are numbered after the tips. For each row, the first column gives the ancestor.
  • edge.length (optional) a numeric vector giving the lengths of the branches given by edge.
  • tip.label a vector of mode character giving the names of the tips; the order of the names in this vector corresponds to the (positive) number in edge.
  • Nnode the number of (internal) nodes.
  • node.label (optional) a vector of mode character giving the names of the nodes.
  • root.edge (optional) a numeric value giving the length of the branch at the root if it exists.

DBI

seqinr

Bioconductor

• course materials
  • 2010 introductory course

Examples

UMLS

library(RMySQL)
library(DBI)
drv = dbDriver("MySQL")
con = dbConnect(drv,username='',password='',dbname='')
que = "select CUI1,CUI2,COF from MRCOC where SAB=\'CCPSS\'"
ccpss = dbSendQuery(con,statement=que)
ccpss.all = unlist(fetch(ccpss,n=-1))
ccpss.all = matrix(ccpss.all,length(ccpss.all)/3,3)
colnames(ccpss.all) = c('CUI1','CUI2','COF')
tmp = unique(c(ccpss.all[,1],ccpss.all[,2]))
cof = matrix(NA,10877,10877)
dimnames(cof) = list(tmp,tmp)
apply(ccpss.all,1,function(x) { cof[x[1],x[2]]=x[3] })
cof[1:10,1:10]

Links

• News article about R in NYTimes