9

次のようなRのサンプルコードがあります。

library(igraph)
rm(list=ls())
dat=read.csv(file.choose(),header=TRUE,row.names=1,check.names=T) # read .csv file
m=as.matrix(dat)
net=graph.adjacency(adjmatrix=m,mode="undirected",weighted=TRUE,diag=FALSE)

ここで、次のデータを含む csv ファイルを入力として使用しました。

    23732   23778   23824   23871   58009   58098   58256
23732   0   8   0   1   0   10  0
23778   8   0   1   15  0   1   0
23824   0   1   0   0   0   0   0
23871   1   15  0   0   1   5   0
58009   0   0   0   1   0   7   0
58098   10  1   0   5   7   0   1
58256   0   0   0   0   0   1   0

この後、次のコマンドを使用して重量値を確認しました。

E(net)$weight

予想される出力は次のようになります。

> E(net)$weight
 [1]  8  1 10  1 15  1  1  5  7  1

しかし、私は奇妙な値を取得しています(そして毎回異なります):

> E(net)$weight
 [1] 2.121996e-314 2.121996e-313 1.697597e-313  1.291034e-57 1.273197e-312 5.092790e-313 2.121996e-314 2.121996e-314 6.320627e-316 2.121996e-314 1.273197e-312 2.121996e-313
[13] 8.026755e-316  9.734900e-72 1.273197e-312 8.027076e-316 6.320491e-316 8.190221e-316 5.092790e-313  1.968065e-62 6.358638e-316

どこで何が間違っているのかわかりませんか? 正しい期待される結果を得るために私を助けてください。

ありがとう、ニティン

4

2 に答える 2

2

問題は、行列要素のデータ型が原因のようです。graph.adjacencyタイプの要素が必要ですnumeric。バグかどうかはわかりません。

あなたがした後、

m <- as.matrix(dat)

モードを次のように設定しますnumeric。

mode(m) <- "numeric"

そして、次のようにします。

net <- graph.adjacency(m, mode = "undirected", weighted = TRUE, diag = FALSE)
> E(net)$weight
[1]  8  1 10  1 15  1  1  5  7  1
于 2013-01-15T07:06:33.200 に答える