User:Timothee Flutre/Notebook/Postdoc/2011/11/07: Difference between revisions
From OpenWetWare
(→About R: add links helping in package dev) |
(→About R: add swirl) |
||
Line 12: | Line 12: | ||
* '''Documentation''': | * '''Documentation''': | ||
** try it [https://www.codeschool.com/courses/try-r online] | ** try it [https://www.codeschool.com/courses/try-r online] | ||
** learn how to use it with [http://swirlstats.com/ swirl] | |||
** official introductory [http://cran.r-project.org/doc/manuals/R-intro.html manual] | ** official introductory [http://cran.r-project.org/doc/manuals/R-intro.html manual] | ||
** well-organized [http://www.statmethods.net/ how-to] | ** well-organized [http://www.statmethods.net/ how-to] |
Revision as of 11:08, 6 February 2014
Project name | <html><img src="/images/9/94/Report.png" border="0" /></html> Main project page <html><img src="/images/c/c3/Resultset_previous.png" border="0" /></html>Previous entry<html> </html>Next entry<html><img src="/images/5/5c/Resultset_next.png" border="0" /></html> |
About R
S <- 3 # nb of subgroups V <- 7 # nb of observations z <- matrix(c(0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,1,1,1,0,0), nrow=V, ncol=S, byrow=TRUE) myheatmap <- function(z, out.file="") { def.par <- par(no.readonly=TRUE) par(mar=c(4,5,3,2), font=2, font.axis=2, font.lab=2, cex=1.5, lwd=2) if (out.file != "") pdf(out.file) layout(mat=cbind(1, 2), width=c(7,1)) # plot + legend mycol <- rev(heat.colors(4)) image(x=1:NCOL(z), y=1:NROW(z), z=t(z), xlim=0.5+c(0,NCOL(z)), ylim=0.5+c(0,NROW(z)), xlab="", ylab="Observations sorted by cluster", main="Custom heatmap", axes=FALSE, col=mycol) axis(1, 1:NCOL(z), labels=paste("subgroup", 1:NCOL(z)), tick=0) par(mar=c(0,0,0,0)) plot.new() legend("center", legend=sprintf("%.2f", seq(from=min(z), to=max(z), length.out=5)[-1]), fill=mycol, border=mycol, bty="n") if (out.file != "") dev.off() par(def.par) } myheatmap(mydata.sort) |