/* je génère une table de données*/

data ma_base ;
format date DDMMYY10. ;
date="2MAR1999"d ;
do i=1 to 1000 ;
date=date+1 ;
sp500=rannor(-1) ;
output ;
end ;
drop i ;
run;

/*Calcul du rendement*/

data ma_base ; set ma_base ;

annee=year(date) ; /*création de la colone année*/
mois=month(date) ; /*création mois, mais ça vous l'avez déja*/

sp500_1=lag1(sp500) ; /*serie retardée*/

rdt=(sp500-sp500_1)/sp500_1 ; /*calcul du rendrement (si mes souvenirs sont bons...)*/ 

drop sp500_1 ;

run;

/*Calcul du rendement moyen par mois annee*/

proc sort data=ma_base ;by annee mois;run;
proc means data=ma_base noprint ;
by annee mois ;
var rdt ;
output out=moy_rdt_y_m(drop=_type_ _freq_) mean=rdt_m ;
run;

/*petit graphique*/

data moy_rdt_y_m ; set moy_rdt_y_m ;
if mois < 10 then do ;
time=compress(annee!!"/"!!"0"!!mois) ;
end ;
else do ;
time=compress(annee!!"/"!!mois) ;
end ;
run;

ods html style=journal ;
goption reset=all ;
goption device=javaimg ;
symbol i=join v=dote ;
axis1 label=("Rendement moyen") ;
proc gplot data=moy_rdt_y_m ;
plot rdt_m*time / vaxis=axis1 ;
run;
quit ;
ods html close ;