Discussion :

[zacharie.moulin]

Bonjour,

Je souhaite faire un tableau croisé me calculant le total de surfaces d'un état A vers un état B. Mon fichier .dbf en entrée (issu d'un shapefile), contient 3 champs : l'état initial, l'état final et la surface. Mes individus sont des polygones (étant donné que ça vient d'un SIG).

Exemple :

Etat initial	Etat final	Surface
Etat A	Etat A	30
Etat B	Etat A	10
Etat A	Etat A	20
Etat B	Etat A	20
Etat B	Etat B	40

Je souhaite obtenir le tableau suivant :

initial \ final	Etat A	Etat B
Etat A	50	0
Etat B	30	40

Le but au final est de faire un graphe de transition d'un état à l'autre, comme on peut le voir sur ce document page 8.

Je précise que je n'ai pas nécessairement le même nombre d'états initiaux et finaux, et donc par conséquent certains états initiaux ne se retrouvent pas en finaux et inversement.

J'espère avoir été assez clair dans l'expression de mon besoin. Merci de votre aide !

[mgdondon]

Bonjour,

Une proposition de code avec la librairie reshape et la fonction cast(). C'est l'équivalent d'un tableau croisé dynamique :

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
> toto <- data.frame(initial=c("A","B","A","B","B"),
+                    final=c("A","A","A","A","B"),
+                    value=c(30,10,20,20,40))
> toto
  initial final value
1       A     A    30
2       B     A    10
3       A     A    20
4       B     A    20
5       B     B    40
> library(reshape)
> cast(toto,initial~final,sum)
  initial  A  B
1       A 50  0
2       B 30 40

Cordialement,

[zacharie.moulin]

Désolé pour le retard de réponse, j'ai eu d'autres affaires entre temps.
Merci beaucoup c'est exactement ce que je cherchais

Zacharie

[zacharie.moulin]

Bonjour,

je ré-ouvre la discussion, puisque j'ai un petit problème.

La fonction cast marche très bien, lorsque ma dernière colonne est celle utilisée pour le value. Est-ce forcément la dernière colonne ou peut-on spécifier le nom de la colonne à utiliser? En effet j'ai rajouté une colonne "annee" à mon data frame, et je voudrais lui spécifier d'utiliser la colonne "surface".

J'ai donc colonnes dans ;on data frame : initial, final, surface et anneeJ'ai essayé avec le code suivant : cast(toto,initial~final,sum, surface=guess_value(toto))d'après ce que j'avais pu lire dans l'aide de la fonction.
J'ai un message d'erreur "invalid 'type' (character) of argument" pourtant ma colonne est bien numérique.

Peut-être pouvez vous m'aider mgdondon,
Merci beaucoup !

[mgdondon]

Bonjour,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
> toto <- data.frame(initial=c("A","B","A","B","B"),
+                    final=c("A","A","A","A","B"),
+                    surface=c(30,10,20,20,40),
+                    annee=c(1980,1990,1994,1982,1990))
> toto
  initial final surface annee
1       A     A      30  1980
2       B     A      10  1990
3       A     A      20  1994
4       B     A      20  1982
5       B     B      40  1990

Code :

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

1
2
3
4
5
> library(reshape)
> cast(toto,initial~final,value="surface",sum)
  initial  A  B
1       A 50  0
2       B 30 40

Cordialement,

[zacharie.moulin]

Oui, tout simplement, les ""...

Merci beaucoup, je clos à nouveau le sujet