Tenho um grafo
e calculei a distribuição dos graus e o grau como se segue:
dd <- degree_distribution(graph)
d <- degree(graph)
A partir disso, estimei a Power Law
, para ver se a minha distribuição segue a "lei de Potência" :
degree = 1:max(d)
probability = dd[-1]
# Exclui zeros, pois log de 0 é infinito!
nonzero.position = which(probability != 0)
probability = probability[nonzero.position]
degree = degree[nonzero.position]
reg = lm(log(probability) ~ log(degree))
cozf = coef(reg)
#Estima a power law com base nos valores
power.law.fit = function(x) exp(cozf[[1]] + cozf[[2]] * log(x))
A partir disso, plotei os pontos e a power law
usando o ggplot2
.
Resultando na imagem a seguir:
df <- data.frame(x = degree, y = probability)
print(
ggplot(df, aes(x,y,colour="Distribuição"))+
geom_point(shape = 4) +
stat_function(fun = power.law.fit, geom = "line", aes(colour="Power Law"))+
labs(title = "Grafo", subtitle = "Distribuição dos Graus",
x="K", y="P(k)", colour="Legenda")+
scale_color_brewer(palette="Dark2")
)
Como é possível observar, a minha distribuição não segue a Power Law
! Gostaria de estimar a distribuição de Poisson
e plotar no mesmo gráfico.
Mesmo não tendo certeza que a minha distribuição não segue (ou segue) a Poisson
, gostaria de plotar junto com a Power Law
. Não tenho a mínima ideia de como estimar essa distribuição (Poisson
) a partir dos dados, e calcular o grau médio.
Alguém poderia me ajudar? Grato.
- O grafo usado para calcular a distribuição e o grau é muito grande (700 mil vértices), por isso não coloquei os dados dos grafos. A explicação da resposta pode ser baseada em um grafo qualquer.
r
dizR é um ambiente e linguagem de programação de código aberto para computação estatística...
r
, que é uma linguagem de programação voltada para estatística, ainda que a pergunta seja mais voltada para estatística. Na maioria das faculdades de programação nós temos que aprenderestatística
para poder lidar com certos casos envolvendo programação. Eu entendi o que você disse, mas mesmo assim não acredito que fuja do escopo. Agora sobre o Cross Validated sim, realmente seria melhor respondida lá, acho isso poderia ter sido comentado de primeira.