similar to: [R-pkgs] vcdExtra 0.5-0 is released to CRAN

For the vcdExtra package, I'm exploring methods to generate and fit collections of glm models, and handling lists of such model objects, of class "glmlist". The simplest example is fitting all k-way models, from k=0 (null model) to the model with the highest-order interaction. I'm having trouble writing a function, Kway (below) to do what is done in the example below >

I received the following email note re: the vcdExtra package > A vcd update has shown that packages TIMP and vcdExtra are not > declaring their dependence on colorspace/MASS: see > > http://cran.r-project.org/web/checks/check_results_vcdExtra.html But, I can't see what to do to avoid this, nor understand what has changed in R devel. Sure enough, CRAN now reports errors in

I'm checking a new release of vcdExtra via win builder with R-devel 2014-10-29 r66897 and have run into a Warning I don't know how to fix. I have one S3 generic, mosaic3d() that uses rgl and don't want to have rgl always loaded via Depends:. Instead, the mosaic3d.default() method includes if (!require(rgl)) stop("rgl is required") This always worked in the past, but

In checking my vcdExtra package, the following NOTE newly appeared (R-Forge, using R version 3.0.1 Patched (2013-08-20 r63635)) Package in Depends field not imported from: ?gnm? These packages needs to imported from for the case when this namespace is loaded but not attached. In the DESCRIPTION file, I have Depends: R (>= 2.10), vcd, gnm (>= 1.0.3) In NAMESPACE: # we are a vcd

Context: I have two or more rgl-based views of a given data set, perhaps fitting different models, or showing different things across views. I want to be able to hand-rotate, zoom, scale one view to something I like, and then show the other views with matching viewpoints and scaling. so that one could flip back/forth among graphs and see only the relevant differences. In 2D, all this usually

I'd like to fit and explore a collection of hierarchical loglinear models that might range from the independence model, ~ 1 + 2 + 3 + 4 to the saturated model, ~ 1 * 2 * 3 * 4 I can use add1 starting with a baseline model or drop1 starting with the saturated model, but I can't see how to get the model formulas or terms in each model as a *list* that I can work with further. Consider

Dear all, I have a problem to make figures with two columns in package vcd. Here an example code I take from "\library\vcd\html\plot.loglm.html" What I need, I want to make two figures in one plot. How could I do that. I have tried with layout(rbind(c(1, 1, 2, 2))) but the same result, two plot. Best wishes, Muhammad Subianto library(vcd) oldpar <- par(mfrow=c(1, 2)) ## mosaic

helpeRs, I have a nice looking mosaic plot in an article to be published soon. Sadly, the published version will be in black and white and so ruin the advantage of the default shading scheme of tiles. What would readers suggest as an alternative shading scheme? If I have a black-and-white shading scheme graduated according to suitable cutoffs I won't be able to tell positive from

Dear useRs, a completely revised version of the `vcd' ("Visualizing Categorical Data") package is now available from CRAN. This major revision includes the following enhancements: * grid-based: The package is now entirely based on `grid', the new R graphics system, thus exploiting its unique functionalities. Powered by grid, it is now possible, e.g., to simply compose complex

Dear useRs, a completely revised version of the `vcd' ("Visualizing Categorical Data") package is now available from CRAN. This major revision includes the following enhancements: * grid-based: The package is now entirely based on `grid', the new R graphics system, thus exploiting its unique functionalities. Powered by grid, it is now possible, e.g., to simply compose complex

I have been tinkering with mosaicplot() and friends as a way of learning R. As part of this, I've written a pair.table() method for mosaic matrices, and would like to extend mosaicplot to work with loglin and logln (MASS) objects. I'm using R 1.4.0 on Win 98. I've been trying to figure out the formula interface, and think there's a bug, but not sure how to find it, yet alone fix

I posted on this topic to r-help, but never got a sufficient answer, so I'm reposting here. [Env: R 2.11.1, Win Xp, using Eclipse/StatET] In a package I'm working on, I want to create as.matrix() and as.array() methods for a particular kind of object (log odds ratios). These are returned in a loddsratio object as the $coefficients component, a vector, but really reflect an underlying

Hello: What is the most efficient way to change the plotted variable names in mosaic plots in the vcd package? Should one do a separate contingency table first, change the dimension names there and then pass that to mosaic? Or is there a way to do it simply within mosaic. I was thinking something like: mosaic(~var1+var2, labelling_args=list(varnames=c('newvar1', 'newvar2')) Simon

The NEWS for R 2.3.0 states that "R CMD check works for packages whose package name is different from the directory name in which it is located." However that hasn't been my experience. I ran R CMD check on package sources located in a directory with the same name as the package and it worked as expected. Then I renamed the directory and tried again. The first attempt got stuck

Hi Milan, I expect I may be able to do something about the way the terms are evaluated, to ensure the evaluation is done in the gnm namespace (while still ensuring the variables can be found!). In the meantime, I think the following will work: Mult <- gnm::Mult f <- Freq ~ Eye + Hair + Mult(Eye, Hair) gnm::gnm(f, family=poisson, data=dat) Hope that helps, Heather On Wed, Apr 29, 2015,

Hi! Some time ago, I replaced calls to library() with calls to requireNamespace() in my package logmult, in order to follow the new CRAN policies. But I just noticed it broke jackknife/bootstrap using several workers via package parallel. The reason is that I'm running model replicates on the workers, and the formula includes non-standard terms like Mult() which are provided by gnm. If gnm

Hi all I want to construct a generalised nonlinear model (binomial family) using gnm, of the form: Response = a + b variable1 + c variable2 + d variable3 - d b variable4 - d c variable5, with the parameters b, c, and d appearing more than once. Hence, I think I need to use a custom nonlin function with gnm. One of my predictor variables is categorical, so I have created a dummy variable for

[Env: R 2.12.2, WinXp] In a vignette for the vcdExtra package, I had a text file, tv.dat under data/, that I used in the vignette as <<tv1,results=verbatim>>= tv.data<-read.table(system.file("data","tv.dat",package="vcdExtra")) head(tv.data,5) @ I was told that this now generates a warning for non-Rdata files in R CMD check. But I'm now

The null loglinear model is an intercept-only model for log frequency, log(f) = \mu For a one-way table the test of the null model is the same as the chisq.test. This can be fit using loglin(), but I don't think there is any way to specify this using MASS::loglm > t1<- margin.table(Titanic,1) > t1 Class 1st 2nd 3rd Crew 325 285 706 885 > loglin(t1, NULL) 0

loglm() is a port of an original written for S-PLUS. The fact that it carries only a reference to the data frame is neither intentional nor unintentional, but an unnoticed side-effect. I can see advantages both ways. (I'm not so sure, either, that what you say is standard behaviour for model fitting functions really is so universal.) It would not be too hard to come up with a version that