topic can't assign an Array in a Dirichlet prior, but SAS requires a vector for Dirichlet in Statistical Procedures

can't assign an Array in a Dirichlet prior, but SAS requires a vector for Dirichlet

csfcgua — Mon, 06 Jul 2015 03:13:33 GMT

t Hello,

I put a Model in Proc MCMC, and set one of the parameters as Dirichlet prior:

Proc MCMC Data=test ...... ;

Parms Y X;

Array X[3] x1 x2 x3;

Array P[3] p1 p2 p3;

Model A1 ~ normal( Y + x1, sd=1);

Model A2 ~ normal( Y + x2, sd=1);

Model A3 ~ normal( Y + x3, sd=1);

Prior Y ~ beta(....);

Prior X ~ Dirichlet(P);

Run;

SAS kept telling me that "The name X is already defined as a variable and cannot be redeclared as an array by the ARRAY statement"

But if I delete the Array statement, SAS showed ERROR: The DIRICHLET distribution must be applied to an array.

This is very frustrating. Can somebody please tell me how to fix this problem? Thanks.

Ksharp — Mon, 06 Jul 2015 13:19:40 GMT

Then use another name

Proc MCMC Data=test ...... ;

Parms Y X;

Array XX[3] x1 x2 x3;

Array P[3] p1 p2 p3;

Model A1 ~ normal( Y + x1, sd=1);

Model A2 ~ normal( Y + x2, sd=1);

Model A3 ~ normal( Y + x3, sd=1);

Prior Y ~ beta(....);

Prior X ~ Dirichlet(P);

Run;

csfcgua — Wed, 08 Jul 2015 02:24:46 GMT

Thanks.

But both X and P in dirichlet must be declared as array.

charlie26 — Tue, 24 Jan 2017 15:55:45 GMT

I ran into the same problem. Any solution? Thanks in advance!