Гамма-случайные числа, MATLAB и Python

У меня есть эта строка кода в программе MATLAB:

x(:,i) = gamrnd(a(i),1,dim,1)

Мне было интересно, как я могу написать эту же строку на Python. Я думаю, что эквивалентное утверждение:

gamma.rvs(a, size=1000) 

Однако это продолжает давать мне ошибку индекса.

Вот мой полный код для этой части:

x = np.array([])
for i in range(N-1):
    # generates dim random variables
    x[:, i] = gamma.rvs(a[i], dim-1)    # generates dim random variables
                                            # with gamma distribution

Спасибо за помощь!