Dear all
If I run the model and get the estimated parameter a11. Then I want to
use the estimated parameter to calculate the elasticities by using the
formula e11=a11/mw1-1. What I have done is using the command of linear.
Hypothesis.
> formulas1=dWfresh~dlnPfresh+dlnPfrozen+dlnPsmoke+dlnQP+cosL1+sinL1+cosL2
> formulas2=dWsmoke~dlnPfresh+dlnPfrozen+dlnPsmoke+dlnQP+cosL1+sinL1+cosL2
> formulas12=list(formulas1,formulas2)
> > salmon1=systemfit(formulas12, data=Fsalmon2, method="SUR")
> summary (salmon1)
>MWfresh<-mean (Wfresh2, data=Fsalmon)
>MWfrozen<-mean (Wfrozen2, data=Fsalmon)
>MWsmoke<-mean (Wsmoke2, data=Fsalmon)
>restricte11=c("eq1_dlnPfresh/MWfresh-1=0")
>linear.hypothesis(salmonsdlnf1r, restricte11, test="Chisq")
It does not work. However, if I replace MWfresh with the pure numerical
number
restricte11=c("eq1_dlnPfresh/0.42-1=0")
It works.
Can anyone help me to work it out.?
Thanks
Jinghua Xie
University of Tromsø, Norway
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
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.
