setwd("~/Mestrado/Ciencias de Dados em Psicometria/Aula 25092018")
con<-url("http://www.labape.com.br/rprimi/ds/enem_2015.RData") 
library(tidyverse)
## Warning: package 'tidyverse' was built under R version 3.5.1
## -- Attaching packages -------------------------------------------------------------------------- tidyverse 1.2.1 --
## v ggplot2 2.2.1     v purrr   0.2.4
## v tibble  1.4.2     v dplyr   0.7.6
## v tidyr   0.8.0     v stringr 1.3.0
## v readr   1.1.1     v forcats 0.3.0
## Warning: package 'dplyr' was built under R version 3.5.1
## -- Conflicts ----------------------------------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library(readxl)
bd_tcfi <- read_excel("bd_tcfi.xlsx")

bd_tcfi <- bd_tcfi %>% select(c(2, 4, 8:17))

1º Passo

Gather: cria uma coluna chamada “item” e reescreve o banco de dados considerando esta nova coluna

bd_tcfi <- bd_tcfi %>% gather(At3.01:At3.10, key = "item", value = "scr") %>% arrange(ID, RaterN, item)

2º Passo

Spread: Organiza a tabela a partir do nome do avaliador

bd_tcfi <- bd_tcfi %>% spread(key = "RaterN", value = "scr" )

3º Fazer a análise de concordância

Realiza a análise de concordância dos juízes

library(irr)
## Loading required package: lpSolve
agree(bd_tcfi[ , 3:8])
##  Percentage agreement (Tolerance=0)
## 
##  Subjects = 0 
##    Raters = 6 
##   %-agree = NaN
library(psych)
## 
## Attaching package: 'psych'
## The following objects are masked from 'package:ggplot2':
## 
##     %+%, alpha
corr.test(bd_tcfi[, 3:8])
## Warning in sqrt(n - 2): NaNs produzidos
## Warning in corr.test(bd_tcfi[, 3:8]): Number of subjects must be greater
## than 3 to find confidence intervals.
## Warning in sqrt(n[lower.tri(n)] - 3): NaNs produzidos
## Call:corr.test(x = bd_tcfi[, 3:8])
## Correlation matrix 
##           alexandre andre carla everson jurandir marco
## alexandre      1.00    NA    NA    0.46       NA  0.31
## andre            NA  1.00  0.46      NA       NA  0.31
## carla            NA  0.46  1.00      NA       NA  0.47
## everson        0.46    NA    NA    1.00       NA  0.45
## jurandir         NA    NA    NA      NA     1.00  0.81
## marco          0.31  0.31  0.47    0.45     0.81  1.00
## Sample Size 
##           alexandre andre carla everson jurandir marco
## alexandre       662     0     0     661        0   449
## andre             0   707   705       0        0   460
## carla             0   705   713       0        0   461
## everson         661     0     0     679        0   465
## jurandir          0     0     0       0      540   509
## marco           449   460   461     465      509  1628
## Probability values (Entries above the diagonal are adjusted for multiple tests.) 
##           alexandre andre carla everson jurandir marco
## alexandre         0    NA    NA       0       NA     0
## andre            NA     0     0      NA       NA     0
## carla            NA     0     0      NA       NA     0
## everson           0    NA    NA       0       NA     0
## jurandir         NA    NA    NA      NA        0     0
## marco             0     0     0       0        0     0
## 
##  To see confidence intervals of the correlations, print with the short=FALSE option