Bonjour, Je ne pense pas que le problème viennent de df, mais de density : hist(sim,prob=T,xlim=c(0,2),breaks=c(seq(0,2,le=30),max(sim))) curve(df(x,df1,df2),add=TRUE,col=2) La on voit bien que les simulations suivent une loi de Fisher. Je pens...
Bonjour, Oui pour chaque axe il suffit de diviser la valeur propre par la somme des valeurs propres et ça te donne l'inertie expliqué par l'axe : library(ade4) data(banque) banque.acm <- dudi.acm(banque,scannf=F,nf=20) banque.acm$eig/sum(banque.acm$eig)*100 # voir aus...
Bonjour, La différence que Nico observe provient de control1.2 et de la modalité qui est prise en compte dans le modèle. Normalement il devrait prendre en compte pour l'interaction le fait d'être à la fois 2 de control1 et 2 de control2 donc la modalité 2.2 or dans le modèle de Nico ce n'est pas cet...
Bonjour, Je crois qu'ici on est en droit de se poser quelques questions tout de même : quel importance attribué à une analyse de variance sur 12 valeurs avec 3 facteurs qui prennent respectivement 2 ou 3 modalités, et plus que 6 df pour estimer la variance résiduelle ? Quelle signification ce SCEres...
Bonjour, Premièrement je te conseille de regarder le lien suivant : https://www.biostat.envt.fr/wiki/index.php?title=Analyse_de_variance&printable=yes#Conditions_d.27utilisation_d.27une_ANOVA Et tu verras que quand tu emploies la synthaxe facteur1*facteur2 tu ne considères pas l'effet du facteur...
Bonjour, Pour faire des manips recursiuive je vois deux possibilités (qui ne sont surement pas les seules) : 1) tu utilises un while et tant que ta condition n'est pas vérifié tu continues 2) tu utilises une fonction qui fait appel a elle même si tu dois continuer. En terme d'optimisation je ne sais...
Re, le as.is par défaut est FALSE donc je ne l'ai pas modifié ct juste pour attirer l'attention sur ce paramètre. Sinon ce que j'ai modifié c'est le dec="." au lieu de dec=",". Si tu lui dis que le séparateur décimal est une virgule il va traduire tes 0.258 en carcactère puisque ...
Bonjour, Pour te répondre Elodie. Le problème de ton apparté vient du dernier ";" sur la ligne, comme c'est un séparateur il attend une valeur derrière n'en trouve pas et retourne donc un NA. Essaie ça : essai <- read.table("outpg3.fasta",header=F,sep=";",dec="...
Bonjour, Essaye peut-être de mettre a jour la librairie lattice parce que en reprenant exactement ton code, j'ai bien obtenu ce qui tu cherchais à obtenir. J'ai bien des points séparés avec des couleurs différentes ou des droites qui ne se superposent pas. Sous quel forme sont stocké tes différents ...
Bonjour, Pour enlever les points : col.p="transparent". Après il se peut aussi que tes classes n'aient pas de sens et que de ce fait les ellipses se superposent. Soit c'est un problème graphique soit ça n'a pas de sens. Pour ce qui est des valeurs de l'ACP Nico te propose d'utiliser les co...
Re, tu n'as pas de 3ème axe car ayant mis scannf=FALSE tu ne choisis pas le nombre d'axe et par défaut nf=2 donc que les deux premiers axes. Sois tu vires scannf=FALSE et après tu rentres 3 quand il te pose la question soit tu gardes scannf=F, mais derrière tu rajoutes nf=3. Pour ce qui est du trois...
Bonjour, Je pense que le problème vient du facteur que tu as rentré dans discrimin. Quel est la tête de T$Palmier ? Es-tu sur qu'il ne prend que trois modalités (à priori ici je dirais qu'il en prend 10 soit une par palmier). De plus je ne sais pas ce que tu peux avoir comme information sur des palm...
Bonsoir, Moi je te répondrai ceci : ce forum n'est pas à priori un forum pour débattre de la recherche française ou des états d'âme des uns et des autres. Crois moi je sais ce que c'est d'avoir un encadrement pourri pendant un stage important, pour autant tout les gens ne sont pas à mettre dans le m...
Bonjour, Je ne connais pa SAS alors je n'ai pas d'idée précise concernant la tête des sorties de SAS mais je pense que ça peut t'aider : x <- rnorm(200) y <- gl(2,100) test <- kruskal.test(x~y) bob <- cbind(test$data,test$statistic,test$p.value) dimnames(bob) ...