USMeatConsump: U.S. Meat Consumption Data
In micEconAids: Demand Analysis with the Almost Ideal Demand System (AIDS)

Description Usage Format Source Examples

The USMeatConsump data set contains quarterly retail prices and consumption quantities for four meat product categories: beef, pork, chicken, and turkey. The data period ranges from the first quarter of 1975 to the third quarter of 1999. Hence, there are 99 observations.

1	data(USMeatConsump)

This data frame contains the following columns:

year: Year.
qtr: Quarter of the year.
t: Time trend.
pop: Population [million].
cpi: Consumer price index.
total_exp: Total per capita expenditure.
meat_exp: Per capita expenditure on meat.
beef_q: Per capita consumption of beef [pound].
pork_q: Per capita consumption of pork[pound].
chick_q: Per capita consumption of chicken [pound].
turkey_q: Per capita consumption of turkey [pound].
beef_p: Retail price of beef [cents / pound].
pork_p: Retail price of pork [cents / pound].
chick_p: Retail price of chicken [cents / pound].
turkey_p: Retail price of turkey [cents / pound].
beef_w: Expenditure share of beef (in meat).
pork_w: Expenditure share of pork (in meat).
chick_w: Expenditure share of chicken (in meat).
turkey_w: Expenditure share of turkey (in meat).

SAS, SAS/ETS Examples: Estimating an Almost Ideal Demand System Model, https://support.sas.com/rnd/app/ets/examples/aids/index.htm.

   ## replicating the LA-AIDS estimation of the SAS example
   # loading data set
   data( USMeatConsump )

   # adding shifter variables for modeling seasonal effects
   USMeatConsump$co1 <- cos( 1 / 2 * 3.14159 * USMeatConsump$t )
   USMeatConsump$si1 <- sin( 1 / 2 * 3.14159 * USMeatConsump$t )

   # Scaling prices by their means
   USMeatConsump$beef_pm <- USMeatConsump$beef_p / mean( USMeatConsump$beef_p )
   USMeatConsump$pork_pm <- USMeatConsump$pork_p / mean( USMeatConsump$pork_p )
   USMeatConsump$chick_pm <- USMeatConsump$chick_p / mean( USMeatConsump$chick_p )
   USMeatConsump$turkey_pm <- USMeatConsump$turkey_p / mean( USMeatConsump$turkey_p )

   # Estimation of the model
   meatModel <- aidsEst( c( "beef_pm", "pork_pm", "chick_pm", "turkey_pm" ),
      c( "beef_w", "pork_w", "chick_w", "turkey_w" ),
      "meat_exp", shifterNames = c( "co1", "si1", "t" ),
      priceIndex ="S", data = USMeatConsump, maxiter=1000 )
   summary( meatModel )

Loading required package: lmtest
Loading required package: zoo

Attaching package: 'zoo'

The following objects are masked from 'package:base':

    as.Date, as.Date.numeric

Loading required package: micEcon

If you have questions, suggestions, or comments regarding one of the 'micEcon' packages, please use a forum or 'tracker' at micEcon's R-Forge site:
https://r-forge.r-project.org/projects/micecon/

Demand analysis with the Almost Ideal Demand System (AIDS)
Estimation Method: Linear Approximation (LA) with Stone Index (S)
Estimated Coefficients:
             Estimate  Std. Error  t value  Pr(>|t|)    
alpha 1    3.9790e-01  3.9850e-01   0.9985  0.318915    
alpha 2    2.4096e-01  2.8614e-01   0.8421  0.400470    
alpha 3    5.7960e-01  1.7562e-01   3.3004  0.001092 ** 
alpha 4   -2.1846e-01  2.9296e-01  -0.7457  0.456477    
beta 1     2.3424e-02  4.4072e-02   0.5315  0.595501    
beta 2     3.3858e-03  3.1644e-02   0.1070  0.914870    
beta 3    -5.3550e-02  1.9440e-02  -2.7546  0.006267 ** 
beta 4     2.6740e-02  3.2459e-02   0.8238  0.410767    
gamma 1 1  4.1175e-02  1.2616e-02   3.2638  0.001238 ** 
gamma 1 2 -4.0660e-03  8.0822e-03  -0.5031  0.615313    
gamma 1 3 -4.9260e-02  5.3807e-03  -9.1550 < 2.2e-16 ***
gamma 1 4  1.2151e-02  9.5605e-03   1.2710  0.204806    
gamma 2 1 -4.0660e-03  8.0822e-03  -0.5031  0.615313    
gamma 2 2  4.4733e-02  8.7628e-03   5.1048 6.176e-07 ***
gamma 2 3 -3.0740e-02  5.9886e-03  -5.1331 5.387e-07 ***
gamma 2 4 -9.9267e-03  7.5019e-03  -1.3232  0.186855    
gamma 3 1 -4.9260e-02  5.3807e-03  -9.1550 < 2.2e-16 ***
gamma 3 2 -3.0740e-02  5.9886e-03  -5.1331 5.387e-07 ***
gamma 3 3  1.0139e-01  7.5646e-03  13.4029 < 2.2e-16 ***
gamma 3 4 -2.1388e-02  6.9571e-03  -3.0742  0.002322 ** 
gamma 4 1  1.2151e-02  9.5605e-03   1.2710  0.204806    
gamma 4 2 -9.9267e-03  7.5019e-03  -1.3232  0.186855    
gamma 4 3 -2.1388e-02  6.9571e-03  -3.0742  0.002322 ** 
gamma 4 4  1.9163e-02  1.2331e-02   1.5541  0.121316    
delta 1 1 -1.6311e-02  1.6215e-03 -10.0587 < 2.2e-16 ***
delta 1 2 -2.2825e-03  1.5621e-03  -1.4611  0.145122    
delta 1 3 -1.4775e-03  4.8256e-05 -30.6175 < 2.2e-16 ***
delta 2 1  9.6112e-03  1.1803e-03   8.1427 1.342e-14 ***
delta 2 2  6.4541e-03  1.1305e-03   5.7090 2.932e-08 ***
delta 2 3  2.7848e-04  3.5980e-05   7.7397 1.890e-13 ***
delta 3 1 -5.6697e-03  6.9619e-04  -8.1439 1.332e-14 ***
delta 3 2 -8.3559e-04  6.6021e-04  -1.2656  0.206707    
delta 3 3  8.9276e-04  2.6218e-05  34.0517 < 2.2e-16 ***
delta 4 1  1.2369e-02  1.1458e-03  10.7952 < 2.2e-16 ***
delta 4 2 -3.3360e-03  1.0925e-03  -3.0537  0.002481 ** 
delta 4 3  3.0626e-04  4.6136e-05   6.6380 1.680e-10 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
R-squared Values of expenditure shares:
    beef_w     pork_w    chick_w   turkey_w 
 0.9542558  0.7775530  0.9783787 -1.3797828 
R-squared Values of quantities:
  q_beef_w   q_pork_w  q_chick_w q_turkey_w 
-1.0106297  0.2530835 -2.6273590 -0.9436939