search for: vgam

The symptom triggering this email is that an S4 summary method sometimes refuses to be invoked, even when a package is explicitly loaded, if the first load of the package is implicit. It may or may not be specific to 'summary' methods and/or the 'VGAM' package. I've sent to R-devel because (i) it looks like some kind of bug to me, but I'm not sure; (ii) it's not something I personally need any help with; and (iii) it seems a bit specialized for R-help. Here's the case notes. I have an object 'nf1' of S4 class 'vg...

Hello, I am trying to build an R package for Windows from an R package for Linux. There seems to be a "windres" problem somewhere. Here's what I get. C:\Temp>Rcmd build --binary VGAM > mm make[2]: windres: Command not found make[2]: *** [VGAM_res.o] Error 127 make[1]: *** [src/VGAM.dll] Error 2 make: *** [pkg-VGAM] Error 2 In more details, it goes like the following. C:\Temp>Rcmd build --binary VGAM * checking for file 'VGAM/DESCRIPTION' ... OK make: Enteri...

I just discovered that when the VGAM package (not on CRAN) is loaded, glm() doesn't work. This is because VGAM defines a family function() which gets found by glm() in place of the family function from stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I though...

Hello R-Users, I have recently run into several problems using vgam() in the VGAM package. I am hoping someone might have some solutions... Briefly, I have been trying to fit GAM models for zero-altered negative binomial models. 1. When fitting smoothed parameters (e.g. s(X, df=2)) changing the degrees-of-freedom has no effect on the level of smoothing (e.g....

Dear R-help gurus (and T.Yee, the VGAM maintainer) - I've been banging my head against the keyboard for too long now, hopefully someone can pick up on the errors of my ways... I am trying to use optim to fit a zero-inflated negative binomial distribution. No matter what I try I can't get optim to recognize my initial paramete...

Dear, When i multiply the y of a regression by 10, I would expect that the coefficient would be multiply by 10 and the z value to stay constant. Here some reproducible code to support the case. *Ex 1* library(mvtnorm) library(VGAM) set.seed(1) x=rmvnorm(1000,sigma=matrix(c(1,0.75,0.75,1),2,2)) summary(vglm("y~x",family=studentt2,data=data.frame("y"=x[,1],"x"=x[,2]))) summary(vglm("y~x",family=studentt2,data=data.frame("y"=x[,1]*10,"x"=x[,2]))) summary(vglm("y~x...

Hello R users, I have a puzzle with the VGAM package, on my first excursion into generalized additive models, in that this very nice package seems to want to do either more or less than what I want. Precisely, I have a 4-component outcome, y, and am fitting multinomial logistic regression with one predictor x. What I would like to find out i...

Hi! R users I would like to ask you where could we find the VGAM Package. I don't find it in the list of packages. Thak you for your help Lassana KOITA Etudes de S?curit? et d'Exploitation a?roportuaires / Safety Study & Statistical analysis Service Technique de l'Aviation Civile (STAC) / Civil Aviation Technical Department Direction G?n?rale...

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é...

Dear Prof. Thomas Yee I$B!G(Bm very interested in your R program VGAM. I tried below your data: # Nonparametric proportional odds model data(pneumo)pneumo = transform(pneumo, let=log(exposure.time))vgam(cbind(normal,mild,severe) ~ s(let), cumulative(par=TRUE), pneumo) However, the results by Version of VGAM are different; ----------The result by Version 0.7-7...

I would like to turn off all the messages during library(aPackage) or require(aPackage) I tried different commands: invisible, capture.output, sink but none of them is working. For example, loading VGAM, gives a lot of unnecessary messages: > library(VGAM) Attaching package: 'VGAM' The following object(s) are masked from package:splines : bs The following object(s) are masked from package:splines : ns The following object(s) are masked...

Dear useRs, upon request, the VGAM package (currently version 0.7-1) has been officially released on CRAN (the package has been at my website http://www.stat.auckland.ac.nz/~yee/VGAM for a number of years now). VGAM implements a general framework for several classes of regression models using iteratively reweighted least squares (I...

...trying to install a package by name 'netmodels' keep on getting the following error : > install.packages("netmodels") Warning in install.packages("netmodels") : argument 'lib' is missing: using '/usr/lib64/R/library' also installing the dependency ‘VGAM’ trying URL 'http://cran.stat.ucla.edu/src/contrib/VGAM_0.7-10.tar.gz' Content type 'application/x-tar' length 1608488 bytes (1.5 Mb) opened URL ================================================== downloaded 1.5 Mb trying URL 'http://cran.stat.ucla.edu/src/contrib/netmodels_0.2...

Hello everyone, using the VGAM package and the following code library(VGAM) bp1 <- vglm(cbind(daten$anzahl_b, daten$deckung_b) ~ ., binom2.rho, data=daten1) summary(bp1) coef(bp1, matrix=TRUE) produced this error message: error in object$coefficients : $ operator not defined for this S4 class I am bit confu...

Hi Folks, I wonderif someone who is familiar with the details of vglm in the VGAM package can assist me. I'm new to using it, and there doesn;t seem much in the documentation that's relevant to the question below. Say I have a vector x of 0/1 responses and another vector y of 0/1 responses, these in fact being a bivariate set of 0/1 responses equivalent to cbind(x,y)....

...English_United Kingdom.1252;LC_CTYPE=English_United Kingdom.1252;LC_MONETARY=English_United Kingdom.1252;LC_NUMERIC=C;LC_TIME=English_United Kingdom.1252 attached base packages: [1] stats4 splines stats graphics grDevices utils datasets methods base other attached packages: [1] VGAM_0.7-5 xlsReadWrite_1.3.2 -------------------------------- I'm having some problems with a Poisson-lognormal density (mass?) function. VGAM has the dpolono function, but that doesn't work for x-values over 170, and I need to go to *much* bigger numbers. It fails first because of ga...

...nto the integrand "feature" mentioned in the note. I resolved it by limiting the range of integration as shown here: -------------------------------------------------- function (x, meanlog = 0, sdlog = 1, sdlim=6, rel.tol = .Machine$double.eps^0.5,...) { #+ # K. Jewell Feb 2008 # Based on VGAM::dpolono, modified to cover wider range of x # Some simplification, but major change is to avoid integration problems by # limiting range of integration to sdlim (default = 6) standard deviations each side # of mean of Poisson or lognormal. Also increase default precision of integration. # For...

Anyone know how to get p-values for the t-values from the coefficients produced in vglm? Attached is the code and output ? see comment added to output to show where I need p-values + print(paste("********** Using VGAM function gamma2 **********")) + modl2<- vglm(MidPoint~Count,gamma2,data=modl.subset,trace=TRUE,crit="c") + print(coef(modl2,matrix=TRUE)) + print(summary(modl2)) [1] "********** Using VGAM function gamma2 **********" VGLM linear loop 1 : coefficient...

sorry...its the same post but this time with subject dear all: I want try the vgam function in VGAM library, but I'm stucked since I can't define the weights. This is because I dont understand well the concept of "matrix-band" (even after read the reference card). I was wondering if someone have some tutorial examples of how to define weights on a matrix form,...

Hi all: I want to fit a multinomial logit model with VGAM package, however I cannot find a way to check the residuals since residuals(my_model) and resid(my_model) does not work. Any suggestions? Thanks in advance. Alex Brito [[alternative HTML version deleted]]