search for: cendoya

Hi R-Helpers: I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows 98),and I get the following error: > pd5 = GLMM(nplant~sitio+ + fert+ + remo+ + sitio:fert+ + remo:sitio+ + remo:fert+ + remo:fert:sitio + data=datos, + family=binomial, +

...tualice 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 Biometría Profesor Adjunto Cátedra de Estadística y Diseño Facultad de Ciencias Agrarias Universidad Nacional de Mar del Plata ______________________________ ___________________________________________________________________________ Aviso: ===== El contenido del presente e-mail y s...

Hi R-Helpers: I?m trying to fit the same linear model to a bunch of variables in a data frame, so I was trying to adapt the codes John Fox, Spencer Graves and Peter Dalgaard proposed and discused yesterday on this e-mail list: for (y in df[, 3:5]) { mod = lm(y ~ Trt*Dose, data = x, contrasts = list(Trt = contr.sum, Dose = contr.sum)) Anova(mod, type = "III") } ## by John Fox or for

Estimados me ha surgido una duda con un simple stem and leaf el mismo es el siguiente: > stem.leaf(Datos$ej1.21) 1 | 2: represents 0.12 leaf unit: 0.01 n: 40 LO: 0.72 0.85 3 10 | 9 11 | 4 12 | 4 5 13 | 7 7 14 | 710 (*) 9 15 | 18 16 16 | 3447899 (6) 17 | 045599 (***) 18 18 | 2568 14 19 | 23710 (**) 10 20 | 389 7 21

Buen día a todos el problema que tengo es el siguiente: tengo una matriz 0 0 0 1 5 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 1 6 0 0 0 0 0 0 0 0 0 1 4 0 1 0 7 2 9 0 1 necesito hacer el primer conteo de números hasta que encuentre el primer cero. para la fila 1: tiene 1 y 5 son dos elementos fila 2: solo tiene 1 elemento fila 3: tiene 1 fila 4: 1 y 6 dos elementos fila 5: no tiene elementos distintos de

...gt; "Re: Contents of R-help-es digest...". Además, por favor, incluya en > la respuesta sólo aquellas partes del mensaje a las que está > respondiendo. > > > Asuntos del día: > >   1. Re: Votaciones concurso logotipo (Usuario R) >   2. Re: Paquete lsmeans (Gabriela Cendoya) >   3. Leer datos de Unicode (Usuario R) >   4. Re: Leer datos de Unicode (Carlos J. Gil Bellosta ) >   5. Re: Leer datos de Unicode (Juan Jos?é Vidal Agust?ín) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 21 Ap...

Hi: Does R have a function as gsorth is SAS, that perform a the Gram-Schmidt orthonormal factorization of the m ?n matrix A, where m is greater than or equal to n? That is, the GSORTH subroutine in SAS computes the column-orthonormal m ?n matrix P and the upper triangular n ?n matrix T such that A = P*T. or any other version of Gram-Schmidt orthonormal factorization? I search the help, but I

Hi: I was trying to download the package MCMCglmm and it give me this message: Aviso: dependency ?tensorA? is not available probando la URL 'http://www.vps.fmvz.usp.br/CRAN/bin/windows64/contrib/2.11/MCMCglmm_2.06.zip' Content type 'application/zip' length 8988896 bytes (8.6 Mb) URL abierta downloaded 8.6 Mb package 'MCMCglmm' successfully unpacked and MD5 sums checked my

Hi: I'm using function lmer from package lme4, and I get this message: " There were 12 warnings (use warnings() to see them)" So I checked them: Warnings 1 to 11 said: 1: optim returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: "LMEoptimize<-"(`*tmp*`, value = structure(list(maxIter = 50, ... and Warning 12 said: 12: IRLS iterations for glmm did

Estimados Escribo para consultar sobre el uso de modelos mixtos anidados. Los datos que estoy analizando provienen de censos de malezas en cuatro tipos de paisajes de la región pampeana, en los que seleccioné al azar igual número de lotes agrícolas cultivados con tres cultivos (maíz, soja y trigo-soja). En cada lote censé el número de especies de malezas en tres posiciones: el alambrado, el borde

Hi! I have a table representing both absolute and relative frequency, for example (code to get example data under the signature): Italy Germany absolute 100 105 relative 40.51 41.18 How can I print a different number of decimal digits? I try to transform to as.character, but cells result aligned to left and I don't like this solution. At the

Hi, I did the analysis of variance of a split-plot and the effect of treatment was significant. I would like compare treatment means using Tukey. I can't extract the mean square to apply HSD.test to use in agricolae package. anava = aov(ganhos ~ Blocos + Trat*Supl + Error(Blocos/Trat)) names(anava) summary(anava) require(agricolae) HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05)

...ry in the previos replay I forgot monday' scale so weight should be: weight <- c(0.91,1,1,1,1,1.21,1.22) names(weight) <- c("mon","tue","wed","thu","fri","sat","sun") ----- Original Message ----- From: "Gabriela Cendoya" <gcendoya at balcarce.inta.gov.ar> To: "Robin Williams" <> Cc: <> Sent: Wednesday, July 23, 2008 4:14 PM Subject: Re: [R] Using if, else statements > Hi Robin: > I think you can avoid the loops doing this: > > my.df<-data.frame(d.o.w...

Hi: I think I found an error in plot.lm with the option which=5, of course I can be wrong , as usually happen, but I had work on it for a while and show it to some other people that work with R, and so far I don't see what I can be interpreting wrong. I also worked over the plot.lm's code and change some lines to get what I call "the right plot", if any body is

Suppose I have the following data: tmp <- data.frame(var1 = sample(c(0:10), 3, replace = TRUE), var2 = sample(c(0:10), 3, replace = TRUE), var3 = sample(c(0:10), 3, replace = TRUE)) I can run the following double loop and yield what I want in the end (rr1) as: library(statmod) Q <- 2 b <- runif(3) qq <- gauss.quad.prob(Q, dist = 'normal', mu = 0, sigma=1) rr1 <- matrix(0,

When I import a simple dataset, run LDA, and then try to use the model to forecast out of sample data, I get a forecast for the training set not the out of sample set. Others have posted this question, but I do not see the answers to their posts. Here is some sample data: Date Names v1 v2 v3 c1 1/31/2009 Name1 0.714472361 0.902552278 0.783353694 a 1/31/2009 Name2 0.512158919 0.770451596

Hola lista, tengo una base de datos muy grande de un datalogger. Consiste en un vector con 0's y 1's, y tengo que ver cuando hay 600 o más unos seguidos. Se me ocurrió hacerlo con un loop for. Pero tarda demasiado. También intente usar which para que seleccione solo los 1 para empezar a sumar pero no hay gran diferencia. A alguien se le ocurre alguna solución para hacerlo más rápido? Acá

Hola, mi pregunta es la siguiente: Tengo un modelo genético aditivo, donde los individuos estan emparentados. Quiero que la matriz de parentesco sea la matriz de coeficientes de las correlaciones entre genotipos. En SAS, con PROC MIxed podia agregarla como un archivo y usarla en la sentencia random...No encuentro como hacerlo en R, alguna idea?? En resumen, puedo imponerle a los efectos

I have just started looking at R, and are getting more and more irritated at myself for not having done that before. However, one of the things I have not found in the documentation is some way of preparing output from R for convenient formatting into something like MS Word. An example: If you use summary(lm(....)) you get nice output. However, if you try to paste that output into the word

Podrías mandar un summary de tu conjunto de datos o un str para ver por que no puede estimar esos efectos, también sería útil el resultado de with(tusdatos, table(Landscape,Crop,Position)) Gabriela