我正在绘制图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")
此代码绘制此scatterplot ,
这与我的需求非常接近,但有一些问题:
1)x标签是错误的,因为它们不代表度数但直方图断裂。
2)轴条很杂乱。如何删除空的?
3)如何绘制回归线?我尝试使用abline
和lm(histdata$mids~histdata$count)
,但没有任何内容被绘制。
感谢任何提示!
更新:这个情节可能是完全错误的。见http://www.hpl.hp.com/research/idl/papers/ranking/ranking.html