similar to: restrictions for attribute access <Watchdog: Virus checked>

Dear R devels, let's say I have three packages "pkg1", "pkg2" and "pkg3" which all contain new S4 classes and methods. Where "pkg3" depends on "pkg2" and "pkg2" depends on "pkg1". Moreover, all three packages have namespaces. 1) I use ".onLoad <- function(lib, pkg) require(methods)". Do I also have to

The current coxme code has functions that depend on bdsmatrix and others that depend on Matrix, both those pacakges define S4 methods for diag. When loaded, the message appears: replacing previous import ?diag? when loading ?Matrix? Questions: 1. Do I need to worry about this? If so, what can I do about it? I suppose I could add an importFrom directive, but it will be a pain unless there

I have a question about what to do in the following situation (please bear with the setup): Package A defines an S4 generic 'foo' and as well as S4 methods for 'foo' and has exportMethods("foo") in its NAMESPACE file. Package B defines another method for 'foo' for class "bar" and has importFrom(A, "foo") exportMethods("foo")

It is the old story of defined behaviour and expected outcomes. Hard to change now. So I would suggest you do something like this in your ~/.Rprofile: R> smry <- function(...) summary(..., digits=6) R> smry(155555L) Min. 1st Qu. Median Mean 3rd Qu. Max. 155555 155555 155555 155555 155555 155555 R> Maybe call it Summary() instead. Dirk --

I have encountered an issue which I have been unable to resolve, involving an S3 generic (print) being declared S4 generic in a package, and the method being exported. This all works fine - the problem occurs when I try to import the method to another package. Here is the bit that works fine. ------------- #the .r file for package bar

Dear list, I have another (again possibly boneheaded) puzzle about importing, again encapsulated in a nearly trivial package. (The package is posted at <http://www.math.mcmaster.ca/bolker/misc/coefsumtest_0.001.tar.gz>.) The package consists (only) of the following S3 method definitions: coeftab <- function(object, ...) UseMethod("coeftab",object) coeftab.default <-

Hello, Two lapply questions (system info and sample data below): 1) Why does the first form of command1 add the name of y _after_ the str() output rather than before as does the second (preferred) form? # command1 version1 invisible(lapply(ls(pattern='bn'), function(y) cat(y, "\n", str(get(y)), "\n") )) # command1 version2 (preferred output)

Hi folks, I'm not sure if this is a bug or not so I thought I'd check here first. I came across it while working on an update to my package where I try to get AICtab slot from the summary of an lmer object ( summary(my_lmer)@AICtab ). The attached contains a minimal example, where the code below will work if you load the f() function by sourcing the package's R code but will throw an

I am using the AICmodavg package and using R version 2.15.1. I need help with code that is instead being read as text. Below is a subset of code... I actually have 12 models, but I am trying to get this to work for 2 below right now. Everything 'appears' to work through the line starting with Modnames. After that the code starting with aictab and beyond is recognized as text and not

Hola buenas noches, tengo un problema bastante gordo. ¿A alguno le ha dejado de funcionar las funciones pamer.fnc y mcp.fnc con la nueva versión de R? La semana pasada formatee el ordenador y ahora scripts antiguos no funcionan. La cuestión es que me precupa que no funcione el ejemplo de tutorial del autor. Os dejo un script que debería de funcionar y no lo hace

Hi R community! I've fitted three mixed-effects regression models to a thousand bootstrap samples (case-resampling regression) using the lme4 package in a custom-built for-loop. The only output I saved were the inferential statistics for my fixed and random effects. I did not save any output related to the performance to the machine learning algorithm used to fit the models (REML=FALSE).

Using the Cement hardening data in Anderson (2008) Model Based Inference in the Life Sciences. A Primer on Evidence, and working with the best model which is lm ( y ~ x1 + x2, data = cement ) the AIC value from R is model <- lm ( formula = y ~ x1 + x2 , data = cement ) AIC ( model ) 64.312 which can be converted to AICc by adding the bias

I'm having problems using Microsoft Access 2000 mdb files on a Samba share. The same database worked very well on a Windows server. I had hoped Samba would be a drop in replacement for a Windows server but apparently Access needs special Samba configuration. In Access I sometimes get Run-time error '3035': System resource exceeded. I don't see any errors in the log file when

Dear R users I have a question related to the interpretation of results based on GAMMs using Simon Woods package gamm4. I have repeated measurements (hours24) of subjects (vpnr) and one factor with three levels (pred). The outcome (dv) is binary. In the first model I'd like to test for differences among factor levels (main effects only): gamm.11<-gamm4(dv ~ pred +s(hours24), random = ~

I want to call summary on a mer object (from lme4) within my package but I can't seem to get the namespace to import the necessary method. I've simplified my package to this one function: --------------------------- ss <- function(m) { summary(m) } --------------------------- And my namespace file looks like this, where I've attempted to follow the instructions in "Writing

Hello, after upgrading from 3.0.8 to 3.0.14a I started have problems in my network. Environment: - clients with W2KSP4 - Linux Debian Sarge server with Samba acting as PDC. There was no other changes on server configuration. Problems: 1) Only some users cannot reach their home share automatically after login. I am using "login drive = J:" to auto-setup connection after login. 2) Only

>>>>> Dirk Eddelbuettel <edd at debian.org> >>>>> on Fri, 19 Aug 2016 11:40:05 -0500 writes: > It is the old story of defined behaviour and expected outcomes. Hard to > change now. yes... not impossible though... see below > So I would suggest you do something like this in your ~/.Rprofile: R> smry <- function(...)

I've been trying with no success to model mixtures of Gamma distributions using the package flexmix (see examples below). Can anyone help me get it to model better? Thanks very much. -Ben ## ## Please help me get flexmix to correctly model mixtures of ## Gamma distributions. See examples below. ## library('flexmix') ## ## Plot a histogram of dat and the Gamma mixture model given

Dear all, I have the following question. I have to fit the hierarchical model for the hypothesis concern the individual-level effects by controlling for the individual -level attributes and national-level contextual effects on individuals by using R. O have to obtain the estimates of the impact of the second-level (national: GDP per capita) effects on individuals ( in this instance the impact

Dear list, I am using optim() function to MLE ~55 parameters, but it is very slow to converge (~ 25 min), whereas I can do the same in ~1 sec. using ADMB, and ~10 sec using MS EXCEL Solver. Are there any tricks to speed up? Are there better optimization functions? Thanks Toshihide "Hamachan" Hamazaki, $B_@:j=S=((JPhD Alaska Department of Fish and Game: