similar to: Pregunta sencilla

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

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, +

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

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

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

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)

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

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á

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

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

Dear Dr. Bates, Thank you very much for your response. I had consulted the algorithm described in Pinheiro and Bates. However, what I don't understand (among other things) is why my two parameters appear to be estimated at different grouping levels (based on the DF values). Affect this different values of DF at the estimates parameters? The estimates fixed effects were get at the same level of

Hello, I am creating a plot/image using different data and a couple fit lines (see attached image). In the legend, I want the Default and Exponential symbol to be a line. I am using the pch command, I tried to use "-" to represent a line but does not work so I currently have set as a "1". Any thoughts or suggestions would be greatly appreciated. Below are the commands

Hi all! Is there a function that provides the VarCov matrix for a nlme objects? How I can extract the matrix for a nlme model fitted? I would appreciate any guidance Regards Lic. Gabriela Escati Pe?aloza Biolog?a y Manejo de Recursos Acu?ticos Centro Nacional Patag?nico(CENPAT). CONICET Bvd. Brown s/n?. (U9120ACV)Pto. Madryn Chubut Argentina Tel: 54-2965/451301/451024/451375/45401 (Int:277)

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 have great problems with my work in R. I look for to model the growth of fish. I have "Longitudinal data", a serie of repeated measures for each individual. Using the corresponding packages "nlme" in R. I treat to fit to the data different growth functions, wich were entered by me. Unfortunately for no it was arrived at the convergence, several error messages appeared. I

Hello, I'm new to boostrapping and I'd need some help to understand the error message that pops up when I run my script. I have a data.frame with 73 lines and 21 column. I am running a stepwise regression to find the best model using the R function "step". I apply bootstrapping to obtain model coefficients. This is my script: # "datare80" is the name of the

Hi, I' m having a hard time understanding the computation of degrees of freedom when runing nlme() on the following model: > formula(my data.gd) dLt ~ Lt | ID TasavB<- function(Lt, Linf, K) (K*(Linf-Lt)) my model.nlme <- nlme (dLt ~ TasavB(Lt, Linf, K), data = my data.gd, fixed = list(Linf ~ 1, K ~ 1), start = list(fixed = c(70, 0.4)), na.action= na.include,

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

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

Dear users I got this msg accessing s share: server:~# smbclient '\\localhost\Adm1' -U gabriela Added interface ip=192.168.0.254 bcast=192.168.0.255 nmask=255.255.255.0 Password: Domain=[EQUATORIAL] OS=[Unix] Server=[Samba 2.0.5a] smb: \> put lixo2 ERRDOS - ERRnoaccess (Access denied.) opening remote file \lixo2 The smb.conf says [Adm1] comment = Adm-Financeiro path