Le forum du Master ESA économétrie et statistique appliquée - Université d'Orléans

Vous n'êtes pas identifié.

Annonce

Vous êtes sur le forum du master ESA !

Le site du master ESA - description de la formation, notes de cours, contacts... vient de déménager !!!

Venez visiter notre nouveau site : www.master-esa.fr

#1 12-03-2007 23:18:04

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Pour les M1 : maximisation de la log vraisemblance (cas general)

solution du problème de vendredi soir : (comme ca on pourra passer à autre chose cette semaine! )
pour ceux qui veulent chercher encore la solution de leur côté, qu'ils ne lisent pas la suite de ce message !! ( je ne veux pas me prendre des coups de trousses dans la geule en cours ! )
S'il y a d'autres methodes je prend ! big_smile et faites pas les radins! vous faites tourner un peu!!
bon let's go :
libname db 'E:\Cours\EVQ\td';
%let listvar= x1......x50; ** c'est tout ce que l'utilisateur doit entrer pour que le prog lui sorte les rersultats, la dernière var a mettre est la var explicative; au passage pas besoin de crée un col pour la constante;
proc iml;
use input;
read all var{&listvar} into x;

nobs=nrow(x);
nbrcol=ncol(x);
start loglikj(theta) global(x);
nobs=nrow(x);
nbrcol=ncol(x);
ll=0;
do i=1 to nobs;
bidul=0;
do j=1 to nbrcol-1;
bidulj=-theta[j+1]*x[i,j];
bidul=bidul+bidulj;
end;
score=bidul-theta[1];
pi=1/(1+exp(score));
lli=x[i,nbrcol]*log(pi)+(1-x[i,nbrcol])*log(1-pi);
ll=ll+lli;
end;
return(ll);
finish loglikj;

options = {1,2};
parmin=j(1,nbrcol,0);
call nlpnra(rc,thetahat,"loglikj",parmin,options);
call nlpfdd(f,g,hess,"loglikj",thetahat);
param=t(thetahat);                             
v=inv(hess);                                   
stderrors=sqrt(abs(vecdiag(v)));               
WaldChiSquare=(param/stderrors)##2;             
Pvalue=1-cdf('CHISQUARE',WaldChiSquare,1);     
print,"Résultats de l'estimation par Maximum de Vraisemblance",,
param stderrors WaldChiSquare Pvalue;
quit;

S'il y a des trucs que vous voulez que j'exlique, demandez !!

Hors ligne

 

#2 13-03-2007 07:55:40

gang
Member
Lieu: Orleans
Date d'inscription: 05-03-2007
Messages: 324

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

c'est quoi ts ces trucs infectes..... ggrrrrrrrrrrrrrrr
je me demande des fois si on bossait tt ca l'année passée lol!!!!!

Hors ligne

 

#3 13-03-2007 14:31:11

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

ah bon....? on-m'aurait menti ? yikes
ca sert à qlq chose au moins ?

Hors ligne

 

#4 13-03-2007 14:45:08

gang
Member
Lieu: Orleans
Date d'inscription: 05-03-2007
Messages: 324

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

bien sur rassure toi ca sert à bcp de choses et bcp trop de choses heureusement!!!!
mais bon le forum est relativement calme aujourd'hui.... :-(
ok moi je me retrouve au boulot je reviens faire un tour à ma prochaine pause café (et oui je suis tombé ds une boite ou la pause café est sacrée lol.....)

Hors ligne

 

#5 13-03-2007 15:06:03

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

lol t'en a de la pause café dis-donc !!

Hors ligne

 

#6 13-03-2007 15:09:17

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

c vrai que c trop calme d'habitude ce forum ! on est plus 3 survivants sur toutes les promos ! tiens je suis plus new member ! je suis member ! awesome !!!

Hors ligne

 

#7 15-03-2007 02:36:43

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

On en découvre des trouc shous ShASh! bon y'en qui ont pas reussi a faire tourner mon prog ! La g la version finale ! alors la si vous arrivez pas a la faire tourner je me tire un balle mad!!!!!
libname db 'E:\Cours\EVQ\td'; **n'oubliez pas de changer ca svp !! ;
%macro maxvrai(base=,listvar=);
proc iml;
use &base; ** encore plus subtile !! ;
read all var{&listvar} into x;
nobs=nrow(x);
nbrcol=ncol(x);
start loglikj(theta) global(x);
nobs=nrow(x);
nbrcol=ncol(x);
ll=0;
do i=1 to nobs;
bidul=0;
do j=1 to nbrcol-1;
bidulj=-theta[j+1]*x[i,j];
bidul=bidul+bidulj;
end;
score=bidul-theta[1];
pi=1/(1+exp(score));
lli=x[i,nbrcol]*log(pi)+(1-x[i,nbrcol])*log(1-pi);
ll=ll+lli;
end;
return(ll);
finish loglikj;
options = {1,2};
parmin=j(1,nbrcol,0);
call nlpnra(rc,thetahat,"loglikj",parmin,options);
call nlpfdd(f,g,hess,"loglikj",thetahat);
param=t(thetahat);                             
v=inv(hess);                                   
stderrors=sqrt(abs(vecdiag(v)));               
WaldChiSquare=(param/stderrors)##2;             
Pvalue=1-cdf('CHISQUARE',WaldChiSquare,1);     
print,"Résultats de l'estimation par Maximum de Vraisemblance",,
param stderrors WaldChiSquare Pvalue;
quit;
%mend;
%maxvrai(base=db.base,listvar=duree montcred ancienn incident); * C la ou vous rentrez votre base et vos variable !! bien sur l'expliquée en dernier pour l'instant c la seule alternative;

Hors ligne

 

#8 15-03-2007 15:09:53

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

c juste un test !

Code:

libname db 'E:\Cours\EVQ\td'; [color=green]**n'oubliez pas de changer ca svp !! ;[/color]
%macro maxvrai(base=,listvar=);
proc iml;
use &base; [color=green]** encore plus subtile !! ;[/color]
read all var{&listvar} into x;
nobs=nrow(x);
nbrcol=ncol(x);
start loglikj(theta) global(x);
nobs=nrow(x);
nbrcol=ncol(x);
ll=0;
do i=1 to nobs;
bidul=0;
do j=1 to nbrcol-1;
bidulj=-theta[j+1]*x[i,j];
bidul=bidul+bidulj;
end;
score=bidul-theta[1];
pi=1/(1+exp(score));
lli=x[i,nbrcol]*log(pi)+(1-x[i,nbrcol])*log(1-pi);
ll=ll+lli;
end;
return(ll);
finish loglikj;
options = {1,2};
parmin=j(1,nbrcol,0);
call nlpnra(rc,thetahat,"loglikj",parmin,options);
call nlpfdd(f,g,hess,"loglikj",thetahat);
param=t(thetahat);                             
v=inv(hess);                                   
stderrors=sqrt(abs(vecdiag(v)));               
WaldChiSquare=(param/stderrors)##2;             
Pvalue=1-cdf('CHISQUARE',WaldChiSquare,1);     
print,"Résultats de l'estimation par Maximum de Vraisemblance",,
param stderrors WaldChiSquare Pvalue;
quit;
%mend;
%maxvrai(base=db.base,listvar=duree montcred ancienn incident); [color=green]* C la ou vous rentrez votre base et vos variable !! 
                                                                      bien sur l'expliquée en dernier pour l'instant c la seule 
                                                                        alternative;[/color]

Hors ligne

 

#9 15-03-2007 15:12:01

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

arf ! elles passent pas les couleurs ! c embetant ca !

Hors ligne

 

#10 15-03-2007 16:19:02

esa_sr
Administrator
Date d'inscription: 21-02-2007
Messages: 5898
Site web

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

ben oui, il ne lit pas les BBcode de couleur à l'intérieur des du BBcode code...
on ne peut pas tout avoir !
cordialement
SR

Hors ligne

 

#11 18-09-2007 14:08:52

rim
membre extérieur
Date d'inscription: 18-09-2007
Messages: 1

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

bonjour à tous je suis une etudiante en master prévision et j ai un problème de modèle
application du modèle GARCH sur le logiciel GAUSS
prière de m aider
j ensuis reconnaissante
merci

Hors ligne

 

#12 19-09-2007 14:26:06

esa_sr
Administrator
Date d'inscription: 21-02-2007
Messages: 5898
Site web

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

sous gauss ? vous n'aurez pas beaucoup d'aide par ici
(on ne pratique pas)

bonne chance

SR

Hors ligne

 

#13 16-09-2010 00:41:52

mirra
membre extérieur
Date d'inscription: 15-09-2010
Messages: 1

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

salut c'est mirra je voulai que vous m'aidiez à trouver un théme d'ingeniorat ,je voulai combiner
l'estimation, la recherche opérationnelle et la prévision dans le domaine de la télécommunication ou du petrol..... 
je serai reconnaissante merci

Hors ligne

 

#14 16-09-2010 07:25:14

chiraz
Member
Date d'inscription: 18-12-2008
Messages: 104

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

recherche opérationnelle+estimation +prévision=réseau de neurones

Hors ligne

 

#15 16-09-2010 15:01:56

sspdiddy
Member
Date d'inscription: 29-04-2009
Messages: 146

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

tu aimes les reseaux de neurones chiraz smile

Hors ligne

 

#16 16-09-2010 15:56:01

chiraz
Member
Date d'inscription: 18-12-2008
Messages: 104

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

non en ce moment j'arrive à les détester

Hors ligne

 

#17 16-09-2010 15:57:37

chiraz
Member
Date d'inscription: 18-12-2008
Messages: 104

Re: Pour les M1 : maximisation de la log vraisemblance (cas general)

mais si tu dit à un cancérologue que tu as mal à la tête il va penser à un cancer du cerveau

Hors ligne

 

Pied de page des forums

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson

[ Generated in 0.027 seconds, 6 queries executed ]