Collider Example

February 24, 2011
By

(This article was originally published at Adventures in R, and syndicated at StatsBlogs.)

library (MASS)
covar<-mvrnorm(250, c(0, 0), matrix(c(1, 0.00, 0.00, 1), 2, 2))
mydata<-data.frame(covar)
names(mydata)<-c("sat", "mot")
mydata$admin<- mydata$sat + mydata$mot
mydata$admin2<- ifelse (mydata$admin >=quantile(mydata$admin, .85), "pass", "fail")
 
library(ggplot2)
qplot(sat,mot, data = mydata, color = admin2)

Created by Pretty R at inside-R.org




Please comment on the article here: Adventures in R

Tags:

Subscribe

Email:

  Subscribe