Deybi Morales
2010-Jul-21 18:44 UTC
[R-es] Resp. (Solución)-cómo aplicar restricción lineal
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