08 novembre
Par
Le
Dans graphique
Avec 3 Commentaires

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-08b2<-0.2652811b3<- 1986.235s<-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 axesplot.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égendesaxis(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 graphiquetitle(main= »force of infection »)#légende de l’ordonnée n°2mtext(« rs »,side=4,line=2.5)#on termine le graphiquebox()Read More →

06 novembre
Par
Le
Dans base indispensable
Avec 1 Commentaire

Quand on travaille sur un grand jeu de données, on peut voir à quoi il ressemble en rentrant son nom directement dans la console : data(iris) iris Mais pour de gros objets ce n’est pas du tout pratique… et souvent votre ordinateur peut avoir du  mal à tout vous afficher. Utilisez plutôt : data(iris) View(iris) # pensez bien au V majuscule Cela ouvre une nouvelle fenêtre avec l’intégralité du jeu de données. Et vous pouvez facilement vous balader dedans. Par contre si seul l’aspect général du tableau vous intéresse, vous pouvez utiliser la fonction head : data(iris) head(iris) Elle vous affiche les 5 premières lignesRead More →

05 novembre
Par
Le
Dans algorithmique, base indispensable
Avec 1 Commentaire

i<-0 while (i<10){ print(i) i<-i+1 } print(« on sort de la boucle ») While va réaliser ce qui est écrit entre les accolades {} tant que ce qui est dans les parenthèses () est vrai. Ce code peut donc se traduire de la façon suivante : i vaut 0 tant que (i est inférieur à 0){ afficher i augmenter i de 1 } Dans cet exemple au moment où i vaudra 10, on sortira de la boucle.Read More →

05 novembre
Par
Le
Dans algorithmique, base indispensable
Avec 11 Commentaires

for ( i in 1:10) { print(i) } Cette commande peut se traduire par : Pour (i allant de 1 à 10) { affiche i} Il faut noter que les parenthèses () servent à définir la variable et les valeurs qu’elle va prendre successivement à chaque tour de boucle. Les accolades {} servent à délimiter les actions à effectuer pour chacune des valeurs prises par la variable. IMPORTANT : R n’aime pas vraiment les boucles for, il est beaucoup plus efficace d’utiliser apply. Tout particulièrement pour les très grandes et longues boucles, apply fait cela en une fraction de seconde… alors que for peut mettreRead More →

04 novembre
Par
Le
Dans configuration de R, message d'erreur
Avec 0 Commentaire

Ce message apparaît lorsque l’on essaye de lancer un script initialement écrit sous windows sur une station linux (et inversement). Il s’agit d’un souci d’encodage. il est possible de corriger cela grâce au script suivant (remplacer latin1 par l’encodage d’origine du fichier) options(encoding= »latin1″) Ou alors options(encoding= »utf-8″)Read More →

04 novembre
Par
Le
Dans Transformation de données
Avec 0 Commentaire

MonVecteur=as.numeric(as.character(MonVecteur)) Cette commande permet de changer la classe de MonVector, reconnue comme « factor » par R mais codé à l’aide de chiffre, en classe « numeric ». Le passage par la classe « character » évite les changements des dits chiffres de MonVecteur souvent observés lors de l’utilisation de la fonction as.numeric(). exemple: vec<-as.factor(c(« 3″, »7″, »11″, »13 »)) vecas.numeric(vec) #ne donne pas ce que l’on veutas.numeric(as.character(vec)) # permet d’obtenir le résultat voulu.Read More →