MultBar : Advanced multiple barplot with SEM

October 26, 2009
By
MultBar : Advanced multiple barplot with SEM

Producing this kind of graphs (below) in R can be a pain in the a*s. Here is a simple code that requires that data are presented in lists (see the example below). multbar <- function(list.of.lists,...,condnames=0,pal=colorRampPalette(c('grey','cornsilk')),seriesnames=0,legendpos='topleft',legh=TRUE,do.pty='s') { par(pty=do.pty,mgp=c(1.9,0.8,0),oma=c(0,0,0,0),mar=c(4,3,2,1),bg='transparent',bty='o',tck=0.02,yaxs='i') NofList <- length(list.of.lists) NofSubList <- length(list.of.lists]) if(condnames==0){condnames=c(1:NofList)} if(seriesnames==0){seriesnames=c(1:NofSubList)} dim.mat.treat <- (NofList+1)*NofSubList pos <- c(1:dim.mat.treat) pos <- matrix(pos,nrow=(NofList+1)) nbreaks

Read more »

R 2.10.0 is Out!

October 26, 2009
By
R 2.10.0 is Out!

The new R 2.10.0 is out! Get it from here. If you like take a look at these posts for some miscellaneous advices to make the upgrade easier. Fell free to contribute with suggestions about how to upgrade your R installation.

Read more »

Free statistics e-books for download

October 25, 2009
By

This post will eventually grow to hold a wide list of books on statistics (e-books, pdf books and so on) that are available for free download.  But for now we’ll start off with just one several books: The Elements of Statistical Learning written by Trevor Hastie, Robert Tibshirani and Jerome Friedman. you can legally download

Read more »

Example 7.16: assess robustness of permutation test to violations of exchangeability assumption

October 24, 2009
By
Example 7.16: assess robustness of permutation test to violations of exchangeability assumption

Permutation tests (section 2.4.3) are a form of resampling based inference that can be used to compare two groups. A simple univariate two-group permutation test requires that the group labels for the observations are exchangeable under the null hypothesis of equal distributions, but allows relaxation of specific distributional assumptions required by parametric procedures such as the t-test (

Read more »

A primer to Sweave with LyX

October 24, 2009
By

This is a step-by-step guide to 'Sweave' a document in R with Lyx on Ubuntu (9.04).System: Ubuntu 9.04, LyX 1.6.2, R 2.11.0 (Development version, r50178)Install LyX with Synaptic if it is not installed yet. (Current version 1.6.2)Start LyX, go to menu ...

Read more »

Data Mining and R

October 24, 2009
By
Data Mining and R

This post lists a few data mining resources in R. I also provide a few observations on the distinction between data mining, data analysis, and statistics as it pertains to the analysis work that I do in psychology.Online ResourcesThe classic book The E...

Read more »

RSPerl : Using R from within Perl

October 23, 2009
By

Some things I write in perl some in R, sometime I use perl to write R and run R. One thing that I find very useful is the functionality of RSPerl which enables you to call R functions from within perl and on perl variables. It can also call perl from R, though I have no idea why...

Read more »

RSPerl : Using R from within Perl

October 23, 2009
By

Some things I write in perl some in R, sometime I use perl to write R and run R. One thing that I find very useful is the functionality of RSPerl which enables you to call R functions from within perl and on perl variables. It can also call perl from R, though I have no idea why...

Read more »

Scivews-K got updated

October 23, 2009
By
Scivews-K got updated

http://www.sciviews.org/SciViews-K/index.htmlThe problem is, I cannot get it to work, on either Windows or Linux.

Read more »

From ORD Sessions to R-Forge in 12 hours with RProtoBuf

October 22, 2009
By

Yesterday, via in invitation from fellow Chicago-area Google Summer of Code mentor Borja Sotomayor, I attended the Second ORD Sessions. These are happening at the HQ of Inventable where a couple of technologists and Open Source geeks from the Chicagoland area get together and riff on code for a few hours after work over some pizza and beer. Sounded good, and I...

Read more »

From ORD Sessions to R-Forge in 12 hours with RProtoBuf

October 22, 2009
By

Yesterday, via in invitation from fellow Chicago-area Google Summer of Code mentor Borja Sotomayor, I attended the Second ORD Sessions. These are happening at the HQ of Inventable where a couple of technologists and Open Source geeks from the Chicagola...

Read more »

From ORD Sessions to R-Forge in 12 hours with RProtoBuf

October 22, 2009
By

Yesterday, via in invitation from fellow Chicago-area Google Summer of Code mentor Borja Sotomayor, I attended the Second ORD Sessions. These are happening at the HQ of Inventable where a couple of technologists and Open Source geeks from the Chicagoland area get together and riff on code for a few hours after work over some pizza and beer. Sounded good, and I...

Read more »

ISO week

October 22, 2009
By

I am working with a model that produces estimates of snow water equivalent through time. Because I deal with large spatial extents, I decided to have the model produce weekly averages. The problem with this is knowing which file to access for a given d...

Read more »

ISO week

October 22, 2009
By

I am working with a model that produces estimates of snow water equivalent through time. Because I deal with large spatial extents, I decided to have the model produce weekly averages. The problem with this is knowing which file to access for a given d...

Read more »

Leap years

October 22, 2009
By

A quick function that when provided a numeric vector of years returns a boolean vector where TRUE == Leap year.is.leapyear=function(year){ #http://en.wikipedia.org/wiki/Leap_year return(((year %% 4 == 0) & (year %% 100 != 0)) | (year %% 400 == 0))}

Read more »

Leap years

October 22, 2009
By

A quick function that when provided a numeric vector of years returns a boolean vector where TRUE == Leap year.is.leapyear=function(year){ #http://en.wikipedia.org/wiki/Leap_year return(((year %% 4 == 0) & (year %% 100 != 0)) | (year %% 400 == 0))}

Read more »

Aggreate electoral targeting with R

October 22, 2009
By
Aggreate electoral targeting with R

Aggregate electoral targeting is the process of determining the likely hood of a citizen choosing to vote turnout), and which candidate that person is most likely to vote for (partisan bias) in a given race using historical turnout and partisan bias. The output from this targeting allows a campaign to project the likely number of voters, what percentages their...

Read more »

ggplot2: Two Color XY-Area Combo Chart

October 21, 2009
By
ggplot2: Two Color XY-Area Combo Chart

David@Work blog shows how to fill in the area between two crossing lines in an Excel chart. This post was also published as a guest-post on PTS blog. Let’s try to replicate this graph in ggplot2. First, load ggplot2 and generate the data frame to be used in the example (I am using a slightly

Read more »

InferenceR looks interesting

October 21, 2009
By
InferenceR looks interesting

I obtained a free academic copy of the InferenceR, which is a R IDE for Windows. I has several interesting features such as a visual debugger. It is definitely worth some further exploration. Of course, the problem (and the only problem for now) is tha...

Read more »

[, [[, $: R accessors explained

October 21, 2009
By
[, [[, $: R accessors explained

For more than ten years, I have been teaching R both formally and informally. One thing that I find often trips up students is the use of R’s accessors and mutators. ( For those readers not from a formal computer science background, an accessor is a method for accessing data in an object

Read more »

Lots of maps using R

October 21, 2009
By
Lots of maps using R

http://r-spatial.sourceforge.net/gallery/#fig01.RVery helpful!

Read more »

AlgoCompSynth Reading List

October 21, 2009
By

Gerhard Niehaus Algorithmic Composition: Paradigms of Automated Music Generation by Gerhard Nierhaus Iannis Xenakis Formalized Music by Iannis Xenakis Music and Architecture: Architectural Projects, Texts, and Realizations by Iannis Xenakis Harry Partc...

Read more »

The geese are migrating, and so have I

October 20, 2009
By

topsyWidgetPreload({ "url": "http%3A%2F%2Fborasky-research.net%2F2009%2F10%2F20%2Fthe-geese-are-migrating-and-so-have-i%2F", "style": "big", "title": "The geese are migrating, and so have I" }); As someone wiser than I once said, “It’s fall – time to roast footballs and play touch marshmallow.” Or something like that. It’s the time of the year when geese migrate. From somewhere to somewhere else..Somewhere south of where they started.

Read more »

Editing while running R within Emacs

October 20, 2009
By
Editing while running R within Emacs

Simply press "Ctrl + g" will return the control back to the user.

Read more »

Running JAGS via R2jags

October 20, 2009
By
Running JAGS via R2jags

I finally figure out how to run WinBUGS examples in JAGS via R2jags. JAGS is preferred over WinBUGS for its cross-platform compatibility and performance (it was written in C++).

Read more »

Some Interesting Packages in R: swfDevice, RGoogleTrends, FlashMXML, SVGAnnotation

October 19, 2009
By

I love R because there are always exciting new packages which can be far beyond your imagination. Here I’d like to introduce a couple of packages that look really awesome: 1. swfDevice: R graphics device for SWF output (by Cameron Bracken) This package is still at a pre-alpha stage but you can see a sketch now in

Read more »

RPostgreSQL 0.1-6

October 19, 2009
By

A quick bug fix released of RPostgreSQL was uploaded to CRAN earlier today; it updates last week's 0.1-5 release. An editing error left one of the helper functions broken. I also updated a few of the test scripts to also allow for optional environmen...

Read more »

RPostgreSQL 0.1-6

October 19, 2009
By

A quick bug fix released of RPostgreSQL was uploaded to CRAN earlier today; it updates last week's 0.1-5 release. An editing error left one of the helper functions broken. I also updated a few of the test scripts to also allow for optional environment ...

Read more »

RPostgreSQL 0.1-6

October 19, 2009
By

A quick bug fix released of RPostgreSQL was uploaded to CRAN earlier today; it updates last week's 0.1-5 release. An editing error left one of the helper functions broken. I also updated a few of the test scripts to also allow for optional environmen...

Read more »