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
|
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.000 loops, 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´s a 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 |
