我有一个邻接矩阵，另一个表示两个节点是否共享属性。将其视为同质测试。我们要测试形成连接的可能性是否取决于两个节点具有共同属性的事实。 现在，使用R和SNA包，我运行了一个关联，并且通过QAP测试来检验其重要性：
g <- array(dim=c(2,nrow(x),nrow(x)))
g[1,,] <- x
g[2,,] <- y
q.12 <- qaptest(g, gcor, reps = 2000, g1=1, g2=2, diag=FALSE)
相关系数为0.7479487，p值为0
QAP Test Results
Estimated p-values:
p(f(perm) >= f(d)): 0
p(f(perm) <= f(d)): 1
然后我将logit拟合到该数据
nl <- netlogit(y, x, mode="digraph", diag=FALSE, nullhyp="qap", reps=2000)
，但其系数不显着。那怎么可能？
Network Logit Model
Coefficients:
Estimate Exp(b) Pr(<=b) Pr(>=b) Pr(>=|b|)
(intercept) -2.940634 5.283224e-02 0.000 1.000 0.00
x1 21.506702 2.188981e+09 0.519 0.481 0.97
Goodness of Fit Statistics:
Null deviance: 17234.41 on 12432 degrees of freedom
Residual deviance: 4617.118 on 12430 degrees of freedom
Chi-Squared test of fit improvement:
12617.29 on 2 degrees of freedom, p-value 0
AIC: 4621.118 BIC: 4635.974
Pseudo-R^2 Measures:
(Dn-Dr)/(Dn-Dr+dfn): 0.5036986
(Dn-Dr)/Dn: 0.7320989