Discussion :

[mathieu_r]

Bonjour,


J'ai un petit problème de calcul de matrice à vous confier.

dans le cadre d'une analyse géographique multi-critères, je dispose de données diverses: la pente, l'altitude, ainsi que la densité de population sur un secteur donné.

A partir de cela, je dois définir des aires d'aménagement optimales en appliquant des poids à chaque variable de manière relative.

Ainsi, la pente sera deux fois plus importante à prendre en compte par exemple que l'altitude, la population 4 fois plus que la pente.

Cela doit aboutir à une matrice dont les dimensions en colonnes et en lignes sont identiques et dont la diagonale égale 1. La diagonale constitue une sorte de miroir, un axe symétrique de part et d'autre de laquelle on trouve l'inverse des valeurs d'un côté.

Sachant que la matrice de poids relatifs a été établie de cette façon:


Comment aboutir rapidement à ceci (fonction déjà existante?)


Pour l'instant ma méthode est assez la laborieuse:
je définis la matrice de poids mat1

Code :

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

mat1[2,1]=2

etc...
Puis je trouve la réciproque de la matrice de laquelle j'inverse les valeurs

Code :

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

mat2=1/t(mat2)

je remplace les valeurs nulles des matrices par 0 pour pouvoir effectuer la somme de la matrice et de sa réciproque

Code :

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

1
2
3
mat[sapply(mat, is.na)]=0
mat2[sapply(mat2, is.na)]=0
mat3=mat+mat2

Je définis la diagonale égale à 1:

Code :

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

diag(mat3)=1

Assez laborieux tout de même: je pense vraiment qu'il existe une fonction pour ça...

[lilly74]

Bonjour,
je vous propose:

Code :

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

1
2
3
4
n=2; a=2
M <- matrix(rep(c(a^(0:n), 0), n+1), nrow=n+1, ncol=n+1)
M[3,2] <- 0.5
M[upper.tri(M)] <- 1/t(M[lower.tri(M)])

Je n'ai pas trouvé mieux...

Bonne journée

[mathieu_r]

La formule magique consistait donc à utiliser upper.tri et lower.tri.
C'est nickel..

Sinon, une question qui n'a rien à voir.

Pourquoi utilise-t-on généralement le signe <- et pas = ?

[ced]

Pourquoi utilise-t-on généralement le signe <- et pas = ?

Pour des raisons "historiques", et notamment de compatibilité avec le langage S et S-Plus.

En outre, les 2 opérateurs ("=" et "<-") n'ont pas tout-à-fait le même comportement. Le "<-" peut être utilisé dans toutes les expressions, alors que l'usage du "=" est plus restreint.

Exemples :

Code :

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

1
2
3
4
> median(x = 1:10)
[1] 5.5
> x
Erreur : objet 'x' introuvable

Alors que :

Code :

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

1
2
3
4
> median(x <- 1:10)
[1] 5.5
> x
 [1]  1  2  3  4  5  6  7  8  9 10

[mathieu_r]

ok ça marche