1. Grâce à ce formulaire vous pouvez proposer une astuce ou un script sur R.
  2. Votre script doit pouvoir être lancé en l'état. Veuillez penser à :
    • inclure le chargement des "library" nécessaires
    • construire un petit jeu de données si besoin est.
    • commenter les lignes de codes pour en faciliter la compréhension.
  3. Le titre de votre script ou astuce doit être clair et explicite.
  4. Pensez à mettre votre code entre les balises [R] et [/R]; Pour cela, vous pouvez utiliser le bouton

Vous pouvez utiliser vos comptes Facebook, twitter ou google pour vous identifer (google est compatible yahoo, openID...)
L'ideal étant de vous connecter si vous avez un compte utilisateur, ou faire une demande de compte utilisateur si vous n'en avez pas encore.
Créer un compte va vous permettre de pouvoir éditer vos codes et de mettre en avant votre site internet.
Sinon vous pouvez soumettre anonymement en remplissant les champs ci-après.







Choisissez les catégories correspondantes à votre Code:

  • algorithmique
  • Analyse de survie
  • base indispensable
  • bayésien
  • configuration de R
  • exportation de données
  • fonctions utiles
  • graphique
  • importation de données
  • inférence
  • manipulation de données
  • message d'erreur
  • modélisation
  • Non classé
  • optimisation
  • planification
  • programmer avec R
  • regression linéaire
  • Test
  • tidyverse
  • Transformation de données


Comment tracer un graphique à deux ordonnées dans R ?

3
dans graphique
Il existe plusieurs façons de faire un graphique avec deux ordonnées. En voici une qui utilise les outils graphiques de base

 
 # Données d'exemple (peu importe...)
times<-seq(0,3000)
p<- 0.002197451 * exp(- 0.0009076665 *times)
b1<- 7.376812e-08
b2<-0.2652811
b3<- 1986.235
s<-b1*exp(-0.5*(log(times/b3)/b2)^2)
# On ouvre une nouvelle fenêtre plot.new()
# On choisit les paramètres de la fenêtre, voir ?par, ici mar correspond aux marges par(mar=c(5,4,3,4))
# On met le premier graphique en définissant les limites des axes
plot.new()
plot.window(xlim=c(0,3000),ylim=c(0,0.0022))
 
lines(p~times,type='l',col='burlywood1',lwd=3)
# on ajoute les axes et leurs légendes
axis(1)
axis(2)
title(xlab="time")
title(ylab="rp")
# On superpose le graphique avec un axe des ordonnées différent qui sera à droite du graphique (axis(4))
plot.window(xlim=c(0,3000),ylim=c(0,8e-08))
lines(s~times,type='l',col='burlywood3',lwd=3)
axis(4)
#titre du graphique
title(main="force of infection")
#légende de l'ordonnée n°2
mtext("rs",side=4,line=2.5)
#on termine le graphique
box() 
 
 

Proposé par Melen Leclerc.

Ce script vous à rendu service? remerciez l'auteur en votant ici:
- ça ne sert à rien -- c\'est interessant - (score de +11 sur 11 votes)
Loading ... Loading ...

3 Comments

  1. Céline
    Posté le 13 février 2014 a 10 h 43 min | Permalink

    Bonjour,merci beaucoup pour ce script, seulement quand je le teste, je n’arrive pas à avoir la légende de l’ordonnée 2. Est-ce que cela peut venir d’un package absent? Merci

  2. Vincent
    Posté le 13 février 2014 a 10 h 50 min | Permalink

    Bonjour Celine,

    le script tourne ( je viens de le retester). qu’entendez vous par legende? s’il s ‘agit des numeros sur les axes ils sont présent, s ‘il s ‘agit d ‘une « vrai » legende . voir ?legend

  3. Céline
    Posté le 13 février 2014 a 11 h 26 min | Permalink

    bonjour Vincent,en fait c’est cette partie qui ne marche pas #légende de l’ordonnée n°2 mtext(« rs »,side=4,line=2.5)mais en mettant line = 0.5, j’arrive à voir mon texte.je pense que c’est un problème de marge de mon garphique!!en attendant merci beaucoup, et excusez moi pour le dérangement !Céline

Poster un commentaire


Votre email ne seras jamais communiqué champs requis désigné par une *

*
*


cinq × = 10

Formation logiciel R