similar to: nlme

To all: I currently download the R binaries for Redhat 7.x Linux. There is considerable turmoil in the vendors of Linux. Redhat apparently is changing it's business model to paid versions. This might motivate my department to use a different vendor of Linux. Is there anything predictable about which vendors/versions of Linux will have R binaries in the future? Thanks, Nathan PS I

Running under Redhat 7.3 Linux I have installed R Version 1.8.0 (2003-10-08) from a binary download using the rpm manager. Using packageStatus(), I am attempted to upgrade several packages. This is unsuccessful with the following error messages: * Installing *source* package 'MASS' ... ** libs gcc -I/usr/lib/R/include -I/usr/local/include -D__NO_MATH_INLINES -mieee-fp -fPIC -O2

Hi, I using the boot package 1.2-20 on R 2.0.1. My statistics function estimates 6 parameters. In a small percentage of resampled data sets my statistics function doesn't produce an estimate for one parameter and the boot function stops with an error. I can write an ifelse(exists('parameter.estimate'), parameter.estimate, NA) statement within the statistic function to substitute

Hi All, I'm new at programming in R. Some functions create objects for which extractor functions are written to pull out some partial result. I wish to extract partial results from functions without extractor functions: for example, to pass a vector of results to another function. Even after looking at V&R S Programming I don't see a general approach. Specifically, how would I

> Hello, > > I'm running R 1.8.1 on a Mac G4 with Dual 800 MHz CPU, 1.5 GB RAM, and > sufficient disk space. > > I have loaded the boot package version 1.2-15. > > I defined a function to pass to the statistic argument which estimates > a vector (length = 6) of parameter estimates; none are variances. > > I defined a function to pass to the ran.gen argument

Hi All, I have a four panel xyplot. I wish to plot each point as an open or filled circle depending on the value of an indicator variable. I assume I need to use panel.superpose(), but I can't figure out the syntax from lattice documentation. Running R 2.1 under Mac OS X 10.4.2. Any suggestions would be appreciated. Nathan Nathan Leon Pace, MD, MStat University of Utah Salt Lake

Hi, Has anyone written code for isotonic regression on ordered rectangular grids? Nathan Nathan Leon Pace, MD, MStat University of Utah Salt Lake City, UT 84132 Office: 801.581.6393 Fax: 801.581.4367 Cell: 801.558.3987 Pager: 801.291.9019 Home: 801.467.2925 [[alternative text/enriched version deleted]]

This is a question of clarification. IN 2009 Higgins, Thompson and Spiegelhalter (J R Statist Soc A 172:137-159) gave WinBUGs code to get credible intervals from random effects meta analysis for the prediction interval of a new study. It appears that the predict.rma function creates approximate credible intervals (pending a function revision by the author) for that purpose. Is my assumption

I'm trying to set up a mixed model to solve using lme. It will have 3 fixed effects, two random effects and two interaction terms. I've been reading Pinheiro's and Bates's book on the nmle library, but find the part about display functions to be unclear. When creating a groupedData object from a data.frame, you need to enter a function of the form: response ~primary|grouping

we have CENTOS 5.5 on X86. I tried to create a "raw partition" (NOT FS) on a disk and it continue to show "ext3". How can I get ride of it? === procedures===== # parted /dev/sde GNU Parted 1.8.1 Using /dev/sde Welcome to GNU Parted! Type 'help' to view a list of commands. (parted) p Model: DELL PERC

Hi, I've got a linear mixed model created using lmer: A6mlm <- lmer(Score ~ division + (1|school), data=Age6m) (To those of you to whom this model looks familiar, thanks for your patience with this & my other questions.) Anyway, I was trying this to look at the significance of my fixed effects: A6post <- mcmcsamp(A6mlm, 50000) library(coda) HPDinterval(A6post) ..but I got this

I am running R 2.1.1 on a Mac g5 under Mac OS 10.4.2. I have an xyplot with a single conditioning variable (8 levels) . Here is the code for the conditioning variable used in the formula argument of xyplot: factor( drugauthoryear, levels = c( 'bupicapogna1999', 'levobenhamou2003', 'ropicapogna1999',

Hi, I was wondering: I've got a dataset where I've got student 'project's nested within 'school's, and 'division' (elementary, junior, or senior) at the student project level. (Division is at the student level and not nested within schools because some students are registered as juniors & others as seniors within the same school.) So schools are random,

I have a function (see below) that does some bootstrapping (I am happy to expand offline why I could use existing functions.) I put my results into and empty matrix and add a row of results with each iteration. My problem is i am a new user to R and I don't understand data frames, matrices, elements, and vectors well. What I would like is to have a data frame I can manipulate outside of the

In a generalized linear model with k covariates, there are 2(kth power) - 1 possible models (excluding interactions). Awhile ago a posting to R-help suggested Model Selection and Multimodel Inference, 2nd ed, by Burnham and Anderson as a good source for understanding model selection. They recommend (page 71) computing AIC differences over all candidate models in the set of possible models. After

After compiling the kernel 4.9RC1 dmesg gives the following error lpt0: Interrupt-driven port ppi0: <Parallel I/O> on ppbus0 ed1: <Realtek Plug & Play Ethernet Card> at port 0x220-0x23f irq 5 on isa0 ed1: address 52:54:40:28:83:9a, type NE2000 (16 bit) ad0: 4112MB <WDC AC24300L> [8912/15/63] at ata0-master WDMA2 Mounting root from ufs:/dev/ad0s1a module_register: module

Dear Friends, I'm trying to learn to use parallel computation using snow & rpvm Now, I have two boxes: 1) Pentium 4, pvm 3.4.5, R 2.2.0 with FreeBSD 5.4 (box uffbsd); 2) PowerPC, pvm 3.4.5, R 2.1.1 with Debian Linux latest stable (box powerpclinux); I'm using the instructions in (found in snow.pdf) http://www.stat.uiowa. edu/~luke/R/cluster/cluster.html. Now ** from the box

Full_Name: Jelle Goeman Version: 2.0.0 OS: windows Submission from: (NULL) (145.88.209.33) I get a strange error: > is.vector(1:10) Error: recursive default argument reference What's recursive about is.vector? Kind regards, Jelle Goeman My R: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2

Hello R experts, Is the a function to preform a Hosmer Lemshaw goodness of fit? I am running R 2.8. Thanks, AlonBen-Ari Department of Anesthesiology and CCM Univ. Pittsburgh [[alternative HTML version deleted]]

Dear R Masters, I'm an anesthesiology resident trying to make his way through basic statistics. Recently I have been confronted with longitudinal data in a treatment vs. control analysis. My dataframe is in the form of: subj | group | baseline | time | outcome (long) or subj | group | baseline | time1 |...| time6 | (wide) The measured variable is a continuous one. The null hypothesis in