################################################################################
#    **************************** R companion for ************************** 
#
# Kenny, D. A., Kashy, D. A., & Cook, W. L. (2006). Dyadic data analysis. 
# New York: Guilford Press.
#
# written by Avi Kluger: avik@savion.huji.ac.il
#
#                              CHAPTER 3 -- Table 3.3
# Data for Table 3.3 is not provided in the book's site.  Alternative dataset
# is provided here.  These data (Lehmann, 2018) are ratings of humility by 
# listeners and speakers (within dyad variable) in a control group of a larger
# experiment.
################################################################################
rm(list = ls())                               # Clean the Global Environment
cat ("\014")                                  # Clean the R console
if (is.null(dev.list()) == FALSE) dev.off()   # Clean Plots

# Example data DIFFERENT FROM THE BOOK based on Lehmann study of humility

urlRemote   <- "https://raw.githubusercontent.com/"
pathGithub  <- "avi-kluger/RCompanion4DDABook/master/Chapter%203/"
fileName    <- "Lehmann%20Humility%20data.csv"

humility_df <- read.csv(paste0(urlRemote, pathGithub, fileName))

if (!require('tidyverse')) install.packages('tidyverse'); library('tidyverse')
humility_df %>% summarize(m.l = mean(Humility.listener),
                         sd.l = sd(Humility.listener),
                          m.s = mean(Humility.Speaker),
                         sd.s = sd(Humility.Speaker))

# ICC based on Pearson's r

r   <- round(cor(humility_df[, 
                         c("Humility.listener", "Humility.Speaker")]), 2)[1, 2]
r
fit <- (t.test(humility_df$Humility.listener, 
       humility_df$Humility.Speaker, 
       paired=TRUE, 
       conf.level=0.95))
dD  <- (as.numeric(fit$statistic) * sqrt(2))/sqrt(nrow(humility_df))
d   <- dD * sqrt(1 - r)
dD
d