我正在绘制图形的顶点度数(入射边数)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