1 Star 0 Fork 7

夏汉林/nwcommands

forked from 连享会/nwcommands 
加入 Gitee
与超过 1200万 开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)
免费加入
文件
该仓库未声明开源许可证文件(LICENSE),使用请关注具体项目描述及其代码上游依赖。
克隆/下载
nwkatz.ado 1.34 KB
一键复制 编辑 原始数据 按行查看 历史
ThomasGrund 提交于 2015-09-18 19:13 . v1.5.1
capture program drop nwkatz
program nwkatz
version 9
syntax [anything(name=netname)] , alpha(real) [ GENerate(string) unconnected(integer 0)]
_nwsyntax `netname', max(9999)
_nwsetobs
if `networks' > 1 {
local k = 1
}
if "`generate'" == "" {
local generate = "_katz"
}
local generate_all ""
qui foreach netname_temp in `netname' {
nwname `netname_temp'
local directed = r(directed)
mata: katz = J(`r(nodes)', `r(nodes)', `alpha')
tempname geo dist
nwgeodesic `netname_temp', nosym name(`geo') unconnected(`unconnected')
nwtomata `geo', mat(`dist')
mata: katz = katz :^ `dist'
mata: katz_out = (colsum(katz))'
mata: katz_in = (rowsum(katz))
if "`directed'" == "true"{
capture drop `generate'_out`k'
capture drop `generate'_in`k'
nwtostata, mat(katz_out) gen(`generate'_out`k')
nwtostata, mat(katz_in) gen(`generate'_in`k')
local generate_all "`generate_all' `generate'_out`k' `generate'_in`k'"
}
else {
capture drop `generate'`k'
nwtostata, mat(katz_out) gen(`generate'`k')
local generate_all "`generate_all' `generate'`k'"
}
capture nwdrop `geo'
local k = `k' + 1
}
mata: st_rclear()
di "{hline 40}"
di "{txt} Network name: {res}`netname'"
di "{hline 40}"
di "{txt} Katz centrality"
sum `generate_all'
end
*! v1.5.0 __ 17 Sep 2015 __ 13:09:53
*! v1.5.1 __ 17 Sep 2015 __ 14:54:23
Loading...
马建仓 AI 助手
尝试更多
代码解读
代码找茬
代码优化
1
https://gitee.com/xiahanlin/nwcommands.git
git@gitee.com:xiahanlin/nwcommands.git
xiahanlin
nwcommands
nwcommands
master

搜索帮助