similar to: Building R fails

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 <-

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 <-

Full_Name: Jerome Asselin Version: 1.6.2 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) In the following example, the line type lty=2 does not show properly across the entire line. x <- c(seq(0,.5,.001),seq(.6,1,.1)) y <- rep(1,length(x)) plot(x,y,type="s",lty=2) Sincerely, Jerome Asselin

I will be giving a course in survival analysis using R (of course!) for people who know nothing about the subject (including R), but know basic statistics. I'm looking for a suitable course book. Therneau & Grambsch (2000) is an excellent book, but too much for this course. I need somthing more elementary. I have a vague memory saying that such books exist, but I cannot find any for the

Hi R-Help! My question: I have lifetime/failure data of machines with different stress levels and i think an weibull/extreme value distribution would fit this data. So I did: model1 <- survreg(Surv(lfailure)~stress,data=steel,dist="extreme") (where lfailure=log(failure)) Now I would like to do a likelihood ratio test to test the hypothesis H0: location parameters of the

This happened to me with R-1.7.0, Linux RH9: ------------------------------------------------------------------ > q() Save workspace image? [y/n/c]: y Warning messages: 1: namespaces may not be available when loading 2: names in persistent strings are currently ignored gb at lasker:~/R/test$ R R : Copyright 2003, The R Development Core Team Version 1.7.0 (2003-04-16) R is free software

Is it possible to let a function be known under two names without having two identical copies of the function body? --- G?ran Brostr?m tel: +46 90 786 5223 Department of Statistics fax: +46 90 786 6614 Ume? University http://www.stat.umu.se/egna/gb/ SE-90187 Ume?, Sweden e-mail: gb at stat.umu.se

Hi, I was wondering if someone can help me interpret the results of running weibreg. I run the following and get the following R output. > weibreg(Surv(time, censor)~covar) fit$fail = 0 Call: weibreg(formula = Surv(time, censor)~covar) Covariate Mean Coef Rel.Risk L-R p Wald p covar 319.880 -0.002 0.998 0.000 log(scale) 0.000 8.239

I?m new to R so maybe this issue has been asked before and I still could not read the complete set of past messages sent to the list. I found a weird behabiour that I will explain with a simple example. Lets consider the following block of commands: > x <- diag(c(1,4,10)) > x [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 4 0 [3,] 0 0 10 > invx <- x^-1 > invx

I'm looking for a Windows distribution of TeX that works with R, after a few years' absence from Windows. On Duncan Murdoch's Rtools page fptex is still recommended, but it turns out that fptex is "defunct" as of May 2005, see http://www.metz.supelec.fr/~popineau/xemtex-7.html So, what is suggested? TUG (tug.org) recommends something called proTeXt, which is said to be

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).

G=F6ran, At 11:04 07/12/01 +0100, G=F6ran Brostr=F6m wrote: >On Wed, 5 Dec 2001, G=F6ran Brostr=F6m wrote: > >[...] >=20 >> My real problem is how to create a data frame in a sequentially growing >> manner, when I know the final size (no of cases). I want to avoid to >> call 'rbind' many times, and instead create an 'empty' data frame in >> one

The release of R-1.9.1 is scheduled for Monday, June 21. Automatic generation of daily alpha releases should start Monday, June 7 and switch to beta status on Monday, June 14. It would be good if package maintainers could get any planned changes done as soon as possible, and test their packages carefully against the alpha/beta releases. -- O__ ---- Peter Dalgaard Blegdamsvej 3

I use the recommendations in "R coding standards", i.e., I put ;;; C (add-hook 'c-mode-hook (lambda () (c-set-style "bsd"))) ;;; ESS (add-hook 'ess-mode-hook (lambda () (ess-set-style 'C++ 'quiet) (add-hook 'local-write-file-hooks (lambda

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 >

I'm still working to add namespace support to some of my packages. I've removed the 'Depends' line from the DESCRIPTION file, and created an appropriate NAMESPACE files. Strangely, whenever I use 'importFrom(package, function)' R CMD check generates "Namespace dependencies not required" warnings . Without the import statements, the warning does not occur, but

Dear all, I have data from 1970 to 1990 for people above age 50. Now I want to calculate survival curves by age starting at age 50 using the Kaplan Meier Estimator. The problem I have is that there are already people in 1970 who are older than 50 years. I guess this is called delayed entry or left truncation (?). I thought the code would be: roland <- survfit(Surv(time=age.enter,

For anyone who might be interested, I now have a home page (finally) on which can be found the R code, data, and output for all the examples and some of the exercises in the second edition of my Models for Repeated Measurements (being printed, out in June or July), as well as the latest version of my R libraries: http://www.luc.ac.be/censtat/ Jim

For anyone who might be interested, I now have a home page (finally) on which can be found the R code, data, and output for all the examples and some of the exercises in the second edition of my Models for Repeated Measurements (being printed, out in June or July), as well as the latest version of my R libraries: http://www.luc.ac.be/censtat/ Jim

Brian D Ripley writes: Interactive rotation of perspective plots would be a great boon, as would shading options for perspective plots, including contouring on such plots. I am envisioning a division into static and dynamic plots. The static version will concentrate on rendering as nice a picture as possible - anti-aliasing (if supported) arbitrary text rotation etc, whereas the dynamic