Как создать список ребер из матрицы в R?
Отношение выражается в виде матрицы x как это:
A B C D
A 0 2 1 1
B 2 0 1 0
C 1 1 0 1
D 1 0 1 0
Записи относятся к числу соединений, которые они имеют. Может ли кто-нибудь показать мне, как записать его в качестве крайнего списка? Спасибо! Я собираюсь написать это как крайний список, например:
A B
A B
A C
A D
B C
Но прав ли этот список границ, если я хочу создать сетевой сюжет? Спасибо!
3 ответа
С использованием igraph пакет:
x <- matrix(c(0,2,1,1,2,0,1,0,1,1,0,1,1,0,1,0), 4, 4)
rownames(x) <- colnames(x) <- LETTERS[1:4]
library(igraph)
g <- graph.adjacency(x)
get.edgelist(g)
# [,1] [,2]
# [1,] "A" "B"
# [2,] "A" "B"
# [3,] "A" "C"
# [4,] "A" "D"
# [5,] "B" "A"
# [6,] "B" "A"
# [7,] "B" "C"
# [8,] "C" "A"
# [9,] "C" "B"
# [10,] "C" "D"
# [11,] "D" "A"
# [12,] "D" "C"
Я также рекомендую вам потратить некоторое время на чтение igraph документацию по http://igraph.sourceforge.net/index.html поскольку многие из ваших недавних вопросов относятся к простому случаю.
(В качестве бонуса, plot(g) ответит на другой ваш вопрос Как построить отношения в R?)
С помощью melt в reshape2, а затем удалите вес ==0. если нет необходимости печатать вес. просто удали это.
x
sample1 sample2 sample3 sample4
feature1 0 2 1 1
feature2 2 0 1 0
feature3 1 1 0 1
feature4 1 0 1 0
melt(x)
Var1 Var2 value
1 feature1 sample1 0
2 feature2 sample1 2
3 feature3 sample1 1
4 feature4 sample1 1
5 feature1 sample2 2
Попробуй это
M <- matrix( c(0,2,1,1,2,0,1,0,1,1,0,1,1,0,1,0), 4, 4, dimnames=list(c("A","B","C","D"), c("A","B","C","D")))
eList <- NULL
for ( i in 1:nrow(M) ){
for ( j in 1:ncol(M)) {
eList <- c(eList, rep(paste(dimnames(M)[[1]][i], dimnames(M)[[2]][j] ), M[i,j]))
}
}
Выход
> eList
[1] "A B" "A B" "A C" "A D" "B A" "B A" "B C" "C A" "C B" "C D" "D A" "D C"