search for: r46864

...39;, ylab='Ratio',fit=c("linear","quadratic")) ) # Option 2 require(scatterplot3d) with(res, scatterplot3d(h4,ypred,h11, xlab='Chr4',zlab='Chr11',ylab='Ratio') ) # ------------- # Session info # ------------- R version 2.8.0 Patched (2008-11-08 r46864) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] tcltk stats graphics grDevices utils datasets methods base other att...

Could any one help ? I start to learn the glmnet package. I tried with the example in the manual. x=matrix(rnorm(100*20),100,20) y=rnorm(100) fit1=glmnet(x,y) When I tried to fit the model, I received the error message: Error in validObject(.Object) : invalid class "dgCMatrix" object: row indices are not sorted within columns Thank you very much!

My two matrices are roughly the sizes of m1 and m2. I tried using two apply and cor.test to compute the correlation p.values. More than an hour, and the codes are still running. Please help to make it more efficient. m1 <- matrix(rnorm(100000), ncol=100) m2 <- matrix(rnorm(10000000), ncol=100) cor.pvalues <- apply(m1, 1, function(x) { apply(m2, 1, function(y) { cor.test(x,y)$p.value

...ystem.time(apply(combs,1,function(x) > rcorr(m1[,x[1]],m2[,x[2]])$P[2])->pvalues2) > # user system elapsed > # 7.00 0.00 7.02 > > > HTH, > > Jorge > > > # ------------- Session Info ---------------------------- > R version 2.8.0 Patched (2008-11-08 r46864) > i386-pc-mingw32 > > locale: > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.1252;LC_MONETARY=English_United > States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 > > attached base packages: > [1] stats graphics grDevices utils...