similar to: Difference between two correlation matrices

Hi I'm seeing some "odd" behaviour with cbind(). My code is: > cat <- read.table("cogs_category.txt", sep="\t", header=TRUE, quote=NULL, colClasses="character") > colnames(cat) [1] "Code" "Description" > is.factor(cat$Code) [1] FALSE > is.factor(cat$Description) [1] FALSE > is.factor(rainbow(nrow(cat))) [1]

Hello list, I have two vectors x and x2: x=runif(10) x2=runif(10) and one vectors with their names : my.names=c("x","x2") I would like to cbind these two vectors using their names contained in the vector my.names. I can create a string with comma ncomma=paste(my.names,collapse=",") and now, I just need a function to transform this string into a adequate

I am not sure a previous e-mail reached the list (no mail aknowledgement from R-boundle etc.). The question was how to write polygon or segment coordinates into a shapefile set or any other ArcGIS supported format. The library shapefiles seems to do something but the documentation is a bit beyond of my mind.... and I cannot get the meaning of the functions write**** and its application to the case

Hello, suppose I have a list with matrices: a=list(x1=matrix(rnorm(10),5,2),x2=matrix(rnorm(10),5,2),x3=matrix(rnorm(10),5,2)) I want to compute for all combination of xi and xj (x1,x2 x1,x3 and x2,x3) a value. This value is given for the pair x1,x2 by trace(x1%*%t(x1)%*%x2%*%t(x2)) / trace(x1%*%t(x1))*trace(x2%*%t(x2)) I know that product matrices t(xi)%*%xi can be obtained by:

Hello, I have a question that is not directly related to R ... but I try to do it in R ;-) : I would like to generate a matrix Q satisfying (for a given Z, X and W) the two following conditions: t(Q)%*%Q=Z (1) XQ=W (2) where: Q is m rows and r columns X is p rows and m columns D is p rows and r columns C is r rows and r columns with m>p,r e.g: m=6, p=2 r=3

Hello list, Sorry, these questions are not directly linked to R. If I consider an indefinte real matrix, I would like to know if the symmetry of the matrix is sufficient to say that their eigenvectors are real ? And what is the conditions to ensure that eigenvectors are real in the case of an asymmetric matrix (if some conditions exist)? Thanks in Advance, St?phane DRAY

Hi all, I am trying to compare two distance matrices with R. I would like to create a XY plot of these matrices and do some linear regression on it. But, I am a bit new to R, so i have a few questions (I searched in the documentation with no success). The first problem is loading a distance matrix into R. This matrix is the output of a the Phylip program Protdist and lookes like this: 5

Hello I have a data base consisting of soil parameters, and tree species densities. The vcd (visualizing categorical data) package includes the ternaryplot function which plots the gravitation center of 3 prameters. I'd like to find a similar function for use with continuous data rather than categorical. Does anyone know of a function suitable for this objective? Thanks Steve

Dear R-users How can I get the eigenvalues out of an lda analysis? thanks a lot christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

Hi, I am looking for some code for Moran's I. Has anyone previously done this? I have been unable to find it in the search engines. James

Dear All, I am using the scatter.dudi finction in the 'ade4' package to produce correspondence analysis (nice) plots. I do not seem to figure out how to plot the raw coordinates only -- or column coordinates only. I would appreciate any help in doing that. Here is the example I am following -- from the package. data(banque) > banque.acm <- dudi.acm(banque, scann = FALSE, nf =

Does R have a function to calculate Bray-Curtis distance measures, which is probably one of the most frequently used and recommended dissimilarity measures in ecology? It isn't mentioned in dist().

Bonjour à tous, Je débute dans l'utilisation de R, et suis à la recherche de toute piste me permettant l'introduction à ce logiciel. J'ai pour le moment la fiche d'Emmanuel Paradis: "R pour les débutants", mais j'aimerais savoir où je peut trouver d'autres documentations notamment pour le traitement statistique de données (Anova, modèle linéaire et modèle

Hello I have some kind of problem with R. I want to do a pca (I work on R and ade4 library) but R doesn't wan to realize my pca because of my data: I have some binary variable (presence/absence of birds, hole, vegetals...) and metric variable (cliff 's size...) and each time I have an error message because of the mix of the 2 type of variable. So, what can I do? Famille SALMON

Is there a command to get cluster criterion for the cluster methods? SAS has its criterion, but I prefer to do it in R. If there is not a command is there code to produce criteria to choose the number of clusters? Adrian Katschke Statistics Grad Student University of Nebraska-Lincoln [[alternative HTML version deleted]]

Numeric data that is part of a mixed type data frame is converted into character. How can I tell apply to maintain the original class of a column and not convert it into character. I would like to do this of the vector and not inside the apply function individually over each element. Consider the following two scenarios, in the second column 'a' maintained its class while it lost its

Hi, i'm having trouble using Statis in ADE4 package. I want to study a matrix with 18x414...could any body help me?!?!' Urgent Thanks Miguel Ribeiro Portugal [[alternative HTML version deleted]]

Is it possible to make disappear the id numbers from scatter.dudi (mc analysis) ? a <- as.factor(c(1, 2, 3, 2, 1)) b <- as.factor(c(3, 2, 3, 1, 1)) x <- as.factor(c(1, 2, 2, 1, 3)) y <- as.factor(c(2, 2, 3, 1, 1)) dat <- data.frame(a=a, b=b,x=x,y=y) summary(dat) dat require(ade4) dat.acm <- dudi.acm(dat, scann = FALSE, nf = 2) scatter.dudi(dat.acm) Thank you very much !

Does anyone know R is to be cited in a publication? I've looked everywhere on the website and cannot find this. Thanks, Heather -- Heather Maughan Department of Ecology and Evolutionary Biology Biosciences West 310 University of Arizona Tucson, AZ 85701 Phone: 520-626-5108 Fax: 520-621-9190 hmaughan at u.arizona.edu

Hello list, I have used scan function to import data into R. I have done some analysis and find strange results. I have found my problem : when importing data with scan, this can slightly modify the data : > write(c(0.251,3.399,-0.481,0.266),"essai.txt") > scan("essai.txt") Read 4 items [1] 0.251 3.399 -0.481 0.266 > print(scan("essai.txt"),17) Read