similar to: boot package

In the non linear mixed effects package a groupedData object can be created to facilitate modeling. The gD object includes a formula of the form 'response variable' ~ 'primary covariate' | 'grouping factor'. In experiments creating response surfaces there are 2 or more primary covariates. Is there any way to use the groupedData() function to include 2 primary

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]]

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

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

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

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

> 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'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

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',

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

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

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,

Hi, LLVMers, has someone implemented bitwidth analysis for LLVM? I was looking for something similar to the bitwise compiler described in "Bidwidth analysis with application to silicon compilation, by Mark Stephenson, Jonathan Babb and Saman Amarasinghe" e.g.: http://portal.acm.org/citation.cfm?id=349299.349317 all the best, Fernando

I'm doing some analyses of historical data from France in 1830 on 'moral statistics' that I'd like to show on a map. I've done most of my analyses in SAS, but a few things would work better in R. To do this, I have to adjust the modern map, library(maps) map('france') to adjust for changes in departments (86 in 1830, to 97 now). I've read the documentation

Hello, This is a pretty simple question, but after spending quite a bit of time looking at "Hmisc" and using Google, I can't find the answer. If I use stat_summary(fun.data="mean_cl_boot") in ggplot to generate 95% confidence intervals, how many bootstrap iterations are preformed by default? Can this be changed? I would at least like to be able to report the number of

Overview ======== This module implements the Center for Internet Security (CIS) Security Configuration Benchmark for Red Hat Enterprise Linux 6 v.1.1.0 (avilable at http://benchmarks.cisecurity.org). Each scored control has been implemented as a class or a custom fact. Installation ============ Please either: - Clone git repo from https://github.com/arildjensen/cis-puppet - Run "puppet

I have a matrix containing means and CIs (lower and upper in two columns, so three columns for every data point) for several points. I have to build a graph of these means accompained by the CIs (as wiskers). No problems with making the graph of means, but I don't know how to introduce CIs. Can anybody advise? -- View this message in context:

One variable contains values (1.30 - one hour and thirty minutes, 1.2 (which is supposed to be 1.20 - one hour and twenty minutes)). I would like to convert to a minute variable so 1.2 is converted to 80 minutes. How?

I have a system running rsync in daemon mode, the version information: > rsync --version rsync version 2.6.8 protocol version 29 On another system I have: > authoritative# rsync --version rsync version 3.0.7 protocol version 30 I have the following problem when trying to fetch files from the system running the old daemon: # /usr/local/bin/rsync --protocol=29 -vvvvvvvv

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