similar to: subsetting a data.frame to the 'unique' of a column

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

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

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

Hi, attach(data) if (TREUE[TREUE ==11] && EMPFEHL[EMPFEHL < 11]) { TREUE <- EMPFEHL } what's wrong ( TREUE ==EMPFEHL works not,too )? many thanks & regards,Christian -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

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

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

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

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

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

Hello I would like to compare the results obtained with a classical non parametric proportionnal hazard model with a parametric proportionnal hazard model using a Weibull. How can we obtain the equivalence of the parameters using coxph(non parametric model) and survreg(parametric model) ? Thanks Virginie

Hi, guys, I need to simulate Pareto distribution. But I found 'rpareto' didn't exist in R. And it seems that Pareto distribution don't have mathematical relationships with other distributions. What can I do? Thanks a lot. Ni --------------------------------- [[alternative HTML version deleted]]

I'm using Redhat linux 7 and R 1.2.2. When trying to run a program I get the following error: /usr/lib/R/bin/R.bin: error while loading shared libraries: /usr/lib/R/bin/R.bin: symbol __sysconf, version GLIBC_2.2 not defined in file libc.so.6 with link time reference. What does this mean and how can I fix it? Thank you ===== Janine Porter Paladin Data Systems 206-910-8370

Hi All, We're interested in finding other users of R in pharmaceutical companies. Please reply if you are or know of someone who is. -Greg Gregory R. Warnes, Ph.D. Pfizer Global Research and Development LEGAL NOTICE Unless expressly stated otherwise, this message is confidential and may be privileged. It is intended for the addressee(s) only. Access to this E-mail by anyone else is

Hello, A small problem I can't solve > p <- 0.0001 > p [1] 1e-04 How can I force the printout of p to 0.0001? I have tried 'format', 'round', 'signif', 'print' in different combinations without success. Fredrik Lundgren

Dear R experts, How should lambda and gamma (with std.errors) be calculated for a weibull model with age as an independent predictor? I have assumed that this can be done with survreg with e. g. (summary(survreg(Surv(time, status) ~ age, dist = 'weibull')) ) and predict.survreg with e.g. (predict(model, se.fit = T, newdata = data.frame(age = seq(50, 80, 5)) but unfortunately I'm

Dear R users, I have recently started using R. I need to execute a FORTRAN file within R. Does it allow to call executable (or FORTRAN code) file within R ? Thank you very much in advance. Kind Regards Seyit Ali KAYIS ___________________________________________________________________________ THE FOLLOWING WARNINGS ATTACH TO THIS COMMUNICATION This material is intended

I'm having trouble with save and load in R 1.4.0 (on RH 7.2) I originally noticed this with a large R object, but it happens with a tiny one as well: > a <- c(1,2,3) > save(a,file="test.RData") > a1 <- load("test.RData") > print(a1) NULL The file "test.RData" is created, and it does have something in it. When I tried saving a