Hello everybody!
In this my first message I feel I must first congratulate (and thank)
Allin and Riccardo for the support to this wonderful adventure that
Gretl has become; and all the rest of 'open hands' that gives breath to
the sail of this swift boat.
I install Gretl 1.9.12 in a desktop (Fujitsu Siemens, Win XPPro 32-bits
SP3, Intel Pentium 4550, 1 Gb RAM) and in a laptop (ASUS, Win XPPro
32-bits SP3, Intel Mobile Core 2 2Ghz, 2 Gb RAM) and couldn't replicate
Manny´s problem nor Hèlio Guillerme's.
I expect not to add noise in this search for a solution to Manny´s
question because my inexperience (or ignorance), but I have also a crash
with my desktop Windows XP trying to 'pseudo-simulate' a normal
bivariate distribution with 3000 observations and to program Mardia´s test.
When I run this script in my desktop for 200.000 loops (modifying 3th
line), it always crashes. I suppose that´s a problem of memory enough
cause in laptop it runs ok 1.000.000loops, and Manny´s question seems a
problem of graphics specification but Gretl never crashes until I
install all the Function Files from server. The installation of some
function file could affect the graphics behaviour of Model menu option?
Sorry if it´sa stupid idea.
The script was:
nulldata 3000
ProbMardOm=0
loop i=1..100000 --quiet
matrix Matriz1=mnormal(3000,2)
if det(mcorr(Matriz1))>=0.99999999999999
Datos=(Matriz1-ones(3000,1)*meanc(Matriz1))*(inv(mcov(Matriz1).*I(2))).^0.5
series X=Datos[,1]
series err=Datos[,2]
matrix Chol=cholesky(inv(mcov(Datos)))
A=vec(Datos*Chol*Chol'*Datos').^3
Asimetr=A'*ones(9000000,1)/(3000^2)
K=diag(Datos*Chol*Chol'*Datos').^2
Kurtose=K'*ones(3000,1)/3000
MardiaOm=(3000/6)*Asimetr+3000*((Kurtose-8)/8)^2
Prob=pvalue(X,5,MardiaOm)
if Prob>ProbMardOm
ProbMardOm=Prob
X2=X
err2=err
endif
endif
endloop
print ProbMardOm
corr X2 err2
... but please don´t waste your time in solve my problem. Don´t
worthwhile. And yes!, you can laugh if you want, :-D :-D :-D (suppose it
could be done in a more efficient manner)... Maybe if I learn someday
how hansl can be used and their commands power.
Many thanks Allin and 'Jack'
Juan C. Estévez
Universidade de Santiago de Compostela (Spain)
P.D. Oh! Can I contribute to Gretl translation? Who can I send a
translation when I finished it?
I have 1/3 of gretl.po translated to galician (the 3th spanish
'cooficial' language like basque and catalan).
Sorry for my English, I guarantee my galician is a little better. :-D