similar to: 'summary' with logicals (PR#2629)

In the first line, use the dist function, found in library mva, to get the distance between each pair of rows. From this calculate an incidence matrix for which element i,j is true if row i in dat equals row j in dat (and false elsewhere). In the second line, for each row calculate the indices of the matching rows and take the minimum of those as the key. incid <-

NAME Douglas Bates EMAIL bates@stat.wisc.edu VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Debian 2.1 CC/FC/MAKE egcs/g77/make NAME Martyn Plummer EMAIL plummer@iarc.fr VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Redhat 5.1 CC/FC/MAKE gcc/egcs-g77/make NAME Göran Broström EMAIL gb@stat.umu.se VERSION 0.63.3 PLATFORM

I wanted to add some text to a plot and got (R-1.6.0, Linux): > text(x = c(1, 4), y = 5, labels = x) Error in text.default(x = c(1766, 1895), y = 5, labels = x) : Object "x" not found With the default value of 'labels': > text(x = c(1, 2), y = 5, labels = seq(along = x)) Error in seq(along = x) : Object "x" not found A scoping bug? :) But >

On Thu, 18 Apr 2002 ripley@stats.ox.ac.uk wrote: > On Thu, 18 Apr 2002, Martin Maechler wrote: > > > >>>>> "Jason" == Jason Turner <jasont@indigoindustrial.co.nz> writes: > > > > Jason> strptime() mysteriously adds a day to a date, unless the year > > Jason> is specified. Tested on: > > Jason> Linux (RedHat

SunOS fluke 5.8 Generic_108528-07 sun4u sparc SUNW,Ultra-Enterprise R-1.4.0 Running from the command line R CMD check chlib results in the log file 00check.log: * using log directory `/home/woodstock/hoffmann/R/Sources/chlib.Rcheck' * checking for file `chlib/DESCRIPTION' ... OK * checking package directory ... OK * checking for sufficient/correct file permissions ... WARNING *

Hi, the link to the R code for repeated measurement analyses of J Lindsey is unfortunately not working. I am desperate for repeated measurements in R; could you please help me out. Sincerely, Dr. G. Stoet -- Dr. Gijsbert Stoet email: stoet@thalamus.wustl.edu Web: http://eye-hand.wustl.edu/lab/people/stoet.html Phone: (314)7474095 Fax: (314)7474370

On Thu, 18 Apr 2002, Martin Maechler wrote: > >>>>> "Jason" == Jason Turner <jasont@indigoindustrial.co.nz> writes: > > Jason> strptime() mysteriously adds a day to a date, unless the year > Jason> is specified. Tested on: > Jason> Linux (RedHat 6.0) - R version 1.4.1 and R-devel. > Jason> Windows - R version 1.4.1 >

On Thu, 19 Sep 2002 10:43:01 -0400 (EDT), you wrote in message <200209191443.KAA17404@falmouth.bwh.harvard.edu>: >Hi Duncan -- I am curious as to whether you are building/distributing >R 1.6.0 for windows? There's nothing urgent about it, I have one. >The question is how to identify a distribution URL if one exists. I will be building and distributing the final release.

I want the indices i for which x[i] < 0 (say): > x <- c(1, -1, 3, 3, -2) > where.negative(x) [1] 2 5 Surely where.negative is something simple, but how? G?ran -------------------------------------------------------------- G?ran Brostr?m Department of Statistics tel: +46 90 786-5223 Ume? University fax: +46 90 786-6614 S-90187 Ume?, Sweden

I find 'How do I convert factors to numeric?' in the FAQ (7.12), but not the other way around. Trivial maybe, but > codes(factor(c(2, 10)) [1] 2 1 > codes(factor(c(2, 9))) [1] 1 2 How do I get the levels attached to the codes in numeric order? G?ran --- G?ran Brostr?m tel: +46 90 786 5223 Department of Statistics fax: +46 90 786 6614 Ume? University

I need to order a long vector of integers with rather few unique values. This is very slow: > x <- sample(rep(c(1:10), 50000)) > system.time(ord <- order(x)) [1] 189.18 0.09 190.48 0.00 0.00 But with no ties > y <- sample(500000) > system.time(ord1 <- order(y)) [1] 1.18 0.00 1.18 0.00 0.00 it is very fast! This gave me the following idea: Since I don't care

I've rolled up R-1.6.0.tgz a short while ago. This is a major upgrade, with several new features. A few of the changes generate compatibility issues. You may wish to consult http://developer.r-project.org/160update.txt. We no longer ship the set of recommended packages separately, but bundle them up with the main distribution. Thus there is now only a single file to fetch and (barring

I've rolled up R-1.6.0.tgz a short while ago. This is a major upgrade, with several new features. A few of the changes generate compatibility issues. You may wish to consult http://developer.r-project.org/160update.txt. We no longer ship the set of recommended packages separately, but bundle them up with the main distribution. Thus there is now only a single file to fetch and (barring

I have a problem with some functions handling dates, in packages 'date' and 'survival' (they seem to be identical). For instance, from the documentation, -------------------- mdy.date {survival} R Documentation Convert to Julian Dates Description Given a month, day, and year, returns the number of days since January 1, 1960. Usage mdy.date(month, day, year, nineteen = TRUE,

R 1.6.0 has gone into final feature freeze as of today. As a new feature, we'll make interim beta versions available via ftp://cran.us.r-project.org/pub/R/src/base alias http://cran.us.r-project.org/src/base (filename R-1.6.0beta_*.tar.gz, where * is the creation date). If you want to help ensure that the final 1.6.0 works satisfactorily on *your* platform, it might be a good idea to

On Wed, Aug 1, 2012 at 3:08 PM, <fra.meucci@hotmail.it> wrote: > Dear Göran Broström, > I am trying to use AFTREG function for R to estimate a loglogistic > survival function, including time dependent covariates. > Actually, my Subset includes some partial events; the idea is to model > this kind of events using something similar to “weights” in the SURVREG > function.

Only three entries thus far. Please keep them coming in. -k NAME Martyn Plummer EMAIL plummer@iarc.fr VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Redhat 5.1 CC/FC/MAKE gcc/egcs-g77/make NAME Douglas Bates EMAIL bates@stat.wisc.edu VERSION 0.63.3 PLATFORM i386-unknown-linux SYSTEM Debian 2.1 CC/FC/MAKE egcs/g77/make NAME Thomas

Dear List, After including cluster() option the coxreg (from eha package) produces results slightly different than that of coxph (from survival) in the following time-dependent treatment effect calculation (example is used just to make the point). Will appreciate any explaination / comment. cheers, Ehsan ############################ require(survival) require(eha) data(heart) # create weights

I have a data frame containing children, with variables 'year' = birth year, and 'm.id' = mother's id number. Let's assume that all the births of each mother is represented in the data frame. Now I want to create a subset of this data frame containing all children, whose mother's first birth was in the year 1816 or later. This seems to work: mid <-

When reading an spss file, I get the following message: Warning message: ../totmorH.sav: Unrecognized record type 7, subtype 13 encountered in system file. What does it mean? Do I have to worry? The result looks as expected. I do not have (easy) access to SPSS, so I cannot check exactly that I got what I should have. Thanks, G??ran And, I'm using R-2.0.1 on Debian testing (from source).