similar to: An artifact of base being namespace

I have some questions about S4 style methods and generics. First of all, is there any way of using default values for arguments in the generics/methods? My own experiments show that such arguments are always ignored. The green book seems to be silent on this matter. The second question is about using a non-trivial function body for generics. Page 351 of the green book gives an example of this

I've noticed quite a few messages relating to large data sets bedeviling R users, and having just had to program my way through one that actually caused a "Bus error" when I tried to read it in, I'd like to ask two questions. 1) Are there any facilities for aggregation of data in R? ( I admit that this will not do much for the large data set problem immediately) 2) Is there any

Dear List, RMySQL is old and wouldn't install on R 2.0.1. How should the source be rewritten to be accepted by R 2.0.1. With thanks Fredrik Lundgren

Dear List, While trying to define a customised CITATION file for a package, following R-exts, I realised that if I use only one 'citEntry' I got both a text description and a BibTex entry for the package, as for the 'nlme' package: -------------------------- citation('nlme') To cite package 'nlme' in publications use: Jose Pinheiro, Douglas Bates, Saikat

If .../<pkg>/man/ contains a filename beginning with [, such as [.datetimes, Rdconvlib.pl cannot convert that help file into latex at INSTALL time. [Using --no-latex is a workaround.] The error given is: # R INSTALL ssCatalogues : help >>> Building/Updating help pages for package `ssCatalogues' Formats: text html latex example /^[.datetimes$/: unmatched [] in regexp at

Hello, I am interested in using nlme to model repeated measurements, but I don't seem to get good CIs. With the code below I tried to generate data sets according to the model given by equations (1.4) and (1.5) on pages 7 and 8 of Pinheiro and Bates 2000 (having chosen values for beta, sigma.b and sigma similar to those estimated in the text). For each data set I used lme() to fit a model,

Forgive me if this isn't relevant, but: are you stripping off the output you need (e.g., parameter estimates) and throwing away what you don't (e.g., residuals)? Reid Huntsinger -----Original Message----- From: Saikat DebRoy [mailto:saikat at stat.wisc.edu] Sent: Wednesday, March 21, 2001 12:15 PM To: Han Lai Cc: R-help Subject: Re: [R] Is there any help for "Cannot allocate vector

Saikat DebRoy is working on a revised "foreign" package (currently in the src/contrib/Devel area on CRAN) that will include the ability to read the SPSS save format. Since we don't use SPSS we don't have easy access to a sample data set and have been using PSPP (see http://pspp.stat.wisc.edu/) to generate test data. It would be more realistic to have a small, redistributable

Did not seem to reach Saikat DebRoy, this might be the forum. Bendix -----Original Message----- From: BXC (Bendix Carstensen) Sent: 24. november 2008 15:00 To: 'saikat at stat.wisc.edu'; 'stvjc at channing.harvard.edu' Subject: read.ssd It's always annoyed me that that read.ssd crashed on datasets with long variable names, but the other day a collegue of mine pointed out to

Dear list, i am a student of psychology and have to do a multilevelanalysis on some data. About that i have one general and one specific question. This is what i have copied from the help-file on lme: data(bdf) fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf, random = ~ IQ.ver.cen | schoolNR) summary(fm) after summary(fm) i get the following error:

> library(foreign) > da<-read.dta(file.choose()) > da startdat starttim enddate endtime days hoursmin secused 1 2005-01-11 2 2005-12-15 20.19 NA NA 9 > attributes(da) $datalabel [1] "Example of use of date and time functions" $time.stamp [1] "15 伿伄伓侢佋伮 2005 20" $names [1] "startdat" "starttim"

Full_Name: W.B.Kloke Version: 1.8.0 OS: FreeBSD-4.7 Submission from: (NULL) (195.253.22.63) I found a parameter evaluation inconsistency in NLME package. I tried to use gls() inside a function, and I wanted use gls() for different subsets of a data frame: prgls <- function(name){ gls( log10(Y)~(cond-1)+(cond-1):t ,pr,subset=subject==name)} Applying this function with a string as parameter

Dear all, Given a LME model (following the notation of Pinheiro and Bates 2000) y_i = X_i*beta + Z_i*b_i + e_i, is it possible to extract the variance-covariance matrix for the estimated beta_i hat and b_i hat from the lme fitted object? The reason for needing this is because I want to have interval prediction on the predicted values (at level = 0:1). The "predict.lme" seems to

Hi, I need to write a function that would look something like this: S <- function(b=betas){ expression(b[1] * f(b[2] * x * f(b[3] * x * f(...b[n-1] * x * f(b[n] * x)))...) } Where n is the number of element in b. Further I need to be able to evaluate S at some x numerically of course and I need to use "deriv" and produce dS/dx such that I can evaluate it also at some x. I

Hello- I understand that it's convention, when comparing two models using the anova function anova(model1, model2), to put the more "complicated" (for want of a better word) model as the second model. However, I'm using lme in the nlme package and I've found that the order of the models actually gives opposite results. I'm not sure if this is supposed to be the case

Dear all, I would like to (i) produce boxplot graphs with axis in logarithm in base 10 and (ii) showing the values on the axis in 10^exponent format rather than 10E+exponent. To illustrate with an example, I have some widely spread data that I chart plot using boxplot() [figure on the left]; the log="y" option of boxplot() I obtained the natural logarithm conversion of the data and

Dear All, Using:- R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 and Jose Pinheiro, Douglas Bates, Saikat DebRoy and Deepayan Sarkar (2005). nlme: Linear and nonlinear mixed effects models. R package version 3.1-65. on a WINDOWS 2000 machine I am trying to run the scripts from the Mixed Models book and am running into some

Even though the FDA has no policies at all that limit our choices of statistical software, there is one defacto standard in place: reliance of the SAS transport file format for data submission (even though this format is deficient for this purpose, e.g., it does not even document value labels or units of measurement in a self-contained way). Because of the widespread use of SAS transport files in

Dear all, I try to load the nlme library for Linear and nonlinear mixed effects models by Jose Pinheiro, Douglas Bates and Saikat DebRoy but I get the following error message: > library(nls) > library(nlme) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "C:\Program\R\rw1010/library/nlme/libs/nlme.dll": LoadLibrary failure: Det

Recent changes to read.spss() in the foreign package return a dataframe containing additional attributes. For example, >TEMP<-read.spss(choose.files(), to.data.frame=T,use.value.labels=F) > str(TEMP) `data.frame': 780 obs. of 8 variables: $ EXPOS01: atomic 1 1 2 1 2 3 2 4 2 1 ... ..- attr(*, "value.labels")= Named num 5 4 3 2 1 .. ..- attr(*,