similar to: silent loading of packages

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

Hi, I have a S3 package with namespace called "myS3Pkg". Inside my package I would like to create a S4 class which extends (adds 2 slots) another S4 class from some other package. The class should be created in "myPkg" environment (and not global environment). Using: setClass("myS4class", representation("otherS4class", mydata = "numeric"),

I have a question about "deparse" function in R What is the reason that "deparse" use an argument like "width.cutoff" ? Why the maximum cutoff is 500? I was manipulating an R formula and used "deparse". Since the length of user's formula was greater then 500, my code didnt work. thanks Johan johan Faux <johanfaux@yahoo.com> wrote: I have a

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 parameters. I think the problem is that dnbinom allows either

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 is, is there a single nonlinear function f(x) which acts

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

I would like R CMD build to run some R code which does some stuff and save the result as a file in /inst/docs folder. Is there any way of doing this. Thank you. Johan ____________________________________________________________________________________ We won't tell. Get more on shows you hate to love [[alternative HTML version deleted]]

Hello , a<-c(1) for(i in 2:length(a)) do.something with a[[i]] I get : Error in a[[i]] : subscript out of bounds Am I missing something here? Doesnt R check the value of i inside "for" and if the condition is not tru, dont do anything ???? thanks, johan --------------------------------- Got holiday prints? See all the ways to get quality

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.

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 thought namespaces should protect us from this happening? Kjetil --

Dear all , I am 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

Hi, I am trying to fit a line in the log plot of my networks degree distribution to show that it is a power-law distribution. I am using the following commands. However, I am not able to see the fitted line. Any comments to help? I am using following packages: igraph, splines,base,VGAM, netmodels. g is my network, d is the degree of nodes in the network, and dd is the degree distribution d

Hi Folks, after some days of intensive google-work I didn't find a solution for joining Windows 10 Clients (release version) to Samba 3.6.6 LDAP based NT Domain. I know, that I should switch to Samba4 but we got hundreds of scripts, services, devices and other stuff (our coffeemachine estimates the preferred strength for each user ... mission critical ;) ) that use the ldap service for

Hi, I was going through the examples in smartpred. It seems there's an unstated dependency on the fitted() function in package flexmix. n = 20 set.seed(86) x = sort(runif(n)) y = sort(runif(n)) library(splines) fit = lm(y ~ ns(x, df=5)) plot(x, y) lines(x, fitted(fit)) # won't work w/o prior loading of the flexmix package. newx = seq(0, 1, len=n) points(newx, predict(fit,

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))

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 (IRLS). The key ideas are Fisher scoring, generalized linear and

I've been running a glmer logit on a very large data set (600k obs). Running on a 10% subset works correctly, but for the complete data set, R completes apparently without error, but does not display the results. Given these jobs take about 200 hours, it's very hard to make progress by trial and error. I append the code and the sample and complete output. As is apparent, I upgraded R

I believe there is a coding error in the first part of predict.mlm in the splines package, but perhaps someone could explain the logic to me if I'm wrong. I don't see how the second if (missing(newdata)) could ever be true. (I would show the code but I'm using email on a different machine than R today.) Paul Gilbert

Dear R developers: I am trying to fit a PPO model using vglm from the library VGAM, and get an error while executing the code. Here is the data, code, and error: Data: first row is the column names. a = age, and 1,2,3, 4 and 5 are condition grades. a 1 2 3 4 5 1 1 0 0 0 0 2 84 2 7 10 2 3 16 0 6 6 2 4 13 0 3 4 0 5 0 0 0 1 0 Library(VGAM)

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))