similar to: R INSTALL -l

Using R (version 0.63.3) for MS windows, I try the following command > shell(paste("cd ",getenv("RHOME"),sep="")) which replies with an error message Too many parameters - FILES\RW0633 which appears to suggest that the space if the path name is causing difficulties to the cd command. getenv returns > getenv("RHOME") RHOME

Hello everyone, (This is my second question posted today on the R list). I am carrying out a competing risks analysis using the cuminc function...this takes the form: cuminc(ftime,fstatus,group) In my study, fstatus has 3 different causes of failure (1,2,3) there are also censored cases (0). "group" has two levels (0 and 1). I therefore have 6 different cumulative incidence curves:

Hi, for example one could do it this way: v <- summary(fit)$coefficients[,4] the coefficient attribute is a matrix, and with the "4" you refere to the pvalue (at least in lm - don't know if summary(glm) produces sligthely different output). to skip the intercept (1st row): v <- summary(glmfit)$coefficients[-1,4] hope this helps, Arne -- Arne Muller, Ph.D.

Updating from centos 5.4 to 5.5 using an update rather than a rebuild, I get the following complaint from yum:- lvm2-2.02.56-8.el5_5.1.x86_64 from setup has depsolving problems --> Missing Dependency: device-mapper >= 1.02.39-1 is needed by package lvm2-2.02.56-8.el5_5.1.x86_64 (setup) This is from a local mirror, so checking the upstream feed in case the mirroring is incomplete I see

Dear List Apologies for such a simple question: I have a vector of 738 elements, coded with values between 1 and 16 (but not containing 7, 10, 11 or 13) and wish to recode value 14 to 1, 4 to 2, 1 to 3, 2 to 4 and all other values to 5. I've been trying to use the replace function (in base) and %in%, but without success. Many thanks Roy

Hello I'm very new to R (which is excellent), so apologies if this has already been raised. In the spdep package, I'm trying to undertake an autoregressive mixed model using the lagsarlm function. This is working fine, but there does not appear to be a method of including an explanatory variable without it automatically being included as a lagged term. I'm after something along the

I apologise if this type of question has been answered before, however how do I count the number of TRUE's in the upper half (non-symmetric) matrix of logicals. Thank you Ross Darnell University of Newcastle -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I am trying to adjust the sigma symbol down in the y direction with the following command but R recycles the sigma and puts another at -0.1 in the x axis. Any suggestions would be appreciated. mtext(expression(sigma),1,outer=TRUE,adj=c(0.5,-0.1)) Regards Ross Darnell -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I am trying to place a legend onto an image but the legend to clipped/nonexitent outside the boundary of the image. I am wanting to place the legend where the main title normally would appear. I would appreciate someone telling me why I cannot do this. (Which part of the help file I haven't read most likely) Regards Ross Darnell

I think I have a data.frame called "statlab". > is.data.frame(statlab) [1] TRUE and I can perform summary, dim, ... but when I try > attach(statlab) Error in attach(statlab) : attach only works for lists and data frames I am using > version _ platform i686-unknown-linux arch i686 os linux system i686,

I would appreciate help translating the following lme model to an lmer function. lme(lognrms ~ Group*Rotation*muscle*side*support*arms, random=~1|Subject/Stratum2/rep, data=Data) Many thanks Ross Darnell r.darnell at uq.edu.au

I have generated a graph on the X11 device, but when I > dev.copy(postscript) postscript 3 > dev.off() X11 2 > the Rplots.ps file is generated but contains a prologue only. Can anyone suggest what I am doing wrong? Just to anticipate the upgrade reply ... > version _ platform i686-pc-linux-gnu arch i686 os linux-gnu

I am relatively new to R, so this may be a much simpler question than it seems to me. I am trying to create a function that includes two summations, and can't figure out how. I am attaching the equation as a pdf file. This function will then be optimized over a chosen range of values, but right now I just need help with the function. Thank you. Zack Darnell Duke University -- Zack

Dear R-help I have two separate experiments, one a repeated-measures design, the other a split-plot. In a standard ANOVA I have usually undertaken a multiple-comparison test on a significant factor with e.g TukeyHSD, but as I understand it such a test is inappropriate for repeated measures or split-plot designs. Is it therefore sensible to use Helmert contrasts for either of these designs?

I have fitted a logistic regression model > failed.lr2$call lrm(formula = failed ~ Age + task2 + Age:task2, data = time.long, na.action = na.omit) using the Design package functions and would like to generate a nomogram from this model. the datadist information is generated and stored in > ddist time.long$Age time.long$task2 Low:effect 45

Hello everyone, These are some questions about the 'pec' function in R. These questions deal with prediction error curves and their derivation. Prediction error curves are documented in, for example, "Efron-type measures of prediction error for survival analysis" by Gerds and Schumacher. I have detailed some syntax that I have used at the bottom of this email. The associated

I would like to create a page of two graphs (2 rows by 1 col) and then draw vertical lines (segments?) on both graphs from the minimum values to the corresponding maximum value. So I have tried # > y <- rnorm(3000) > par(mfrow=c(2,1)) > plot(y,type="l") > plot(cumsum(y),type="l") > segments(1000,min(cumsum(y)),1000,max(cumsum(y))) > par(mfg=c(1,1)) >

I would like to define the order of the levels of a factor. The relevel function would work but since I have 20 levels I would prefer to declare the order explicitly. Using a smaller example levels(oldfactor) "b1" "b2" "r1" "r2" nufactor <- order(oldfactor,order=c("b1","r1","b2","r2")) # my fabricated function

Maybe a useful addition to the predict functions would be to return the values of the predictor variables. It just (unless there are problems) requires an extra line. I have inserted an example below. "predict.glm" <- function (object, newdata = NULL, type = c("link", "response", "terms"), se.fit = FALSE,

Dear R developers, Running a R CMD check under R2.2 gives me the error messages (appended to this e-mail) while the with R2.1 the check on the same package directory runs fine. To my knowledge the configuration of R2.2 and R2.1 on my machine are identical. * checking S3 generic/method consistency ... WARNING Error: .First.lib failed for 'SBMLodeSolveR' Call sequence: 2: