similar to: contour labelling [was "Re: image ()"]

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

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

I am using rw0641 with Windows 98. To list just the influential repetitiones that result from "influence.measures", I am using the input result <- lm(y~x) and the code from the example in the help for "influence.measures" INFLM <- function(result){ inflm <- influence.measures(result) which(apply(inflm$is.inf,1,any)) } It works fine up to now with the

I'm using R for calculations and plots in a statistical paper. Is there a canonical reference to use in the bibliography? (Or is it not common to acknowledge computer programs in journal papers?) -- Bj?rn-Helge Mevik <bhm at math.uio.no> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

> Hi, > I recently downloaded R on window and running Rgui on NT. I have never use > R before and trying to learn R by following An introduction to R. > I generated a table of 900x5 in excel and saved as tab delimited txt file. > I then do read.table but I keeping getting the following message > > data <- read.table("ML4mm25.txt", header=T, sep=" ")

The lm() function in R seems to handle the inversion of singular X'X matrices (where there is collinearity between regression inputs) in a way where one of the inputs is dropped and this also seems to be the default behavior in SAS (please let me know if i'm wrong about this). In some other packages (i.e. octave ols() function) the pseudo inverse is computed where singular values less

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

Hi all, I've got at really basic Q for you all: How do I produce a 3D plot in R? Like the one in the screenshot on the R-project wab-page. What is the command? /Fredrik -- ----------------------------------------------------------------- Fredrik Karlsson, Research Student Department of Philosophy and Linguistics Ume? University S-901 87 UME? SWEDEN Tel: +46 90 786 56 84 Web:

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

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

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

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

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

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

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

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