Bonjour, Voici une définition tout droit sortie de l'un de mes cours de stats sur les GLM, néanmoins elle s'applique à tout les modèles je pense : La vraisemblance de la valeur d ’un paramètre est la probabilité des données si le paramètre a cette valeur. C ’est donc une quantité qui mesure « l ’acc...
Re, Pour mettre en liste les résultats tu peux faire comme ça avec a=-2 et b=5 : x <- NULL for (i in 1:3){ u <- combn(c(1,2,3),i) x <- c(x, apply(u,2,function(w) paste("c(",-2,",",paste(w,collapse=","),",...
Bonjour, Je n'ai pas tout saisi ce que tu comptes faire, mais tu peux générer ta séquence de cette façon : x <- NULL for (i in 1:3){ u <- combn(c(1,2,3),i) x <- c(x, apply(u,2,function(w) paste("[a,",paste(w,collapse=","),&...
Bonjour, Le problème vient surement d'une ou de plusieurs valeurs manquantes dans ton tableau données, des NA sous R. Pour le savoir tu peux faire : any(is.na(données)), si as pour résultat TRUE alors c'est que quelque part dans ton tableau tu as une valeur manquante d'ou le message d'ereur. Pour sa...
Re, pardon, ted2 n'est pas une fonction mais le nom de la matrice. J'avais l'exemple sur mon ordi avec ted2 comme nom de matrice et ensuite j'avais voulu harmonisé avec le nom de ta matrice à toi. Une fonction étant suivie de parentèse alors que les crochet servent à l'indexation. Avec correction ça...
Bonjour, Voici deux solutions qui je pense te donneront le résultat que tu attends : GG3 <- cbind(c(rep(0,p-1),1),rbind(diag(1,p-1,p-1),rep(0,p-1))) # ou GG3.2 <- diag(1,p) GG3.2 <- ted2[,c(p,1:(p-1))] all.equal(GG3,GG3.2...
Re, Le problème ne vient pas je pense de cl$cluster. Il vient de acp$li. Dans mon exemple acp$li faisait référence à une acp faite avec la fonction dudi.pca de ade4. Avec princomp il te faut utiliser acp$scores. Les deux étant identiques au signe près. Personnellement je suis plus habitué à ade4 que...
En fait ça doit venir de la dernière ligne de ton fichier qui doit être une ligne de NA. donnees[115400,] pour voir. Le all... te sert à savoir si toutes tes données ne sont pas des valeurs manquantes et on voit bien que ce n'est pas le cas.
ça veut bien dire que tu as une valeur de NA dans chacune de tes colonnes. Pour savoir laquel tu peux faire : lapply(donnees,function(x) which(is.na(x))). Sinon pour vérfier l'ensemble : all(!is.na(donnees)).
Bonjour, Une des premières choses à faire quand on se lance dans ce genre d'analyse est de comprendre un minimum la "philosophie" de l'analyse et ce qu'elle fait. Je t'invite à regarder ce lien qui te donnera pas mal d'explications quand au sortie de la cca sous ade4. http://pbil.univ-lyon...
Bonjour, Je t'invite à lire le bas de la page 38 et le début de la page 39 de : Practical Regression and Anova using R, Julian J. Faraway. Tu trouveras le lien dans la rubrique liens. Cette corrélation y est expliquée. Je t'invite aussi à lire : COLINÉARITÉ ET RÉGRESSION LINÉAIRE, Thierry FOUCART qy...