cna error

Using cna for community network analysis , I am getting error in cmdscale.

Following steps were followed:

pdb<-read.pdb("../b_Firstframe.pdb")
inds<-atom.select(pdb, elety="CA")
cij <-read.table("lmi.txt")
cij<-data.matrix(cij)
net_cutoff0.5 <- cna(cij, cutoff.cij=0.5)
Warning messages:
1: In (function (graph, weights = E(graph)$weight, directed = TRUE, :
At community.c:460 :Membership vector will be selected based on the lowest modularity score.
2: In (function (graph, weights = E(graph)$weight, directed = TRUE, :
At community.c:467 :Modularity calculation with weighted edge betweenness community detection might not make sense -- modularity treats edge weights as similarities while edge betwenness treats them as distances
3: In cna.dccm(cij, cutoff.cij = 0.5) :
The number of communities is larger than the number of unique
'colors' provided as input. Colors will be recycled
plot(net_cutoff0.5, pdb)
Obtaining layout from PDB structure
Error in cmdscale(dist(cent), k = k) : NA values not allowed in 'd'
print(net_cutoff0.5)

Call:
cna.dccm(cij = cij, cutoff.cij = 0.5)

Structure:

NETWORK NODES#: 314 EDGES#: 22310
COMMUNITY NODES#: 119 EDGES#: 1904
attr: network, communities, community.network,
community.cij, cij, call

‌

sessionInfo()
R version 3.2.3 (2015-12-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.5 LTS

locale:
[1] LC_CTYPE=en_IN LC_NUMERIC=C LC_TIME=en_IN
[4] LC_COLLATE=en_IN LC_MONETARY=en_IN LC_MESSAGES=en_IN
[7] LC_PAPER=en_IN LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C LC_MEASUREMENT=en_IN LC_IDENTIFICATION=C

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] bio3d_2.3-4

loaded via a namespace (and not attached):
[1] magrittr_1.5 parallel_3.2.3 igraph_1.2.2 Rcpp_0.12.19
[5] grid_3.2.3 pkgconfig_2.0.2

‌

‌

Comments (2)