Discussion :

[marlene2110]

bonjour,
étant habituée au logiciel R je ne sais meme pas faire un calcul tout simple sous SAS...
J'ai récupéré les réisdus de mon modèle et j'aimerai ensuite calculer le RMSE (racine de la moyenne des résidus carré) mais je sais pas comment l 'écrire...

voici le code que j'ai commencé a écrire, si qn pouvait m'éclairer...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
proc glm data=data1;
class group;
model ch=group;
output out=residus r=resid;
lsmeans group;
run;
 
/*proc print data=residus;
run;*/
 
data data1;
RMSE = sqrt((sum(resid^2)/75));
run;
proc print data=RMSE;
run;

MERCI!!!!!!!!!!!!!!!!!

[Mag35]

Bonjour,

Pour mettre au carré la syntaxe est **

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
proc glm DATA=data1;
class GROUP;
model ch=GROUP;
output out=residus r=resid;
lsmeans GROUP;
run;
 
/*proc print data=residus;
run;*/
 
DATA data1;
set residus;
RMSE = sqrt((sum(resid**2)/75));
run;
proc print DATA=RMSE;
run;

J'espère que ça va fonctionner.

Magali

[marlene2110]

non ca ne marche pas ms en rajoutant une ligne c un peu mieux :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
DATA data1;
set residus;
RMSE = sqrt((sum(resid**2)/75));
put RMSE=;
run;

ms le pbm c est que g 75 obs ds ma table et donc 75 RMSE qui sortent je comprend pas...

[olivier.decourt]

L'étape DATA ne fonctionne qu'en lisant une ligne à la fois dans les données. Les calculs sont donc effectués ligne par ligne. Pas de souci pour la mise au carré, mais la somme n'est pas verticale. La fonction SUM dans l'étape DATA sert à additionner les valeurs de plusieurs variables dans une même ligne.
Ce sera plus simple en SQL.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
PROC SQL ;
SELECT sqrt((mean(resid**2))) AS RMSE FROM work.residus ;
QUIT ;

Et si tu veux mettre ça dans une table, il faut mettre

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

CREATE TABLE work.res1 AS

devant SELECT.

[marlene2110]

ok merci!!!