Hola. La verdad es que yo suelo usar el Likelihood ratio Test, mediante la función Anova del paquete car, pero no está de más saber que existe, pero por si alguna vez me piden este ya sé dónde mirar.
Saludos El 12/05/14 18:21, Lorena Tudela Marco escribió: > Gracias Carlos, muy interesante el enlace, le echo un vistazo! ;) > > > > > > > El 12 de mayo de 2014, 14:04, Carlos Ortega > <[email protected]>escribió: > >> Hola Lorena, >> >> Mira el ejemplo que aparece en esta página (aplicado sobre la regresión >> logÃstica) sobre cómo se aplica el test de Wald y el significado de >> "Terms": >> >> http://www.ats.ucla.edu/stat/r/dae/logit.htm >> >> Saludos, >> Carlos Ortega >> www.qualityexcellence.es >> >> >> El 12 de mayo de 2014, 12:57, Lorena Tudela Marco < >> [email protected]> escribió: >> >> Buenos dÃas, >>> Gracias Carlos, siguiendo el ejemplo que comentas, esto es lo que he >>> introducido en el Scrip de RStudio: >>> >>> *library(xlsx)* >>> *library(xlsxjars)* >>> *library(rJava)* >>> *library(aod)* >>> >>> *R<-read.csv("2002.CSV", sep=";", dec=",", header=T)* >>> *attach(R)* >>> >>> *group<-gl(2,670,1340,labels= c("AVE", "Log.Imports.Value.in.1000.USD"))* >>> *weight<-c(AVE,Log.Imports.Value.in.1000.USD)* >>> *TUN<-lm(weight~group)* >>> >>> *wald.test(b = coef(TUN), Sigma = vcov(TUN), Terms=2)* >>> >>> Obtengo esto en la consola de RStudio: >>> >>> *Wald test:* >>> *----------* >>> *Chi-squared test:* >>> *X2 = 6070.7, df = 1, P(> X2) = 0.0* >>> >>> Pvalue=0?? creo que estoy cometiendo algún error...Adjunto al email la >>> base de datos en excel, por si puede servir de ayuda! ;) >>> >>> ¿Por otro lado no entiendo el termino de "terms"? (coloco "2" siguiendo >>> el modelo, pero no se si es correcto) >>> >>> Gracias por vuestra ayuda de antemano, >>> >>> Un abrazo, que tengáis buen dÃa! >>> >>> Lorena >>> >>> >>> >>> >>> >>> >>> El 10 de mayo de 2014, 22:28, Carlos Ortega >>> <[email protected]>escribió: >>> >>> Hola, >>>> Mira la ayuda del wald.test() para entender lo que te piden... >>>> Tienes que proporcionar los coeficientes del modelo, la matriz de >>>> covarianzas y los términos del modelo sobre los que quieres aplicar el >>>> test. >>>> >>>> Este es un ejemplo: >>>> >>>> require(graphics) >>>> ## Annette Dobson (1990) "An Introduction to Generalized Linear Models". >>>> ## Page 9: Plant Weight Data. >>>> ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) >>>> trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) >>>> group <- gl(2, 10, 20, labels = c("Ctl","Trt")) >>>> weight <- c(ctl, trt) >>>> lm.D9 <- lm(weight ~ group) >>>> >>>> wald.test(b = coef(lm.D9), Sigma = vcov(lm.D9), Terms=2) >>>> >>>> Saludos, >>>> Carlos Ortega >>>> www.qualityexcellence.es >>>> >>>> >>>> >>>> >>>> El 10 de mayo de 2014, 19:43, Lorena Tudela Marco < >>>> [email protected]> escribió: >>>> >>>>> Hola a todos y todas, >>>>> >>>>> Gracias por vuestro apoyo en cantidad de preguntas anteriores, de nuevo >>>>> os >>>>> escribo para compartir una duda: >>>>> >>>>> Estoy trabajando con un modelo bien sencillo, es una regresión simple, >>>>> pero >>>>> me gustarÃa comprobar la significación estadÃstica de cada uno de los >>>>> coeficientes de regresión en el modelo. La idea es hacer un contraste de >>>>> hipótesis. >>>>> >>>>> Me he descargado el paquete "aod", para poder utilizar el test de Wald, >>>>> pero mi duda es: *¿Que comandos tengo que utilizar para introducir el >>>>> test >>>>> de wald en el modelo?* >>>>> >>>>> >>>>> Esto es lo que estado probando pero la consola no muestra los >>>>> estadÃsticos >>>>> de wald, solo resuelve el modelo. >>>>> >>>>> *modeloTUN<-lm(as.numeric(AVE)~ Log.Imports.Value.in.1000.USD + Tariff >>>>> + M >>>>> )* >>>>> *summary(modeloTUN, Wald=TRUE)* >>>>> >>>>> >>>>> Muchas gracias, seguimos avanzando! >>>>> >>>>> Un abrazo, >>>>> >>>>> Lorena >>>>> >>>>> [[alternative HTML version deleted]] >>>>> >>>>> >>>>> _______________________________________________ >>>>> R-help-es mailing list >>>>> [email protected] >>>>> https://stat.ethz.ch/mailman/listinfo/r-help-es >>>>> >>>>> >>>> >>>> -- >>>> Saludos, >>>> Carlos Ortega >>>> www.qualityexcellence.es >>>> >>> >> >> -- >> Saludos, >> Carlos Ortega >> www.qualityexcellence.es >> > [[alternative HTML version deleted]] > > > > _______________________________________________ > R-help-es mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es [[alternative HTML version deleted]]
_______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
