search for: burrard

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

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

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

...[4,] 0.1522791 -0.4205906 [5,] 0.1522791 -0.4205906 [6,] 0.1522791 -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

...tial) likelihood. Any comment on this issue would be much appreciated. Sincerely, Jerome Asselin R 1.6.2 on Red Hat Linux 7.2 Package: survival Version: 2.9-7 -- 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

...od="ML")) DF[i] <- ifelse(attr(tryerror,"class")[1]=="try-error",NA,fit.nlme$fixDF$X) } names(DF) <- 1:nseeds DF ########### -- 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

Hi, I'm trying to construct a legend which has four lines of text and associated symbols. The first two symbols need to be normal lines which vary only in colour. The second two symbols should have filled boxes. How do I suppress the fill boxes in the first two lines? J.

Hi, is there any R-package providing a Passing Bablok regression procedure? This Question seems to have occurred once or twice before, but I never saw an answer to the previous inquiries. Best regards Benjamin Otto ====================================== Benjamin Otto University Hospital Hamburg-Eppendorf Institute For Clinical Chemistry Martinistr. 52 D-20246 Hamburg

...BOUNDS IS ZERO #Error in survreg(time59 ~ ecog.ps + rx, data = newovarian, # dist = "lognormal") : # Invalid survival times for this distribution -- 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

...kernel image vmlinuz... - stops after receiving only a few blocks over tftp (system hangs) - next step would be to transfer the initrd but it doesn't get here Thanks, Mike --------------------- Mike Bethune, SCSA, SCNA, LCP, LCA Senior UNIX Engineer Fusepoint Managed Services Inc. 2323 - 595 Burrard Street P.O. Box 49336, Vancouver, BC, V7X 1L4 T. 604.687-7757 x303 F. 604.687.7761 Email:mike.bethune at fusepoint.com Web : www.fusepoint.com

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

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 ;

Hello, I have a beginner's Q: if i want to know all the functions provided by a package, what is command for that? in another word, Is there a command to list all the commands available in a packege? thanks a lot, yan

...314 choose(9,n11)*choose(15-9,13-n11)*theta^n11/ sum(choose(9,n11)*choose(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

Hi, I'm facing this problem quite a lot, so it seems worthwhile to check to see what the most efficient solution is. I've two vectors x (values ordered) and y. I've ranges x < x0, x0 <= x < x1, x1 <= x < x2, x2 <= x < x3, x > xn and want to construct a subvector yprime of y which consists of the first/last value of y whose x values are in the range. For

...() fixes this problem. I don't know if this correction would create 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