Today I have been able to get 2 more templates working: The ternary plot and the correlation ellipses. They were already “working”, but just the initial script. which can be found here: ternary plot: http://addictedtor.free.fr/graphiques/RG...

Today I have been able to get 2 more templates working: The ternary plot and the correlation ellipses. They were already “working”, but just the initial script. which can be found here: ternary plot: http://addictedtor.free.fr/graphiques/RG...

After 2 restarts of my PC and 3 times restarting Rserve and Knime, the red text of doom (look at the bottom of Day 21) was finally gone. So I could start on my correlation ellipses. this wasn’t so hard. I added a panellistbox to my template, whic...

After 2 restarts of my PC and 3 times restarting Rserve and Knime, the red text of doom (look at the bottom of Day 21) was finally gone. So I could start on my correlation ellipses. this wasn’t so hard. I added a panellistbox to my template, whic...

For my second lecture today, I need to plot a likelihood surface for a basic two-component mixture with only the means unknown: here is the R code to speed up things llsurf=function(trumyn=2.,wayt=.3,var2=1.,ssiz=500){ # draws the log-likelihood surface and a random sample sd2=sqrt(var2) parti=(runif(ssiz)>wayt) sampl=(1-parti)*rnorm(ssiz)+parti*(trumyn+sd2*rnorm(ssiz)) mu2=mu1=seq(min(sampl),max(sampl),.1) mo1=mu1%*%t(rep(1,length(mu2))) mo2=(rep(1,length(mu2)))%*%t(mu2) ca1=-0.5*mo1*mo1 ca2=-0.5*mo2*mo2 like=.1*(ca1+ca2) # log prior N(0,10) for

The R Core Group has released the latest update to open-source R, R 2.13.0 is now available. You can download the new source distribution from the R Project website, and I see from the Download R tool at inside-R.org that binaries for Windows and Mac are already available at my local CRAN mirror (and, shortly, yours too). The (long)...

So everyone's excited about the new R 2.13 release because of the compiler package.Apparently it is easy to get a 3x speed increase by simply compiling a function.Doing a lot of the MCMC stuff, I am particularly excited about speed in R. I just compile...

So everyone's excited about the new R 2.13 release because of the compiler package.Apparently it is easy to get a 3x speed increase by simply compiling a function.Doing a lot of the MCMC stuff, I am particularly excited about speed in R. I just compile...

A brand new 0.9.4 release of Rcpp is now on CRAN and Debian. This version contains an improvement to loading and initialization of Rcpp modules, a bug fix for vectors of factors, another build issue fix as well as (per common practice with JSS) ...

It's been a year since I began keeping a web log. This post presents some thoughts related to the experience. Blogging is Sharing Ideas Blogging is online self-publishing. There is no faster way to share your ideas so broadly. Last year at the useR! conference (in Gaithersburg, MD, just a few months after I started

A maintenance release 0.2.3 of RProtoBuf is now on CRAN. RProtoBuf provides GNU R bindings for the Google Protobuf data encoding library used and released by Google. The NEWS file entry follows below: 0.2.3 2011-04-12 o Protect UINT64...

The recent opening of the Heritage Health Prize both represents a milestone and raises a cautionary flag. On the one hand, crowdsourced analytics prizes have never tackled anything so noble (not to discount predicting movie ratings), but on the other … Continue reading →

This is the last lesson of the R Videotutorial for spatial statistics.It is all about cokriging in gstat. For this lesson I used the meuse dataset, available within gstat, for the references to this dataset take a look at the script.The videotutorial i...

This is the last lesson of the R Videotutorial for spatial statistics.It is all about cokriging in gstat. For this lesson I used the meuse dataset, available within gstat, for the references to this dataset take a look at the script.The videotutorial i...

The latest edition of the Journal of Statistical Software is out, with plenty of interesting articles for R users. A must-read is Hadley Wickham's article on "The Split-Apply-Combine Strategy for Data Analysis", which makes a compelling argument for the use of the plyr package to partition datasets and apply sub-group analyses. Also, anyone who hasn't yet purchased a copy...

R 2.13 introduces an option to specify a custom graphics device in an Sweave code chunk. This is really cool and allows you to use tikzDevice output like pgfSweave does. In an pinch, say when you dont have access to any non-core packages, you can use tikzDevice output with the regular Sweave driver (RweaveLatex) more

I'm working on a project using next-gen sequencing to fine-map a genetic association in a gene region. Now that I've sequenced the region in a small sample, I'm picking SNPs to genotype in a larger sample. When designing the genotyping assay the lab wi...

R 2.13.0 will be released tomorrow. It contains a new package by R Core member Luke Tierney: compiler. What a simple and straightforward name, and something that Luke has been working on for several years. The NEWS file says o Package compiler is now provided as a standard package. See ...

R 2.13.0 will be released tomorrow. It contains a new package by R Core member Luke Tierney: compiler. What a simple and straightforward name, and something that Luke has been working on for several years. The NEWS file says o Package compiler is now provided as a standard package. See ...

Recently I started to build a new package for R containing Wilcox’ collection of functions for robust statistics. Wilcox provides 700+ functions for robust statistics, including: robust correlations (e.g. percentage bend correlation) robust measures of location and mean differences (e.g. … Continue reading →

Hi,in these days I'm creating lots and lots of box plot with ggplot2.The look of them is really good and you can change every bit of code so that you can customize the plot completely.Here is the code I'm using with a test data file to try it:BoxPlot.z...