Buenos días, *Muchas gracias, todas las aportaciones han sido bien útiles.* Las he tenido en cuenta y he pasado los datos con el R, siguiendo el siguiente comando: *modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8 + d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20 +d21 + Tariff*d1 + Tariff*d2 + Tariff*d10)* *summary(modeloRTUN2)* Siendo: AVE. = Variable dependiente (explicada) Tariff = Variable dependiente (explicativa) d1.....d22= los diferentes productos (dimensión producto: para comprobar si el tipo producto tiene efecto sobre AVE.) Tariff * d1 = efecto fijo del producto y tariff En resumen me gustaría comprobar si existe trade off entre las dos variables NTM y Tariff y conocer que productos son significativos (conocer que productos tienen presentan mas protección en frontera). Me surge ahora otra duda ...que diferencia hay en colocar "Tariff*d2" y "Tariff:d2", en ambos casos es la multiplicación de las dos variables? Adjunto la base de datos en formato excel al email, para que sea mas fácil entenderme. Un abrazo mas y* muchas gracias, que tengáis un buen día.* Realmente todos los consejos han sido aportaciones muy valiosas para que pudiese seguir avanzando, os lo agradezco! Seguimos hablando ;) Lorena El 13 de enero de 2014, 9:46, José Trujillo <trujillo en unex.es> escribió:> Si AVE es la variable dependiente y Tarrif "LA" variable independiente: > > ¿Quién son d1 y d1T? > > Como dice Carlos, as.factor dentro de la fórmula no es buena idea, es > más seguro que la variable d1T sea convertida en factor previamente. > > Saludos. > > > El 12/01/14 13:53, Lorena Tudela Marco escribió: > > Buenos días, > > > > Me gustaría aplicar una regresión múltiple a los datos con los que > trabajo pero > > no se como introducir los datos en R. He probado introducir el siguiente > > comando: > > > > modeloM<-lm(AVE.~ d1 + as.factor(d1T)*Tariff) > > summary(modeloM) > > > > Pero me da el siguiente error: > > > >> summary(modeloM) > > Error in if (attr(z$terms, "intercept")) sum((f - mean(f))^2) else > sum(f^2) > > : > > argument is of length zero > > > > Mi intención es calcular la siguiente regresión: Ave = cte + > > d1*coeficiente1 + (d1T*Tariff)*coeficiente 2 > > Siendo: AVE, la variable dependiente y Tarrif, la variable independiente. > > > > Imagino que no estoy introduciendo correctamente el comando inicial del > > modelo y por eso me da error ¿Alguien esta trabajando en temas similares > o > > podría orientarme ? > > > > Muchas gracias, cualquier aportación sera bienvenida! > > > > Un abrazo, que tengáis un buen día. > > > > Lorena > > > > [[alternative HTML version deleted]] > > > > > > > > _______________________________________________ > > R-help-es mailing list > > R-help-es en r-project.org > > https://stat.ethz.ch/mailman/listinfo/r-help-es > > > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20140114/45f0658f/attachment-0001.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: TUNISIA.docx Type: application/vnd.openxmlformats-officedocument.wordprocessingml.document Size: 22094 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20140114/45f0658f/attachment-0001.bin>
Buenos días, *Muchas gracias, todas las aportaciones han sido bien útiles.* Las he tenido en cuenta y he pasado los datos con el R, siguiendo el siguiente comando: *modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8 + d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20 +d21 + Tariff*d1 + Tariff*d2 + Tariff*d10)* *summary(modeloRTUN2)* Siendo: AVE. = Variable dependiente (explicada) Tariff = Variable dependiente (explicativa) d1.....d22= los diferentes productos (dimensión producto: para comprobar si el tipo producto tiene efecto sobre AVE.) Tariff * d1 = efecto fijo del producto y tariff En resumen me gustaría comprobar si existe trade off entre las dos variables NTM y Tariff y conocer que productos son significativos (conocer que productos tienen presentan mas protección en frontera). Me surge ahora otra duda ...que diferencia hay en colocar "Tariff*d2" y "Tariff:d2", en ambos casos es la multiplicación de las dos variables? Adjunto la base de datos en formato excel al email, para que sea mas fácil entenderme. Un abrazo mas y* muchas gracias, que tengáis un buen día.* Realmente todos los consejos han sido aportaciones muy valiosas para que pudiese seguir avanzando, os lo agradezco! Seguimos hablando ;) Lorena El 14 de enero de 2014, 12:24, Lorena Tudela Marco < lorenatudelamarco en gmail.com> escribió:> Buenos días, > > *Muchas gracias, todas las aportaciones han sido bien útiles.* Las he > tenido en cuenta y he pasado los datos con el R, siguiendo el siguiente > comando: > > *modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8 + > d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20 +d21 + Tariff*d1 > + Tariff*d2 + Tariff*d10)* > > *summary(modeloRTUN2)* > > Siendo: > > AVE. = Variable dependiente (explicada) > Tariff = Variable dependiente (explicativa) > d1.....d22= los diferentes productos (dimensión producto: para comprobar > si el tipo producto tiene efecto sobre AVE.) > Tariff * d1 = efecto fijo del producto y tariff > > En resumen me gustaría comprobar si existe trade off entre las dos > variables NTM y Tariff y conocer que productos son significativos (conocer > que productos tienen presentan mas protección en frontera). > > Me surge ahora otra duda ...que diferencia hay en colocar "Tariff*d2" y > "Tariff:d2", en ambos casos es la multiplicación de las dos variables? > Adjunto la base de datos en formato excel al email, para que sea mas fácil > entenderme. > > Un abrazo mas y* muchas gracias, que tengáis un buen día.* > > Realmente todos los consejos han sido aportaciones muy valiosas para que > pudiese seguir avanzando, os lo agradezco! > Seguimos hablando ;) > > Lorena > > > El 13 de enero de 2014, 9:46, José Trujillo <trujillo en unex.es> escribió: > > Si AVE es la variable dependiente y Tarrif "LA" variable independiente: >> >> ¿Quién son d1 y d1T? >> >> Como dice Carlos, as.factor dentro de la fórmula no es buena idea, es >> más seguro que la variable d1T sea convertida en factor previamente. >> >> Saludos. >> >> >> El 12/01/14 13:53, Lorena Tudela Marco escribió: >> > Buenos días, >> > >> > Me gustaría aplicar una regresión múltiple a los datos con los que >> trabajo pero >> > no se como introducir los datos en R. He probado introducir el >> siguiente >> > comando: >> > >> > modeloM<-lm(AVE.~ d1 + as.factor(d1T)*Tariff) >> > summary(modeloM) >> > >> > Pero me da el siguiente error: >> > >> >> summary(modeloM) >> > Error in if (attr(z$terms, "intercept")) sum((f - mean(f))^2) else >> sum(f^2) >> > : >> > argument is of length zero >> > >> > Mi intención es calcular la siguiente regresión: Ave = cte + >> > d1*coeficiente1 + (d1T*Tariff)*coeficiente 2 >> > Siendo: AVE, la variable dependiente y Tarrif, la variable >> independiente. >> > >> > Imagino que no estoy introduciendo correctamente el comando inicial del >> > modelo y por eso me da error ¿Alguien esta trabajando en temas >> similares o >> > podría orientarme ? >> > >> > Muchas gracias, cualquier aportación sera bienvenida! >> > >> > Un abrazo, que tengáis un buen día. >> > >> > Lorena >> > >> > [[alternative HTML version deleted]] >> > >> > >> > >> > _______________________________________________ >> > R-help-es mailing list >> > R-help-es en r-project.org >> > https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> >> [[alternative HTML version deleted]] >> >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> >------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20140114/8760ae1b/attachment-0001.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: TUN.xlsx Type: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet Size: 133372 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20140114/8760ae1b/attachment-0001.bin>
Hola,
No sé exactamente lo que quieres hacer pero aparentemente se
trataría de un modelo de análisis de la covarianza. Entiendo que las
variables di son niveles de un factor (e.g. producto) y yo lo
codificaría en una única variable (además, en los datos que envías,
muchas de estas son siempre 0, por lo que no tiene sentido incluirlas -
no te salen resultados "extraños" con el modelo que utilizaste?). Se
entiende que el efecto lineal de la variable Tariff (explicativa?)
depende del nivel del factor. Tienes que mirar la formula
correspondiente al modelo buscado. Por ejemplo en cuanto a lo que
comentas, como dice en la ayuda: "The |*| operator denotes factor
crossing: |a*b| interpreted as |a+b+a:b"|, y R no trata igual a una
variable numérica que a un factor. Igual "AVE~producto*Tariff" es
justo
lo que buscas...
En resumen, mi recomendación es que estudies un poco más el tema
para tenerlo claro. Mirar un libro de modelos lineales con R puede estar
bien. Por ejemplo:
2005 - Linear Models with R - Faraway
Otros que tengo por aquí que te pueden valer:
2002 - An R and S Plus Companion to Applied Regression - J. Fox - SAGE
2012 - An R Companion to Linear Statistical Models - C. Hay-Jahans
Un saludo,
Rubén Fdez Casal
El 14/01/2014 12:28, Lorena Tudela Marco escribió:> Buenos días,
>
> *Muchas gracias, todas las aportaciones han sido bien útiles.* Las he
> tenido en cuenta y he pasado los datos con el R, siguiendo el
> siguiente comando:
>
> /modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8
> + d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20 +d21 +
> Tariff*d1 + Tariff*d2 + Tariff*d10)/
>
> /summary(modeloRTUN2)/
>
>
> Siendo:
>
> AVE. = Variable dependiente (explicada)
> Tariff = Variable dependiente (explicativa)
> d1.....d22= los diferentes productos (dimensión producto: para
> comprobar si el tipo producto tiene efecto sobre AVE.)
> Tariff * d1 = efecto fijo del producto y tariff
>
> En resumen me gustaría comprobar si existe trade off entre las dos
> variables NTM y Tariff y conocer que productos son significativos
> (conocer que productos tienen presentan mas protección en frontera).
>
> Me surge ahora otra duda ...que diferencia hay en colocar
"Tariff*d2"
> y "Tariff:d2", en ambos casos es la multiplicación de las dos
variables?
> Adjunto la base de datos en formato excel al email, para que sea mas
> fácil entenderme.
>
> Un abrazo mas y* muchas gracias, que tengáis un buen día.*
>
> Realmente todos los consejos han sido aportaciones muy valiosas para
> que pudiese seguir avanzando, os lo agradezco!
> Seguimos hablando ;)
>
> Lorena
>
>
> El 14 de enero de 2014, 12:24, Lorena Tudela Marco
> <lorenatudelamarco@gmail.com
<mailto:lorenatudelamarco@gmail.com>>
> escribió:
>
> Buenos días,
>
> *Muchas gracias, todas las aportaciones han sido bien útiles.* Las
> he tenido en cuenta y he pasado los datos con el R, siguiendo el
> siguiente comando:
>
> /modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 +
> d8 + d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20
> +d21 + Tariff*d1 + Tariff*d2 + Tariff*d10)/
>
> /summary(modeloRTUN2)/
>
>
> Siendo:
>
> AVE. = Variable dependiente (explicada)
> Tariff = Variable dependiente (explicativa)
> d1.....d22= los diferentes productos (dimensión producto: para
> comprobar si el tipo producto tiene efecto sobre AVE.)
> Tariff * d1 = efecto fijo del producto y tariff
>
> En resumen me gustaría comprobar si existe trade off entre las dos
> variables NTM y Tariff y conocer que productos son significativos
> (conocer que productos tienen presentan mas protección en frontera).
>
> Me surge ahora otra duda ...que diferencia hay en colocar
> "Tariff*d2" y "Tariff:d2", en ambos casos es la
multiplicación de
> las dos variables?
> Adjunto la base de datos en formato excel al email, para que sea
> mas fácil entenderme.
>
> Un abrazo mas y*muchas gracias, que tengáis un buen día.*
>
> Realmente todos los consejos han sido aportaciones muy valiosas
> para que pudiese seguir avanzando, os lo agradezco!
> Seguimos hablando ;)
>
> Lorena
>
>
> El 13 de enero de 2014, 9:46, José Trujillo <trujillo@unex.es
> <mailto:trujillo@unex.es>> escribió:
>
> Si AVE es la variable dependiente y Tarrif "LA" variable
> independiente:
>
> ¿Quién son d1 y d1T?
>
> Como dice Carlos, as.factor dentro de la fórmula no es buena
> idea, es
> más seguro que la variable d1T sea convertida en factor
> previamente.
>
> Saludos.
>
>
> El 12/01/14 13:53, Lorena Tudela Marco escribió:
> > Buenos días,
> >
> > Me gustaría aplicar una regresión múltiple a los datos con
> los que trabajo pero
> > no se como introducir los datos en R. He probado introducir
> el siguiente
> > comando:
> >
> > modeloM<-lm(AVE.~ d1 + as.factor(d1T)*Tariff)
> > summary(modeloM)
> >
> > Pero me da el siguiente error:
> >
> >> summary(modeloM)
> > Error in if (attr(z$terms, "intercept")) sum((f -
> mean(f))^2) else sum(f^2)
> > :
> > argument is of length zero
> >
> > Mi intención es calcular la siguiente regresión: Ave = cte +
> > d1*coeficiente1 + (d1T*Tariff)*coeficiente 2
> > Siendo: AVE, la variable dependiente y Tarrif, la variable
> independiente.
> >
> > Imagino que no estoy introduciendo correctamente el comando
> inicial del
> > modelo y por eso me da error ¿Alguien esta trabajando en
> temas similares o
> > podría orientarme ?
> >
> > Muchas gracias, cualquier aportación sera bienvenida!
> >
> > Un abrazo, que tengáis un buen día.
> >
> > Lorena
> >
> > [[alternative HTML version deleted]]
> >
> >
> >
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es@r-project.org <mailto:R-help-es@r-project.org>
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
> [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org <mailto:R-help-es@r-project.org>
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]