Je souhaite réaliser un histogramme empilé à partir du tableau de données suivant:
Code : Tout sélectionner
> head(data)
codsta etat.scenario1 etat.scenario2 etat.scenario3 morphologie
1 6177970 bon bon bon faible
2 6153060 bon bon bon faible
3 6005550 moyen moyen moyen fort
4 6168200 bon tres_bon bon moyen
5 6300600 médiocre moyen moyen fort
6 6016940 bon bon mauvais fort
J'ai d'abord calculé les pourcentages de stations dans les 5 classes d'état selon la morphologie, et ce séparément pour les 3 scénarios étudiés.
Code : Tout sélectionner
prop1<-prop.table(table(data$etat.scenario1,data$morphologie),2)
prop1
faible fort moyen
bon 0.29383886 0.30603448 0.28240741
mauvais 0.00000000 0.07758621 0.03240741
mediocre 0.04739336 0.22413793 0.11574074
moyen 0.12796209 0.24568966 0.26388889
tres_bon 0.53080569 0.14655172 0.30555556
prop2<-prop.table(table(data$etat.scenario2,data$morphologie),2)
prop3<-prop.table(table(data$etat.scenario3,data$morphologie),2)
barplot(prop1,legend=T)
barplot(prop2)
barplot(prop3)
Comment faire pour représenter sur le même graphique ces trois barplot? C'est à dire que pour chacune des modalités de ma variable "morphologie", j'aimerai afficher 3 barres empilées cote à cote, correspondant aux 3 scenarios d'évaluation de l'état.
J'ai essayé de passer par ggplot mais je n'arrive pas à mettre en forme mes données tout en calculant les pourcentages.
Merci pour votre aide.