グラフの頂点度 (入射エッジの数) をプロットしていますg。
deg <- degree(g, v=V(g), mode = c("in"), loops = TRUE)
histdata <- hist( deg, breaks=1000, plot=FALSE )
plot(histdata$count, log="xy", type="p", col="blue", bg = "blue", pch=20,
xlim=c(1,max(deg)),
ylim=c(1,max(histdata$count)),
ylab="Frequency", xlab="Degree")
このコードは、この散布図をプロットします。
これは私が必要とするものに非常に近いですが、いくつかの問題があります:
1) x ラベルは度数を表していないが、ヒストグラムが壊れているため、間違っています。
2) 軸棒が乱雑です。空のものを削除するにはどうすればよいですか?
3)回帰直線をプロットするにはどうすればよいですか? 私は試してみましablineたlm(histdata$mids~histdata$count)が、何もプロットされません。
ヒントをありがとう！
更新: このプロットはおそらく間違っています。http://www.hpl.hp.com/research/idl/papers/ranking/ranking.htmlを参照してください。