[R-es] Resp. (Solución)-cómo aplicar restricción lineal

Deybi Morales morales.economia en gmail.com
Mie Jul 21 20:44:09 CEST 2010


He podido aplicar la restricción lineal a los coeficientes sin muchas
complicaciones aplicando el paquete systemfit.

Este es un ejemplo de como he seguido el proceso para aplicar la
restricción a los coeficientes de un Mínimo Cuadrado Ordinario.

#Paquete necesario
#aplique si no lo tiene instalado
install.packages("systemfit")
#correr paquete systemfit
library(systemfit)

#regresión OLS
probando <-systemfit(p1q1 ~ p1+m+p2+p3, method = "OLS")
summary(probando)
print(probando)

#aplicar la restricción a los coeficientes
restriccion <- "eq1_(Intercept)+ eq1_p1 + eq1_m + eq1_p2 + eq1_p3 = 1"
probando <-systemfit(p1q1 ~ p1+m+p2+p3, method = "OLS",
restrict.matrix =+restriccion)
summary(probando)

#comprobando si se cumple la restricción
print(probando)
b1<-as.numeric(probando$coefficients[1])
b2<-as.numeric(probando$coefficients[2])
b3<-as.numeric(probando$coefficients[3])
b4<-as.numeric(probando$coefficients[4])
b5<-as.numeric(probando$coefficients[5])
sumadecoeficientes<-b1+b2+b3+b4+b5
sumadecoeficientes

Muchas gracias a los señores  Olivier Nuñez y Juan José Gibaja
Martíns, por responder a mi ayuda. Tomé muy en cuenta sus
recomendaciones, espero que también les sea útil el proceso que seguí
con el paquete systemfit.

Deybi Morales León

morales.economia en gmail.com
http://moraleseconomia.blogspot.com
http://www.facebook.com/pages/Economia-Aplicada-Por-una-aficion-la-economia/113850775323574



Más información sobre la lista de distribución R-help-es