Eu estou tendo problemas com a programação R. De uma matriz, como posso encontrar a mediana de cada linha excluir a diagonal?
Ex: matriz 4x4
0 1 2 3
1 0 1 2
2 1 0 1
3 2 1 0
Eu quero encontrar a mediana de cada linha excluindo a diagonal (neste ex, o diag = 0)
Eu tentei:
diag(A) <- NA
mean(A, na.rm = TRUE) # doesn"t work
apply(A, 1, median) # it works but the calculation including the diagonal.
Respostas:
2 para resposta № 1Tente isto:
A <- matrix(c(0,1,2,3,1,0,1,2,2,1,0,1,3,2,1,0),nrow=4)
sapply(1:4, function(x) median(A[x,-x]))
[1] 2 1 1 2
1 para resposta № 2
Eu não tenho certeza porque você está usando mean()
; por que você não tenta median()
? Aqui vai:
A <- read.table(text="0 1 2 3
1 0 1 2
2 1 0 1
3 2 1 0", header=F)
A <- as.matrix(A)
diag(A) <- NA
A
V1 V2 V3 V4
[1,] NA 1 2 3
[2,] 1 NA 1 2
[3,] 2 1 NA 1
[4,] 3 2 1 NA
myFun <- function(x){ median(x, na.rm=T)}
apply(A, 1, FUN=myFun)
[1] 2 1 1 2