Groupe des utilisateurs du logiciel R

En fait je me demande si la question posée ne serait pas un peu plus pertinente s'il y avait un effet de la variable etud? Non, au contraire : le fait qu'etud soit (faussement) significatif dans un cas et pas dans l'autre accentue même le problème. Mais l'analyse sortie par défaut est incorrecte, q...

[Pour info ce message sera prochainement édité et déplacé dans Fiches]

Bonjour, Il y a eu récemment une question sur la syntaxe aov : http://forums.cirad.fr/logiciel-R/viewtopic.php?t=587 qui me conduit à la question suivante. Supposons un essai avec 2 facteurs de contrôle control1 et control2 , 1 facteur étudié etud , et une variable x . Le jeu de données fictif suiva...

Bonjour, Considérons le graphe suivant : df <- data.frame(x=rep(1:20,each=12),y=rep(1:12,20)) windows(6,4) par(mar=c(3,3,2,2)) symbols(df,squares=rep(0.5,nrow(df)),inches=F,asp=1) title(paste(nrow(df),"point...

Vous pouvez aussi poster dans "Questions en cours" un message décrivant votre document, par exemple dans ce style : http://forums.cirad.fr/logiciel-R/viewtopic.php?t=35 Le message sera déplacé par les modérateurs dans "Fiches" et il apparaîtra ainsi avec votre nom.

Il faut sortir pdf() et dev.off() de la fonction : graphes <- function(x, namesvar){ dens=density(x) ... print(qqmath(resid, xlab="quantiles theoriques", ylab="quantiles experimentaux")) } pdf("plot.pdf",onefile=T) graph...

Voir aussi la fonction chull du package graphics.

Bonjour,
Merci pour ce lien, cependant il est hors sujet car il ne concerne pas le logiciel R.

Une petite précision sur l'utilisation du forum : de plus en plus d'intervenants postent dans un même sujet une succession de messages dont chacun complète ou rectifie ses propres messages précédents. Plutôt que d'utiliser cette méthode, chacun a la possibilité (avec la touche "éditer" en ...

J'ai essayé ceci et je ne comprends pas les résultats de ces commandes : > data3 <- data[col1==1] > data3 col1 col2 1 1 ab 2 1 cd 3 2 ef 4 2 gh C'est très logique. Un data.frame est un cas particulier de "list" dont les éléments sont les colonnes. Donc tout index simple (à une seule dimen...

Bonjour à tous, J'interviens en tant qu'administrateur du forum (mais pas, loin s'en faut, le principal animateur!). Il faut tout d'abord expliquer l'origine du forum. Il est issu d'un groupe de travail sur les logiciels statistiques, interne à notre organisme [le Cirad : Centre de coopération inter...

Bonjour, Pour calculer la différence entre éléments successifs mieux vaut utiliser la fonction diff. Par exemple on peut efficacement remplacer w <- sapply(1:(100000-1),function(x){FicBase[x,1]-FicBase[x+1,1]-FicBase[x,2]+FicBase[x+1,2]}) par w <- diff(FicBase[,...

Francois, est ce que ce n'est pas plutot if (length(x)>1) panel.histogram(x,...) a la place de if (length(x)>0) panel.histogram(x,...) Non car les NA ne sont pas transmis, pour s'en persuader on peut ajouter un "print" dans l'appel à "histogram" histogram(~ x | z, data = myd...

Si on veut les 3 histogrammes superposés on peut faire comme ça : mydata <- data.frame( x = rnorm(100) + rep(c(0, 1), each = 50), z = rep(c("A", "B"), each = 50)) mydata2 <- rbind(data.frame(x=NA,z="Total"),myd...

Oui mais ça ne marche pas si les dates ne pas les mêmes car maree$date, étant de classe "factor", est mal convertie Donc : maree$date <- as.Date(as.character(maree$date), format = "%d/%m/%Y") en effet : > as.Date(maree$date, format = "%d/%m/%Y") ...