PERMA Positive Emotion

setwd("/Users/levibrackman/Desktop/Adult_study")
data_pre <- read.csv("pre.csv")
data_post <- read.csv("post.csv")

table(data_pre$GROUP, dnn = "group")

## group
##  0  1 
## 35 32

data <- merge(data_pre, data_post, by = "ID", all = TRUE)

library(psych)
library(lattice)
# Only Positive Emotion questions
data$T1PERMA_POS_EMO <- apply(data[, c("PERMA3.x", "PERMA7.x", "PERMA16.x")], 
    1, mean, na.rm = TRUE)
data$T2PERMA_POS_EMO <- apply(data[, c("PERMA3.y", "PERMA7.y", "PERMA16.y")], 
    1, mean, na.rm = TRUE)

plot(data$T1PERMA_POS_EMO, data$T2PERMA_POS_EMO, ylab = "Pre", xlab = "Post", 
    main = "POS EMO")

# pre test plots
bwplot(GROUP.x ~ T1PERMA_POS_EMO, ylab = "GROUP", xlab = "PERMA_POS_EMO", main = "Pre test", 
    data = data)

# post test plots
bwplot(GROUP.x ~ T2PERMA_POS_EMO, ylab = "Group", xlab = "PERMA", main = "Post test", 
    data = data)

# Pre test
t.test(T1PERMA_POS_EMO ~ GROUP.x, data = data)

## 
##  Welch Two Sample t-test
## 
## data:  T1PERMA_POS_EMO by GROUP.x
## t = -0.3279, df = 62.73, p-value = 0.7441
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.9566  0.6870
## sample estimates:
## mean in group 0 mean in group 1 
##           6.771           6.906

t.test(T2PERMA_POS_EMO ~ GROUP.x, data = data)

## 
##  Welch Two Sample t-test
## 
## data:  T2PERMA_POS_EMO by GROUP.x
## t = -1.243, df = 44.15, p-value = 0.2205
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -1.5155  0.3592
## sample estimates:
## mean in group 0 mean in group 1 
##           7.295           7.873

# Ancova, Model for MLQ
PERMA_POS_EMO_ANCOVA <- lm(T2PERMA_POS_EMO ~ as.factor(GROUP.x) + T1PERMA_POS_EMO, 
    data = data)
# check assumptions visually
plot(PERMA_POS_EMO_ANCOVA)

# see results
summary(PERMA_POS_EMO_ANCOVA)

## 
## Call:
## lm(formula = T2PERMA_POS_EMO ~ as.factor(GROUP.x) + T1PERMA_POS_EMO, 
##     data = data)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -4.034 -0.458  0.005  0.434  1.707 
## 
## Coefficients:
##                     Estimate Std. Error t value Pr(>|t|)    
## (Intercept)           1.3925     0.6413    2.17    0.035 *  
## as.factor(GROUP.x)1   0.4466     0.2777    1.61    0.115    
## T1PERMA_POS_EMO       0.8447     0.0879    9.61  2.2e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.945 on 44 degrees of freedom
##   (22 observations deleted due to missingness)
## Multiple R-squared:  0.687,  Adjusted R-squared:  0.673 
## F-statistic: 48.4 on 2 and 44 DF,  p-value: 7.72e-12