fditraglia

#Function to draw a single Chi-squared random variable with degrees of freedom equal to df
my.rchisq <- function(df){
  
  #Draw df independent standard normals
  normal.sims <- rnorm(df)
  
  #Square them and sum the result
  chi.sims <- sum(normal.sims^2)

fditraglia

ones <- rep(1, 10)
ones
cumsum(ones)
cumsum(1:10)



n <- 10000
sims <- rnorm(n, mean = 0, sd = 10)
running.mean <- cumsum(sims)/(1:n)

fditraglia

%  I learned how to make these flowcharts by looking at some very helpful examples online:
%
%    http://www.texample.net/tikz/examples/simple-flow-chart/
%    http://www.texample.net/tikz/examples/borrowers-and-lenders/
%    http://www.texample.net/tikz/examples/entity-relationship-diagram/
%
%  For a long list of tikz examples, see http://www.texample.net/tikz/examples/		

\documentclass{beamer}  
\usepackage{amssymb, amsmath, amsthm} 

fditraglia

pnorm(10, mean = 12, sd = 2)



pnorm(250, mean = 300, sd = 10)



n <- 950:961
width <- 2 * qnorm(1 - 0.01/2) * 6/sqrt(n)

fditraglia

group <- c('treatment', 'control', 'refused')
n.children <- 1000 * c(200, 200, 340)
n.polio <- c(57, 142, 157)
rate <- n.polio/n.children
polio.data <- data.frame(group, n.children, n.polio, rate)
polio.data



treatment <- subset(polio.data, group == 'treatment')

fditraglia

coin.power <- function(a, p, n){
  c <-  qnorm(1 - a/2)  
  mu <- sqrt(n) * (2 * p - 1)   
  sigma <- sqrt(4 * p * (1 - p))
  
  less.than <- pnorm( -c, mean = mu, sd = sigma)
  greater.than <- 1 - pnorm(c, mean = mu, sd = sigma)
  power <- less.than + greater.than
  return(power)
}

fditraglia

data.url <- "http://www.ditraglia.com/econ103/survey_clean.csv"
survey <- read.csv(data.url)
anchoring <- survey[,c("rand.num", "africa.percent")]
lo <- subset(anchoring, rand.num == "10")$africa.percent
hi <- subset(anchoring, rand.num == "65")$africa.percent
lo <- lo[!is.na(lo)]
hi <- hi[!is.na(hi)]
SE.lo <- sd(lo)/sqrt(length(lo))
SE.hi <- sd(hi)/sqrt(length(hi))
SE <- sqrt(SE.hi^2 + SE.lo^2)

fditraglia

install.packages("arm")
library('arm')

data <- read.csv("http://www.ditraglia.com/econ103/child_test_data.csv")
attach(data)


head(data)

fditraglia

data.url <- 'http://www.ditraglia.com/econ103/ex_13_5.csv'
election <- read.csv(data.url)
lm(y ~ x1 + x2, data = election)



library(arm)
data.url <- "http://www.ditraglia.com/econ103/child_test_data.csv"
data <- read.csv(data.url)
attach(data)

fditraglia

<snippet>
	<content><![CDATA[
\documentclass[12pt]{article}
\usepackage[margin=1.5in]{geometry}
\usepackage{todonotes}
\usepackage{amssymb,amsmath,amsthm,graphicx}
\usepackage{geometry}
\usepackage{rotating}
\usepackage{amssymb, amsmath, amsthm, graphicx} 
\usepackage{multirow}