similar to: Duda glmer

Buen dia!! Pues me encuentro trabajando con un conjunto de datos simulados que se ajustan a un modelo Ar(1) y pues queria saber si existe un comando para realizarle una prueba de homocedasticidad, pues la prueba que hay en el R es la bartlett.test pero pues no estoy muy seguro para usarla. Gracias [[alternative HTML version deleted]]

Hi, I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and \beta_1, this can be achieved by solving the following three equations: n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) - \sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} 1/(psihat(i)) - (\alpha+1) \sum\limits_{i=1}^{n} ( 1 / (x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} (

Hi Listers, seems I have forgotten some basics re offset in glm: data: counts (y) from locations off different size (area), explanatory variable: x Model: y ~ x+offset(area) Predictions (pred) using Poisson errors plot(x,y) and points(x,pred) gives neat "line" of estimated values. However, for ease of understanding graphs are better using plot(x,y/area). Question: How to display

Hola, ¿qué tal? Cierto, cierto, había un error en el código que publiqué. Pero el diagnóstico es parecido. Cuando los datos se generan con el coeficiente de x2 igual a 7, los coeficientes estimados tienen una distribución extraña, bimodal (aparentemente), en lugar de _normalmente_ distribuida alrededor del 7 como se espera. Supongo que depende del número de casos en que x2 = 1 e y = 0. Un

In the logistic regression model, there is no residual log (pi/(1-pi)) = beta_0 + beta_1*X_1 + ..... But glm model will return residuals What is that? How to understand this? Can we put some residual in the logistic regression model by replacing pi with pi' (the estimated pi)? log (pi'/(1-pi')) = beta_0 + beta_1*X_1 + .....+ ei Thanks! [[alternative HTML version deleted]]

Hi all Assume I have a data set xx; Group: 1=group1 ?, 2=group2 IQ: ?1= High, 0 =low fit <- glm(IQ ~group, data = xx, family = binomial()) summary(fit) Results ?????? ????????????Estimate Std. Error z value Pr(>|z|) (Intercept) -2.55456??? 0.210 -12.273? < 5e-16 *** group????????? 0.36180 ?????0.076?? 3.952 ????5.24e-05 *** the odd ratio = exp(0.36180 )= 1.435912 My question

Hi all I know about producing a minimal example to show my problem. But I'm having trouble producing a minimal example that displays this behaviour, so please bear with me to begin with. Observe: I create an array called model.mat. Some details on this: > str(model.mat) num [1:18, 1:4] -0.170 -0.304 -2.617 2.025 -1.610 ... - attr(*, "dimnames")=List of 2 ..$ : chr

Estimados Les consulto por redes neuronales, hay diversos artículos como los siguientes (el último tienen un error actualmente). Pero mi pregunta va un poco por otro lado. http://www.r-bloggers.com/build-your-own-neural-network-classifier-in-r/ http://www.r-bloggers.com/classification-using-neural-net-in-r/ Básicamente se puede calcular un valor, por ejemplo doblar 2,4 grados a la derecha, luego 1

Hola, Esta es una forma de hacerlo... Mira que lo primero que he modificado es el fichero "x.csv" para sustituir los espacios en los nombres por "_". Y también he quitado los acentos y las eñes... He utilizado el paquete RNNS y la función "mlp()" para ajustar la red. #------------------------------------------- > x <- read.csv("x.csv",

Por favor, estoy aprendiendo R para aplicarlo exclusivamente en modelos jerárquicos lineales o multinivel. Todo lo que me puedan ayudar les quedaré muy agradecido. Cordialmente, Jairo P.D.: Presento mis disculpas por el correo anterior [[alternative HTML version deleted]]

Hi R users! I am trying to learn some multilevel analysis, but unfortunately i am now very confused. The reason: http://www.ats.ucla.edu/stat/hlm/seminars/hlm_mlm/mlm_hlm_seminar.htm http://www.ats.ucla.edu/stat/sas/seminars/sas_mlm/mlm_sas_seminar.htm and MlmSoftRev. pdf from mlmRev package. >From what i see, the first two links seem to declare the level one variable as a random part (i

Hi Dr. Snow, I am a graduate student working on analyzing data for my thesis and came across your post on an R forum: The default link function for the glm poisson family is a log link, which means that it is fitting the model: log(mu) ~ b0 + b1 * x But the data that you generate is based on a linear link. Therefore your glm analysis does not match with how the data was generated

Hola, ¿qué tal? Tengo una pregunta de esta

Dear List: My question is more statistical than R oriented (although it originates from my work with nlme). I know statistical questions are occasionally posted, so I hope my question is relevant to the list as I cannot turn up a solution anywhere else. I will frame it in the context of an R related issue. To illustrate the problem, consider student achievement test score data with multiple

Gracias Jorge. No entiendo bien; la variable objetivo es ya factor. El árbol me la predice bien, como factor, también. Es al ir construyendo el vector que lo anota con un nº, según de cuál de las 4 categorías se trate. Quoting Jorge I Velez <jorgeivanvelez en gmail.com>: > Estimado Manuel, > > Debes definir ecsta como factor usando, por ejemplo, > > factor(ecsta,

Buenas, Estoy intentando hacer un trabajo con una compañera, a la que le pedí que me lo explicara en unas lineas. Esas lineas son las que siguen: Quisieramos analizar algunos comportamientos humanos a lo largo del tiempo; esto es, encontrar patrones de en algunos comportamientos del habla y la comunicación: vocalizaciones, gestos (divididos en categorías), mirada (en categorías). Conocemos un

Buenas tardes, tengo un loop que hace un árbol de clasificación cada vez y va creando un vector con una predicción que hace. Son 4 categorías (pongamos a, b, c y d), pero en vez de ir añadiendo la categoría predicha me añade al vector el nº (del 1 al 4) al que corresponde esa categoría. Supongo que se puede hacer que añada la categoría (la letra), pero no sé cómo. preds <- {} for (i

--- El lun, 12/9/11, Cristalina <pa100cia77@yahoo.es> escribió: De: Cristalina <pa100cia77@yahoo.es> Asunto: Re: [R-es] Cosinor Analysis Para: "Carlos Ortega" <coforfe@gmail.com> Fecha: lunes, 12 de septiembre, 2011 08:43 Hola,   Carlos, muchas gracias.   El método empleado en http://tolstoy.newcastle.edu.au/R/e6/help/09/01/0626.html (el url que se referencia

Bueno, creo que no contesté tu pregunta. Con training <- data[-i, ] crea una df llamada training, sin la muestra i, que después utiliza para entrenar el algoritmo. Quoting Javier Marcuzzi <javier.ruben.marcuzzi en gmail.com>: > Estimado Manuel Mendoza > > Con sus datos y a modo de curiosidad, ¿que pasa en training <- data[-i, ]? > > Javier Rubén Marcuzzi >

hola a todos y todas, Estoy haciendo un estudio sobre las mareas y tengo dos variables que quiero relacionar, la altura del mar y el tipo de marea viva, muerta o intermedia. Unos datos simulados podrían ser: > datos <- data.frame(v1=sin(1:50), v2= rep(c("a","b","c"), each = 5, > len = 50)) Ahora mi pregunta, si dibujo la altura de la marea sería: