similar to: xyplot problem

To enhance my understanding, and that of my students, I have a question about cox.zph in the survival package. If I have correctly gleaned the high-level point from the 1994 Biometrika paper of Grambsch and Therneau, it looks to me like cox.zph provides a mechanism to test for a simple trend in plots of a function of time, g(t) versus the scaled schoenfeld residuals and it also provides some

I downloaded R-Patched today (to see if another problem I want to ask about is still present or if its just me - as per the posting guide). I ran "tools/rsync-recommended" successfully. I then ran "configure --enable-R-shlib" successfully. Then make stops with the following error. gcc -shared -L/usr/local/lib -o tools.so text.o init.o Rmd5.o md5.o -L../../../../lib -lR

Hi friends, I'm beginning in R and I have simple question. I have this piece of my program and how you see, that's ok (whit > num<- 0.002) num<-0.002 # ok, but not when I change whit num<-0... ? factor1<-1; while(1) { if (num*factor1<1) factor1<-factor1*10 else { print("out ok!!"); break; } } [1] "out ok!!" but when I change (whit >

I seem to remember reading somewhere about some style guide regarding R the number of comment characters (#) prior to the comment meaning something. Anyone know to what I'm referring? Where? ---------------------------------------------------------- SIGSIG -- signature too long (core dumped)

> From: Duncan Murdoch <murdoch at stats.uwo.ca> > To: ramasamy at cancer.org.uk > Date: Fri, 10 Mar 2006 08:33:09 -0500 > Subject: Re: [R] To improve my understanding of workspaces > Other than Emacs, I use the same work habits as Adai. An advantage of > this workflow is that almost everything is stored in text format, so it > is easy to compare different versions to

Dear all, I apologize if this is a FAQ (seems a bit like one, but I didn't find anything). I'm looking for an easy way to cut one value out of a vector and shorten the vector accordingly. Something like: x <- c(1, 1, 0, 6, 2) throwaway(x[3]) which will return x = 1 1 6 2, with length(x) = 4. I know one could do this by hand, but then one would have to create a second vector y in a

Dear All, I am wondering on how to use the abline.xyplot with xyplot such that I will have different vertical lines for each panel. More sepcifically, suppose that the xyplot generates 4 panels defined by the combination of two binary variables: X_1 and X_2. i.e. xyplot(Y ~ Z | X_1*X_2, data = df) I want something like: abline(v = 5) if X_1=0 and X_2 = 0 abline(v =

I just upgraded to version 2.1.1 (from 2.0.1) today. > R.version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 2 minor 1.1 year 2005 month 06 day 20 language R I am using SuSE 9.2 and did the upgrade using rpm -U with the RPM available on CRAN. After upgrading r-base, I ran update.packages(). Four previously

Hi, I am using Sweave for chapters in my thesis that contain results. In the beginning of each chapter, I use this to load libraries I need. <<init,echo=FALSE,quiet=TRUE>>= library(gplots) library(Hmisc) library(e1071) @ What I want is, of course, to supress messages written by this code, but what I get in the end is X-init.tex with this the contents below. How do I really supress

Dear All, I wonder if there is any R function to generate a sequence of vectors from existing ones. For example: x 1<- c(1,2,3) x2 <- c(4,5) x3 <- c(6,7,8) The desired output is a list of all 3*2*3 = 18 possible combinations of elements of x1,x2 and x3. One element for example is (1,4,6). Many thanks in advance, Bernard

Dear All I am totally new to R and I would like to know whether R is able and appropriate to illustrate to my students the Central Limit Theorem, using for instance 100 independent variables with uniform distribution and showing that their sum is a variable with an approximated normal distribution. Thanks in advance, Paul

Greetings all: OK, this is bugging the @#@%* out of me. I know the answer is simple and straightforward but for the life of me I cannot find it in the documentation, in the archives, or in my notes (because I know I've encountered this in the past). My problem is: I have a data frame with columns A, B, C, D, and E. A, B, and E are factors and C and D are numeric. I need a new data frame with

R-experts: A quick question, please. >From a lab exp, I got 12 positives out of 50. To get 90% CI for this , I think binom.test might be the one to be used. Is there a better way or function to calculate this? > binom.test(x=12, n=50, p=12/50, conf.level = 0.90) Exact binomial test data: 12 and 50 number of successes = 12, number of trials = 50, p-value = 1 alternative

Is there a way to do text plots in R? I'd like to do some simple XY plots in R with the output in text (ascii). For example, with gnuplot I can do the following: echo 'set terminal dumb ; plot sin(x)' | gnuplot To generate a simple sin wave. Since I connect to a remote Linux machine using SSH, being able to generate a rough idea of what a plot will look like in text would be

Greetings, I have generated several ROC curves and would like to compare the AUCs. The data are cross sectional and the outcomes are binary. I am testing which of several models provide the best discrimination. Would it be most appropriate to report AUC with 95% CI's? I have been looking in to the "net reclassification improvement" (see below for reference) but thus far I can only

I'm sure this is simple enough, but an R site search on my subject terms did suggest a solution. I have a numeric vector with many values that I wish to create a factor from having only a few levels. Here is a toy example. > x <- 1:10 > x <-

Hello. I am posting this on behalf of a student in my class. I have been introducing them to R. I suspect that the problem I need help with is not so much an R issue, but a MAC config issue, but I don't have enough MAC experience to know what to do. Anyway, on to the problem. We are using Frank Harrell's Hmisc and Design packages. This student also wishes to use the latex formatting

Dear All, I wonderif there is a simple way to draw a regression line in the xyplot: more specifically, let: age <- c(20:30, 31:40 ) age.cut <- cut(age, breaks = 2 ) y<- rnorm(20) x <- rnorm(20,4,1) xyplot(y ~ x| age.cut, xlab="x", ylab="y") How to draw (in the plot given by xyplot) the two regression lines (y ~ x) corresponding to the two category

I have a dataset at a hospital level (as opposed to the patient level) that contains number of patients experiencing events (call this number y), and the number of patients eligible for such events (call this number n). I am trying to model logit(y/n) = XBeta. In SAS this can be done in PROC LOGISTIC or GENMOD with a model statement such as: model y/n = <predictors>;. Can this be done

I tried to use write.foreign() to export to SAS this morning and got an error. When I looked at the code for writeForeignSAS() I saw this line: dfn < -df which I think should be dfn <- df So, I tried to run update.packages() to see if there was an updated version and got the following result. > update.packages(c("foreign")) Warning message: In list.files(lib) :