similar to: Possible Bug in Effects Package

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 a tod en s, estoy trabajando con modelos de regresión mixta con el paquete "nlme", concretamente con la función "lme". Una vez que hago mi modelo y el summary del mismo obtengo por ejemplo la significación (valor p) de la variable independiente, así como el valor del AIC para comparar mi modelo con otro similar, pero no sé cómo obtener un valor del grado de ajuste de la

Buenas tardes a todos. Mi problema es el siguiente: Calculo una muestra de unos índices (s) y a partir de esa muestra elijo aquellos elementos de una tabla de datos que tengo (en este caso XP4) que cumplen la condición de que el índice de XYP sea el de la muestra. Ahora comienza mi problema, la única forma que he encontrado para que me R me tome todos los elementos que necesito es creando un

Hola:   Tengo una serie de funciones que he integrado formando un paquete, y dichas funciones manejan una lista global en la cual se alojan los parámetros de configuración del análisis, pero también los data.frame de resultados numéricos.   La pregunta es la siguiente: en cada análisis algunos data.frame de resultados se actualizarán pero otros no (si el usuario no activó ese análisis en

Hola a tod en s, mi compañero y yo intentamos ver la correlación de nuestros datos mediante regresiones logísticas. Trabajamos con proporciones (1 variable dependiente y 1 independiente) mediante modelos mixtos (los datos están agrupados porque hay pseudoreplicación). Hemos usado el paquete "lme4" y la función "lmer". Encontramos "overdispersion" en el resultado

Hola amigos de R en Español, Estoy eliguiendo unas variables desde una dataframe, pero cuando creo un objeto nuevo ( con cbind) las nueva tabla (que es una matriz) se pierden los nombres de los factores. Aparecen numeros (1,2,3..). Alguien me puede decir si existe un comando para extraer las variables a un nuevo dataframe, pero sin cambiar el nombre, ni la codificacion ? Gracias, José

Dear All, probably a simple enough solution but don;t seem to be able to get my head around it...example based on a publicly available data set: mydata <- read.csv("https://stats.idre.ucla.edu/stat/data/binary.csv") mylogit <- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial") library(effects) plot(allEffects(mylogit) ? ?

Hi! I try to explain the efffect of (1) forest where i took samples's soils (* Lugar*: categorical variable with three levels), (2) nitrogen addition treatments (*Tra*: categorical variable with two levels) on total carbon concentration's soil samples (*C: *continue* *variable) during four months of sampling (*Time:* categorical and ordered variable with four levels). I fitted the

Hola a todos. En linux uso rkward o gedit con el plugin Rgedit, pero me gustaría saber qué editores se pueden usar en windows. He probado Tinn-R , ¿qué usáis vosotros? Gracias

Dear R users, I have some trouble using the allEffects() function to compute and display effect plots for a linear model. My data is quite simple, it concerns effects of 3 treatments on the tumoral volume of mice. vTum codes for the qualitative initial volume, from small to big, temps is the time in month since beginning of treatment, and S?rie codes for the batch. Data is unbalanced. >

Dear everybody! Need help with graphics. I am runnig a simple lm and then using allEffects from 'effects' package: require(effects) model<-lm(Y~A+B, data=mydataframe) I am trying to build (for each predictor - A and then B) a plot of means on Y. I was successful doing it like this - in one swoop: ml.eff<-allEffects(ml1, se=F) plot(ml.eff,ylab="Title of Y") Is it

Thanks in advance for any help. I am attempting to add points to a plot using the allEffects command in the effects package. When I try to add the points I get the following error message: Error in plot.xy(xy.coords(x, y), type = type, ...) : plot.new has not been called yet Strangely, using the code I've pasted below this has worked for me in the past however figuring out what has

Hello, I am not having success in a simple task. Using the effects package, I would like to add reference lines at probability values of 0.1 – 0.6 on a plot of the effects. The plot command works, but following up with an abline command produces the message “plot .new has not been called yet”, and of course the reference lines were not added. Looking through past R help lists, there was a

Buenos dias a todos. Tengo este Barchart con dos tipos de datos. Desembarcado y Muestreado. Como los datos de Desembarcado son mucho mayores que los de Muestreado, estos ultimos datos apenas si se aprecian en la grafica. Lo que se me ocurre es cambiar la escala de este eje. ¿alguien sabria como hacerlo? O alguna otra sugerencia claro. Muchisimas gracias a todos los que me estan ayudando Jose Luis

Dispongo de un experimento en el que cinco tratamientos ha sido aplicados a cinco grupos de voluntarios. En cada grupo había tres personas y a cada persona se le tomaron 3 medidas, En total dispongo de 45 medidas, pero evidentemente no son independientes entre sí. Si no tomo en cuenta que las medidas de la misma persona son más parecidas entre sí (bloques anidados) estaría incurriendo en

Dear R-gurus, I have an interpretation problem regarding lme models. I am currently working on dog locomotion, particularly on some variation factors. I try to figure out which limb out of 2 generated more dispersed data. I record a value called Peak, around 20 times for each limb with a record. I repeat the records during a single day, and on several days. I tried to build two models, one

I have been using Mandriva 2008.0 regularly since last Jan so I'm not very good. I have Tomb Raider Anniversary (TRA )and although i had to remove all the graphic options everytime it loads i.e. vsync shadow reflection and water effects, i could play it and i was happy. After some time i had to stop because i had a problem with video codecs and i had to install/remove a lot of things to get

I have a process where I am creating a effects plot similar to the cowles effect example. I would like to add the point estimates to the effects plot, can someone show me the correct syntax. I have included the "R" effects example, so you can show me the correct syntax. Thanks mod.cowles <- glm(volunteer ~ sex + neuroticism*extraversion, data=Cowles, family=binomial)

Dear all, I was wondering if anyone could help solve a problem of a missing interaction effect!! I carried out a 2 x 2 factorial experiment to see if eggs from 2 different locations (Origin = 1 or 2) had different hatching success under 2 different incubation schedules (Treat = 1 or 2). Six eggs were taken from 10 females (random = Female) at each location and split between the treatments,

Hi, I'm running R 2.11.1 on Windows 32-bit. I fitted a glm to a "large" data set with no problems (i.e., nothing is wrong with the created glm object). However, when I try to use the "effects" package to plot the effects I get the error message shown below. I get the same message if I try to plot a single effect. Any idea what could be wrong here? glm.eff <-