<div dir="ltr">Buenos días, <div><br></div><div>Gracias Carlos, siguiendo el ejemplo que comentas, esto es lo que he introducido en el Scrip de RStudio:</div><div><br></div><div><div><i style="background-color:rgb(255,255,255)"><font color="#073763">library(xlsx)</font></i></div>
<div><i style="background-color:rgb(255,255,255)"><font color="#073763">library(xlsxjars)</font></i></div><div><i style="background-color:rgb(255,255,255)"><font color="#073763">library(rJava)</font></i></div><div><i style="background-color:rgb(255,255,255)"><font color="#073763">library(aod)</font></i></div>
<div><i style="background-color:rgb(255,255,255)"><font color="#073763"><br></font></i></div><div><i style="background-color:rgb(255,255,255)"><font color="#073763">R<-read.csv("2002.CSV", sep=";", dec=",", header=T)</font></i></div>
<div><i style="background-color:rgb(255,255,255)"><font color="#073763">attach(R)</font></i></div><div><span style="background-color:rgb(255,255,255)"><font color="#073763"><br></font></span></div><div><i style="background-color:rgb(255,255,255)"><font color="#073763">group<-gl(2,670,1340,labels= c("AVE", "Log.Imports.Value.in.1000.USD"))</font></i></div>
<div><i style="background-color:rgb(255,255,255)"><font color="#073763">weight<-c(AVE,Log.Imports.Value.in.1000.USD)</font></i></div><div><i style="background-color:rgb(255,255,255)"><font color="#073763">TUN<-lm(weight~group)</font></i></div>
<div><font color="#073763"><i><br></i></font></div><div><font color="#073763"><i>wald.test(b = coef(TUN), Sigma = vcov(TUN), Terms=2)</i></font></div><div><i><br></i></div><div>Obtengo esto en la consola de RStudio:</div>
</div><div><i><br></i></div><div><div><i><font color="#073763">Wald test:</font></i></div><div><i><font color="#073763">----------</font></i></div><div><i><font color="#073763">Chi-squared test:</font></i></div><div><i><font color="#073763">X2 = 6070.7, df = 1, P(> X2) = 0.0</font></i></div>
<div><span style="color:rgb(7,55,99)"><br></span></div><div><font color="#000000">Pvalue=0?? creo que estoy cometiendo algún error...</font>Adjunto al email la base de datos en excel, por si puede servir de ayuda! ;)</div>
<div><br></div><div><font color="#000000">¿Por otro lado no entiendo el termino de "terms"? (coloco "2" siguiendo el modelo, pero no se si es correcto)</font></div><div><font color="#000000"><br></font></div>
<div><font color="#000000">Gracias por vuestra ayuda de antemano,</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">Un abrazo, que tengáis buen día!</font></div><div><font color="#000000"><br>
</font></div><div><font color="#000000">Lorena</font></div><div><i><font color="#073763"><br></font></i></div><div><br></div><div><br></div><div><br></div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">El 10 de mayo de 2014, 22:28, Carlos Ortega <span dir="ltr"><<a href="mailto:cof@qualityexcellence.es" target="_blank">cof@qualityexcellence.es</a>></span> escribió:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div style="font-family:verdana,sans-serif">
Hola,<br><br></div><div style="font-family:verdana,sans-serif">Mira la ayuda del wald.test() para entender lo que te piden...<br>
</div><div style="font-family:verdana,sans-serif">Tienes que proporcionar los coeficientes del modelo, la matriz de covarianzas y los términos del modelo sobre los que quieres aplicar el test.<br><br>
</div><div style="font-family:verdana,sans-serif">Este es un ejemplo:<br></div><div style="font-family:verdana,sans-serif"><br>require(graphics)<br>## Annette Dobson (1990) "An Introduction to Generalized Linear Models".<br>
## Page 9: Plant Weight Data.<br>ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)<br>trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)<br>group <- gl(2, 10, 20, labels = c("Ctl","Trt"))<br>
weight <- c(ctl, trt)<br>lm.D9 <- lm(weight ~ group)<br><br>wald.test(b = coef(lm.D9), Sigma = vcov(lm.D9), Terms=2)<br><br></div><div style="font-family:verdana,sans-serif">Saludos,<br>Carlos Ortega<br>
<a href="http://www.qualityexcellence.es" target="_blank">www.qualityexcellence.es</a><br></div><div style="font-family:verdana,sans-serif"><br><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
El 10 de mayo de 2014, 19:43, Lorena Tudela Marco <span dir="ltr"><<a href="mailto:lorenatudelamarco@gmail.com" target="_blank">lorenatudelamarco@gmail.com</a>></span> escribió:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div class="">
Hola a todos y todas,<br>
<br>
Gracias por vuestro apoyo en cantidad de preguntas anteriores, de nuevo os<br>
escribo para compartir una duda:<br>
<br>
Estoy trabajando con un modelo bien sencillo, es una regresión simple, pero<br>
me gustaría comprobar la significación estadística de cada uno de los<br>
coeficientes de regresión en el modelo. La idea es hacer un contraste de<br>
hipótesis.<br>
<br>
Me he descargado el paquete "aod", para poder utilizar el test de Wald,<br></div>
pero mi duda es: *¿Que comandos tengo que utilizar para introducir el test<br>
de wald en el modelo?*<div class=""><br>
<br>
Esto es lo que estado probando pero la consola no muestra los estadísticos<br>
de wald, solo resuelve el modelo.<br>
<br></div>
*modeloTUN<-lm(as.numeric(AVE)~ Log.Imports.Value.in.1000.USD + Tariff + M<br>
)*<br>
*summary(modeloTUN, Wald=TRUE)*<div class=""><br>
<br>
Muchas gracias, seguimos avanzando!<br>
<br>
Un abrazo,<br>
<br>
Lorena<br>
<br></div>
[[alternative HTML version deleted]]<br>
<br>
<br>_______________________________________________<br>
R-help-es mailing list<br>
<a href="mailto:R-help-es@r-project.org" target="_blank">R-help-es@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-help-es" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-help-es</a><br>
<br></blockquote></div><span class=""><font color="#888888"><br><br clear="all"><br>-- <br><span style="font-family:verdana,sans-serif">Saludos,</span><br style="font-family:verdana,sans-serif">
<span style="font-family:verdana,sans-serif">Carlos Ortega</span><br style="font-family:verdana,sans-serif">
<span style="font-family:verdana,sans-serif"><a href="http://www.qualityexcellence.es" target="_blank">www.qualityexcellence.es</a></span>
</font></span></div>
</blockquote></div><br></div></div>