Re, Il peut peut être s'agir d'un problème de classe de tes objets d'un côté tu as un array mais de l'autre de quoi s'agit-il ? D'une matrice, d'un data frame ?, d'une liste ? S'il s'agit d'un data frame il te faudra d'abord le passé en matrice pour que ta structute soit conservée. ted <- array(...
Bonjour, L'erreur vient du fait que lorsque tu fais eval(parse(text=paste("sst",i,"_3D",sep="")))[,,j] il te renvoie bien le tableau mais pas l'objet en question donc tu ne peux le modifier l'objet de cette manière. Je peux me tromper mais c'est comme si tu tapais rnorm...
Bonsoir, Il est difficile de répondre à ta question sans avoir une petite idée de tes objets, pourrais tu mettre un exemple que l'on puisse regarder. Une question qui ne mange pas de pain est-ce que tes j et tes i sont bien définis ? Ne manque t'il pas des "" pour entourer datName ? Il sem...
Re, Tu peux tout exporter dans un fichier texte que tu pourras ouvrir dans excel ensuite avec ce script (ya pas mal de trucs qui ne servent qu'a la mise en forme dans le ficher texte) : for (i in 1:length(res)){ if (i<2){ sink("resultat.txt",append=FAL...
Re, Peux tu juste mettre les 10 premières lignes de ton tableau par exemple. # ne connaissant pas vraiment tes données je te propose ça # pour le mois de juillet, data etant le tablmeau de données plot(data$dateJ[format(data$dateJ,"%m")=="07"],data$donne[format(da...
Une question concernant tes données comment les as-tu incoporés dans R, s'agit-il de nombre, de chaine de caractères ou encore de date au sens format date de R ?
Bonjour, Désolé de ne pas avoir pu répondre plus tôt à ta question. Le apply du promier post te permet de récupérer le pourcentage de chaque modalité de chaque réponse. Le résultat est généré sous forme de liste car toutes tes questions n'ont pas le même nombre de réponses donc si tu veux savoir com...
Re, Tu peux aussi utiliser ça : co <- combn(ncol(don)-1,2) # toutes les combinaisons de variables 2 à 2 res <- list() for (i in 1:ncol(co)){ auxi <- table(don[don$fin!=0,co[,i]]) class(auxi) <- "matrix" auxi <- auxi[apply(aux...
Bonjour, Si tu veux en rester a ton test de chi deux tu peux le faire de la manière suivante : p = y/sum(x) # en prenant pour hypothèse que les y sont les valeurs théoriques je recalcule les probabilités associés aux valeurs théoriques chisq.test(x,p=p) Chi-squared test for given pro...
Bonjour, Je te propose le code suivant : apply(don[don$fin!=0,],2,function(x) table(x)/sum(!is.na(x))*100) # et ça te donne ça : $X1_1 x 1 2 3 4 5 4.761905 4.761905 2.380952 14.285714 73.809524 $X1_1_1 x 1 3 85.71429 14.28571 $X1_2 x 1 2 4 27.5 55.0 17.5 $X1_2...
Bonjour, tu peux très bien crée un facteur US contre le reste du monde avec la fonction ifelse : pays <- sample(c("US","FRA","ESP","BRE"),20,rep=T) pays [1] "US" "ESP" "US" "ESP" "FRA" "BRE...
Re, tu peux essayer la méthode de Tukey pour comparer des moyennes deux à deux : y <- rnorm(100) x <- gl(4,25) aov1 <- aov(y~x) TukeyHSD(aov1) Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = y ~ x) $x diff lwr upr p adj ...
Re, Ce n'est pas si "simple" que ça, puisque si tu as suffisement de données la normalité n'est pas le facteur le plus important, contrairement à l'homogénéité des variances, donc si celle-ci est avérée tu peux utiliser des tests paramétriques et revenir à des tests de comparaison deux a d...
Sébastien Boutry a écrit :j'ai utilisé un test de Kruskal Wallis car je dois utiliser un test non paramétrique pour différencier mes groupes qui sont indépendants
C'est un peu plus clair mais qu'est-ce qui vous oblige à utiliser un test non paramétrique ? problème d'homogénéité de variance ? etc ...
Il est difficile de répondre à ce genre de questions sans savoir ce que tu cherches à tester (quelle est la question que tu te poses) ni sans avoir une petite idée des données. Par exemple pourquoi avoir utiliser un test de Kruskal Wallis ? etc.