Learn R Programming

micEcon (version 0.2-4)

aidsTestConsist: Consistency Test of the AIDS

Description

Test whether the specification of the AIDS is consistent with microeconomic demand theory (i.e. utility maximisation).

Usage

aidsTestConsist( pNames, wNames, xtName, data, coef = NULL,
      alpha0 = ifelse( is.null( coef$alpha0 ), 0, coef$alpha0 ) )

Arguments

pNames
a vector of strings containing the names of the prices.
wNames
a vector of strings containing the names of the expenditure shares.
xtName
a string containing the variable name of total expenditure.
data
a data frame containing the data.
coef
a list containing the coefficients alpha, beta and gamma.
alpha0
coefficient $\alpha_0$ of the translog price index.

Value

  • a list containing following elements:
  • monotonya logical vector indicating whether the monotony condition is fulfilled at each observation.
  • mPercentpercent of observations where the monotony condition is fulfilled.
  • concavitya logical vector indicating whether the concavity condition is fulfilled at each observation.
  • cPercentpercent of observations where the concavity condition is fulfilled.
  • cMatricesa list of the 'C' matrices for each observation to check for concavity (see Deaton and Muellbauer, 1980b, p.76 ).

References

Deaton, A.S. and J. Muellbauer (1980a) An Almost Ideal Demand System. American Economic Review, 70, p. 312-326.

Deaton, A.S. and J. Muellbauer (1980b) Economics and Consumer Behavior, Cambridge University Press, Cambridge.

See Also

aidsEst, aidsEla

Examples

Run this code
data( Blanciforti86 )
   # Data on food consumption are available only for the first 32 years
   Blanciforti86 <- Blanciforti86[ 1:32, ]

   pNames <- c( "pFood1", "pFood2", "pFood3", "pFood4" )
   wNames <- c( "wFood1", "wFood2", "wFood3", "wFood4" )

   estResult <- aidsEst( pNames, wNames, "xFood",
      data = Blanciforti86, method = "IL:L" )
   tc <- aidsTestConsist( pNames, wNames, "xFood", Blanciforti86,
      coef = estResult$coef )
   tc$mPercent   # great!
   tc$cPercent   # Oh, that is bad!

Run the code above in your browser using DataLab