50行ぐらいでRとBioCを理解した気になる講義です。 RとBioconductorの講義資料を置いておきます。 Ustはこちら。http://www.ustream.tv/channel/mntv R demo(graphics) # R 52.0/(1.63*1.63) w h w/h^2 22*h^2 # function a b c a[1]/a[2]^2 bmi bmi(a) # quit and session q() # matrix and apply wh wh[,1]/wh[,2]^2 summary(wh) ?summary example(summary) # push q key example(summary) apply(wh,1,bmi) >>> BioC # BioC library(affy) e is(e) head( exprs(e) ) write.table(exprs(e), "e.txt", sep="\t", eol="\n", quote=F, row.names=T, col.names=F) # interesting genes library('genefilter') # template matching e.cor e.cor exprs(e[e.cor$`100001_at`$indices]) library(annotate) # T-test t.filter t.filterfun which.t sum(which.t) e.t e.t # clustering and drawing heatmap pdf("heatmap.pdf", paper="a4") heatmap(exprs(e.t)) dev.off() >>> 詳しくは以下の本を参照してください。 {{amazon '4501622601