Generating random inverse gamma vector

I have a problem with generating a random inverse gamma vector. Consider this example:
a = 100 * ones(10, 1);
b = 1 ./ seqa(90, 1, 10);
c = 1 ./ rndGamma(10, 1, a, 1./b);


However, Gauss gives the error message:

G0618 : Third input, shape argument, must be a scalar

Why? According to the documentation, Gauss allows for vector shape and scale parameters.

2 Answers

0

This is an error in the documentation. Currently, rndGamma only supports a scalar for the shape and scale inputs.

link

admin

62

0

Isn't that a bug? Other distributions, such as the beta, accept vector parameter inputs. It would be great to fix it because using a loop is not very efficient.

link

sebszab

24