similar to: Como imponer una matriz de correlaciones en un modelo mixto?

> > Hola, a todos! Estoy ajustando un modelo componentes de varianza utilizando el paquete lme4, quiero extraer luego los valores estimados de varianza de cada factor aleatorio de este objeto lmer y guardarlo en un vector. Existe la función VarCorr, Esto es lo que hago: Mi ejemplo: tengo 3 factores NDVIaño; Máx.IH; TV # Ajusto el modelos lmer(peligro100 ~ 1 + (1|NDVIaño) + (1|Máx.IH) +

Hola a todos, quisiera saber como hacer para cambiar la escala del eje x, por ejemplo en un gráfico de puntos. Saludos _________________________________________________________________ Windows Live Messenger GRATIS: lo que faltaba en tu BlackBerry http://www.messengerentublackberry.com?ocid=WL_BB_LandPage_TagLine [[alternative HTML version deleted]]

Esta pregunta nos la hicimos muchos, alguien aproxima una respuesta a través de esta publicación: http://www.scielo.org.ar/pdf/ecoaus/v18n2/v18n2a07.pdf Saludos cordiales Teresa -- Ing. Agr. M.Sc.Teresa Boca Instituto de Clima y Agua INTA Castelar. Las Cabañas y Los Reseros s/n° 1712 Castelar Pcia de Buenos Aires, Argentina Cel: 01115-3572-9584 Tel: (54-11)-4621-0125/1684 Fax:(54-11)-4621-5663

Dear list, I'm having trouble to see how my elements on a 4 dimensional array are listed. For example, I generated the following array: junk.melt=melt(occ.data,id.var=c("Especie", "Site", "Rep", "Año"), measure.var="Pres") y=cast(junk.melt, Site ~ Rep ~ Especie ~ Año) Now, I want to be able to look at how my species (Especie) are listed, in

I need to generate output files as .csv file in a loop: let's say, I want to generate a sequence of files according to the loop number "i". the first file genearted should be file1.csv the second should be file2.csv I use: write.csv(temp, "file$i$.csv") in the loop but it did not work. can anyone give suggestions? thanks a lot. Jack

I know this seems like a very easy question (and maye it is) but I've been trying to simulate nested data and been unsucessful so far.. I want to simulate a varying number of species within a group; and then create an array to store the results of my for-loop. For example: groups<-3 species$groups<-as.integer(runif(groups,1,10)) #species per functional group ###create arrays to store

Muy buenas, quiero calcular correlaciones de Pearson entre dos variables (a,b) teniendo en cuenta una tercera (c). Para ello estoy usando una función llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce una función similar en alguna librería de R? Por otro lado, para ver si me cuadraban los resultados,

List, I want to reshape my data, but I'm not sure how to do it... it might be a simple task, but don't know which package does this. "occ.data" (see below) is how my original data are arranged, and I know that with melt() I can reshape it like "y" (see below). However, I just want to build a matrix like the "y" matrix, but with only 2 dimensions. Something

Hello ! I have a dataframe with 6 variables (A1,A2,B1,B2,C1,C2) and 1 factor (F). I would like to produce a graph consisting of 3 boxplots sets, one for every two variables (i.e A1 &A2) by the factor (F). I was looking around and I cannot figure it out, any suggestions? Best Regards, Gabriel [[alternative HTML version deleted]]

Dear list, I'm working with a large data set, where I grouped several species in one group (guild). Then I reshaped my data as shown below. Now, I just want to have "Rep" only as 1 or 0. I'm not being able to change the values of rep>=1 to 1... tried many things and I'm not being successful! > melting=melt(occ.data,id.var=c("guild", "Site",

I've been trying to reshape this database but haven't succeed at it. I tried using loops but can't get it right. I just want to reshape my database from this matrix, to the one below, with only one column of data. Year Route Point Sp1 Sp2 Sp3 2004 123 123-1 0 1 0 2004 123 123-2 0 1 1 2004 123 123-10 1 1 0 Year Route Point 2004 123 123-1 Sp1 0 2004 123 123-2 Sp1 0 2004 123

Hola a todos, Sólo por curiosidad, sabe alguien cómo enviar correos electrónicos desde R? Estoy haciendo algunas simulaciones en diferentes computadores y me gustaría tener la posibilidad de enviar un mensaje a mi correo electrónico cuando cada una termine. En R-help hubo una pregunta similar pero requiere el uso de Outlook y operar "manualmente" algunas cosas, lo cual no es

con as.numeric tampoco me van los cálculos y no conozco otra manera, pues son alrededor de 40 variables las que tiene la tabla. To: valeagui en outlook.es; cof en qualityexcellence.es; r-help-es en r-project.org From: javier.ruben.marcuzzi en gmail.com Subject: RE: [R-es] Problema con la lectura de datos Date: Tue, 10 Nov 2015 16:43:05 -0300 Estimada Valentina as.numeric(...) Mira mi correo

I fit the model fit<-lm(thresh~cond*Ne) where thresh is the reponse cond is a factor with levels a, b, and c Ne is a continuous indep var I think of this full model as having three lines: thresh as a function of Ne for each condition. Thus we have slopea, slopeb, slopec, inta, intb, intc. lm output my params ------------------------- (Intercept) inta condb intb - inta condc

Hola, Estoy tratando de leer una base de datos: tengo 39 columnas, de las cuales 38 son variables y una es el nombre de las empresas. Por lo que escribo el siguiente codigo: Variables <- read.csv("C:/Users/usuario/Documents/variables2.csv", header=TRUE, sep=";", comment.char="" , colClasses=c(Empresas="character",

Hola, sin indicar los tipos de variables si los lee, pero el problema está en que necesito calcular una matriz de correlaciones y genera el siguiente error: Error in cor(Variables) : 'x' must be numeric Date: Tue, 10 Nov 2015 20:35:40 +0100 Subject: Re: [R-es] Problema con la lectura de datos From: cof en qualityexcellence.es To: valeagui en outlook.es CC: r-help-es en r-project.org

Estimadas y estimados miembros de la comunidad. Junto con saludar, les escribo este mensaje porque estoy interesado en desarrollar un taller de estadística descriptiva con R para estudiantes de secundaria. El objetivo es que a través de un conjunto de datos puedan calcular medidas de tendencia central, dispersión, calcular correlaciones y hacer gráficos. Actualmente, he realizado un taller

Valentina, Si entiendo bien la primer columna de Variables es el nombre de Empresas (un character), no uses esta columna en la función cor, sólo las columnas numéricas. cor(Variables[ , -which( names(Variables) == "Empresas")]) Si no es ese el problema envía str(Variables) para entender mejor la estructura de los datos, y si puedes dput(Variables). Suerte, Daniel Merino El 10 de

Hi, I have a client connected to the ''net through 3 ISP''s. Have set up a Linux box to do routing and load sharing for the 3 connections. A fourth interface is connected to the LAN with private IP addresses. Am using iptables to SNAT traffic to the appropriate IP depending on the interface the packet gets routed onto. The setup looks something like this: Interface IP

Hola todos: acabo de actualizar R a R 2.10.1, y cuando actualice Rcmdr (Versión 1.5-4) me cargo la versión en inglés, no me acuerdo que hice antes para tener Rcmdr en español (trabajo con windows XP). alguien podría recordarme lo que hay que hacer? Gracias Gabriela. ______________________________ Lic. María Gabriela Cendoya Magíster en