search for: ajustment

Dear All, I'd like to fit a "kind" of logistic model to small data-set using nonlinear least-squares regression. A transcript of R-script are reproduced below. Estimated B and T (the model's coeff, herein B=-8,50 and T=5,46) seem appropriate (at least visually) but are quite diff from those obtained w/ SPSS (Levenberg-Marquardt): B=-19,56 and T=2,37. Am I doing something wrong in

Oliver gracias por le pista! Es un poco mas complicado, el código tiene varios pasos incluidos calculo de nuevas variables con los betas del glm y otras vueltas. Necesitaría una sentencia que tome cada lote haga los cálculos por separado y luego los une en el dataframe de salida El 30 de mayo de 2016, 10:02, Olivier Nuñez <onunez en unex.es> escribió: > Mira lapply .... > > Si L=

Buenos días tengo una consulta sobre un código rutinario que quiero ejecutar. Tengo 120 Lotes con un ID diferente (del C001 al C120) a cada uno de los cuales debo efectuares un GLM y con los parámetros obtenidos calcular otras variables para luego hacer gráficas, razón por la cual no puedo procesarlos todos juntos. Intente hacer una sentencia repeat() o con for() para que me ejecute de forma

Hi, I make a glm analysis with binomial errors and its OK. Now I need plot the ajusted logist line in my points. How do I make it? | * * * * -> I need plot an ajusted | * * line in this points. | * * | * | * * | * * |* * * * |_________________________ Exist an abline like function for ajust others formulas? Exist in R

Package: xen-utils-3.2-1 Version: 3.2.0-5 Severity: important virtinst Debian package is pointing the pygrub path to /usr/bin/ and hvmloader to /usr/lib/xen/boot/. I suggest ajust the path creating a symbolic link of /usr/lib/xen-3.2-1/bin/pygrub to /usr/bin/pygrub and /usr/lib/xen-3.2-1/boot/hvmloader to /usr/lib/xen/boot/hvmloader. In the last case, need ajust the directory

Carlos, capaz no me explique bien, va de nuevo... Tengo un dataFrame con varias variables que describen Lotes (127), he creado un código que calcula nuevas variables, re codifica otras y ordenas otras tantas. Es muy laborioso correr ese código 127 veces por lo cual quería hacer un sentencia (del tipo if() o repeat()), para obtener al final todos los lotes procesados. Saludos y gracias por tan

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150127/8c94355b/attachment.html>

buenas   Por favor me pueden decir como ajusto unos datos con las distribuciones poisson generalizada poisso-pascalgeneralizada conway-maxwell-poisson   no se en que paquette estan   Gracias [[alternative HTML version deleted]]

Hola: tengo un pequeño problemita que es más bien estadístico que específico de R, pero de todos modos espero me puedan dar una mano. Quiero encontrar una forma de reportar la Bondad de ajuste de un modelo No Lineal. Estuve buscando y leyendo y parece que tal cosa no es posible o al menos no es fiable, desde el punto de vista estadístico. El tema es que quiero mandar un trabajo a una

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150128/9e23caec/attachment.html>

Buenos días para cada una/o. Acabo de suscribirme. Tengo una inquietud que vi tratada por partes (y usando distintos paquetes) en las consultas anteriores. ¿Cuál es la mejor forma de ajustar distribuciones datos a distribuciones de probabilidad teóricas? ¿Cuáles son los paquetes para hacerlo en forma ágil? Gracias por su atención. Saludos. César Escalante C.

Buenas tardes, Quisiera obtener el intervalo de confianza (y también intervalos de predicción) para los valores predichos en un modelo nls. ¿Hay alguna manera que no sea por ggplot2 (me interesaría obtener el valor listado -además de en el gráfico-) o por bootstrap? Os copio el código del ajuste del modelo y predicción para los 3 días siguientes: *#Ajuste del modelo* model = nls(formula =

Muchas gracias por la respuesta. He mirado y los coeficientes no son altos pero sí tengo una gran cantidad de ceros en la variable dependiente (más del 90%). Sin embargo, al incluir otro tipo de variables independientes no me da ese aviso, dejando la misma variable dependiente. ¿Cómo podría utilizar stan/rstan de forma sencilla para diagnosticar el modelo? Muchas gracias El Lun, 9 de Abril de

Buenas a todos. Estoy realizando unos ajustes a una serie de valores mensuales maximos con ismev. Los resultados son muy buenos. Por definicion, los valores que obtengo del ajuste de la funcion de distribucion GEV, me dan valores de periodos de retorno medios para los niveles de retorno que estudio. El problema es que estos valores son los medios esperables para esos periodos de retorno y yo, lo

Hola a todos estoy aplicando la función polinómica de Hossfeld [1], y algunos otros modelos no lineales para tratar de ajustarlos a un grupo de datos forestales,   [1] Y= b*t*exp(c)/(t*exp(c)+a) Al colocar la función en R con parámetros estimados, me devuelve los siguiente: ## model1 <- nls(ho ~ (b*edad*exp(c)/(edad*exp(c)+a)), data=nigra,     start=list(a=0.005,b=0.08,c=-0.00006),

Hola, con vistas a docencia, me gustaría saber si hay alguna librería que haga el test de la chi cuadrado para variables continuas, que no sean la normal. Lo encontré para variables discretas y para la normal. También, obviamente, se puede hacer un pequeño programa que lo calcule "como a mano", pero a los alumnos les obliga a entrar en detalles de programación... :-/ La idea sería

¿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

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, Estoy estimando un modelo binomial negativo de ceros inflados (ZINB) utilizando el comando zeroinfl() del paquete pscl. Al ejecutarlo me da el siguiente aviso: Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred ¿Sabéis que significa y si puedo usar el modelo aún con ese aviso? ¿Los coeficientes son fiables? Muchas gracias, Miriam

Colegas: Tengo una pregunta rara, y ofrezco mis disculpas si es una tonteria. Estoy tratando de hacer unos graficos con datos recolectados cada 24 horas, con varias replicas por cada periodo de tiempo. La mejor representacion que he encontrado es hacer boxplots haciendo el tiempo una variable categorica. Sin embargo, me gustaria tener una linea de regresion, o ajuste usando splines o LOESS a