R/cohen_kappa_mri.R

In rafaelsommer1/emocempr: Functions for analyzing EMOCEMP data

#' Takes a cleaned data frame from MRI study and calculate cohen kappa for selected varibles
#'@param datafile a data frame resulted from clean_brain, clean_orbit or clean_spinal
#'@param varnames a string vector with variable names from the datasets
#'@return a list with kappa statistics for each selected variable
#'@export
#'@import dplyr
#'@importFrom irr kappa2



cohen_kappa_mri <- function (datafile, varnames) {
        #-----------------------------------------------------------------#
        list_k <- list()
        for (i in varnames){
        testdata_BKC <- datafile %>%
                mutate(bkc = factor(datafile[[i]])) %>%
                select(id_participante,avaliador,bkc) %>%
                filter(avaliador == "BKC")
        testdata_RP <- datafile %>%
                mutate(rp = datafile[[i]]) %>%
                select(id_participante,avaliador,rp) %>%
                filter(avaliador == "RP")
        #-----------------------------------------------------------------#
        data4kappa <- merge(testdata_BKC,testdata_RP,by = "id_participante")
        data4kappa <- select(data4kappa, bkc, rp)
        #-----------------------------------------------------------------#
        list_k[[i]] <- kappa2(data4kappa)

        }
        return(list_k)
}