Skip to content

Instantly share code, notes, and snippets.

@thomasjensen
thomasjensen / tikzggplot.r
Created December 7, 2011 10:56
Use tkizDevice with ggplot2
library(tikzDevice)
library(ggplot2)
y <- exp(seq(1,10,.1))
x <- 1:length(y)
data <- data.frame(x = x, y = y)
tikz(file = "test.tex")
plot <- ggplot(data, aes(x = x, y = y)) + geom_line()
@thomasjensen
thomasjensen / epOutliers.r
Created December 6, 2011 22:32
Outliers in the EP
#read the ggplot2 library
library(ggplot2)
#read the data
mepDat <- read.csv("ep6.csv")
#plot the density for the questions variable
plot <- ggplot(mepDat, aes(x = Questions))
plot <- plot + geom_density()
plot <- plot + geom_segment(aes(x = median(mepDat$Questions), xend = median(mepDat$Questions), y = 0, yend = .031), linetype = 2)
@thomasjensen
thomasjensen / simulated_effect.r
Created December 6, 2011 11:29
Simulate confidence intervals
library(MASS)
library(ggplot2)
#generate the fake data
covvar <- matrix(c(1,.6,.6,.6,1,.6,.6,.6,1),nr = 3)
means <- c(20,25,30)
data <- as.data.frame(mvrnorm(100,means,covvar))
#dichotomize the dependent variable
data$V1 <- ifelse(data$V1 > 20,1,0)
@thomasjensen
thomasjensen / interaction.r
Created December 5, 2011 23:43
Interpret Interaction Terms
#generate data
library(MASS)
covvar <- matrix(c(1,.6,.6,.6,1,.6,.6,.6,1),nr = 3)
means <- c(20,25,30)
data <- as.data.frame(mvrnorm(100,means,covvar))
#dichotomize the dependent variable
data$V1 <- ifelse(data$V1 > 20,1,0)
@thomasjensen
thomasjensen / simple.r
Created December 5, 2011 15:59
simpleR
x <- c(1,2,3)
plot(x)