similar to: Retaining regularly used add-ons

Hi all, I'm using R 1.5.1 on Windows 2000. The following snippet of code doesn't seem to do anything - no error is reported, and there is no name change. names(myFrame[,c(1:3)]) <- c("name1", "name2", "name3") This code however works nicely: names(myFrame)[c(1:3)] <- c("name1", "name2", "name3") Can anyone suggest why

Hello All, Specifying 'type = "response"' when using predict() on a model fit using glm(...,family="binomial") returns fitted probabilities. Is it possible to get the same from a model object fit using glmmPQL() ? Thanks in advance, Rob _____________________________________________________ Rob Keefe Lab: (208) 885-5165 M.S. student

Hi All, Is it possible for system.time() to measure the time it takes for a machine to evaluate more than one R expression? For example, # This I can do: > system.time(x <- rnorm(100000)) [1] 0.07 0.00 0.13 0.00 0.00 # But this I can't: > system.time(x <- rnorm(100000); new <- sample(x, 100000, replace=T)) Error: syntax error # Nor this: > system.time(x <-

I'm running a simulation that uses nls on each run. Sometimes, not unexpectedly, nls doesn't return parameter estimates (e.g., "singular gradient" error). Of course when this happens the simulation (in a loop) stops and no results are recorded. What I would like is for the simulation to simply record that nls couldn't obtain estimates on a particular run due to whatever

HI I'm a student in chemical engineering, and i have to implement an algoritm about FIVE PARAMETERS INTERPOLATION for a calibration curve (dose, optical density) y = a + (c - a) /(1+ e[-b(x-m]) where x = ln(analyte dose + 1) y = the optical absorbance data a = the curves top asymptote b = the slope of the curve c = the curves bottom asymptote m = the curve X intercept Have you never seen

Dear R community, In the excellent nlme package the default diagnostic plot graphs the innermost residuals against innermost fitted values. I recently fit a mixed-effects model in which there was a very clear positive linear trend in this plot. I inferred that this trend occurred because my fixed effect was a two-level factor, and my random effect was a 12-level factor. The negative residuals

I don't get the same result, do you have a package loaded that would change the default behavior (such as Hmisc)? > list(grp.1) [[1]] [1] 1 2 Levels: 1 2 > list(grp.1[mask]) [[1]] [1] 1 Levels: 1 2 > library(Hmisc) <<snip>> > list(grp.1[mask]) [[1]] [1] 1 Levels: 1 --Matt > version _ platform i386-pc-mingw32 arch i386 os

Hello, I have been searching for ways to obtain these for combinations of fixed factors and levels other than the 'baseline' group (contrasts coded all 0's) from a mixed-effects model in lme. I've modelled the continuous variable y as a function of a continuous covariate x, and fixed factors A, B, and C. The fixed factors have two levels each and I'd like to know whether

I would appreciate some thoughts on using the bootstrap functions in the library "bootstrap" to estimate confidence intervals of ICC values calculated in lme. In lme, the ICC is calculated as tau/(tau+sigma-squared). So, for instance the ICC in the following example is 0.116: > tmod<-lme(CINISMO~1,random=~1|IDGRUP,data=TDAT) > VarCorr(tmod) IDGRUP = pdLogChol(1)

I have two lists: xa <- list( X=c(1,2,3), Y=c(4,5,6), Z=c(7,8,9) ) xb <- with( barley, tapply( X=seq(1:nrow(barley)), INDEX=site , FUN=function(z)yield[z])) I can convert xa to a dataframe easily with: as.data.frame(xa) But if i try the same with xb I get: as.data.frame(xb) Error in as.data.frame.default(xb) : can't coerce array into a data.frame What

Dear list, I'm trying to download a text file directly from the internet using the RCurl package and the command getURL. Duncan Lang graciously helped me solve the first step in this problem using the following command: ################# txtfile <- getURL('ftp://ftp.wcc.nrcs.usda.gov/data/snow/snow_course/table/history/idaho/13e19.txt', ftp.use.epsv = FALSE) #################

Hi, I just installed the odesolve package and ran the check command on it. It failed trying to execute library(gnlm). Sure enough, there is no gnlm library on my system and I could not find it on the CRAN archive either. Am I missing anything obvious or is gnlm some private library that somehow found its way into the example section of odesolve? Thanks in advance, Andy

Greetings all, I'm trying to figure out how to calculate the inverse CDF (i.e. a quantile) for a non-central F distribution. I could put together a quick numerical solver routine using the CDF, but I wonder if there's a function that I've missed that would be more efficient? Thank-you, Andrew Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885

Hi everyone, I want to create my own distance measure, other than 'euclidean' or 'manhatan', to use in cluster pckgs. To do this I think that I need to change dist(), in mva pckg, or daisy(), in cluster pckg. (or is there a cleaver way ?) But this functions are in fact things like: .Fortran( "daisy", ... ) or .C("dist",...). I tried unsuccessfully to find

Pursing my earlier question, when I tried loading Lindsey's gnlm, I got a message Loading required package: rmutil Warning message: There is no package called 'rmutil' in: library(package, character.only = TRUE, logical = TRUE, warn.conflicts = warn.conflicts, According to the R documentation http://finzi.psych.upenn.edu/R/doc/html/packages.html rmutil is in the standard

Why does stop("we are done") print "Error in eval.with.vis(expr, envir, enclos) :" ? It would seem to me that a plain stop() is not an error, and that it would make more sense to have an error() function that is different from a stop(). Is there a rationale here that I am missing? sincerely, /iaw

Top Ten Female Turn Ons When it comes to sexual turn ons, men and women are very different. It is not a myth that women are attracted to a man who is comfortable with his emotions. http://2site.com/epqjag The Truth About Creative Baby Shower Thank You Poems Where are the cutest baby shower thank you poems? Don''t you know the best thank you poems lie somewhere deep in your

Dear list, I am trying to download a text file from an ftp site using download.file(). I used the following code: url <- "ftp://ftp.wcc.nrcs.usda.gov/data/snow/snow_course/table/history/idaho/13e19.txt" dest <- "c:/test/downloadtest.txt" download.file(url, dest) I get this error message, indicating that R is unable to open the connection: trying URL

There are some rather nice archive replacements for Mailman that provide full searching, etc. I can provide more info if wanted. Geoff

Dear R Help-Team, If I start R and type in search() the result is: [1] ".GlobalEnv" "package:ctest" "Autoloads" "package:base" I often need an Add-on-package. I can load it with "> library(package)" or include a require-command in the .First-function. The .First-function is fine, if I always use the same working directory. But I