T2 Interrater Reliability

library(irr)

library(tidyverse)
library(knitr)
library(stargazer)
library(dplyr)
library(tibble)
library(Hmisc)
library(psych)

#setwd
setwd("C:/Users/C00252837/OneDrive - University of Louisiana Lafayette/PEACE/DATA_T1_T2")

#------------Read in Data
data<-read.csv("t2 behavioral final.csv", header = T)
#names(data)
#names(select(data, 1:20))

FFSF (face-to-face still face) session

Coders: C1(EN), C2(IB), C3(MW), C4(FW)

FFSF1childpos

4 coders

#FFSF1childpos_4 coders
# [51] "FFSF1childpos_C1"                [52] "FFSF1childpos_C2"            [53] "FFSF1childpos_C3"                [54] "FFSF1childpos_C4"
FFSF1childpos_4 <- select(data, 51:54)
names(FFSF1childpos_4)

## [1] "FFSF1childpos_C1" "FFSF1childpos_C2" "FFSF1childpos_C3" "FFSF1childpos_C4"

#unique(FFSF1childpos_4$FFSF1childpos_C1)
#str(FFSF1childpos_4)

#Exclude missing values
FFSF1childpos <- na.omit(FFSF1childpos_4)
#str(FFSF1childpos_4)

#ICC for C1, C2, C3, and C4
icc_result <- icc(FFSF1childpos_4)
icc_value <- icc_result$value
icc_value

## [1] 0.4905009

3 coders

#FFSF1childpos_3 coders
# [51] "FFSF1childpos_C1"                [52] "FFSF1childpos_C2"            [53] "FFSF1childpos_C3"                
FFSF1childpos_3 <- select(data, 51:53)
names(FFSF1childpos_3)

## [1] "FFSF1childpos_C1" "FFSF1childpos_C2" "FFSF1childpos_C3"

#unique(FFSF1childpos_3$FFSF1childpos_C1)
#str(FFSF1childpos_3)

#Exclude missing values
FFSF1childpos_3 <- na.omit(FFSF1childpos_3)
#str(FFSF1childpos_3)

#ICC for C1, C2, and C3
icc_result <- icc(FFSF1childpos_3)
icc_value <- icc_result$value
icc_value

## [1] 0.764638

FFSF1childneg

4 coders

#FFSF1childneg_4 coders
# [55] "FFSF1childneg_C1"                 "FFSF1childneg_C2"              [57] "FFSF1childneg_C3"                 "FFSF1childneg_C4"
FFSF1childneg_4 <- select(data, 55:58)
names(FFSF1childneg_4)

## [1] "FFSF1childneg_C1" "FFSF1childneg_C2" "FFSF1childneg_C3" "FFSF1childneg_C4"

#unique(FFSF1childneg_4$FFSF1childneg_C1)
#str(FFSF1childneg_4)

#Exclude missing values
FFSF1childneg_4 <- na.omit(FFSF1childneg_4)
#str(FFSF1childneg_4)

#ICC for C1, C2, C3, and C4
icc_result <- icc(FFSF1childneg_4)
icc_value <- icc_result$value
icc_value

## [1] 0.6078431

3 coders

#FFSF1childneg_3 coders
# [55] "FFSF1childneg_C1"                 "FFSF1childneg_C2"              [57] "FFSF1childneg_C3"                
FFSF1childneg_3 <- select(data, 55:57)
names(FFSF1childneg_3)

## [1] "FFSF1childneg_C1" "FFSF1childneg_C2" "FFSF1childneg_C3"

#unique(FFSF1childneg_3$FFSF1childneg_C1)
#str(FFSF1childneg_3)

#Exclude missing values
FFSF1childneg_3 <- na.omit(FFSF1childneg_3)
#str(FFSF1childneg_3)

#ICC for C1, C2, and C3
icc_result <- icc(FFSF1childneg_3)
icc_value <- icc_result$value
icc_value

## [1] 0.9220779