Denver Justice 2014 QTR 1

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

denverjustice<-read.csv("/Volumes/TOSHIBA EXT/Dropbox/Schools Study Data/Denver Justice Data and R Scripts/denver_justice_Aug-Oct14_complete.csv", header=T, fill=T)

Meaning in Life Questianir

library(psych)
denverjustice$meanMLQ <- apply(denverjustice[, c("MLQ_1","MLQ_2","MLQ_3","MLQ_4","MLQ_5","MLQ_6", "MLQ_7", "MLQ_8", "MLQ_10")], 1, mean, na.rm = TRUE)
t.test(denverjustice$meanMLQ ~ Time, data = denverjustice, paired = TRUE)

## 
##  Paired t-test
## 
## data:  denverjustice$meanMLQ by Time
## t = 0.383, df = 9, p-value = 0.7106
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.7632  1.0743
## sample estimates:
## mean of the differences 
##                  0.1556

boxplot(denverjustice$meanMLQ ~ Time, data=denverjustice, notch=F, col=(c("gold","blue")), main="MLQ All", xlab="Meaning and Purpose")

``` MLQ Searching

denverjustice$MLQS <- apply(denverjustice[, c("MLQ_2" ,"MLQ_3", "MLQ_7", "MLQ_8", "MLQ_10")], 1, mean, na.rm = TRUE)
t.test(denverjustice$MLQS ~ Time, data = denverjustice, paired = TRUE)

## 
##  Paired t-test
## 
## data:  denverjustice$MLQS by Time
## t = 0.979, df = 9, p-value = 0.3532
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.4194  1.0594
## sample estimates:
## mean of the differences 
##                    0.32

plot(denverjustice$MLQS)

barplot(denverjustice$MLQS)

boxplot(denverjustice$MLQS ~ Time, data=denverjustice, notch=F, col=(c("gold","blue")), main="MLQ", xlab="Searching Meaning and Purpose")

Positive Questions from MLQ

denverjustice$MLQP <- apply(denverjustice[, c("MLQ_1" ,"MLQ_4", "MLQ_5", "MLQ_6")], 1, mean, na.rm = TRUE)
t.test(denverjustice$MLQP ~ Time, data = denverjustice, paired = TRUE)

## 
##  Paired t-test
## 
## data:  denverjustice$MLQP by Time
## t = -0.0952, df = 9, p-value = 0.9263
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -1.239  1.139
## sample estimates:
## mean of the differences 
##                   -0.05

plot(denverjustice$MLQP)

barplot(denverjustice$MLQP)

boxplot(denverjustice$MLQP ~ Time, data=denverjustice, notch=TRUE, col=(c("gold","blue")), main="MLQ", xlab="Meaning and Purpose")

## Warning: some notches went outside hinges ('box'): maybe set notch=FALSE

Culture of Purpose Survey

denverjustice$meanCPS <- apply(denverjustice[, c("CPS_1","CPS_2","CPS_3","CPS_4","CPS_5","CPS_6", "CPS_7", "CPS_8", "CPS_9", "CPS_10")], 1, mean, na.rm = TRUE)
t.test(denverjustice$meanCPS ~ Time, data = denverjustice, paired = TRUE)

## 
##  Paired t-test
## 
## data:  denverjustice$meanCPS by Time
## t = 1.24, df = 9, p-value = 0.2462
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.6343  2.1743
## sample estimates:
## mean of the differences 
##                    0.77

barplot(denverjustice$meanCPS)

boxplot(denverjustice$meanCPS ~ Time, data=denverjustice, notch=TRUE, col=(c("gold","blue")), main="Culture of Purpose", xlab="Culture of Purpose")

## Warning: some notches went outside hinges ('box'): maybe set notch=FALSE

Happiness Questions

denverjustice$meanHAPPI <- apply(denverjustice[, c("HAPPI.1_1",  "HAPPI.2_1",  "HAPPI.3_1", "HAPPI.4_1")], 1, mean, na.rm = TRUE)
t.test(denverjustice$meanHAPPI ~ Time, data = denverjustice, paired = TRUE)

## 
##  Paired t-test
## 
## data:  denverjustice$meanHAPPI by Time
## t = -0.4605, df = 9, p-value = 0.6561
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -1.1824  0.7824
## sample estimates:
## mean of the differences 
##                    -0.2

boxplot(denverjustice$meanHAPPI ~ Time, data=denverjustice, notch=TRUE, col=(c("gold","blue")), main="Happiness", xlab="Happiness")

## Warning: some notches went outside hinges ('box'): maybe set notch=FALSE

Grit

denverjustice$meanRES <- apply(denverjustice[, c("RES.1",  "RES.2", "RES.3")], 1, mean, na.rm = TRUE)
t.test(denverjustice$meanRES ~ Time, data = denverjustice)

## 
##  Welch Two Sample t-test
## 
## data:  denverjustice$meanRES by Time
## t = 0.5912, df = 17.75, p-value = 0.5618
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -1.449  2.582
## sample estimates:
## mean in group 1 mean in group 2 
##           8.000           7.433

boxplot(denverjustice$meanRES ~ Time, data=denverjustice, notch=TRUE, col=(c("gold","blue")), main="Grit", xlab="GRIT")

## Warning: some notches went outside hinges ('box'): maybe set notch=FALSE

ASDQII Questions

denverjustice$meanASDQII <- apply(denverjustice[, c("ASDQII_1",  "ASDQII_2",    "ASDQII_3", "ASDQII_4", "ASDQII_5", "ASDQII_6", "ASDQII_7", "ASDQII_8", "ASDQII_9", "ASDQII_10",    "ASDQII_11",    "ASDQII_12",    "ASDQII_13",    "ASDQII_14",    "ASDQII_15",    "ASDQII_16",    "ASDQII_17",    "ASDQII_18",    "ASDQII_19",    "ASDQII_20")], 1, mean, na.rm = TRUE)
t.test(denverjustice$meanASDQII ~ Time, data = denverjustice, paired = TRUE)

## 
##  Paired t-test
## 
## data:  denverjustice$meanASDQII by Time
## t = -1.747, df = 9, p-value = 0.1146
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.9753  0.1253
## sample estimates:
## mean of the differences 
##                  -0.425

boxplot(denverjustice$meanASDQII ~ Time, data=denverjustice, notch=F, col=(c("gold","blue")), main="Self Concept Academic", xlab="ASDQII")