install.packages(c("sna","foreign","chron","igraph"))
library(sna)
load(file="0Mat.rda") dim(Mat) #[1] 1067 1067 と出力されたはず 1067人×1067人の行列 m<- Mat[1:50,1:50] #一部のみとりだす m gplot(m) #一部のみを社会ネットワークとして描画 gplot(m,label=rownames(m)) #一部のみを社会ネットワークとして描画
gplot(Mat) #全員のを描いてみる sr<-rowSums(Mat) #行方向の和 sc<-colSums(Mat) #列方向の和 par(mfrow=c(1,3)) hist(sr) #ヒストグラム hist(sc) plot(sr,sc) #散布図 g<-sr+sc m2<-Mat[g>20,g>20] #20回以上やりとりしている者のみに限定 dim(m2) #[1] 258 258 par(mfrow=c(1,1)) gplot(m2,label=rownames(m2)) #一部のみを社会ネットワークとして描画 gplot(m2,label=rownames(m2),mode="circle") #プロットの方法を変える gplot(m2,label=rownames(m2),mode="eigen") # gplot(m2,label=rownames(m2),mode="kamadakawai") # gplot(m2,label=rownames(m2),mode="mds") # gplot(m2,label=rownames(m2),mode="spring") #
gplot(m2,label=rownames(m2),mode="spring",interactive=T) #インタラクティブにもできる
#とりあえず回帰分析
res<-lm(sr~sc) summary(res)
Residuals:
Min 1Q Median 3Q Max
-153.658 -1.242 0.992 2.070 168.374
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.069812 0.452570 -4.573 5.36e-06 ***
sc 1.077901 0.004673 230.645 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 14.22 on 1065 degrees of freedom
Multiple R-squared: 0.9804, Adjusted R-squared: 0.9804
F-statistic: 5.32e+04 on 1 and 1065 DF, p-value: < 2.2e-16