similar to: lme4 y residuales

Buenos días a todxs, Estoy comparando la predicción de los valores (0, 1, 2, 3,.....hasta 13) frente a los observados. Con la idea de comparar el modelo Zero inflated y el Binomial negativo y ver cual presenta mas distancia frente a las predicciones observadas. Para ello introduzco los códigos en la consola: #Modelo ZIM pred<-round(colSums(predict(zeroinfl, type="prob") [,1:14]))

Hola Javier, Si, cuando hablo de valor observado me refiero al valor real en campo y el predicho al que estiman los modelos. Disculpa, que no lo detallase así desde el principio. En mi caso trabajo con dos diferentes: Zero inflated y Binomial Negativo y me gustaría comprobar que diferencia (distancia) existe entre cada uno de ellos y la realidad. Estoy trabajando con los siguientes paquetes:

Hola. Tengo que representar en un mapa, valores de proporciones de extranjeros con tarjetas de residencia en vigor, según nacionalidad y país de procedencia. Los datos son de este estilo Rumania Andalucía 0,200365823 Aragón 0,377192073 Asturias 0,205353012 I. Baleares 0,06592893 Canarias 0,03050312 Cantabria 0,190397393 Castilla-León 0,211670421 Castilla-La Mancha

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 =

I'm trying to create a graph where different cells of a grid (a shapefile) will be painted with a color share scale, where the most easy way is to use gray(). Can I somehow get a vector (gradient) of colors, a vector of colors with other methods but gray()? I'm doing this until now quad_N_sp <-

R-Amigos Mis mejores deseos en estas fiestas y que el año nuevo premie la generosidad con la que los miembros de este grupo comparten su sabiduría y sus trucos. Felicidades -- Sent from my IPAD [[alternative HTML version deleted]]

Hola buenos días hice un gráfico combinado de líneas, puntos y barras en ggplot2, pero no sé cómo puedo poner la leyenda de eso gráfico para que me represente para las líneas con puntos los valores estimados por un modelo y observados. este es mi código: ggplot(MLM,aes(x=Individuo)) + geom_bar(aes(y=Observada), stat = "identity", color= "gray") + geom_line(aes(y=Estimada),

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

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

Buenas tardes a la comunidad R Quiero hacer un consulta sobre una duda que tengo en R. Específicamente es saber como puede en R eliminar de un modelo ARMA los coeficientes no significativos, por ejemplo en eviews yo puedo especificar ar(1) ar(4) ar(5) ar(8) los cuales para mi caso son los significativos. Sin embargo en R debo especificar c( 8,0,0) usando la función Arima del paquete forecast,

Buenos días; Aunque ya llevo un tiempo con LaTeX y R no he conseguido que me funcione la librería tikzdevice. El primer problema está con tikzDictionary que me da como error: cannot create lock file. He estado buscando información en Google pero sin éxito. No he conseguido solucionarlo indicando un archivo en RProfile con options(tikzMetricsDictionary = ''path/archivo''). Cuando

Hola a todos, en primer lugar quería agradecer la ayuda recibida desde el foro con respecto a la creación de una red neuronal. Estoy utilizando el paquete Neuralnet, que me parece que es bastante bueno, pero tengo el problema que soy incapaz de hacer las predicciones del modelo. Sé que se hace con el comando "compute", pero no encuentro ningún ejemplo de cómo hacerlo. ¿Alguien me puede

Dear all, I have been using Speex during last year over PC with very good results. Now I'm trying to move my code to the DSP TI C6713. It's the first time I work with a DSP and with Code Composer. I use speex quite well but I'm not very good at compiling for different platforms. After reading all manuals I tried to compile the speex library with Code Composer and the speex options for

OK, para arboles no, pero para otro tipo de modelos como un modelo lineal si no? O para uqe tipo de modelos tengo que normalizar los datos??? ________________________________ De: Manuel Spínola <mspinola10 en gmail.com> Enviado: viernes, 17 de febrero de 2017 19:06 Para: Jesús Para Fernández Cc: r-help-es en r-project.org Asunto: Re: [R-es] Modelos con datos normalizaddos Hola Jesús, No

Hola Buenas noches   Alguien me podria indicar en que paquete puedo encontras modelos lineales y cuadratico para ajustar una serie de datos   Gracias [[alternative HTML version deleted]]

Saludo estimados compañeros y compañeras Tengo una matriz de datos de 51 filas por 160 columnas, Cada columna es una serie de tiempo. Existe alguna libreria q me calcule la posicion de la fila donde se hallan estos maximos y minimos o los eventos de estas curvas Agradezco la atencion CARLOS ANDRES [[alternative HTML version deleted]]

Hola a todos, Tengo algunas dudas sobre como introducir en un modelo lineal factores con efectos fijos y aleatorios. Mi diseño es el siguiente: Factores Line: 40 líneas de trigo (Fijo) Rep: 3 Bloques (Aleatorio) Year: 2 Años (Aleatorio) Variable dependiente alpha.ug.mg Nota: Adjunto tabla de datos Dentro del diseño, 'Rep' y 'Year' se considera de efectos aleatorios y

Hola amigos, Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en el ejemplo de la ayuda: DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) summary(fm1DNase1) Aquí se está modelando la densidad óptica de un ensayo relacionada de forma no lineal (logística) con (el logaritmo) de la concentración de una proteína.

hola soy nuevo usuario de R y necesito crear un objeto p tal que, p=lambda*integral[f´´(x)^2 dx ] donde "lambda" es uno de los parametros que resultan de la funcion "smooth.spline()" y la integral es sobre la derivada 2 de esa misma funcion... dos cuestiones: 1) como extraigo lambda de los resultados de smooth.spline() para usarlo como objeto cuando lo requiera y 2) como

No es un problema de tm ni de SnowfallC ni de mcapply (por el path utilizas linux, en windows mcapply según el manual no va bien) No defines bien los objetos que pasas. Pasas doc.corpus en lugar de corpus ( o asignas a corpus en lugar de a doc.corpus) . Depura los programas cuando salga un error de objeto, como te pone en el Error que pasas . Temporalmente lo tienes arreglado en