This report is a RMarkdown file describing a project on the topic “Televised Presidential Ads and Voter Turnout” on GitHub. This project is created by the user “skleffingwell” on GitHub. It uses the dataset “KrasnoGreen_JOP_2008.csv” for data analysis and interpretation. This dataset could be downloaded from the project link given below.
main <- read.csv("KrasnoGreen_JOP_2008.csv")
library(lmtest)
## Loading required package: zoo
##
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric
lm.1 <- lm(turn00~senate+prezg,data=main)
bptest(lm.1, turn00~lncpop00+pq,data=main,studentize=F)
##
## Breusch-Pagan test
##
## data: lm.1
## BP = 22.564, df = 2, p-value = 1.26e-05
summary(lm.1)$adj.r.squared
## [1] 0.6035078
summary(lm.1)
##
## Call:
## lm(formula = turn00 ~ senate + prezg, data = main)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.110158 -0.025510 -0.000674 0.026931 0.132518
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.2200101 0.0233354 9.428 2.87e-16 ***
## senate 0.7816509 0.0564396 13.849 < 2e-16 ***
## prezg 0.0025384 0.0008621 2.944 0.00386 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.04201 on 125 degrees of freedom
## Multiple R-squared: 0.6098, Adjusted R-squared: 0.6035
## F-statistic: 97.65 on 2 and 125 DF, p-value: < 2.2e-16
coef(lm.1)[["prezg"]]
## [1] 0.002538412
permuteFn_lm.1<-function(y,a){
shuffled_prezg<-sample(main$prezg)
coef<-coef(lm(y~shuffled_prezg+a))[["shuffled_prezg"]]
return(coef)
}
results_p_lm.1<-replicate(10000,permuteFn_lm.1(y=main$turn00, a=main$senate))
mean(results_p_lm.1)
## [1] 4.864166e-06
plot(density(results_p_lm.1/sd(results_p_lm.1)),xlab="Mean Diffs")
curve(dnorm(x,sd=1,mean=0),from=-4,to=4,col="blue")
characters <- colnames(main)
table_prezg <- table(main$prezg)
logic <- is.character("prezg")
lst <- list(characters, table_prezg, logic)