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:
library(pscl)
library(MASS)
library(AER)
library(VGAM)
library(truncreg)
library(censReg)
library(sampleSelection)
library(ggplot2)
library(boot)
library(aod)
library(lmtest)
library(zoo)
library(nlme)
library(lmtest)
library(boot)
library(spatcounts)
Mi duda es: ¿Que comando podrÃa utilizar para calcular por valores
predichos por el Binomial negativo para los 14 primeros valores de la
variable dependiente?
En el caso del ZIM utilizo: pred<-round(colSums(predict(zeroinfl,
type="prob") [,1:14]))
Gracias por vuestra ayuda! ;)
Lore
El 6 de noviembre de 2014, 12:30, "Marcuzzi, Javier Rubén"
<javier.ruben.marcuzzi en gmail.com
<mailto:javier.ruben.marcuzzi en gmail.com>> escribió:
Estimada Lorena Tudela
Voy a escribir algo por como yo lo aprendà (quizás mal redactado
para un estad?itico), lo observado es lo que se mide a campo y lo
predicho es lo que surge a partir del modelo estadÃstico. Luego
aprendà como realizarlo con R, pero no siempre es posible en forma
fácil, porque hay librerÃas que tienen dentro de sus algoritmos la
forma adecuada para la predicción, en cambio otras no, y hay que
realizarlo a mano.
Por lo cuál es necesario conocer que herramienta dentro de R está
utilizando para sus modelos. Porque, puede ser que ella misma
tenga la solución o se deba escribir en código R, y no siempre se
pueden combinar una librerÃa con la otra.
Javier Marcuzzi
El 06/11/14 a las 06:48, Lorena Tudela Marco escibió:
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]))
#Valores observados realmente
obs<-table(IB$nijt)[1:14]
#Tabla comparativa
rbind( pred, obs)
Y obtengo la siguiente tabla:
0 1 2 3 4 5 6 7 8 9 10 11 12 13
pred 3600 589 349 224 151 105 75 55 42 32 25 20 16 13
obs 3529 743 300 203 135 81 76 44 33 37 30 12 14 13
La duda me surge al intentarlo con el Modelo Binomial
negativo.¿Sabeis que
comando podrÃa introducir para obtener los 13 valores
predichos por el
modelo BN?
Muchas gracias por vuestra ayuda y buen dÃa.
Lorena
[[alternative HTML version deleted]]