Hello everybody, I am a student in econometrics and statistics. With my friend, we have to writre a papre about the ARMA indentification with SAS/IML.
At first, we generate multiple ARMA process samples like that:
proc iml;
/* proc this process X_t=1.3X_t-1 -0.4X_t-2 u_t + 0.8u_t-1 ARMA(1;2)*/
/*N=30*/
phi= {1, -0.7};
theta= {1, 0.3, 0.8};
mu=0;
sigma=1;
seed=-12345;
NumSamples=500;
%let N=30;
Y= j(Numsamples, &N, 0);/* J(nrow,ncol,value):renvoie value 1 par defaut pour tout nrow et ncol)*/
do i= 1 to Numsamples;
ARMA21= armasim(phi, theta, mu, sigma, &N, seed);
Y[i,]= t(ARMA21);
end;
print Y;
SampleID= repeat(T(1:NumSamples),1,&N);
create vrai.ARMA21_30sim var {SampleID Y}; append; close vrai.ARMA21_30sim; QUIT;
But some troubles are coming up: we want to retrieve the sum of squares error to calculate the log likehood, then criterion information but we have no idea to proceed that after many attempts. Indeed, we want to have many samples to test the relevance of criteria information (AIC, BIC, HQ,..).
Then, to make our work easier, we want to know is it possible to create a matrix where we compile when a criteria information let a good estimation or not.
Any suggestions that you might have would be greatly appreciated! Thank you.
