The latest version of oro.dicom (0.2.5) has been released on CRAN. New features include:Added "mosaic" capability when creating 3D arrays from DICOMdicomTable() now accepts single DICOM fileBetter handling of SequenceItem tags when reading in DIC...

The latest version of oro.dicom (0.2.5) has been released on CRAN. New features include:Added "mosaic" capability when creating 3D arrays from DICOMdicomTable() now accepts single DICOM fileBetter handling of SequenceItem tags when reading in DIC...

The latest version of oro.dicom (0.2.5) has been released on CRAN. New features include:Added "mosaic" capability when creating 3D arrays from DICOMdicomTable() now accepts single DICOM fileBetter handling of SequenceItem tags when reading in DIC...

This is the announcement as posted in the mailing list : This is to announce that we plan to release R version 2.11.0 on Thursday, April 22, 2010. Those directly involved should review the generic schedule at http://developer.r-project.org/release-checklist.html The source tarballs will be made available daily (barring build troubles) via http://cran.r-project.org/src/base-prerelease/ For the R Core

Over a month ago, David Smith published a call for people to participate in the “Future of Open Source” Survey. 550 people (and me) took the survey, and today I got an e-mail with the news that the 2010 survey results are analysed and where published in the “Future.Of.Open.Source blog” In the following (38 slides) presentation: I would like...

Hadley Wickham, creator of ggplot2, has posted a 2 hour video on data visualization using R. You can find links to the videos and slides over at Revolutions Blog. Check back here soon. I am working with Hadley to arrange a day-long ggplot2 short cours...

Hadley Wickham (the creator of the popular ggplot2 graphics package for R) has posted video of a 2-hour short course on Visualisation in R at his blip.tv channel. The video is split into four thirty-minute segments: Basic Graphics Displaying Large Data Data manipulation and transformations Polishing your plots for publication The course is peppered with self-guided exercises, for which...

Every Six months or so I like to update the Site Mapfor my blog jeromyanglim.blogspot.com. Thus, if you're new to the blog, The Site Map provides an overview of previous blog content organised by topic. RSS Feed: http://feeds.feedburner.com/jeromyang...

Every Six months or so I like to update the Site Mapfor my blog jeromyanglim.blogspot.com. Thus, if you're new to the blog, The Site Map provides an overview of previous blog content organised by topic. RSS Feed: http://feeds.feedburner.com/jeromyang...

Abilities of R for creating graphics is great, but one thing I always missed is the possibility of creating interactive plots and being able to look at graphs while changing one ore more parameters. I know that there is rggobi, but so far I always ran into problems with flexibility each time I wanted to

“How to change the order of legend labels” is a question that gets asked relatively often on ggplot2 mailing list. A variation of this question is how to change the order of series in stacked bar/lineplots. While these two questions seem to be related, in fact they are separate as the legend is controlled by

Here are some (trivial) R tips in the course Stat 511. I’ll update this post till the semester is over. Formatting R Code Reading code is pain, but the well-formatted code might alleviate the pain a little bit. The function tidy.source() in the animation package can help us format our R code automatically. By default

The shiny new 0.2.2 release of RInside has just been uploaded to CRAN; it should hit mirrors tommorow. Sources are also at my RInside page. RInside is a set of convenience classes to facilitate embedding of R inside of C++ applications. It works parti...

The shiny new 0.2.2 release of RInside has just been uploaded to CRAN; it should hit mirrors tommorow. Sources are also at my RInside page. RInside is a set of convenience classes to facilitate embedding of R inside of C++ applications. It works pa...

Announced this morning on the r-announce mailing list is the impending release of R 2.11.0, scheduled for April 22. As usual, the release soon goes into a beta-test phase, with updated sources to be available on the 22nd and binaries to follow a few days thereafter.

Our FEC report file library FECHell has been updated to 0.1.9. The release includes a half dozen bug fixes and the following new features: Speed improvements – Schedule and field names names are matched by compiled regular expressions instead of brute-force string matching, resulting in a ~25% speed increase for large files. DEF file fixes

Download "Getting Started with the Social Media Analytics Research Toolkit" (pdf, 1.25 megabytes) Download the Social Media Analytics Research Toolkit How The Chart Was Made Whenever a Twitter user posts a tweet, an object is created and entered into t...

So I took a brief break from ecology yesterday to pursue another project, modeling the NCAA bracket. I decided to run a stochastic simulation of the tournament. I know this has nothing to do with ecology, but I believe its a good exercise in modeling. Modeling has very little to do with the specifics of a...

The annual March Madness Half Marathon in Cary took place this morning. This is both one of Chicagoland's 'early races' to start the season as well as the classic Boston preparation due to the hilly course. I have now run this consecutively for six y...

Siemens multi-slice EPI data may be collected as a "mosiac" image; i.e., all slices acquired in a single TR (repitition time) of a functional MRI run are stored in a single DICOM file. The images are stored in an MxN array of images. The function create3D() will try to guess the number of images embedded within the single DICOM...

If you have a plot open and want to add a vertical line to it: abline(v=20) #Add vertical line at x=20

There was a post here about obtaining non-standard p-values for testing the correlation coefficient. The R-library SuppDists deals with this problem efficiently. library(SuppDists) plot(function(x)dPearson(x,N=23,rho=0.7),-1,1,ylim=c(0,10),ylab="density") plot(function(x)dPearson(x,N=23,rho=0),-1,1,add=TRUE,col="steelblue") plot(function(x)dPearson(x,N=23,rho=-.2),-1,1,add=TRUE,col="green") plot(function(x)dPearson(x,N=23,rho=.9),-1,1,add=TRUE,col="red");grid() legend("topleft", col=c("black","steelblue","red","green"),lty=1, legend=c("rho=0.7","rho=0","rho=-.2","rho=.9"))</pre> This is how it looks like, Now, let’s construct a table of critical values for some arbitrary or not significance levels. q=c(.025,.05,.075,.1,.15,.2) xtabs(qPearson(p=q, N=23, rho

My parents and I made plans to visit San Jose and Saratoga on my grandmother’s birthday, March 19, since that is where she grew up. I randomly saw someone tweet about the ACM Data Mining Camp unconference that happened to be the next day, March 20, only a couple of miles from our hotel in Santa Clara. This was...

gm(x) But this requires package heR.Misc so you might as well just use: exp(mean(log(x)))

Inspired by CXO group report, I did a rerun of the same strategy on my data. Easter’s dates can be find at wikipedia. Overall, my results are similar to CXO group’s results. In the graph below, I plotted daily returns on Easter week (Monday to Thursday) from 1982 to 2009. I prefer this way of showing

Readers returning to our blog will know that Win-Vector LLC is fairly “pro-R.” You can take that to mean “in favor or R” or “professionally using R” (both statements are true). Some days we really don’t feel that way. Consider the following snippet of R code where we create a list with a single element Related posts: