| aidsCalc {micEcon} | R Documentation |
Given prices, total expenditure and coefficients this function calculates the demanded quantities and expenditure shares based on the Almost Ideal Demand System.
aidsCalc( priceNames, totExpName, data = NULL, px = "TL", lnp = NULL,
coef = NULL, alpha0 = ifelse( is.null( coef$alpha0 ), 0,
coef$alpha0 ) )
priceNames |
a vector of strings containing the names of the prices. |
totExpName |
a string containing the variable name of total expenditure. |
data |
a data frame containing the data. |
px |
prices index (see aidsPx). |
lnp |
log of the price index (only required if 'px' is not 'TL'). |
coef |
a list containing the coefficients alpha, beta and gamma. |
alpha0 |
coefficient α_0 of the translog price index. |
aidsCalc returns a list with following elements:
shares |
a data frame containing the calculated expenditure shares. |
quantities |
a data frame containing the calculated quantites. |
Arne Henningsen ahenningsen@agric-econ.uni-kiel.de
Deaton, A.S. and J. Muellbauer (1980) An Almost Ideal Demand System. American Economic Review, 70, p. 312-326.
data( Blanciforti86 )
# Data on food consumption are available only for the first 32 years
Blanciforti86 <- Blanciforti86[ 1:32, ]
priceNames <- c( "pFood1", "pFood2", "pFood3", "pFood4" )
shareNames <- c( "wFood1", "wFood2", "wFood3", "wFood4" )
## LA-AIDS
estResult <- aidsEst( priceNames, shareNames, "xFood",
data = Blanciforti86, method = "LA:L" )
lnp <- aidsPx( "L", priceNames, shareNames, Blanciforti86 )
fitted <- aidsCalc( priceNames, "xFood", Blanciforti86,
coef = estResult$coef, lnp = lnp )
fitted$shares # equal to estResult$wFitted
fitted$quant # equal to estResult$qFitted
## AIDS
estResult <- aidsEst( priceNames, shareNames, "xFood",
data = Blanciforti86, method = "IL:L" )
fitted <- aidsCalc( priceNames, "xFood", Blanciforti86,
coef = estResult$coef )
fitted$shares # equal to estResult$wFitted
fitted$quant # equal to estResult$qFitted