Skip to content

Instantly share code, notes, and snippets.

View t-redactyl's full-sized avatar

Jodie Burchell t-redactyl

357 followers · 6 following
  • JetBrains
  • Berlin, Germany
View GitHub Profile
@t-redactyl
t-redactyl / calculating_sem_formula.R
Last active September 27, 2017 05:34
Code associated with blog post: http://t-redactyl.github.io/blog/2015/09/A-Gentle-Introduction-to-the-Standard-Error-of-the-Mean.html
# Defining lambda and n
lambda <- 220
n <- 30
# Calculating SEM
sem <- sqrt(lambda / n)
@t-redactyl
t-redactyl / bootstrap_plot_1.R
Last active September 27, 2017 05:38
Code associated with blog post: http://t-redactyl.github.io/blog/2015/09/a-gentle-introduction-to-bootstrapping.html
require(ggplot2); require(gridExtra)
# Set the colours for the graphs
barfill <- "#4271AE"
barlines <- "#1F3552"
line1 <- "black"
line2 <- "#FF3721"
# Plotting histogram of sample of daily page views
g1 <- ggplot(data=as.data.frame(sample), aes(sample)) +
@t-redactyl
t-redactyl / comparing_small_and_larger_samples.R
Last active September 27, 2017 05:32
Code associated with blog post: http://t-redactyl.github.io/blog/2015/09/representative-sampling.html
set.seed(567)
# Sample of 30 (29 from the Poisson distribution and an outlier of 260)
sample1 <- c(rpois(29, lambda = 220), 260)
# Sample of 10 (9 from the Poisson distribution and an outlier of 260)
sample2 <- c(rpois(9, lambda = 220), 260)
@t-redactyl
t-redactyl / generate_95CI_around_t_value.R
Last active October 22, 2015 01:44
Code associated with blog post: http://t-redactyl.github.io/blog/2015/09/two-group-hypothesis-testing-independent-samples-t-tests.html
# Generate the 95% confidence interval.
lci <- -1 * qt(c(.975), 78)
uci <- qt(c(.975), 78)
@t-redactyl
t-redactyl / loading_in_depression_data.do
Last active October 22, 2015 01:41
Code for blog post: http://t-redactyl.github.io/blog/2015/09/interactive-model-decision-trees-in-stata.html
use http://www.ats.ucla.edu/stat/stata/library/depress, clear
reshape long dep, i(subj)
rename _j time
drop pre
@t-redactyl
t-redactyl / perm_alt_sample_plots.R
Last active October 22, 2015 01:43
Code associated with blog post: http://t-redactyl.github.io/blog/2015/10/two-group-hypothesis-testing-permutation-tests.html
# Load required packages
require(ggplot2); require(gridExtra)
# Set the colours for the graphs
barfill <- "#4271AE"
barlines <- "#1F3552"
line1 <- "black"
line2 <- "#FF3721"
# Plotting histogram of sample 1
@t-redactyl
t-redactyl / av_peds_2.R
Last active September 27, 2017 14:05
Code associated with blog post: http://t-redactyl.github.io/blog/2015/10/avoiding-the-crowds-in-melbourne-cbd-using-data.html
av_peds_2 <- ddply(p.subset, c("date", "collapsed_sensors_2"), summarise,
n_peds = sum(Hourly_Counts))
# Extract weekday versus weekend
av_peds_2$day <- weekdays(av_peds_2$date, abbreviate = FALSE)
av_peds_2$weekend <- ifelse((av_peds_2$day == "Saturday" | av_peds_2$day == "Sunday"),
"Weekend", "Weekday")
av_peds_2$weekend <- as.factor(av_peds_2$weekend)
# Extract time of day
@t-redactyl
t-redactyl / centred_chart.R
Created October 29, 2015 00:23
Code associated with blog post
library(ggplot2); library(gridExtra)
g1 <- ggplot(data=mtcars, aes(x=wt, y=mpg)) +
geom_point(alpha = 0.7, colour = "#0971B2") +
ylab("Miles per gallon") +
ylim(10, 35) +
xlab("Weight (`000 lbs)") +
ggtitle("Untransformed Weight") +
geom_vline(xintercept = 0) +
theme_bw()
@t-redactyl
t-redactyl / cleaning_data.R
Created November 4, 2015 05:52
Code associate with blog post
mtcars$am.f <- as.factor(mtcars$am); levels(mtcars$am.f) <- c("Automatic", "Manual")
mtcars$cyl.f <- as.factor(mtcars$cyl); levels(mtcars$cyl.f) <- c("4 cyl", "6 cyl", "8 cyl")
mtcars$vs.f <- as.factor(mtcars$vs); levels(mtcars$vs.f) <- c("V engine", "Straight engine")
mtcars$gear.f <- as.factor(mtcars$gear); levels(mtcars$gear.f) <- c("3 gears", "4 gears", "5 gears")
mtcars$carb.f <- as.factor(mtcars$carb)
@t-redactyl
t-redactyl / cat_class_1_method.py
Created November 12, 2015 03:02
Code associated with blog post:
def name_print(cat):
'''Print the name of the cat.'''
print "The cat is called %s." % cat.name
name_print(felix)