3

我的问题如下：

我正在使用 R SNA 包进行社交网络分析。可以说，我的出发点是具有以下特征的边缘列表。每一行都包含一个公司名称、他们所涉及的项目的 ID 和进一步的特征，比如项目年份。公司可以参与多个项目，一个项目可以由多个公司的合作组成。例子：

Name   Project   Year
AA     1         2003
AB     1         2003
AB     2         2003
AB     3         2004
AC     2         2003
AC     4         2005

对于网络分析，我需要一个所有公司作为行和列标题的邻接矩阵，我构造如下：

grants.edgelist <- read.csv("00-composed.csv", header = TRUE, sep = ";", quote="\"", dec=",", fill = TRUE, comment.char="")

grants.2mode <-  table(grants.edgelist)   # cross tabulate -> 2-mode sociomatrix

grants.adj <- grants.2mode%*%t(grants.2mode)     # Adjacency matrix as product of the 2-mode sociomatrix`

现在我的问题是：我想在邻接矩阵上运行 netlm 回归，在这里我测试给定年份的网络如何解释下一年的网络。但是，因此我想将 grants.edgelist 子集化为（比如说）2003 年和 2005 年的集合。但是，我发现并非所有公司每年都在项目中，因此相应的邻接矩阵具有不同的行和列。

现在我的问题是：如何获得一个包含行和列标题中所有公司的邻接矩阵，但它们的交集设置为零，期望我想观察的年份。我希望我的意思很清楚。

非常感谢您提前。今天这个问题快把我逼疯了！

最好的祝愿

丹尼尔

4

1 回答 1

1

假设同一家公司有可能在多年内从事同一个项目（如果没有，肯定有更简单的解决方案）。一种方法是构造一个 networkDynamic 对象，然后提取您想要的年份并将它们输入netlm.

library(networkDynamic)

# construct example dataset
firmProj <- matrix(
                    c('AA', 1,  2003,
                      'AB', 1, 2003,
                      'AB', 2, 2003,
                      'AB', 3, 2004,
                      'AC', 2, 2003,
                      'AC', 4, 2005),
                       ncol=3,byrow=TRUE)
colnames(firmProj)<-c('Name',   'Project',   'Year')

# make network encompassing all edges
baseNet<-network(firmProj[,1:2])

# get the ids/names of the vertices
ids<-network.vertex.names(baseNet)

# convert original data to a timed edgelist
tel<-cbind(as.numeric(firmProj[,3]),   # convert years to numeric start time
           as.numeric(firmProj[,3])+1, # convert years to numeric end  time
           match(firmProj[,1],ids),    # match label to network id
           match(firmProj[,2],ids))    # match label to network id

# convert to a networkDynamic object
dynFirmProj<-networkDynamic(baseNet, edge.spells = tel)

# bin static networks from the dynamic one, and convert them into list of adjacency matrices
lapply(
       get.networks(dynFirmProj,start=2003, end = 2006, time.increment = 1),
       as.matrix)

[[1]]
   1 2 3 4 AA AB AC
1  0 0 0 0  0  0  0
2  0 0 0 0  0  0  0
3  0 0 0 0  0  0  0
4  0 0 0 0  0  0  0
AA 1 0 0 0  0  0  0
AB 1 1 0 0  0  0  0
AC 0 1 0 0  0  0  0

[[2]]
   1 2 3 4 AA AB AC
1  0 0 0 0  0  0  0
2  0 0 0 0  0  0  0
3  0 0 0 0  0  0  0
4  0 0 0 0  0  0  0
AA 0 0 0 0  0  0  0
AB 0 0 1 0  0  0  0
AC 0 0 0 0  0  0  0

[[3]]
   1 2 3 4 AA AB AC
1  0 0 0 0  0  0  0
2  0 0 0 0  0  0  0
3  0 0 0 0  0  0  0
4  0 0 0 0  0  0  0
AA 0 0 0 0  0  0  0
AB 0 0 0 0  0  0  0
AC 0 0 0 1  0  0  0

但是，我不确定netlm这是否是最好的方法，因为因变量的二分数据“..由于分析的假设而强烈反对”。但也许我不太理解你的问题。

于 2016-07-22T21:39:34.010 回答