[R] How to bootstrap Kaplan-Miere estimator with 95% envelope

Hi all, Given a typical right-censoring data which contains a time variable and a censoring indicator. How do we bootstrap samples to obtain a 95% envelope for the estimated cumulative hazard function? Thanks in advance. Rui __ [EMAIL PROTECTED]

[R] How to obtain a 95% envelope for the estimated cumulaitve hazard function via bootstrap?

Hi all, Given a typical right-censoring data which contains a time variable and a censoring indicator. How do we bootstrap samples to obtain a 95% envelope for the estimated cumulative hazard function? Rui __ [EMAIL PROTECTED] mailing list

Re: [R] Discrepency between R and MlwiN

On Sun, 12 Sep 2004, Damian Betebenner wrote: When playing around fitting unconditional growth models using R and MlwiN today, I produced two different sets of estimates that I can't reconcile and wondered if anyone here has an idea: The data is two-level repeated measures data with

[R] do.call(dim- , ... )

OK guys another problem. I have a 3D array x with dim(x)=c(a,a,b^2) and I want to rearrange the elements of x to make a matrix y with dimensions c(a*b,a*b). Neither a nor b is known in advance. I want the n-th a*a submatrix of y to be x[,,n] (where 1 = n = b^2). Needless to say, this has gotta

[R] Adding ranks to a repeatedly ragged array

How can I add an extra column containing the rank to a ragged array indexed by more than one grouping factors? E.g. with the barley dataset: How can I to add an additional column ``rank'' containing the rank of the ``yield'' of the different varieties in relation to the indices ``year'' and

Re: [R] permuting dimensions (was do.call(dim- , ... ))

What has this to do with the original subject line? Replacement functions are not intended to be used directly, and certainly not in do.call. See ?aperm, as in xx - x dim(xx) - c(2,2,3,3) xx - aperm(xx, c(1,3,2,4)) dim(xx) - c(6, 6) xx as required. BTW, you have a broken package `magic' on

[R] Re: Variable Importance in pls: R or B? (and in glpls?)

On Sunday 12 September 2004 14:12, Christoph Lehmann wrote: Dear R-users, dear Ron I use pls from the pls.pcr package for classification. Since I need to know which variables are most influential onto the classification performance, what criteria shall I look at: a) B, the array of

Re: [R] Adding ranks to a repeatedly ragged array

What do you mean by `a repeatedly ragged array': you haven't defined it? There is no dataset `barley' in vanilla R. Which one did you mean? There is a data frame in package lattice, but that is neither an array nor ragged. If that is what you want, rank.lists will be a list matrix. The

[R] an integration question

Dear all, I'm stuck on a problem concerning integration..Results from the analytical expression and numerical approximation (as returned by integrate()) do not match. It probably depends on some error of mine, so apologizes for this off-topic question. I'm interested in computing the integral of

[R] R and perl on solaris

Dear all, I am developing a GUI in perl-Tk and I would use R engine to make some analysis and some plots. I have seen that the package RSPerl allows R from perl exchanges but this package is developed on and for intel based computers. When I try to install the package on solaris I have some

Re: [Rd] Re: [R] Sweave echoing comments (again)

On Wed, 08 Sep 2004 17:46:51 +0200, Wolski (W) wrote: Hi! I observed it also. There are cases where it is not desirable. It will be quite helpfull, if possible, to have a parameter that allows one to switch of removing the #comments. The problem is that the parser does not keep

Re: [R] Discrepency between R and MlwiN

Prof Brian Ripley [EMAIL PROTECTED] writes: On Sun, 12 Sep 2004, Damian Betebenner wrote: When playing around fitting unconditional growth models using R and MlwiN today, I produced two different sets of estimates that I can't reconcile and wondered if anyone here has an idea: The

RE: [R] SJava, Client X11

This is *R* help, and you are asking about an *Omegahat* package, and indeed about debugging your local installation. In theory Omegahat has its own mailing lists. But to answer your question: Does anyone success on running the SJava examples with an X11 client console? Yes, I have done so

Re: [R] R and perl on solaris

Patrick, This isn't answering your question, but there are other ways to go about this problem. See http://tolstoy.newcastle.edu.au/R/help/04/05/0952.html for instance Further searching might bring more hits. Sean - Original Message - From: [EMAIL PROTECTED] To: [EMAIL PROTECTED]

RE: [R] Rd files with % (was: permuting dimensions)

I suppose the source for ?%in% in base would be a good place to look: \name{match} \alias{match} \alias{\%in\%} [...] \usage{ match(x, table, nomatch = NA, incomparables = FALSE) x \%in\% table } [...] \details{ \code{\%in\%} is currently defined as \cr \code{\%in\%

[R] Mixture Analysis

Dear all, Does anyone know if the mixture analysis (i.e : sexual dimorphism, age group,…) can be done with R software ? Thanks for your help Dan - és pour dialoguer instantanément avec vos amis.Téléchargez GRATUITEMENT ici !

Re: [R] Spare some CPU cycles for testing lme?

On Mon, 2004-09-13 at 07:40, Frank Samuelson wrote: If anyone has a few extra CPU cycles to spare, I'd appreciate it if you could verify a problem that I have encountered. Run the code below and tell me if it crashes your R before completion. library(lme4) data(bdf) dump-sapply(

[R] lmList for glm

Greetings all, Is there a way to do the equivalent of a lmList (fit separate models for each subject) in the package nlme for a glm? Thank you very much, Peter B. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

[R] Smoothing using the kernel distribution

Hi, Does anyone can tell me if R sofware can be used for the smoothing using the kernel distribution ? I will appreciate Dan - és pour dialoguer instantanément avec vos amis.Téléchargez GRATUITEMENT ici ! [[alternative HTML version

Re: [R] Spare some CPU cycles for testing lme?

Marc Schwartz wrote: On Mon, 2004-09-13 at 07:40, Frank Samuelson wrote: If anyone has a few extra CPU cycles to spare, I'd appreciate it if you could verify a problem that I have encountered. Run the code below and tell me if it crashes your R before completion. library(lme4) data(bdf)

[R] bagplot()

Hello, I saw a little discussion about this in the archives, but it was unclear to me whether someone had submitted a port to R of the Splus bagplot() function. If so, does anyone know where I could get it? Thanks. Best, Matt __ [EMAIL PROTECTED]

Re: [R] Smoothing using the kernel distribution

Since you have now asked this three times: Yes, yes, yes. R comes with a pacakge called KernSmooth, for example. On Mon, 13 Sep 2004, bournery alexandre wrote: Does anyone can tell me if R sofware can be used for the smoothing using the kernel distribution ? -- Brian D. Ripley,

Re: [R] bagplot()

There is no bagplot function in S-PLUS. There is one by Rousseeuw et al for S at http://www.agoras.ua.ac.be/Locdept.htm and a reply in the archives about the problems of porting at: http://maths.newcastle.edu.au/~rking/R/help/03b/4916.html My guess is that like several other cases, R porting

RE: [R] Variable Importance in pls: R or B? (and in glpls?)

Christoph: I noted that there were not a great number of people leaping to reply. One reason, I suspect, is that there's really NO GOOD ANSWER to this question. First, there is a huge literature on this -- it's related to variable selection in regression and shrinkage estimates, but, in general,

Re: [R] Kernel distribution

See ?density or try help.search(kernel) On Mon, 2004-09-13 at 16:05, Alexandre Bournery wrote: Hello, Does anyone could tell me if R sofware can be used for smoothing, using the kernel distribution ? I will appreciate Alex [[alternative HTML version deleted]]

[R] How to show the symbol of Angstrom ?

Also, I am wondering if there is any source where the expressions of many symbols are collected. Thanks you very much!!! Li, Xiang(Sean) __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

RE: [R] bagplot()

There's now a vioplot package on CRAN, if that's what you're looking for. Andy From: Matthew David Sylvester Hello, I saw a little discussion about this in the archives, but it was unclear to me whether someone had submitted a port to R of the Splus bagplot() function. If so, does

[R] R-help

Dear all, Does anyone know how to read files with .dbf extension? Thanks for your time. __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] R-help

Did you do a search at www.r-project.org - R site search? Searching for .dbf there just now exposed a read.shape function in the maptools package. hope this helps. p.s. Did you read the posting guide! http://www.R-project.org/posting-guide.html;? Tips provided there may help

Re: [R] lmList for glm

Have you looked at GLMM in package lme4? hope this helps. spencer graves p.s. Have you read the posting guide! http://www.R-project.org/posting-guide.html;? It might help you get answers to questions like this for yourself quicker than waiting for this list to reply. It may also

[R] calculating memory usage

I am comparing two different algorithms in terms of speed and memory usage. I can calculate the processing time with proc.time() as follows but am not sure how to calculate the memory usage. ptm - proc.time() x - rnorm(100) proc.time() - ptm I would like to be within R itself since

Re: [R] R-help

I believe .dbf files are more commonly DBase files, in which case see package RODBC. On Mon, 13 Sep 2004, Spencer Graves wrote: Did you do a search at www.r-project.org - R site search? Searching for .dbf there just now exposed a read.shape function in the maptools package.

Re: [R] calculating memory usage

On Mon, 13 Sep 2004, Adaikalavan Ramasamy wrote: I am comparing two different algorithms in terms of speed and memory usage. I can calculate the processing time with proc.time() as follows but am not sure how to calculate the memory usage. ptm - proc.time() x - rnorm(100)

[R] generalized eigenvalues - porting code using eig() from Matlab to R

Dear list, I am porting some Matlab routines to R. I am looking for the equivalent function in R for the Matlab function eig(), such that eig(A, B) quoteproduces a diagonal matrix D of generalized eigenvalues and a full matrix V whose columns are the corresponding eigenvectors so that A*V =

Re: [R] generalized eigenvalues - porting code using eig() from Matlab to R

On Mon, 13 Sep 2004, Gavin Simpson wrote: I am porting some Matlab routines to R. I am looking for the equivalent function in R for the Matlab function eig(), such that eig(A, B) quoteproduces a diagonal matrix D of generalized eigenvalues and a full matrix V whose columns are the

Re: [R] .Random.seed in R-devel

On Mon, 13 Sep 2004, Frank Samuelson wrote: I'm running R-2.0.0 (yesterday's snapshot)in its own directory, and everything works great, except: .Random.seed Error: Object .Random.seed not found Does 2.0.0 not use .Random.seed for saving, etc, like it says in the help page? `Like

Re: [R] .Random.seed in R-devel

Frank Samuelson [EMAIL PROTECTED] writes: I'm running R-2.0.0 (yesterday's snapshot)in its own directory, and everything works great, except: .Random.seed Error: Object .Random.seed not found Does 2.0.0 not use .Random.seed for saving, etc, like it says in the help page? Thanks

[R] R in Linux Redhat 9.0

Hi, my name is Bryan and I just installed R on my Redhat 9.0. I imported the key from the net as directed, and tested the integreity of the rpm package- it was fine. I then went ahead and installed the i386.rpm package and the installation was carried to completion. However, now I have a hard

Re: [R] .Random.seed in R-devel

Oh. I guess I had a different definition of when one is required than the help page. :) Thanks. Prof Brian Ripley wrote: `Like it says on the help page'! Initially, there is no seed; a new one is created from the current time when one is required. Hence, different sessions will

Re: [R] R in Linux Redhat 9.0

rpm -ql R should show you your installed files. Did you try typing R? Bryan Tseng wrote: Hi, my name is Bryan and I just installed R on my Redhat 9.0. I imported the key from the net as directed, and tested the integreity of the rpm package- it was fine. I then went ahead and installed the

Re: [R] R in Linux Redhat 9.0

Hi! Just open a terminal and type 'R'. Read the R-Intro will help you a lot. Regards. On Mon, 13 Sep 2004, Bryan Tseng wrote: Hi, my name is Bryan and I just installed R on my Redhat 9.0. I imported the key from the net as directed, and tested the integreity of the rpm package- it was

[R] Measures of association for ordinal data

In a classification problem with ordinal data (classes 1 - 4), I used multidimensional optimization to maximize gamma (Goodman's measure of association) between observations and predictions. This resulted in the following frequency table (rows = observations, columns = predictions):

Re: [R] Spare some CPU cycles for testing lme?

On Mon, Sep 13, 2004 at 08:40:15AM -0400, Frank Samuelson wrote: If anyone has a few extra CPU cycles to spare, I'd appreciate it if you could verify a problem that I have encountered. Run the code below and tell me if it crashes your R before completion. library(lme4) data(bdf)

[R] R CMD SHLIB setup problem...

All, When I try and compile a shared library (on WinXP) I get the following error: E:\data\projR CMD SHLIB toy_dll.c Makevars:1: *** missing separator. Stop. Has someone else had this error and fixed it? This code compiles and works fine on Linux (fedora core 2). Everything was working

Re: [R] How to show the symbol of Angstrom ?

Hi xiang li wrote: Also, I am wondering if there is any source where the expressions of many symbols are collected. Thanks you very much!!! (Assuming you mean draw the angstrom symbol on a plot ...) There are several ways: (i) Specify the character code in octal. Assuming ISO Latin 1 encoding,

Re: [R] R CMD SHLIB setup problem...

On Tue, 14 Sep 2004 10:29:31 +1000, [EMAIL PROTECTED] wrote: All, When I try and compile a shared library (on WinXP) I get the following error: E:\data\projR CMD SHLIB toy_dll.c Makevars:1: *** missing separator. Stop. Has someone else had this error and fixed it? Haven't had the error,

[R] as.integer(TRUE)

The fact that as.integer(TRUE) is 1 (and that for FALSE, it gives zero) is a really nice feature, eg when constructing piecewise functions (for example, as in -x*(x0)+x*(x=0)) and for many other things. Since I haven't found a reference about this, I just wanted to ask whether this is officialy

RE: [R] R CMD SHLIB setup problem...

I've compared my Makevars to someone else's and we both have the same thing: CFLAGS+= -I$(RHOME)/src/gnuwin32/graphapp So I assume this is correct? (And the compiler is definitely GNU make). When I tried to recompile another bit of C code I got the following error: make: *** No rule to make

[R] drawing on axes

Hi, I would like to make certain portions of axis lines thicker (as Tufte suggests). How can I draw on axes? I only need a couple of line segments on the left and bottom one. Thanks Tamas __ [EMAIL PROTECTED] mailing list

[R] Signs of loadings from princomp on Windows

I start a clean session of R 1.9.1 on Windows and I run the following code: library(MASS) data(painters) pca.painters - princomp(painters[ ,1:4]) loadings(pca.painters) Loadings: Comp.1 Comp.2 Comp.3 Comp.4 Composition 0.484 -0.376 0.784 -0.101 Drawing 0.424 0.187 -0.280

RE: [R] drawing on axes

I would like to make certain portions of axis lines thicker (as Tufte suggests). How can I draw on axes? I only need a couple of line segments on the left and bottom one. How about: plot (1, 1, xlim= c (-10, 10), bty= 'n') axis (1, at= c (-10, -5), labels= FALSE, tick= T, lwd= 5, tck= 0)