Dear all
I am working on a paper where I am analysing energy-non energy substitution
possibilities in a manufacturing sector. In this case, I am estimating a
translog cost function jointly with factor share equations under SUR
modeling.
The estimation requires me to impose linear homogeneity restrictions, where
coefficients of price variables are required to sum to one as other
coefficients sum to zero. Below are my R commands for this estimation:
D<-c("q1_lnPk + q1_lnPl + q1_lnPe =1", "q1_lnq + q1_lnPklnPk + q1_lnPllnPl
+ q1_lnPelnPe + q1_lnqlnq + q1_lnPklnPl + q1_lnPklnPe + q1_lnPllnPe +
q1_lnPklnq + q1_lnPllnq + q1_lnPelnq=0")
z<-plm(list(q1=lnTC ~ lnPk + lnPl + lnPe + lnq + lnPklnPk + lnPllnPl +
lnPelnPe + lnqlnq + lnPklnPl + lnPklnPe + lnPllnPe + lnPklnq + lnPllnq +
lnPelnq,
q2= Sl ~ slnpl +slnpe +lnq,
q3= Se ~ slnpl+ slnpe+ lnq), Surpanel, effect = "twoways",
restrict.matrix = D, model="pooling" )
summary(z)
The first line of the commands indicates the matrix of restrictions, the
second line indicates the formula and the third line indicates the summary
of findings.
However, after running the commands, I am getting the following error:
*Error in rep(0, nrow(restrict.matrix)) : invalid 'times' argument*
Please help me in getting a solution to this error.
Regards,
Kigundu Ken
plm
2 messages · KENNETH KIGUNDU, Bert Gunter
Well, you may get lucky and get a reply here, but per the posting guide linked below, which you are **supposed to** read and follow **before** posting: "For questions about functions in standard packages distributed with R (see the FAQ Add-on packages in R <http://cran.r-project.org/doc/FAQ/R-FAQ.html#Add-on-packages-in-R>), ask questions on R-help. If the question relates to a *contributed package* , e.g., one downloaded from CRAN, try contacting the package maintainer first. You can also use find("functionname") and packageDescription("packagename") to find this information. *Only* send such questions to R-help or R-devel if you get no reply or need further assistance. This applies to both requests for help and to bug reports." So did you first attempt to contact said maintainer? If so, please tell us with what results. You also probably should tell us what package plm comes from (several contain this function) -- though I suspect only congniscenti who use the package could help you anyway, and they would know this already. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Dec 22, 2020 at 10:30 AM KENNETH KIGUNDU <kigundukenneth at gmail.com> wrote:
Dear all
I am working on a paper where I am analysing energy-non energy substitution
possibilities in a manufacturing sector. In this case, I am estimating a
translog cost function jointly with factor share equations under SUR
modeling.
The estimation requires me to impose linear homogeneity restrictions, where
coefficients of price variables are required to sum to one as other
coefficients sum to zero. Below are my R commands for this estimation:
D<-c("q1_lnPk + q1_lnPl + q1_lnPe =1", "q1_lnq + q1_lnPklnPk + q1_lnPllnPl
+ q1_lnPelnPe + q1_lnqlnq + q1_lnPklnPl + q1_lnPklnPe + q1_lnPllnPe +
q1_lnPklnq + q1_lnPllnq + q1_lnPelnq=0")
z<-plm(list(q1=lnTC ~ lnPk + lnPl + lnPe + lnq + lnPklnPk + lnPllnPl +
lnPelnPe + lnqlnq + lnPklnPl + lnPklnPe + lnPllnPe + lnPklnq + lnPllnq +
lnPelnq,
q2= Sl ~ slnpl +slnpe +lnq,
q3= Se ~ slnpl+ slnpe+ lnq), Surpanel, effect = "twoways",
restrict.matrix = D, model="pooling" )
summary(z)
The first line of the commands indicates the matrix of restrictions, the
second line indicates the formula and the third line indicates the summary
of findings.
However, after running the commands, I am getting the following error:
*Error in rep(0, nrow(restrict.matrix)) : invalid 'times' argument*
Please help me in getting a solution to this error.
Regards,
Kigundu Ken
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.