Hi, R-users Could someone help me to understand this following error. I'm using vglm function in VGAM package Best regards and thank you for your ehlp ######## mydata <- read.table("Data2_overruns.csv", sep =";", header = T, row.names=NULL) > attach(mydata) > > y <- mydata$cat.event > phase.vol <-mydata$phase.vol > pilote <- mydata$pilote > défail.méca <- mydata$défail.méca > mauv.visib <- mydata$mauv.visibilité > piste.cont <- mydata$piste.contam > vent <- mydata$vent > tiers <- mydata$tiers > y <- ordered(y,levels=c("Inci","Inci G","Acci")) > > dat <- data.frame(y,phase.vol,pilote, défail.méca, + mauv.visibilité, piste.cont, vent, tiers)
> pilote_sum<-apply(tab_cont,1,sum) #somme sur les lignes > Event_sum<-apply(tab_cont,2,sum) # somme sur les colonnes > library(VGAM) Le chargement a nécessité le package : splines Le chargement a nécessité le package : stats4 Attachement du package : 'VGAM' The following object(s) are masked from package:splines : bs, ns The following object(s) are masked from package:stats : glm, lm, poly, predict.glm, predict.lm, predict.mlm The following object(s) are masked from package:base : scale.default > result <- vglm(y ~ phase.vol + pilote + défail.méca + mauv.visib + piste.cont + vent + tiers,data = dat, family = cumulative(parallel = F)) Erreur dans dotC(name = "tapplymat1", mat = as.double(mat), as.integer(nr), : NA/NaN/Inf dans un appel à une fonction externe (argument 1) De plus : Warning messages: 1: using type="numeric" with a factor response will be ignored in: model.response(mf, "numeric") 2: production de NaN in: log(x) 3: fitted values close to 0 or 1 in: Deviance.categorical.data.vgam(mu = mu, y = y, w = w, residuals = residuals, 4: production de NaN in: log(x) 5: fitted values close to 0 or 1 in: Deviance.categorical.data.vgam(mu = mu, y = y, w = w, residuals = residuals, > Lassana KOITA Chargé d'Etudes de Sécurité Aéroportuaire et d'Analyse Statistique / Project Engineer Airport Safety Studies & Statistical analysis Service Technique de l'Aviation Civile (STAC) / Civil Aviation Technical Department Direction Générale de l'Aviation Civile (DGAC) / French Civil Aviation Headquarters Tel: 01 49 56 80 60 Fax: 01 49 56 82 14 E-mail: [EMAIL PROTECTED] http://www.stac.aviation-civile.gouv.fr/ [[alternative HTML version deleted]]
______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.