search for: asselin

Full_Name: Jerome Asselin Version: 1.7.1 OS: Red Hat Linux 7.2 Submission from: (NULL) (24.77.125.119) Setting the parameter na.action=na.omit should remove incomplete records in princomp. However this does not seem to work as expected. See example below. Sincerely, Jerome Asselin data(USArrests) princomp(USArrests, c...

Full_Name: Jerome Asselin Version: 1.8.0 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.177.13) I would not expect a segmentation fault; perhaps an error message. > .Call("log") Segmentation fault This is always reproducable for me. Sincerely, Jerome Asselin

...s like survreg() cannot deal with lower bounds of zero, despite the fact that plnorm(0)==0 and pnorm(-Inf)==0 are well defined. Below is a short example to reproduce the problem. Does anyone know why survreg() must behave that way? Is there an alternate solution to this problem? Sincerely, Jerome Asselin library(survival) data(ovarian) newovarian <- ovarian newovarian$lower59 <- newovarian$futime-59 newovarian$time59 <- Surv(newovarian$lower59,newovarian$futime, event=rep(3,nrow(newovarian)),type="interval") survreg(time59~ecog.ps+rx,data=newovarian,dist="lognormal&quo...

Full_Name: Jerome Asselin Version: 1.6.2 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) Should be an easy fix... Consider the examble below: plot(0,0) legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T) It gives the following trace: > plot(0,0) > legend(0,0,c("Hello!&q...

Full_Name: Jerome Asselin Version: 1.6.1 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) Hello, I get a segmentation fault when I run the following code. I wouldn't expect meaningful results because my response variable contains only missing values. However, I would expect something like a regular err...

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

...oes not have to be fancy!! I have never worked with a package where the legend was not automatic. ---------- Forwarded Message ---------- Subject: Re: [R] legend() with option adj=1 Date: Wed, 21 May 2003 09:19:11 +0200 From: Uwe Ligges <ligges at statistik.uni-dortmund.de> To: Jerome Asselin <jerome at hivnet.ubc.ca> Cc: r-help at stat.math.ethz.ch Jerome Asselin wrote: > Hi there, > > I want to justify to right the text of my legend. Consider this short > reproducable example. > > x <- 1:5 > y1 <- 1/x > y2 <- 2/x > plot(rep(x,2),c(y1,y2),type...

...le of such servers, but I'd like to make sure R would work fine. The "R Installation and Administration" guide notes that there may be some problems with BLAS libraries. That's all I could find about R on AMD64. Please share your experience using R on AMD64. Many thanks, Jerome Asselin

...l". axis <- function (side, at = NULL, labels = TRUE, tick = TRUE, line = NA, pos = NA, outer = FALSE, font = NA, vfont = NULL, lty = par("lty"), lwd = par("lwd"), col = par("col"), ...) { [...] } I use R1.7.0 on Red Hat Linux 7.2. Sincerely, Jerome Asselin -- Jerome Asselin (J?r?me), Statistical Analyst British Columbia Centre for Excellence in HIV/AIDS St. Paul's Hospital, 608 - 1081 Burrard Street Vancouver, British Columbia, CANADA V6Z 1Y6 Email: jerome@hivnet.ubc.ca Phone: 604 806-9112 Fax: 604 806-9044

R version: 1.7.1 OS: Red Hat Linux 7.2 In this example, I would expect an error for the overly long variable name. This is always reproducable for me. > formula(paste("y~",paste(rep("x",50000),collapse=""))) Segmentation fault Sincerely, Jerome Asselin -- Jerome Asselin (Jérôme), Statistical Analyst British Columbia Centre for Excellence in HIV/AIDS St. Paul's Hospital, 608 - 1081 Burrard Street Vancouver, British Columbia, CANADA V6Z 1Y6 Email: jerome@hivnet.ubc.ca Phone: 604 806-9112 Fax: 604 806-9044

...el is equivalent to a FIXED effect model. Could someone explain to me why lme() still gives two standard deviation estimates? I would expect lme() to return either: a) an error or a warning for having an unidentifiable model; b) only one standard deviation estimate. Thank you for your time. Jerome Asselin > library(nlme) > simdat <- data.frame(A=1:4,Y=c(23,43,11,34)) > simdat A Y 1 1 23 2 2 43 3 3 11 4 4 34 > lme(Y~1,data=simdat,random=~1|A) <...snip...> Random effects: Formula: ~1 | A (Intercept) Residual StdDev: 12.96007 4.860027 <...snip...> > summary...

Dear All, I am given a set of files names as: velocity1.txt velocity2.txt and so on. I am sure there must be a way to read them automatically in R. It is really taking me longer to read them than to analyze them. Anybody has a suggestion to help me out with this? Many thanks Lorenzo

...bridge Downing Place Cambridge CB2 3EN, UK Telephone: +44 (0)1223 339776 Mobile: 07929 817546 Fax: +44 (0)1223 355674 E-mail: jdb33 at cam.ac.uk E-mail: james_510 at hotmail.com http://www.geog.cam.ac.uk/ccru/CCRU.html ___________________________________________ On Wed, 10 Sep 2003, Jerome Asselin wrote: > On September 10, 2003 04:03 pm, Kevin S. Van Horn wrote: > > > > Your method looks like a naive reimplementation of integration, and > > won't work so well for distributions that have the great majority of > > the probability mass concentrated in a small f...

Full_Name: Jerome Asselin Version: 1.3.0 OS: Windows 98 Submission from: (NULL) (24.77.112.193) I am having accuracy problems involving the computation of inverse of nonnegative definite matrices with solve(). I also have to compute the Choleski decomposition of matrices. My numerical problems involving solve() made me f...

Full_Name: Jerome Asselin Version: 1.3.0 OS: Windows 98 Submission from: (NULL) (24.77.112.193) I have a hard time to magnify the axis annotations in barplot() Here are some examples: #This one has no effect. barplot(1:3,names.arg=1:3,cex.axis=3,xlab="x",ylab="y") #This one magnifies the x and y label...

Full_Name: Jerome Asselin Version: 1.5.1 OS: linux redhat 7.2 Submission from: (NULL) (142.103.173.179) gsub() return different answers depending on how the input variables were created. Here is an example of code that replicates the problem. The vectors y and yy appear to be the same, but gsub() doesn't return the s...

...5,1.5,c("1,000","1,000,000"),lty=1:2,xjust=1,yjust=1,adj=1) Now, I would like to right-justify the text of the legend. As you can see, the option adj=1 does not give satisfactory results. Is this a bug or is there an easy way that I'm missing? Thanks, Jerome -- Jerome Asselin (J?r?me), Statistical Analyst British Columbia Centre for Excellence in HIV/AIDS St. Paul's Hospital, 608 - 1081 Burrard Street Vancouver, British Columbia, CANADA V6Z 1Y6 Email: jerome at hivnet.ubc.ca Phone: 604 806-9112 Fax: 604 806-9044

Full_Name: Jerome Asselin Version: 1.6.0 OS: RedHat 7.2 Submission from: (NULL) (24.83.203.63) Hello, I am trying to draw a legend on top of survival curves using plot.survfit(). As in the example below, I would like to specify a large interval for the x-axis. I can achieve such result using "xlim". However, a...

Full_Name: Jerome Asselin Version: 1.6.2 OS: redhat linux 7.2 Submission from: (NULL) (142.103.173.179) Bug or feature? Hard to say... But it sure would be nice if table() would label the frequency of NA's as it does for NaN's. Regards, Jerome > table(c(2,NA,1,1,1),exclude=NULL) 1 2 3 1 1 > table(c(2,NA,...

....2 When "covmat" is supplied in princomp(), the output value "center" is all NA's, even though the input matrix was indeed centered. I haven't read anything about this in the help file for princomp(). See code below for an example: pc2$center is all NA's. Jerome Asselin x <- rnorm(6) y <- rnorm(6) X <- cbind(x,y) pc1 <- princomp(X) pc1$scores %*% t(pc1$loadings)-X pc1$center pc2 <- princomp(X,covmat=cov(X)*5/6) pc2$scores %*% t(pc2$loadings)-X pc2$center # All NA's (should be the same as pc1$center) > x <- rnorm(6) > y <- rnorm(6)...