search for: hivnet

...ot;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

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!","Hi!"),pch=1:2,lty=1:2,trace=T) xchar= 0.05178 ;

...-0.4205906 > pc2$center # All NA's (should be the same as pc1$center) [1] NA NA -- 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

...cy!! 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="n",xlab=&...

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, cor = TRUE) #THIS WORKS USArrests[1,3] <- NA princomp(USArrests, cor =

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

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 same answer. It should remove all the blanks when I use the vector y, but it

...se(15-9,13-n11)*theta^n11) #[1] 0.5608724 0.3891316 0.0499960 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

...more problems in other cases. Perhaps there is a better fix. 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

...stant Professor of Biostatistics Biostatistics/Univ. of Washington (Th) Box 357232 206-543-1044 (3286=fax) Center for AIDS Research/HMC/UW (M/F) Box 359931 206-731-3647 (3693=fax) VTN/SCHARP/FHCRC (Tu/W) Box 358080 206-667-7025 (4812=fax) rossini@(biostat.washington.edu|u.washington.edu|hivnet.fhcrc.org) http://www.biostat.washington.edu/~rossini -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "bod...

Hi, Are there some packages which can generate multi-normal, multi-t, etc multivariate sampling? thanks! Best wishes, Peng ******************************* Peng Zhang Department of Biostatistics Harvard School of Public Health 655 Huntington Avenue Boston, Massachusetts 02115 ******************************* I believe I can fly I believe I can touch the sky

I'm doing prediction of the survival cases using gene expression profiles(Affymetrix chips). Can somebody tell me how to use the Cox PH model to select genes and make a prediction of survival? Thanks. Guangchun

I have a data file where each entry is indexed by the time in seconds since epoch (e.g. 1046315697). Is there an easy way to convert this time value into a more friendly time (such as Month-Year) when plotting it? I searched through the manual, mailing lists, and functions like as.POSIXct and strptime, but didn't find what I need. Thanks, Sharad.

Hello, does R have functions for setting and testing bit values? I want to conserve memory for storing presence/absence data for large multiple arrays within a single array, using element values like present[x,y] <- ntharray[x,y]*(2^n) where presence is 1, non-presence is 0 and n is the nth array e.g. 1*(2^0) + 0*(2^1) + 0*(2^2) + 1*(2^3) + 0(2^4) for storing the value 9 for presence in

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

...stant Professor of Biostatistics Biostatistics/Univ. of Washington (Th) Box 357232 206-543-1044 (3286=fax) Center for AIDS Research/HMC/UW (M/F) Box 359931 206-731-3647 (3693=fax) VTN/SCHARP/FHCRC (Tu/W) Box 358080 206-667-7025 (4812=fax) rossini@(biostat.washington.edu|u.washington.edu|hivnet.fhcrc.org) http://www.biostat.washington.edu/~rossini -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "bod...

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 error (not a segmentation fault). library(survival) data <-

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,1,1,1,NaN),exclude=NULL) 1 2 NaN 3 1 1 1 #For sake of

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

...stant Professor of Biostatistics Biostatistics/Univ. of Washington (Th) Box 357232 206-543-1044 (3286=fax) Center for AIDS Research/HMC/UW (M/F) Box 359931 206-731-3647 (3693=fax) VTN/SCHARP/FHCRC (Tu/W) Box 358080 206-667-7025 (4812=fax) rossini@(biostat.washington.edu|u.washington.edu|hivnet.fhcrc.org) http://www.biostat.washington.edu/~rossini -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body...