[R] Extract an invertible submatrix

Dear R helpers, I am looking for an efficient way to extract (any) one of the biggest invertible submatrix. I have a rectangular matrix A (p x n), with rank k <= min(p, n). I would like to get a submatrix (k x k) invertible, or even better, the list of rows and columns of A which would form th

[R] Question on LIMMA analysis with covariates and some missing data

t;Moderate", "Severe", "SevereSmokers") design<-model.matrix(~0 + asthma + age + sex) contrast.matrix<-makeContrasts(Control-Moderate, Control-Severe, Control-SevereSmokers, levels=design) fit<-lmFit(data, design) fit2<-contrasts.fit(fit, contrast.matrix) fit2

[R] Elliptic Fourier Analysis and PCA

Hi, I use packages "momocs" and "ade4" to dertermine morphometry of woodpeckers' cavities entrances. I have few questions about the PCA (non-normalized) (see image enclosed) :   1) *What do show the confidence ellipses ?*  Some papers say that ellipses include 95 % of the data, but it seems

[R] gmp: Error in solve.bigz(B) : System is singular

With a matrix such as C I do not have any problem: library(gmp) C V1 V2 V3 V4 V5 V6 V7 [1,] 1 0 0 0 1 0 0 [2,] 0 1 0 0 0 1 0 [3,] 0 0 1 0 0 0 1 [4,] 0 0 0 1 0 0 0 [5,] 0 0 0 0 1 0 0 [6,] 0 0 0 0 0 1 0 [7,] 0 0 0 0 0 0 1 solve.bigz(C) [

[R] [R-pkgs] sensitivity

n rearranged in order to gain computing memory. Some C code additions will be done later. -- Bertrand Iooss CEA Cadarache DEN/CAD/DER/SESI/LCFR 13108 Saint Paul lez Durance France Laboratoire de Conduite et de Fiabilité des Réacteurs [EMAIL PROTECTED] Tel : (33/0)

[R] Random Forest

Hello Is there exists a package for multivariate random forest, namely for multivariate response data ? It seems to be impossible with the "randomForest" function and I did not find any information about this in the help pages ... Thank you for your help

[R] Random Forest and for multivariate response data

Hello Is there exists a package for multivariate random forest, namely for multivariate response data ? It seems to be impossible with the "randomForest" function and I did not find any information about this in the help pages ... Thank you for your help

[R] Random Forest

Hello Is there exists a package for multivariate random forest, namely for multivariate response data ? It seems to be impossible with the "randomForest" function and I did not find any information about this in the help pages ... Thank you for your help

[R] step and glm : keep all models

Hi everyone, I'm running the following command: > step(mydata.glm,directions="both",trace=T) which returns the model with the lowest AIC. Now I'd like to get the list of all the models (or at least the formulas) that were used in between. I've noticed the "keep" option, but couldnt find any hel