cheese                package:bayesm                R Documentation

_S_l_i_c_e_d _C_h_e_e_s_e _D_a_t_a

_D_e_s_c_r_i_p_t_i_o_n:

     Panel data with sales volume for a package of Borden Sliced Cheese
     as well as a measure of display activity and price.  Weekly data
     aggregated to the "key" account or retailer/market level.

_U_s_a_g_e:

     data(cheese)

_F_o_r_m_a_t:

     A data frame with 5555 observations on the following 4 variables.

     '_R_E_T_A_I_L_E_R' a list of 88 retailers

     '_V_O_L_U_M_E' unit sales

     '_D_I_S_P' a measure of display activity - per cent ACV on display

     '_P_R_I_C_E' in $

_S_o_u_r_c_e:

     Boatwright et al (1999), "Account-Level Modeling for Trade
     Promotion,"  _JASA_ 94, 1063-1073.

_R_e_f_e_r_e_n_c_e_s:

     Chapter 3, _Bayesian Statistics and Marketing_ by Rossi et al. 
      <URL:
     http://faculty.chicagogsb.edu/peter.rossi/research/bsm.html>

_E_x_a_m_p_l_e_s:

     data(cheese)
     cat(" Quantiles of the Variables ",fill=TRUE)
     mat=apply(as.matrix(cheese[,2:4]),2,quantile)
     print(mat)

     ##
     ## example of processing for use with rhierLinearModel
     ##
     if(0)
     {

     retailer=levels(cheese$RETAILER)
     nreg=length(retailer)
     nvar=3
     regdata=NULL
     for (reg in 1:nreg) {
             y=log(cheese$VOLUME[cheese$RETAILER==retailer[reg]])
             iota=c(rep(1,length(y)))
             X=cbind(iota,cheese$DISP[cheese$RETAILER==retailer[reg]],
                     log(cheese$PRICE[cheese$RETAILER==retailer[reg]]))
             regdata[[reg]]=list(y=y,X=X)
     }
     Z=matrix(c(rep(1,nreg)),ncol=1)
     nz=ncol(Z)
     ##
     ## run each individual regression and store results
     ##
     lscoef=matrix(double(nreg*nvar),ncol=nvar)
     for (reg in 1:nreg) {
             coef=lsfit(regdata[[reg]]$X,regdata[[reg]]$y,intercept=FALSE)$coef
             if (var(regdata[[reg]]$X[,2])==0)  { lscoef[reg,1]=coef[1]; lscoef[reg,3]=coef[2]}
             else {lscoef[reg,]=coef }
     }

     R=2000
     Data=list(regdata=regdata,Z=Z)
     Mcmc=list(R=R,keep=1)

     set.seed(66)
     out=rhierLinearModel(Data=Data,Mcmc=Mcmc)

     cat("Summary of Delta Draws",fill=TRUE)
     summary(out$Deltadraw)
     cat("Summary of Vbeta Draws",fill=TRUE)
     summary(out$Vbetadraw)

     if(0){
     #
     # plot hier coefs
     plot(out$betadraw)
     }

     }

