metafor.R

p.disagreement <- function(data) {
  n <- ncol(data)
  m <- matrix(0, n, n)
  rownames(m) <- names(data)
  colnames(m) <- names(data)
  for (i in 1:n) {
    for (j in i:n) {
      m[i, j] <- mean(data[, i] != data[, j])
    }
  }
  m[lower.tri(m)] <- t(m)[lower.tri(m)]
  as.dist(m)
}
metafor <- read.delim("~/Tsvetovye_metafory_gumanitariev.csv")
mat <- p.disagreement(metafor)  #Матрица процент несогласия
fit <- hclust(mat, method = "ward")
plot(fit, hang = -1)