1

我在使用 R 编程时遇到问题。从矩阵中,如何找到每行的中位数排除对角线?

例如:矩阵 4x4

0 1 2 3 

1 0 1 2

2 1 0 1

3 2 1 0

我想找到每行的中位数,不包括对角线(在这个例子中,diag=0)

我试过了:

diag(A) <- NA
mean(A, na.rm = TRUE) # doesn't work

apply(A, 1, median) # it works but the calculation including the diagonal. 
4

2 回答 2

2

尝试这个:

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
于 2013-08-14T23:08:49.797 回答
1

我不确定你为什么使用mean(); 你为什么不试试median()?开始:

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
于 2013-08-15T01:33:13.443 回答