pacman::p_load(tidyverse, ggplot2, irr, readxl, GGally, dplyr)
Import Data
prepost <- read_excel("~/Library/CloudStorage/OneDrive-UniversityofEdinburgh/Global Academy/PhD/Survey Development/StatementIOR_240520_Agreement CN LJ Rob.xlsx",
sheet = "Pre vs Post")
#View(prepost)
domaindiff <- prepost[, c("OriginalDomain", "FinalDomain")] #need IOR on this df
combdiff <- prepost[, c("OriginalCOMB", "FinalCOMB")] #need IOR on this df
domaindiffIOR <- kappam.fleiss(domaindiff)
print(domaindiffIOR)
## Fleiss' Kappa for m Raters
##
## Subjects = 78
## Raters = 2
## Kappa = 0.856
##
## z = 21.9
## p-value = 0
combdiffIOR <- kappam.fleiss(combdiff)
print(combdiffIOR)
## Fleiss' Kappa for m Raters
##
## Subjects = 78
## Raters = 2
## Kappa = 0.88
##
## z = 15.2
## p-value = 0
Compare Susan’s results with the final COM-B
combsusan <- prepost[, c("SusanCOMB", "FinalCOMB")] #need IOR on this df
combsusanIOR <- kappam.fleiss(combsusan)
print(combsusanIOR)
## Fleiss' Kappa for m Raters
##
## Subjects = 77
## Raters = 2
## Kappa = 0.523
##
## z = 9.71
## p-value = 0