aidsEst does an econometric estimation
of the Almost Ideal Demand System (AIDS)aidsEst( priceNames, shareNames, totExpName, data = NULL,
instNames = NULL, shifterNames = NULL,
method = "LA:L", hom = TRUE, sym = TRUE, pxBase,
estMethod = ifelse( is.null( instNames ), "SUR", "3SLS" ),
ILmaxiter = 50, ILtol = 1e-5, alpha0 = 0, restrict.regMat = FALSE, ... )## S3 method for class 'aidsEst':
print( x, ... )
aidsPx).systemfit).systemfit).aidsEst.aidsEst are passed to
systemfit;
additional arguments of print.aidsEst are currently ignored.The 'Iterative Linear Least Squares Estimator' (IL) needs starting values for the (translog) price index. The price index used to calculate the initial price index can be specified in the same way as for the LA-AIDS (e.g. 'IL:L')
aidsEst containing following objects:
systemfit.}
aidsEst.}
Blundell, R. and J.M. Robin (1999) Estimationin Large and Disaggregated Demand Systems: An Estimator for Conditionally Linear Systems. Journal of Applied Econometrics, 14, p. 209-232.
summary.aidsEst, aidsElas,
aidsCalc.
[object Object]
## Repeating the demand analysis of Blanciforti, Green & King (1986) ## Note: Blanciforti, Green & King (1986) use scaled data, ## which leads to slightly different results estResult <- aidsEst( c( "pFood1", "pFood2", "pFood3", "pFood4" ), c( "wFood1", "wFood2", "wFood3", "wFood4" ), "xFood", data = Blanciforti86, method = "LA:SL", maxiter = 100 ) print( estResult ) elas( estResult )
## Estimations with a demand shifter: linear trend priceNames <- c( "pFood1", "pFood2", "pFood3", "pFood4" ) shareNames <- c( "wFood1", "wFood2", "wFood3", "wFood4" ) Blanciforti86$trend <- c( 0:( nrow( Blanciforti86 ) - 1 ) ) estResult <- aidsEst( priceNames, shareNames, "xFood", data = Blanciforti86, shifterNames = "trend" ) print( estResult )
# Estimations with two demand shifters: linear + quadratic trend Blanciforti86$trend2 <- c( 0:( nrow( Blanciforti86 ) - 1 ) )^2 estResult <- aidsEst( priceNames, shareNames, "xFood", data = Blanciforti86, shifterNames = c( "trend", "trend2" ) ) print( estResult )