Mise à jour de ce billet: novembre 2020. => passage au tidyverse pour la manipulation des tables de données et les graphiques. A quoi sert-il? Le test U de Mann-Whitney (aussi appelé test de la somme des rangs de Wilcoxon ou plus simplement test de Wilcoxon) sert à tester l’hypothèse selon laquelle la distribution des données est la même pour deux groupes. La p-value associée à ce test va ainsi.. Read More
Ajuster un modèle de distribution à des données
Ce billet vise à vous montrer comment ajuster un modèle de distribution (par exemple une jolie "cloche" gaussienne) à une distribution empirique (i.e. la distribution de vos données telle que vous l'observez, par exemple en traçant l'histogramme correspondant). Il existe plusieurs méthodes possibles pour faire cela, dont la méthode des moments dont je vais parler ici. Considérons les données disponibles ici x=read.csv(paste(dat.path,"data_vraisemblance.csv",sep=""),sep=";")$x head(x) ## [1] 132.5 194.8 179.0 108.7 112.5.. Read More
Faire des boucles for dans R pour automatiser des tâches
La boucle for est un peu le B-A BA de la programmation sous R, mais (pour ceux qui ne connaissent pas), elle a de quoi vous mettre des étoiles dans les yeux (sans exagérer!). En effet, elle permet de vous éviter certaines tâches répétitives en les automatisant, et ainsi de vous faire gagner beaucoup de temps... Voici comment elle est structurée: for (i in 1:n){ ## tâche à répéter n.. Read More
Calculer une aire entre deux courbes
Dans ce post je montre comment utiliser des fonctions que j'ai écrites pour calculer des aires: entre deux courbes entre une courbe et une droite horizontale Cela peut être utile si vous travaillez sur des profils spatiaux qui évoluent au cours du temps (profil topographique par exemple) et que vous souhaitez calculer un bilan (sédimentaire par exemple) entre ces profils. Données d'exemple et script à télécharger Considérons le jeu de.. Read More
Une fonction magique: tapply
Cette fonction permet de calculer une statistique par niveau d'un facteur. Par exemple, considérons les données accessibles ici: ## Warning: attempt to set 'sep' ignored ## site taille ## 1 site1 3.125 ## 2 site1 1.162 ## 3 site1 4.222 ## 4 site1 2.748 ## 5 site1 3.625 ## 6 site1 3.892 Pour connaître la moyenne (par exemple) de la taille en fonction du site, on peut procéder comme suit:.. Read More