similar to: Snow parLapply

Dear R-helpers, I have two data frames, op and em4: > str(op) 'data.frame': 37 obs. of 5 variables: $ m : num 0.202 0.336 0.122 0.139 0.14 ... $ lln : num 0.798 0.643 0.863 0.835 0.823 ... $ rrn : num 0.789 0.702 0.894 0.895 0.923 ... $ asym2: num 0.177 0.304 0.108 0.187 0.274 ... $ asym3: num 0.0755 0.0975 0.0818 0.0651 0.13 ... > str(rownames(op)) chr

Hola compañeros de la lista. Qué tal. Tengo un análisis de componentes principales, en el que se evalúan aproximadamente 1000 variables. Usando la función dudi.pca e inertia.dudi obtengo una cantidad de información sobre la influencia de las variables sobre los dos componentes principales. Me gustaría saber si existe alguna función que sobre esta información me arrojara la lista de

Hi all I?m trying to tests the significance of loadings from a ordination of 46 variables (caategorical, ordinal and nominal). I used dudi.mix from ade4 for the ordination. A years ago Jari Oksanen wrote this script implementing Peres-Neto et al. 2003 (Ecology) bootstraping method: netoboot <- function (x, permutations=1000, ...) { pcnull <- princomp(x, cor = TRUE, ...) res

Hola compañeros de la lista, qué tal. Estoy haciendo un análisis de componentes principales utilizando las funciones "dudi.pca" (paquete "ade4") y "PCA" (paquete "FactoMineR"). Sucede que al comparar las coordenadas de cada individuo que obtiene cada función, las que corresponden al segundo componente principal tienen idéntica magnitud pero con

Hola! No he podido consultar la doc. del paquete ade4, algo debe estar caído en CRAN ahora mismo. Dos cosas sobre la metodología -aun desconociendo los detalles de cómo lo hace ade4: El output de un PCA, los "pesos" de cada variable en las dimensiones de los componentes se interpretan como correlaciones, a mayor valor absoluto mayor asociación variable-componente. Ahora, como tales

Instalo ade4 correctamente y no me abren los datos como por ejemplo data(bsetal97) ¿Qué piensan de eso? De: r-help-es-bounces en r-project.org [mailto:r-help-es-bounces en r-project.org] En nombre de Francesc Carmona Enviado el: Tuesday, October 01, 2013 7:30 AM Para: r-help-es en r-project.org Asunto: Re: [R-es] Análisis de componentes principales con ade4 y FactoMineR Por definición

Hi R-developers In the package Parallel, the function parLapply(cl, x, f) seems to allow transmission of only one parameter (x) to the function f. Hence in order to compute f(x, y) parallelly, I had to define f(x, y) as f(x) and tried to access y within the function, whereas y was defined outside of f(x). Script: library(parallel) f <- function(x) { z <- 2 * x + .GlobalEnv$y # Try to

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

Hi all, thank you for your answers; i've tried both cca from vegan library, and dudi.coa from ade4 library; one last question: my deal is mainly with contingency tables, like the one i'm posting here acciaieria<-c(.41,.02,.44,.04,.09) laminatoio<-c(.34,.28,.26,.01,.11) fonderia<-c(.48,.05,.34,.08,.05) leghe<-c(.45,.19,.25,.03,.08)

invoking parLapply without a cluster fails to find a previously registered cluster > library(parallel) > setDefaultCluster(makePSOCKcluster(2)) > parLapply(X=1:2, fun=function(...) {}) Error in cut.default(i, breaks) : invalid number of intervals This is because in parLapply length(cl) is determined before defaultCluster(cl) is called. By inspection, this appears to be true of

Dear R users, I have some problems with the parLapply function from the "parallel" package: I use parLapply on a pretty big R object without changing the object within the called function. If I execute parLapply alone, everything works fine. It seems that the object resides only once in the memory. But if I use the same call within another function, the object seems to be multiplied to

Hello everybody, I've got a bit of a problem with parLapply that's left me scratching my head today. I've tried this in R 2.11 and the 23 bit Revolution R Enterprise and gotten the same result, OS in question is Windows XP, the package involved is the snow package. I've got a list of 20 rain/no rain (1/0) situations for these two stations i and j, all the items in this list look

I am slightly confused by the way the standard output is redirected in a R snow cluster environment. I am using parLapply from the snow package to execute a function on my MPI/LAM cluster. How can I redirect standard output (produced using "cat") from this function back to the terminal where I invoked it? I intend to transmit some status information in advance to the final result of the

Hi, The following function use to work, but now it doesn't giving the error "> CallSnow(, 100) Using snow package, asking for 2 nodes 2 slaves are spawned successfully. 0 failed. Error in checkForRemoteErrors(val) : 2 nodes produced errors; first error: no applicable method for 'lapply' applied to an object of class "list" ". Where this is the

Dear R-users I was running into problems with my R code trying to run clh sampling (clhs package) in parallel mode (=on various data sets simultaneously). Here is the code (which I developed with some help:)): ****************************************** library("clhs") library("snow") a <- as.data.frame(replicate(1000, rnorm(20))) b <- as.data.frame(replicate(1000,

Dear R users, How can I extract the residuals from a pcaiv/rda in ade4? In Vegan there is the residuals() function, giving the approximation of the original data from the unconstrained ordination Is there something similar in ade4? Nikos

Hola Argel, Puedes encontrar un par de ejemplos que cómo conseguir una valoración (qualitativa y cuantitativa) de la contribución de cada variable a cada componente principal, en el libro "Exploratory Multivariate Analysis by Example Using R". Este libro está basado en la librería "FactoMineR". Los ejemplos están en el capítulo primero. No he mirado como se haría esto con

> Hi > I'm Ph.D student and I need an R code to compute the chi square diistance between n profile rows in a matrix. > > could you help me please? > Thanks > > Paola > ____________________________________________________________

Hello, I not have attached in this e-mail the zipped list of matrices I am using because it has 1 meg once zipped and anyway we cannot send attached files on r-help mailling list. First, after running the code that is written bellow, I realized that the printout of dudi.pca gives me for both of my matrices $cw length that are unequal to either respective $rank or $eig lengths. Is that

Dear all, (Mac OS X 10.4.11, R 2.6.0) I have a quantitative dataset with a lot of Na?s in it. So many, that it is not possible to delete all rows with NA?s and also not possible, to delete all variables with NA?s. Is there a function for a principal component analysis, that can deal with so many NA?s. Thanks in advance Birgit Birgit Lemcke Institut f?r Systematische Botanik