Построение граничной плотности графа с R - PullRequest
Новая
       13

Построение граничной плотности графа с R

1 голос
/ 30 августа 2011

Я строю степень вершины (количество падающих ребер) графа 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")

Этот код отображает этот график рассеяния , enter image description here

, который очень близок к тому, что мне нужно, но имеет несколько проблем:

1) метки x неверны, поскольку они не представляют градусы, а гистограмма нарушается.

2) стержни осей грязные. Как я могу удалить пустые?

3) как построить линию регрессии ? Я пробовал с abline и lm(histdata$mids~histdata$count), но ничего не строится.

Спасибо за любую подсказку!

ОБНОВЛЕНИЕ : этот график, вероятно, совершенно неверный. См http://www.hpl.hp.com/research/idl/papers/ranking/ranking.html

...