similar to: (no subject)

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

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 ;

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.

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

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

Dear sir, I would ask if there are in R some code to generate a random sample from a mvariate student distribution like that one wich generate the multivariate normal one i mean( rmvnorm(n, mu, sigma) Second question : if R can plot density 3Dcurve I don't mean de histogram but de hole density function(normal for example). I use a windows version of The R software Thank you in advance wiyh

Hi, I'm trying to follow the R-admin instructions for using MKL10 as the external BLAS compiling R-2.6.2 under Linux on a RH EL head node of a cluster. The configure process seems to have problems when it checks for dgemm in the BLAS. I'm using configure as: ./configure CC=icc F77=ifort --with-lapack="$MKL" --with-blas="$MKL" where $MKL is defined as in R-admin

Hi Thanks for replying my question! What really interested me is that the package providing some complex form sampling, such as wishart, multinomial, dirichlet. And others for example conditional beta distribution confining the random variable in the interval (a, b). Since these concept are widely used in the baysian, I wonder whether somebody has already written this package. Thanks! Best

Hello All, Has anyone had success building 64-bit R on Solaris 10 Sparc with the Oracle Studio Compiler suite? Could anyone start give me a clue. I have tried to no avail. Thanks, Nynese Nynese Tinsley, BSEE, MSCIS UNIX Systems Analyst Harvard School of Public Health Center for Biostatistics in AIDS Research 651 Huntington Ave, FXB 614 Boston, MA 02115 617-432-3244 office#

Hello there, I just changed to FreeBSD platform, and want to install R on it. I use FreeBSD 5.0 and install nearly all packages on the machine. When I use ports to install R, (cd /usr/ports/math/R-letter, and then type make) I got the following error information. ../../../../library/methods/libs/methods.so is unchanged dumping R code in package 'methods' Fatal error: The X11 shared

Dear All, Is there any function in R for computing "multivariate chi-square distribution"? How about "multivariate gamma distribution"? I appreciate any comment on this subject. Thank you, Amin Zollanvari PhD student Department of Electrical and Computer Engineering, Texas A&M University, College Station, TX

Good Afternoon R Help! On of my users is seriously thinking about moving to R from Splus. In his testing, he has already found some benefits to it. He has, however run into one issue, calling his Fortran functions in R. Previously, in Splus, he would do the following: 1. Compile via f90 the *.f program 2. link the object file to the Fortran IMSL libraries via the 'ld' command 3.

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.

R version: 1.7.1 OS: Red Hat Linux 7.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)

If the observations are normally distributed and the 2xk design is balanced, theory requires that the tests for interaction and row effects be independent. In my program, appended below, this would translate to cntT (approx)= cntR*cntI/N if all R routines were functioning correctly. They aren't. sim2=function(size,N,p){ cntR=0 cntC=0 cntI=0 cntT=0 cntP=0 for(i in 1:N){

The problem occurs when the sample odds ratio is Inf, such as in the following example. Given the fact that both upper bounds of the two 95% confidence intervals are Inf, I would have expected that the two lower bounds be equal, but they aren't. x <- matrix(c(9,4,0,2),2,2) x # [,1] [,2] #[1,] 9 0 #[2,] 4 2 rbind("two.sided.95CI"=fisher.test(x)$conf.int,

Hi, This is not a HW problem, sadly: I was last in a classroom 30 years ago, and can no longer run off to the instructor :-( I apologize but I cut and paste the wrong snippet earlier and made a typo in doing so, but the result is the same with the more appropriate snippet. require(mvtnorm) require(ggplot2) set.seed(1234) xx <- data.frame(rmvt(100, df = c(13, 13))) v <- ggplot(data = xx,

When I run the gam function as part of an optimization and do the optimization many times using a loop, I'm finding that memory use increases over time (based on simply monitoring top). Below is some example code that involves varying the penalty parameter in gam, trying to find the value that gives exactly 50 edf for a simple smoothing problem. I thought I would post to the list to see if