Discussion :

[Lichen]

bonjour,
je lance une AFM avec FactominR mais j'obtiens un message d'erreur "objet res introuvable", et je n'ai aucune sortie.
Ces mêmes données m'avaient pourtant permis de réaliser une ACM sans problème.
L'AFM que je cherche à réaliser est constituée d'un groupe de 8 variables actives (à 3 modalités chacune) et d'un groupe de 5 variables explicatives (2 modalités chacune) + 1 variable explicative (3 modalités)

J'ai essayé de refaire l'analyse en supprimant le groupe de variables explicatives pour lesquelles il y a qques NA mais ça ne fonctionne pas pour autant.
Je n'ai aucune idée de quoi l'erreur peut provenir
Merci d'avance si vous pouvez m'aider à identifier le problème.

[Supernul]

Je n'ai aucune idée de quoi l'erreur peut provenir

Moi non plus.


Bonjour, ma voiture a un problème, je n'arrive pas à la faire avancer.

C'est en substance ce que vous nous donner comme information. Pour vous aider il nous faut plus comme, par exemple votre code, des données pour reproduire l'erreur si le code n'est pas suffisant, ...

[Lichen]

c'est tout simplement que lorsqu'on n'y connait rien, on n'a aucune mesure de l'ampleur du questionnement....
Je ne pouvais ainsi pas savoir si c'était une erreur banale et courante liée à un petit qqchose de connu ou si cela pouvait provenir de multiples raisons.

Voici le code lancé par FActominR (j'avais demandé une classification derrière l'AFM)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
rokval.MFA<-rokval[, c("d_tsol", "d_delai", "d_desh", "d_rdtbas", "d_rdtirg", "d_recolt", "d_smce", "d_info", "d_hum", "Sw_secucat", "Sw_reucat", 
  "Sw_uvNacat", "Sw_Bvcat", "Sw_autocat", "syst_expl")]
res<-MFA(rokval.MFA, group=c(9, 5, 1), type=c("n", "n", "n"), ncp=5, name.group=c("diff", "val", "syst"), num.group.sup=c(2, 3), graph=FALSE)
res.hcpc<-HCPC(res ,nb.clust=-1,consol=TRUE,min=3,max=10,graph=TRUE)
res.hcpc$data.clust[,ncol(res.hcpc$data.clust),drop=F]
res.hcpc$desc.var
res.hcpc$desc.axes
res.hcpc$desc.ind
summary(res, nb.dec = 3, nbelements=10, nbind = 10, ncp = 3, file="")
 
remove(rokval.MFA)

j'ai supprimé ici les lignes de codes liées à la demande d'écriture des différents résultats dans un fichier csv pour alléger.
Merci de votre aide

[VonVelten]

Bonjour,

Si votre message d'erreur est que l'objet res est introuvable, c'est qu'il n'a pas été créé, et donc que l'AFM ne fonctionne pas.

L'idéal serait d'avoir le message d'erreur après avoir exécuté la commande pour réaliser l'AFM. Au vu de ce qui est écrit, je pencherai pour le fait qu'il n'y a qu'un tableau actif pour l'AFM, or elle en nécessite au moins deux pour moi. Sinon, il faut réaliser une ACP qui est plus spécifique de l’analyse à un tableau, et préciser les données supplémentaires qualitatives et quantitatives colonne par colonne.

Cordialement

VV

[Lichen]

j'ai toutes mes données dans un seul tableau.
En fait je pensais que l'AFM pouvait m'apporter un plus par rapport à l'ACM que j'avais réalisée précédemment car elle avait le mérite de pouvoir représenter mes variables par groupe, ce qui correspond à mon cas de figure : j'ai un groupe de variables actives qui correspond à 6 variables "difficultés rencontrées" et un groupe qui comprend 5 variables explicatives "motivations".
Mais tout est stocké dans un même data qui comprends la totalité de ma centaine de variables . Il faudrait que je les sépare physiquement dans deux tableaux ?

En remontant dans les messages d'erreur je note qu'il y a une erreur supplémentaire (mes excuses je ne l'avais pas vue car elle était toute seule au milieu de multiples lignes 'res introuvable')

[45] ERREUR:
dim(X) doit avoir un longueur positive
[46] ERREUR:
objet 'res' introuvable

[VonVelten]

Le principe de l'AFM est qu'elle permet de séparer un tableau en plusieurs sous-tableaux qui sont analysés chacun de leur côté dans un premier temps, puis analysés globalement en équilibrant leur apport d'information (Pour faire simple).

Pour cela, il faut au moins deux sous-tableaux actifs, ce qui n'est pas ce cas ici. Or une AFM à un sous-tableau actif (plus d'éventuels tableaux supplémentaires) revient à faire une ACP (en grossissant le principe de généralisation de l'ACP en AFM). Vu qu'il n'y a donc ici qu'un sous-tableau actif (car 2 sont indiqués en supplémentaires, num.group.sup=c(2, 3)), il revient donc de faire une ACP directement (et d'indiquer dans ce cas les variables supplémentaires).

L'erreur dim(X) doit avoir un longueur positive vient du fait qu'il n'y a qu'un seul tableau avec des variables actives dans l'AFM ici. Les autres sont toutes supplémentaires et ne participent donc pas à la construction des axes factoriels. L'AFM ne peut donc pas converger (car problème avec les données d'entrées).

Je conseille de faire une ACP et de juste renseigner les variables supplémentaires.

J'espère que c'est assez clair, j'ai essayé de ne pas trop rentrer dans le détail, mais j'étais un peu obligé pour expliquer l'erreur.

[Lichen]

Merci beaucoup de votre réponse. J'avais en définitive mal compris les principes de l'AFM. Dans mon cas il faut je suppose une ACM plutôt qu'une ACP car je n'ai que des variables qualitatives.
Je me permets de poser une question complémentaire liée au code que j'ai déjà présenté :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
rokval.MFA<-rokval[, c("d_tsol", "d_delai", "d_desh", "d_rdtbas", "d_rdtirg", "d_recolt", "d_smce", "d_info", "d_hum", "Sw_secucat", "Sw_reucat", 
  "Sw_uvNacat", "Sw_Bvcat", "Sw_autocat", "syst_expl")]
res<-MFA(rokval.MFA, group=c(9, 5, 1), type=c("n", "n", "n"), ncp=5, name.group=c("diff", "val", "syst"), num.group.sup=c(2, 3), graph=FALSE)
res.hcpc<-HCPC(res ,nb.clust=-1,consol=TRUE,min=3,max=10,graph=TRUE)
res.hcpc$data.clust[,ncol(res.hcpc$data.clust),drop=F]
res.hcpc$desc.var
res.hcpc$desc.axes
res.hcpc$desc.ind
summary(res, nb.dec = 3, nbelements=10, nbind = 10, ncp = 3, file="")
 
remove(rokval.MFA)

Jj'ai essayé différentes méthodes maie n'arrive pas à obtenir les données de res.hcpc dans un csv : y'a t'il moyen de les exporter ?
Encore merci

[VonVelten]

Effectivement les résultats de HCPC ne sont pas prévus pour être sauvegardés directement dans un .csv

Il y a deux possibilités, soit aller chercher les éléments un par un qui sont dans des formats classiques (comme la sortie data.clust) et les enregistrer dans des fichiers séparés. C'est long et fastidieux mais ca donne des fichiers propres.


Autre méthode plus barbare : utiliser la fonction sink(). Elle sauvegarde les éléments affichées dans la console dans un fichier voulu. Le résultat est moins propre mais beaucoup plus rapide :

Code R :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
sink("Sauvegarde_HCPC.csv")
for(i in 1:5){  #vérifier l'indexation maximale pour une ACM
print(res.hcpc[[i]])
}
sink()