swh:1:snp:e94096527620720ab530102fd039e636a91a990c
Tip revision: dcd2fa6
covFunction.R
``````##input: X is a dataframe with
##	1st column as a vector of network id and is named 'id'
##	2nd column has names of the node

if(all(dimnames(X)[[2]] != 'id')){stop('network id is not provided. Use name id to name the first column of X')}
if(all(dimnames(X)[[2]] != 'Node')){stop('Nodes are not provided. Use nameNode for the second column')}

netlist = unique(X\$id)
XX = list()
PP = dim(X)[2] - 2
for(ii in 1:length(netlist)){
XX[[ii]] = array(0, dim = c(nnodes[ii],nnodes[ii],PP))
Xsm = X[which(X\$id == ii),]
KK = unique(Xsm\$Node)
if(length(KK) != nnodes[ii])(warning('Missing data replaced by 0'))
dimnames(XX[[ii]])[[1]] = dimnames(XX[[ii]])[[2]] = nodenames[[ii]]
x1 = which(dimnames(Xsm)[[2]] == 'id')
x2 = which(dimnames(Xsm)[[2]] == 'Node')
for(kk in 1:length(KK)){
k1 = which(nodenames[[ii]] == KK[kk])
XX[[ii]][,k1,] = sapply(1:PP,function(x)rep(Xsm[which(Xsm\$Node==KK[kk]),-c(x1,x2)][x],nnodes[ii]))
}
}
return(XX)
}

getSenderCov = 	function(X,nnodes,nodenames){
if(all(dimnames(X)[[2]] != 'id')){stop('network id is not provided. Use name id to name the first column of X')}
if(all(dimnames(X)[[2]] != 'Node')){stop('Nodes are not provided. Use nameNode for the second column')}

netlist = unique(X\$id)
XX = list()
PP = dim(X)[2] - 2
for(ii in 1:length(netlist)){
XX[[ii]] = array(0, dim = c(nnodes[ii],nnodes[ii],PP))
Xsm = X[which(X\$id == ii),]
KK = unique(Xsm\$Node)
if(length(KK) != nnodes[ii])(warning('Missing data replaced by 0'))
dimnames(XX[[ii]])[[1]] = dimnames(XX[[ii]])[[2]] = nodenames[[ii]]
x1 = which(dimnames(Xsm)[[2]] == 'id')
x2 = which(dimnames(Xsm)[[2]] == 'Node')
for(kk in 1:length(KK)){
k1 = which(nodenames[[ii]] == KK[kk])
XX[[ii]][k1,,] = sapply(1:PP,function(x)rep(Xsm[which(Xsm\$Node==KK[kk]),-c(x1,x2)][x],nnodes[ii]))
}
}
return(XX)
}

getEdgeCov = function(X,nnodes,nodenames){
if(all(dimnames(X)[[2]] != 'id')){stop('network id is not provided. Use name id to name the first column of X')}
if(all(dimnames(X)[[2]] != 'Sender')){stop('Sender nodes are not provided. Use name Sender for the second column')}
netlist = unique(X\$id)
XX = list()
PP = dim(X)[2] - 3
for(ii in 1:length(netlist)){
XX[[ii]] = array(0, dim = c(nnodes[ii],nnodes[ii],PP))
Xsm = X[which(X\$id == ii),]
if(length(KK) != nnodes[ii])(warning('Missing data replaced by 0'))
dimnames(XX[[ii]])[[1]] = dimnames(XX[[ii]])[[2]] = nodenames[[ii]]
KK1 = Xsm\$Sender
x1 = which(dimnames(Xsm)[[2]] == 'id')
x2 = which(dimnames(Xsm)[[2]] == 'Sender')
for(kk in 1:nrow(Xsm)){
k1 = which(nodenames[[ii]] == KK1[kk])
k2 = which(nodenames[[ii]] == KK2[kk])
XX[[ii]][k1,k2,] = as.numeric(Xsm[kk,-c(x1,x2,x3)])
}
}
return(XX)
}

``````