Bonjour à tous,
J'ai récupéré un programme SAS codé majoritairement en SAS IML (que je connais assez peu) et que je ne peux pas modifier à ma guise (pas de passage vers du code non IML notamment).
Je dispose par exemple d'une table de ce genre avec en ligne des âges et en colonne des taux à différentes années :
On souhaiterait modifier, à partir de 2013 seulement (colonne tx13), les taux pour tous les âges supérieurs ou égaux à 18, en les multipliant par une constante.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11 data testIML; input age tx12 tx13 tx14 tx15; cards; 17 0.001 0.002 0.003 0.004 18 0.004 0.005 0.006 0.007 19 0.005 0.008 0.009 0.010 20 0.010 0.011 0.012 0.013 99 0.099 0.014 0.015 0.016 100 1.000 0.017 0.018 0.019 ; run;
Actuellement, le code ressemble en gros à ça :
L'utilisation d'indices i, k, age est compliquée à (re)lire et à maintenir/modifier. Si on a un changement dans la table d'entrée avec l'ajout des âges entre 0 et 17 par exemple, ou encore une multiplication par une autre matrice dans lesquelles on n'a pas tous les mêmes âges, on se retrouve vite avec des i-13, k+3... d'autant plus compliqués à gérer.
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
17 proc iml; use tableinit var _all_; setin tableinit; read all into Minit; %let constante = 3; Mtransformee = Minit; do k=2 to &HORIZON.; /* k=indice des colonnes*/ do age=18 to &MaxAge.; /* age=indice des lignes*/ Mtransformee [age-16,k] = &constante. * Minit[age-16,k]; end; end; create tabletrans from Mtransformee; append from Mtransformee; quit;
Les variables contenues dans la table initiale ("importée" ensuite dans une matrice) étant suffixées avec une année, j'aimerais bien pouvoir faire un code plus lisible qui irait chercher directement la bonne ligne et la bonne colonne à partir de l'année (pour les colonnes) et de l'âge (pour les lignes) et ce en s'affranchissant de la structure de la table initiale.
Pour cela, j'ai trouvé notamment la fonction loc. Mais pour le moment je n'arrive à agir que sur la matrice initiale en "vectorisant" les variables suffixées et en allant chercher les bonnes lignes grâce à la fonction loc :
Quelqu'un aurait-il une idée pour s'en affranchir dans la première partie de la formule (celle où on affectite les valeurs dans la matrice de résultat Mout) ?
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 %macro test0; proc iml; use testIML3; read all var {age %do i=12 %to 15; tx&i. %end;}; read all /*var _all_*/ into MTempIML; Mout=MTempIML; %do an=13 %to 15; do agex=18 to 20; Mout[agex-16,&an.-10] = 2*tx&an.[loc(age=agex)]; end; %end; PRINT MTempIML; PRINT Mout; quit; %mend; %test0;
Parallèlement à ça, est-il possible de gérer en amont une erreur dans le cas où on ne trouve pas une variable (ou un âge), par exemple dans le cas où on change la table initiale avec les données de 2015 à 2018 ?
Merci d'avance
Alqua
Partager