*        VRAI.PRG : CALCUL D'UNE FONCTION DE VRAISEMBLANCE SIMPLE
*
ALL 300
COMPUTE NBEG=1 , NEND=300 , MOY=2 , VARI=1.0
*
*        SIMULATION DES DONNÉES
*        Y SUIT N(MOY,VARI)
*
CLEAR Y
EQUATION 1 Y
# CONSTANT
ASSOCIATE(VARIANCE=VARI) 1
# MOY
SIMULATE 1 NEND-1 NBEG
# 1 Y NBEG
*
LINREG(EQUATION=1) Y NBEG NEND
*
*        RECHERCHE DE LA VRAISEMBLANCE MAXIMALE PAR BALAYAGE
*
COMPUTE %BETA(1)=0.6
ZERO FUNCTION; ZERO PARA;  *  INTIALISATION DE DEUX SERIES
*
DO I=1,12
   COMPUTE %BETA(1)=%BETA(1)+0.2
   COMPUTE FV=0.0
   LINREG(EQUATION=1,CREATE,NOPRINT) Y NBEG NEND
   COMPUTE FV=-(NEND/2)*LOG(2*%PI)-(NEND/2)*LOG(VARI)-(1/(2*VARI))*%RSS
   DISPLAY %BETA(1) FV
   COMPUTE FUNCTION(I) = FV
   COMPUTE PARA(I) = %BETA(1)
END DO I
*
OPEN PLOT C:\T837\GRA\VRAI.RGF
SCATTER(HEADER='FONCTION DE VRAISEMBLANCE',STYLE=LINES,$
        VLABEL='VRAISEMBLANCE log L',HLABEL='BETA') 1
# PARA FUNCTION 1 12
END