cod'erreur 1 dans le sous programme Lapack 'dgesdd'

En voulant réaliser une ACM sur un fichier de données composées de 754 403 lignes et 9 variables toutes 'factor', R me retroune le message d'erreur suivant :

" cod'erreur 1 dans le sous programme Lapack 'dgesdd' "

Pour info, je réalise mon ACM avec le package FactoMineR. Et la commande MCA.

Merci, en espérant une réponse à ce problème. :ccool:

Bonjour,
Est ce que vous pourriez poster votre bout de code et le résultat de la commande suivante

Code:

---

str(Mesdonnes)# Mes données (les données que vous utilisz

---

Merci

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215

data<-read.table("C:/Users/Miles/Desktop/projet sous r/Femmes_frais_mensuels.txt",sep=";",dec=",",header=TRUE)   Region<-as.factor(data$REGION) Surcomp<-as.factor(data$Surcomplémentaire) Type<-as.factor(data$Type) Age<-as.factor(data$Tranche.d.age) Ancienneté<-as.factor(data$Expr1) Code.activité<-as.factor(data$DC...LibellÃ..code.activitÃ.) Regroup.Poste.soins<-as.factor(data$Regroupement.1)   don<-data.frame(Frais.Mensuels,Quantité.Actes,Region,Surcomp,Type,Age,Ancienneté,Code.activité,Regroup.Poste.soins)   row.names(don)<-row.names(data)   options(max.print=100000000)     don$Region <- ifelse( test=as.character(don$Region) %in% c("Aquitaine", "Midi Pyrénées","Languedoc Roussillon", "Provence Alpes Cote d'Azur","Rhône Alpes","Centre","Auvergne","Auvergne ","Limousin"), yes= "CENTRE_SUD", no = as.character(don$Region)) don$Region <- as.factor(don$Region)     don$Region <- ifelse( test=as.character(don$Region) %in% c("Bretagne", "Pays de Loire","Poitou Charentes"), yes= "OUEST", no = as.character(don$Region)) don$Region <- as.factor(don$Region)       don$Region <- ifelse( test=as.character(don$Region) %in% c("Normandie ( Haute )", "Normandie ( Basse )","Ile de France"), yes= "Region Parisienne + Normandie", no = as.character(don$Region)) don$Region <- as.factor(don$Region)     don$Region <- ifelse( test=as.character(don$Region) %in% c("Nord Pas de Calais", "Picardie","DOM TOM ou autres"), yes= "NORD + DOM TOM et autres", no = as.character(don$Region)) don$Region <- as.factor(don$Region)     don$Region <- ifelse( test=as.character(don$Region) %in% c("Alsace", "Franche Comté","Champagne Ardennes", "Lorraine","Bourgogne"), yes= "EST", no = as.character(don$Region)) don$Region <- as.factor(don$Region)     don$Surcomp <- ifelse( test=as.character(don$Surcomp) %in% c("OUI"), yes= "Surcomplémentaire", no = as.character(don$Surcomp)) don$Surcomp <- as.factor(don$Surcomp)   don$Surcomp <- ifelse( test=as.character(don$Surcomp) %in% c("NON"), yes= "Base", no = as.character(don$Surcomp)) don$Surcomp <- as.factor(don$Surcomp)       don$Ancienneté <- ifelse( test=as.character(don$Ancienneté) %in% c("- 1 mois","1 à 3 mois","3 à 6 mois","6 à 9 mois","9 à 12 mois","12 à 15 mois","15 à 18 mois","18 à 21 mois"), yes= "partiellement", no = as.character(don$Ancienneté)) don$Ancienneté <- as.factor(don$Ancienneté)     don$Ancienneté <- ifelse( test=as.character(don$Ancienneté) %in% c("21 à 24 mois"), yes= "totalement", no = as.character(don$Ancienneté)) don$Ancienneté <- as.factor(don$Ancienneté)     don$Type <- ifelse( test=as.character(don$Type) %in% c("Isolé"), yes= "1 PERSONNE", no = as.character(don$Type)) don$Type <- as.factor(don$Type)   don$Type <- ifelse( test=as.character(don$Type) %in% c("Couple","Isolé * 2enfants","Famille * 1enfant","Isolé * 1enfant"), yes= "2 ou 3 PERSONNES", no = as.character(don$Type)) don$Type <- as.factor(don$Type)   don$Type <- ifelse( test=as.character(don$Type) %in% c("Famille * 2enfants","Famille * 3 ou + enfants","Isolé * 3 ou + enfants"), yes= "+ 3 PERSONNES", no = as.character(don$Type)) don$Type <- as.factor(don$Type)     don$Age <- ifelse( test=as.character(don$Age) %in% c("17 à 22","22 à 27"), yes= "17 à 27", no = as.character(don$Age)) don$Age <- as.factor(don$Age)     don$Age <- ifelse( test=as.character(don$Age) %in% c("27 à 32","32 à 37"), yes= "27 à 37", no = as.character(don$Age)) don$Age <- as.factor(don$Age)     don$Age <- ifelse( test=as.character(don$Age) %in% c("37 à 42","42 à 47"), yes= "37 à 47", no = as.character(don$Age)) don$Age <- as.factor(don$Age)     don$Age <- ifelse( test=as.character(don$Age) %in% c("47 à 52","52 à 57"), yes= "47 à 57", no = as.character(don$Age)) don$Age <- as.factor(don$Age)   don$Age <- ifelse( test=as.character(don$Age) %in% c("57 à 62","62 à 67","Plus de 67"), yes= "57 et plus", no = as.character(don$Age)) don$Age <- as.factor(don$Age)       don$Code.activité <- ifelse( test=as.character(don$Code.activité) %in% c("Agroalimentaire","Constr. navales","Industrie des métaux","Ind. du meuble","Ind. mecaniques","Imprimerie","Papier Verre Carton","Automobiles","Autres industries","Textile Habil.","Chimie/Plasturg","Pharmacie","Energie/Envir.","Batiment/TP","Mat. construct."), yes= "SECONDAIRE", no = as.character(don$Code.activité)) don$Code.activité <- as.factor(don$Code.activité)     don$Code.activité <- ifelse( test=as.character(don$Code.activité) %in% c("Banques","Distribution","Secteur Public","Informatique","Sans objet","Transport","Services","V.I.P."), yes= "TERTIAIRE/AUTRES", no = as.character(don$Code.activité)) don$Code.activité <- as.factor(don$Code.activité)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Consultations - visites","Consultations","Majorations","Indemnités de déplacement","Visites"), yes= "CONSULTATIONS", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Lentilles","Montures","Verres"), yes= "OPTIQUE", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Honoraires","Forfait hospitalier","Chambre particulière","Frais de séjour"), yes= "HOSPI.", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Prothèses","Orthodontie","Soins"), yes= "DENTAIRE", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Analyses médicales","Autres soins","Autres","Actes de radiologie","Actes de spécialité","Soins à l'étranger"), yes= "ANALYSES/RADIO ET ACTES SPE.", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Pharmacie"), yes= "PHARMA.", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)     don$Regroup.Poste.soins <- ifelse( test=as.character(don$Regroup.Poste.soins) %in% c("Transport","Auxiliaires médicaux","Divers","Obsèques","Cure thermale","Maternité","Autres prothèses",""), yes= "TRANSPORT/AUX. ET DIVERS", no = as.character(don$Regroup.Poste.soins)) don$Regroup.Poste.soins <- as.factor(don$Regroup.Poste.soins)   > str(don) 'data.frame': 754403 obs. of 9 variables: $ Frais.Mensuels : Factor w/ 3 levels "1","2","3": 2 3 3 1 3 2 2 2 2 1 ... $ Quantité.Actes : Factor w/ 3 levels "1","2","3": 3 1 3 1 2 2 2 3 2 2 ... $ Region : Factor w/ 5 levels "CENTRE_SUD","EST",..: 1 1 1 1 1 1 1 1 1 1 ... $ Surcomp : Factor w/ 2 levels "Base","Surcomplémentaire": 1 1 1 1 1 1 1 1 1 1 ... $ Type : Factor w/ 3 levels "+ 3 PERSONNES",..: 1 1 1 1 1 2 2 2 2 2 ... $ Age : Factor w/ 5 levels "17 à 27","27 à 37",..: 3 3 3 3 3 5 5 5 5 5 ... $ Ancienneté : Factor w/ 3 levels "12 à 15 mois",..: 2 2 2 2 2 3 3 3 3 3 ... $ Code.activité : Factor w/ 2 levels "SECONDAIRE","TERTIAIRE/AUTRES": 2 2 2 2 2 2 2 2 2 2 ... $ Regroup.Poste.soins: Factor w/ 7 levels "ANALYSES/RADIO ET ACTES SPE.",..: 2 5 6 3 5 1 1 1 7 1 ...   >library(FactoMineR)   > ACM<-MCA(don) Erreur dans La.svd(x, nu, nv) : code d'erreur 1 dans le sous-programme Lapack 'dgesdd'

---

Je me permets de te signaler que tous tes blocs ifelse sont faux.

Pour un rappel des regroupements sans if else, voir ici:

http://www.developpez.net/forums/d13...ent-modalites/

Je m'excuse, j'ai joins le mauvais code (ancien script).
Normalement mon ifelse est de la sorte :

Code:

---

1 2 3 4 5 6

REGION <- ifelse( test=as.character(data$Region) %in% c("Aquitaine", "Midi Pyrénées","Languedoc Roussillon", "Provence Alpes Cote d'Azur","Rhône Alpes","Centre","Auvergne","Auvergne ","Limousin"), yes= "CENTRE_SUD", no = as.character(data$Region)) REGION<- as.factor(data$Region)

---

C'est toujours faux.

Citation:

---

the condition has length > 1 and only the first element will be used

---

N'aurais tu pas ça qui revient à chaque fois par hasard quand tu lances tous tes if ? ^^

Autre erreur : REGION<- as.factor(data$Region) après ton if annule carrément tout ce que tu fais dans ton if.

Voir l'ancien sujet que tu avais posté, tout est expliqué pour regrouper des modalités.

Bonjour,
Je pense qu'à un moment donné, il y a des calculs de décomposition en valeurs singulières dans ton code avec la fonction La.svd(x, nu, nv).
Vous pouvez essayer d'utiliser la fonction svd() de base (package base) en changeant l'argument LINPACK = FALSE en LINPACK = TRUE.
Pour ce faire, vous copiez le contenu de la fonction (tapez juste svd sur votre console) et vous le recollez sur votre console en changeant l'argument LINPACK.
Ou sinon, copiez-collez le code ci-dessous sur votre console et recommencez pour voir ce que ça donne.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54

svd <- function (x, nu = min(n, p), nv = min(n, p), LINPACK = TRUE) { x <- as.matrix(x) if (any(!is.finite(x))) stop("infinite or missing values in 'x'") dx <- dim(x) n <- dx[1L] p <- dx[2L] if (!n || !p) stop("0 extent dimensions") if (is.complex(x)) { res <- La.svd(x, nu, nv) return(list(d = res$d, u = if (nu) res$u, v = if (nv) Conj(t(res$vt)))) } if (!is.double(x)) storage.mode(x) <- "double" if (!LINPACK) { res <- La.svd(x, nu, nv) return(list(d = res$d, u = if (nu) res$u, v = if (nv) t(res$vt))) } if (nu == 0L) { job <- 0L u <- double() } else if (nu == n) { job <- 10L u <- matrix(0, n, n) } else if (nu == p) { job <- 20L u <- matrix(0, n, p) } else stop("'nu' must be 0, nrow(x) or ncol(x)") job <- job + if (nv == 0L) 0L else if (nv == p || nv == n) 1L else stop("'nv' must be 0 or ncol(x)") v <- if (job == 0L) double() else matrix(0, p, p) mn <- min(n, p) mm <- min(n + 1L, p) z <- .Fortran("dsvdc", as.double(x), n, n, p, d = double(mm), double(p), u = u, n, v = v, p, double(n), as.integer(job), info = integer(1L), DUP = FALSE, PACKAGE = "base")[c("d", "u", "v", "info")] if (z$info) stop(gettextf("error %d in 'dsvdc'", z$info), domain = NA) z$d <- z$d[1L:mn] if (nv && nv < p) z$v <- z$v[, 1L:nv, drop = FALSE] z[c("d", if (nu) "u", if (nv) "v")] }

---

Pour plus d'explications, vous pourrez parcourir cette discussion.

Le code premier des ifelse est correct selon moi. Je ne vois pas mon erreur, si ce n'est que que je remplace mes anciennes modalités de don$nom_variable par un nouveau regroupement.

Le str() affiche bien ces nouveaux regroupements ... pourtant :roll:

La preuve :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34

> summary(don) Frais.Mensuels Quantité.Actes Region 1:251529 1:257820 CENTRE_SUD :162824 2:251492 2:252746 EST :153102 3:251382 3:243837 NORD + DOM TOM et autres :103767 OUEST :153066 Region Parisienne + Normandie:181644     Surcomp Type Age Base :610085 + 3 PERSONNES :179137 17 à 27 : 76013 Surcomplémentaire:144318 1 PERSONNE :254254 27 à 37 :224793 2 ou 3 PERSONNES:321012 37 à 47 :229606 47 à 57 :183783 57 et plus: 40208     Ancienneté Code.activité 12 à 15 mois : 22834 SECONDAIRE :292483 partiellement:188969 TERTIAIRE/AUTRES:461920 totalement :542600         Regroup.Poste.soins ANALYSES/RADIO ET ACTES SPE.:184519 CONSULTATIONS :158637 DENTAIRE : 72675 HOSPI. : 53054 OPTIQUE : 86755 PHARMA. : 84719 TRANSPORT/AUX. ET DIVERS :114044

---

Comment ex^miquer ces regroupements, si mon code est faux ? Les totaux correspondent bien en faisant la somme des anciens.

En parallèle, je te réponds arm et te remercie pour cette réponse. Je pense que R m'a renvoyé cette erreur due à une incapacité de traitement de la machine. Après avoir mis en variable supp. Surcomp, j'ai réussi à sortir l'ACM.

Si jamais cette erreur, re-survient alors je prendrais ta fonction.

Merci pour tout. Cependant, je ne comprends toujours pas ma première erreur. Puisque les regroupements se font effectivement bien ...