그래프의도 분포에서 샘플링하기

Question

나는 간단하고 어리 석었던 파이썬 문제가 있습니다. 그래프가 주어진다면 그래프 분포도와 같은 확률 변수로부터 샘플을 추출하려고합니다.

이것은 매우 간단해야합니다. 그러나 어쨌든 나는 이것을 계속 엉망으로 다루고있다. 내 코드는 다음과 같습니다

import numpy as np 
import scipy as sp 
import graph_tool.all as gt 

G = gt.random_graph(500, deg_sampler=lambda: np.random.poisson(1), directed=False) 
deg = gt.vertex_hist(G,"total",float_count=False) 

# Extract counts and values 
count = list(deg[0]) 
value = list(deg[1]) 

# Generate vector of probabilities for each node 
p = [float(x)/sum(count) for x in count] 

# Load into a random variable for sampling 
x = sp.stats.rv_discrete(values=(value,p)) 
print x.rvs(1) 

을하지만,이 오류를 반환이 실행시 :

Traceback (most recent call last): 
    File "temp.py", line 16, in <module> 
    x = sp.stats.rv_discrete(values=(value,p)) 
    File "/usr/lib/python2.7/dist-packages/scipy/stats/distributions.py", line 5637, in __init__ 
    self.pk = take(ravel(self.pk),indx, 0) 
    File "/usr/lib/python2.7/dist-packages/numpy/core/fromnumeric.py", line 103, in take 
    return take(indices, axis, out, mode) 
IndexError: index out of range for array 

나는이 왜 모르겠어요. 코드의 경우 내가 대신 쓰기 위 :

x = sp.stats.rv_discrete(values=(range(len(count)),p)) 

그런 다음 코드는 잘 실행하지만 이상한 결과를 제공 -이 분포를 지정한 명확하게 방법을 "0"의 값이 가장되어야한다 공유지. 그러나이 코드는 높은 확률로 "1"을 제공하고 "0"을 반환하지 않으므로 어떤 식 으로든 어떻게 든 이동합니다.

여기에서 무슨 일이 일어나는지 명확히 알 수 있습니까? 어떤 도움이라도 대단히 감사하겠습니다!

Answer 1

x.rvs()에 대한 첫 번째 인수는 loc arg 일 것입니다. x.rvs(1)으로 전화하여 loc=1을 작성하면 모든 값에 1이 추가됩니다.

# Extract counts and values 
count = list(deg[0]) 
value = list(deg[1]) 

# Generate vector of probabilities for each node 
p = [float(x)/sum(count) for x in count] 

으로 :

이

count, value = deg  # automatically unpacks along first axis 
p = count.astype(float)/count.sum() # count is an array, so you can divide all elements at once