R/aidsConsist.R

In micEconAids: Demand Analysis with the Almost Ideal Demand System (AIDS)

aidsConsist <- function( priceNames, totExpName, coef, data,
      priceIndex = "TL", basePrices = NULL, baseShares = NULL,
      shareNames = NULL ) {

   result <- list()

   nGoods <- length( coef$alpha )

   result$addingUp <- all(
      all.equal( sum( coef$alpha ), 1 ) == TRUE,
      all.equal( sum( coef$beta ), 0 ) == TRUE,
      all.equal( colSums( coef$gamma ), rep( 0, nGoods ),
         check.attributes = FALSE ) == TRUE
      )

   result$homogeneity <- all.equal( rowSums( coef$gamma ),
      rep( 0, nGoods ), check.attributes = FALSE ) == TRUE

   result$symmetry <- isSymmetric( coef$gamma, tol = 1e-10,
      check.attributes = FALSE ) == TRUE

   result$mono <- aidsMono( priceNames = priceNames, totExpName = totExpName,
      data = data, coef = coef, priceIndex = priceIndex,
      basePrices = basePrices, baseShares = baseShares )

   if( is.character( priceIndex ) ) {
      if( priceIndex == "TL" && result$symmetry ) {
         result$concav <- aidsConcav( priceNames = priceNames,
            totExpName = totExpName, data = data, coef = coef,
            shareNames = shareNames )
      }
   }

   class( result ) <- "aidsConsist"
   return( result )
}