Primero los paquetes
library("irr")
library("icc")
library("tidyverse")
Dataframe
df <- read_csv("icc_felipe.csv")
Reformatear de wide a long
long_df1 <- df %>%
select(`1`:`12`) %>%
gather("Grupo_antes", "Valor_antes")
long_df2 <- df %>%
select(`a`:`l`) %>%
gather("Grupo_despues", "Valor_despues")
long_df <- bind_cols(long_df1, long_df2)
rm(df, long_df1, long_df2)
glimpse(long_df)
long_df <- long_df %>%
select(Valor_antes, Valor_despues)
Ahora el ICC icc(ratings, model = c(“oneway”, “twoway”), type = c(“consistency”, “agreement”), unit = c(“single”, “average”), r0 = 0, conf.level = 0.95)
icc(long_df, model = "twoway", type = "agreement")
Single Score Intraclass Correlation
Model: twoway
Type : agreement
Subjects = 156
Raters = 2
ICC(A,1) = 1
F-Test, H0: r0 = 0 ; H1: r0 > 0
F(155,142) = 94504 , p = 2.56e-313
95%-Confidence Interval for ICC Population Values:
1 < ICC < 1
LS0tCnRpdGxlOiAiSUNDIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tClByaW1lcm8gbG9zIHBhcXVldGVzCmBgYHtyIHBhcXVldGVzLCBldmFsPUZBTFNFfQpsaWJyYXJ5KCJpcnIiKQpsaWJyYXJ5KCJpY2MiKQpsaWJyYXJ5KCJ0aWR5dmVyc2UiKQoKYGBgCkRhdGFmcmFtZQpgYGB7ciBkYXRhZnJhbWV9CmRmIDwtIHJlYWRfY3N2KCJpY2NfZmVsaXBlLmNzdiIpCmBgYApSZWZvcm1hdGVhciBkZSB3aWRlIGEgbG9uZwpgYGB7cn0KbG9uZ19kZjEgPC0gZGYgJT4lCiAgc2VsZWN0KGAxYDpgMTJgKSAlPiUgCiAgZ2F0aGVyKCJHcnVwb19hbnRlcyIsICJWYWxvcl9hbnRlcyIpCgpsb25nX2RmMiA8LSBkZiAlPiUKICBzZWxlY3QoYGFgOmBsYCkgJT4lIAogIGdhdGhlcigiR3J1cG9fZGVzcHVlcyIsICJWYWxvcl9kZXNwdWVzIikKCgpsb25nX2RmIDwtIGJpbmRfY29scyhsb25nX2RmMSwgbG9uZ19kZjIpCgpybShkZiwgbG9uZ19kZjEsIGxvbmdfZGYyKQoKZ2xpbXBzZShsb25nX2RmKQoKbG9uZ19kZiA8LSBsb25nX2RmICU+JSAKICBzZWxlY3QoVmFsb3JfYW50ZXMsIFZhbG9yX2Rlc3B1ZXMpCmBgYApBaG9yYSBlbCBJQ0MKaWNjKHJhdGluZ3MsIG1vZGVsID0gYygib25ld2F5IiwgInR3b3dheSIpLCAKICAgIHR5cGUgPSBjKCJjb25zaXN0ZW5jeSIsICJhZ3JlZW1lbnQiKSwgCiAgICB1bml0ID0gYygic2luZ2xlIiwgImF2ZXJhZ2UiKSwgcjAgPSAwLCBjb25mLmxldmVsID0gMC45NSkKCmBgYHtyfQppY2MobG9uZ19kZiwgbW9kZWwgPSAidHdvd2F5IiwgdHlwZSA9ICJhZ3JlZW1lbnQiKQogICAgCmBgYA==