calcul sur une matrice

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:
http://data.imagup.com/10/1128420585.gif

Comment aboutir rapidement à ceci (fonction déjà existante?)
http://data.imagup.com/12/1128420633.gif

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

Code:

---

mat1[2,1]=2

---

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

Code:

---

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:

---

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:

---

diag(mat3)=1

---

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

Bonjour,
je vous propose:

Code:

---

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 :)

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 = ?

Citation:

---

Envoyé par mathieu_r

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:

---

1 2 3 4

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

---

Alors que :

Code:

---

1 2 3 4

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

---

ok ça marche