similar to: R and Data Storage

Thompson's Manual to Accompany Agresti's book refers to a package named "repeated". It's not on CRAN from what I can see. I have seen rpm's for it. Where is the best place to download this package? Rick B.

I fit the following model using glmmPQL from MASS: fit.glmmPQL <- glmmPQL(ifelse(class=="Disease",1,0)~age+x1+x2,random=~1|subject,family=binomial) summary(fit.glmmPQL) The response is paired (pairing denoted by subject), although some subjects only have one response. Also, there is a perfect positive correlation between the paired responses. x1 and x2 can and do differ within each

> rab45+ at pitt.edu wrote: >> Jerry, >> >> Here are the two R packages. >> >> Thanks. >> >> Rick > All the fortran codes compile OK so the problem is invoking the compiler > you > want. One approach which is a work around is to rename the gfortran > executable > to somthing else like gfortran-save and create a link called gfortran that

Is there an R package available that would facilitate doing a power/sample size analysis for linear mixed effects models? I have seen the Java applets made available by Russell Length which would seem to be able to handle most any lme, but there is little documentation and it's not clear how the models need to be formulated. Rick B.

I'm estimating models using lme in a for loop. Sometimes the model doesn't converge or there are other problems. This causes the evaluation to stop prematurely. I can't remember the function name that I need to use to allow the loop to continue until the end. Could someone remind me the name of the function? I've tried searching but haven't hit upon the function. Rick B.

The stat.table function in the Epi package won't do standard deviations. It didn't seem that it would be difficult to add an "sd" function to the stat.table function. Following the example for the mean, I set up a similar function for the sd (and included it as an options) but it just won't work. (I tried sending messages to the Epi mailing list after subscribing but my mail

Has anyone experimented with the Network Common Data Format (netcdf) described at http://www.unidata.ucar.edu/packages/netcdf for archiving data? Is this format widely used outside of atmospheric research? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I'm fitting a mixed effects model: fit.1 <- lme(y~x,random=~1|id,data=df) There are two different observations for each id for both x and y. When I use plot(fit.1), there is a strong increasing linear trend in the residuals versus the fitted values (with no outliers). This also happens if I use random=~x|id. Am I specifying something incorrectly? Rick B.

Hi, we have a first draft of R functions reading/writing data to XML files including a rather general DTD ... which borrows heavily from the data types of a certain programming language :-) The basic idea is to create an XML standard for data exchange, together with import/export functions for as many applications as possible. We here will need R, Matlab & Octave for our research program,

I'm using groupedData from nlme. I set up a groupedData data.frame with outer=~group1. When I try to plot with outer=TRUE, I get "subscript out of bounds." This happens most of the time. When it works, I get spaghetti-type plots for comparing groups. But I don't understand why it doesn't usually work. > longa.mod.1.gd <- groupedData(mod1.logit~time|

When I try to run the example from Variogram with an lme object, I get an error (although summary works): R : Copyright 2005, The R Foundation for Statistical Computing Version 2.2.1 (2005-12-20 r36812) ISBN 3-900051-07-0 ... > fm1 <- lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat) Error: couldn't find function "lme" > Variogram(fm1, form = ~ Time | Rat, nint =

Hi folks, I've been dickering around again, wasting the time away.... and today Brian Zisk (chief Vorbis Evangelist among other things) said, "Hey! We need to demo to the world on the 4th of April!" Hooo boy. Time to suck it up and bang it out. The Vorbis libs are not a problem. I'm in process of checking in my new psychoacoustics (they're real pretty!) based on

Given that we occasionally run into problems with comparing Excel results to R results, and other spreadsheet-induced errors, I thought this might be of interest. http://www.nextnewdeal.net/rortybomb/researchers-finally-replicated-reinhart-rogoff-and-there-are-serious-problems The punchline: "If this error turns out to be an actual mistake Reinhart-Rogoff made, well, all I can hope is that

Hi, I was wondering if Statdataml is currently the preferred way to represent statistical data in XML in R. And also if the Statdataml api provides ways to load the XML as a HTTP GET? If so can you give me an easy example of this. Thanks. Cheers, Bryan Rasmussen

Hi, I just had a very quick look at the StatDataML proposal --- nice work! At the risk of showing my ignorance, I want to mention my first impressions. My first impression is that defining datasets in terms of arrays and list is a bit too high a level. What about simpler vectors, scalars? (I know that R/S don't have scalars, but other systems/applications do.) Can we think of a core

Thomas (I'm getting a bit confused about whether this discussion is taking place on r-devel or omega-devel, and perhaps everyone is on both lists) I just posted something on R-devel which I hope addresses many of the points you mentioned, in particular pointing out that I meant encryption for the purpose of authentication, not for the purpose of making something secret, and I was using the

Hi all, I have a very large music collection that I keep on a portable hard drive (to plug in to car USB, carry with when I'm at my office, etc). All but a few dozen files are part of an album, and not a single audio file. This translates to an insane amount of files stored on my hard drive. I am very strict about how I label and file my media (I loathe unlabeled audio tracks), but

I thought I would take a look at the StatDataML package for some work I'm doing, but I receive the following error when starting: > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 0.1 year 2004 month 11

Hi, Has anyone out there written a function to take a data.frame as input and generate XML that conforms to the DTD for ggobi ("ggobi.dtd")? In other words, like a simple version of the writeSDML function in the StatDataML package, but using ggobi.dtd instead of StatDataML.dtd. It looks easy to write such a function to handle data.frames with only numeric data, but a bit of work with

how about -allow-unused-prefixes=needs_review On Mon, Nov 9, 2020 at 3:54 PM David Blaikie <dblaikie at gmail.com> wrote: > I don't feel too strongly about the "=true" bit for legitimate cases - > I'd prefer not to have it, but it's not the end of the world. > I do feel more strongly that if you're going to automatically add it > to cases that might