https://github.com/cran/HLSM
Tip revision: c9af23eaea43a10ec3a4129404c6df71182ea949 authored by Samrachana Adhikari on 18 June 2014, 00:00:00 UTC
version 0.4
version 0.4
Tip revision: c9af23e
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
getReceiverCov = 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)
}
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')}
if(all(dimnames(X)[[2]] != 'Receiver')){stop('Receiver nodes are not provided. Use name Receiver for the third 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),]
KK = unique(c(Xsm$Receiver,Xsm$Sender))
if(length(KK) != nnodes[ii])(warning('Missing data replaced by 0'))
dimnames(XX[[ii]])[[1]] = dimnames(XX[[ii]])[[2]] = nodenames[[ii]]
KK1 = Xsm$Sender
KK2 = Xsm$Receiver
x1 = which(dimnames(Xsm)[[2]] == 'id')
x2 = which(dimnames(Xsm)[[2]] == 'Sender')
x3 = which(dimnames(Xsm)[[2]] == 'Receiver')
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)
}