similar to: Modelo para datos de conteos (Manuel Spínola)

Estimados miembros de la lista. Estoy usando 2 funciones del paquete vcd y obtengo resultados que parecen discrepar entre ambas funciones. > dput(peces) c(26L, 63L, 1L, 21L, 20L, 50L, 0L, 104L, 19L, 0L, 8L, 32L, 0L, 24L, 5L, 2L, 0L, 36L, 0L, 5L, 27L, 1L, 13L, 24L, 15L, 25L, 23L, 5L, 37L, 26L, 29L, 14L, 15L) La función Ord_plot me dice que la mejor distribución para los datos es binomial

Hi. I need to draw a circle whit center (a,b) and radio r. So I use the R code below a<-1.975 # valore x del centro b<-1.215 # valores y del centro r<-1.46 # radio x1<-seq(a-r,a+r,by=0.01); #los valores de x yp<-sqrt(r^2-(x1-a)^2)+b; # los valores y a partir de la ra??z positiva yn<-(-1)*sqrt(r^2-(x1-a)^2)+b; # los valores y a partir de la ra??z negativa x<-c(x1,x1);

Hello, The following piece of code works fine in R.2.10 (ubuntu): switch(distr,         normal    = {if (is.infinite(param["desv"]))                      n <- c(n,"La desv. estándar no puede ser Inf.")                    if (param["desv"]<0)                      n <- c(n,"La desv. estándar no puede ser <0.")                    },        

Buenas tardes, ¿Cómo interpretarías el intercepto que da R en un modelo de ceros inflados? Por un lado en la parte de conteo tengo un intercepto de -4.2 y en la parte de ceros de 102, ambos salen significativos (***). ¿Qué me dirían? Gracias

El problema del módulo es que pierde el signo. En tu caso sale igual porque has invertido el signo del coeficiente en el polinomio (en realidad se me pasó a a mí advertir que el término independiente debe ir con signo negativo): .> polyroot(z=c(0.5,0,0,0,0,1)) [1] 0.7042902+0.5116968i -0.2690149+0.8279428i -0.2690149-0.8279428i [4] 0.7042902-0.5116968i -0.8705506+0.0000000i .> .>

Gracias Emilio. Si, son asimétricos porque se mueven en rangos distintos. De hecho, miden cosas distintas, pero una muestra no puede ser positiva y negativa a la vez. Me planteo multiplicar las positivas por 2 y al editar el mapa corregirlo para que indique los valores reales. Eso funcionaría, pero son muchos mapas y quizás haya una función en ggplot que evite el problema directamente. Manuel El

Hola de nuevo Carlos, he probado a quitar esa variable categórica y me sigue dando el aviso... El Lun, 9 de Abril de 2018, 20:17, Carlos J. Gil Bellosta escribió: > Si, creo que el motivo del warning puede ser ese. Es hipotético, pero > plausible. Sobre todo cuando tienes más de un 90% de ceros. > > El coeficiente de ese nivel para el modelo de la mixtura (ceros vs > binomial >

No sé muy bien, Carlos, como hacerlo, pero no parece más sencillo que lo de multiplicar por 2 los positivos y después editar el mapa. De lo que decías, Emilio, de cambiar el color del punto central creo que haría más confusa la interpretación. El vie., 20 sept. 2019 a las 10:33, Carlos Ortega (<cof en qualityexcellence.es>) escribió: > Hola, > > Una alternativa que suele

Estimados, Les consulto por el siguiente problema. Variable respuesta: Número de picadas de alimentación de los insectos (picadas) Factores: 1) Especie de insecto (especie), 2) Estado de desarrollo de los insectos (edad) y 3) Estado de desarrollo de frutos (estado). Modelo: MLG con distribución binomial negativa como componente aleatorio. Estrategia de selección de modelo: Evaluación del cambio

Manuel: I changed your variable names from x to 'long' and y to 'lat' on the riqueza_out.csv file. The code below should do what you want. Also, since the legend title is kind of long, I broke it down into three lines so you can see more plot area. I am cc'ing the other groups so more people use it if needed. library(rgdal) library(ggplot2) library(sp) library(maptools)

Hola, LLevo buscando desde hace tiempo como hacer el Goodness of fit test en R. Es decir, me explico, intento hacer una cosa parecida que se hace en Minitab, por ejemplo, yo tengo un conjunto de datos, y lo que quiero es sabes que tipo de distibución es, en minitab se hace un histograma para ver si se ajusta bien o no a la campana de Gauss, luego vemos si aproximar la distribución de la muestra

¿Quieres decir que para un nivel de una variable categorica todas las observaciones de la variable respuesta sean ceros? Gracias El Lun, 9 de Abril de 2018, 19:59, Carlos J. Gil Bellosta escribió: > ¿Podría ser que para algún nivel de alguna variable independiente > categórica solo hubiese ceros? En ese caso, casi seguro, aparecería ese > tipo de warning. > > El lun., 9 abr. 2018 a

(Perdón, no puse asunto en el post original) Ah, corremos R 3.0.2 y l alibrería MASS 7.3-29 Estimados todos: Os escribo porque tengo un problema que nos está dejando un poco trastornadillos. Probablemente sea uno de esos que los ve un experto a un kilómetro y dice ¡oh! toca aquí. Pero nos trae de cabeza. El caso es que aquí usamos un script para calcular ciertos parámetros relacionados con el

Hola, Estoy intentando esta convolución con una discretizada Gamma con Binomial negativa pero me da el siguiente error. Alguien sabe como puedo hacer funcional este algoritmo con estos parámetros? > fx<-discretize(pgamma(x,shape=2.4149,scale=5742.2), + method="rounding",from=0, to =100, step=0.5) > Fs <- aggregateDist("recursive", model.freq = "negative

I'd like to announce the availability of the "GWRM" package. It provides a set of functions for fitting Generalized Waring Regression Models. It also includes the dataset and the example of Rodriguez-Avi, J; Conde-Sanchez, A; Saez-Castillo, A.J., Olmo-Jimenez, M. J. and Martinez Rodriguez, A. M.(2009). A generalized Waring regression model for count data. Computational Statistics and

I'd like to announce the availability of the "GWRM" package. It provides a set of functions for fitting Generalized Waring Regression Models. It also includes the dataset and the example of Rodriguez-Avi, J; Conde-Sanchez, A; Saez-Castillo, A.J., Olmo-Jimenez, M. J. and Martinez Rodriguez, A. M.(2009). A generalized Waring regression model for count data. Computational Statistics and

What do you use to administer the PBX accounts? (Please don't tell me LDIF text files). Any luck with Luma? (http://luma.sourceforge.net) __________________________________________________________ Mensagem enviada usando o Webmail da ViaLink ver. 2.7.8

Je vous invite ? faire partie de mon r?seau professionnel en ligne sur le site LinkedIn. Manuel Dahmen Manuel Dahmen Developer chez Ibiiz Tera Veuillez confirmer que vous connaissez Manuel Dahmen https://www.linkedin.com/e/nxxpbc-gmf1uzfh-5t/isd/2681722075/t3f3BniE/EML-invite_guest_snackified_59/ ------------------------------------------------------------------------------- Apr?s vous ?tre

Tengo una consulta, alguien tiene algun script en R para hacer graficos de anomalias estandarizadas (precipitacion, temperatura, etc) y que me permita pintar de rojo anomalias positivas y azules las negativas, ​Muchas gracias de antemano​ -- Luis Alfonso LOPEZ ALVAREZ MSc- Meteorologia - Ing. Sistemas PROFESIONAL ESPECIALIZADO (E) METEOROLOGO - IDEAM e-mail: llopez@ideam.gov.co

Buenas, Tengo unas rectas en las que cada x número de ellas hay algún error en la captura de los datos y se produce un outlier. Os pongo un ejemplo: x<-rnorm(15,12,2) y<-jitter(x,200) y[7]<-20 plot(x,y) En este ejemplo se ve muy claro que el outlier y quitarlo seria facil, pero... ¿como hacerlo de manera automatica? es que hay mil curvas com oesta y puede ser una locura. Gracias