Blog Archives

Plug for RStudio: powerful, free, and easy to use interactive development environment for R

February 28, 2011
By
Plug for RStudio: powerful, free, and easy to use interactive development environment for R

(click for a bigger picture)As a longtime SAS user, one obstacle for me in using R professionally has been figuring out a process for saving and testing code across several work sessions and integrating code composition and execution. There are a coup...

Read more »

Example 8.20: Referencing lists of variables, part 2

January 10, 2011
By
Example 8.20: Referencing lists of variables, part 2

In Example 8.19, we discussed how to refer to a group of variables with sequential names, such as varname1, varname2, varname3. This is trivial in SAS and can be done in R as we showed.It's also sometimes useful to refer to all variables which begin w...

Read more »

Example 8.18: A Monte Carlo experiment

December 13, 2010
By
Example 8.18: A Monte Carlo experiment

In recent weeks, we've explored methods to fit logistic regression models when a state of quasi-complete separation exists. We considered Firth's penalized likelihood approach, exact logistic regression, and Bayesian models using Markov chain Monte Ca...

Read more »

Example 8.17: Logistic regression via MCMC

December 6, 2010
By
Example 8.17: Logistic regression via MCMC

In examples 8.15 and 8.16 we considered Firth logistic regression and exact logistic regression as ways around the problem of separation, often encountered in logistic regression. (Re-cap: Separation happens when all the observations in a category sha...

Read more »

Example 8.16: Exact logistic regression

November 30, 2010
By
Example 8.16: Exact logistic regression

In example 8.15, on Firth logistic regression, we mentioned alternative approaches to separation troubles. Here we demonstrate exact logistic regression. The code for this appears in the book (section 4.1.2) but we don't show an example of it there. ...

Read more »

SAS and R joins SAS-x

November 29, 2010
By
SAS and R joins SAS-x

Tal Galili, organizer of the R-bloggers blog aggregator, has opened a new aggregator for people blogging about SAS. If you're unfamiliar with the concept, an aggregator is a single blog which republishes (with permission, in this case) the entries fro...

Read more »

Example 8.15: Firth logistic regression

November 22, 2010
By
Example 8.15: Firth logistic regression

In logistic regression, when the outcome has low (or high) prevalence, or when there are several interacted categorical predictors, it can happen that for some combination of the predictors, all the observations have the same event status. A similar e...

Read more »

Example 8.13: Bike ride plot, part 2

November 8, 2010
By
Example 8.13: Bike ride plot, part 2

Before explaining how to make and interpret the plot above, Nick and I want to make a plea for questions--it's hard to come up with useful questions to explore each week!As shown in Example 8.12, data from the Cyclemeter app can be used to make interes...

Read more »

Example 8.12: Bike ride plot, part 1

November 1, 2010
By
Example 8.12: Bike ride plot, part 1

The iPhone app Cyclemeter uses the phone's GPS capability to record location and other data, and infer speed, while you ride. I took a ride near my house recently, and downloaded the data. I'd like to examine my route and my speed. A simple plot of ...

Read more »

Example 8.9: Contrasts

October 12, 2010
By
Example 8.9: Contrasts

In example 8.6 we showed how to change the reference category. This is the natural first thought analysts have when their primary comparisons aren't represented in the default output. But our interest might center on a number of comparisons which don...

Read more »