tests/indices.R

In micEconIndex: Price and Quantity Indices

library( "micEconIndex" )
options( digits = 4, warn = 1 )

## function for printIndexing indices
printIndices <- function( what, ... ) {
  for( i in c( "Laspeyres", "Paasche", "Fisher" ) ) {
    cat( "\n", i, "\n" )
    if( what == "p" ) {
      index <- priceIndex( ..., method = i, weights = TRUE )
    } else {
      index <- quantityIndex( ..., method = i, weights = TRUE )
    }
    print( index )
    testRowSums <- rowSums( attributes( index )$weights[ !is.na( index ), ] )
    names( testRowSums ) <- NULL
    if( all.equal( testRowSums,
      rep( 1, sum( !is.na( index ) ) ) ) != TRUE ) {
      cat( "\nrowSums are not equal to one!!!\n\n" )
    }
  }
}

## Missong03E7.7
if( requireNamespace( 'micEcon', quietly = TRUE ) ) {
  data( Missong03E7.7, package = "micEcon" )
  ## price indices for Missong03E7.7
  printIndices( "p",  c( "p.beef", "p.veal", "p.pork" ),
    c( "q.beef", "q.veal", "q.pork" ), 1, Missong03E7.7 )
  
  ## quantity indices for Missong03E7.7
  printIndices( "q",  c( "p.beef", "p.veal", "p.pork" ),
    c( "q.beef", "q.veal", "q.pork" ), 1, Missong03E7.7 )
}

## Bleymueller79E25.1
if( requireNamespace( 'micEcon', quietly = TRUE ) ) {
  data( Bleymueller79E25.1, package = "micEcon" )
  ## price indices for Bleymueller79E25.1
  printIndices( "p",  c( "p.A", "p.B", "p.C", "p.D" ),
    c( "q.A", "q.B", "q.C", "q.D" ),  1, Bleymueller79E25.1 )
  
  ## quantity indices for Bleymueller79E25.1
  printIndices( "q",  c( "p.A", "p.B", "p.C", "p.D" ),
    c("q.A", "q.B", "q.C", "q.D" ), 1, Bleymueller79E25.1 )
}

## Electricity (Christensen & Greene 1976)
if( requireNamespace( 'Ecdat', quietly = TRUE ) ) {
  data( "Electricity", package = "Ecdat" )
  Electricity <- Electricity[ 1:35, ]
  ## preparing electricity data
  pNames <- c( "pl", "pk", "pf" )
  qNames <- c( "ql", "qk", "qf" )
  sNames <- c( "sl", "sk", "sf" )
  for( i in 1:3 ) {
    Electricity[[ qNames[ i ] ]] <- Electricity[[ sNames[ i ] ]] *
      Electricity[[ "cost" ]] / Electricity[[ pNames[ i ] ]]
  }
  allObs <- !is.na( Electricity$pl )
  
  ## price indices for Electricity data
  printIndices( "p",  pNames, qNames, 1, Electricity )
  
  ## quantity indices for Electricity data
  printIndices( "q",  pNames, qNames, 1, Electricity )
  
  ## price indices for Electricity data and base=mean
  printIndices( "p",  pNames, qNames, allObs, Electricity )
  
  ## quantity indices for Electricity data and base=mean
  printIndices( "q",  pNames, qNames, allObs, Electricity )
}

## Electricity data with some NA prices
## manipulating data of Electricity data
if( requireNamespace( 'micEcon', quietly = TRUE ) ) {
  ElectricityNaP <- Electricity
  for( i in 1:3 ) {
    ElectricityNaP[[ pNames[ i ] ]][ c( 2, i * 4, i * 8 ) ] <- NA
  }
  
  ## price indices for Electricity data with some NA prices
  printIndices( "p",  pNames, qNames, 1, ElectricityNaP )
  
  ## quantity indices for Electricity data with some NA prices
  printIndices( "q",  pNames, qNames, 1, ElectricityNaP )
  
  ## price indices for Electricity data with some NA prices and na.rm=TRUE
  printIndices( "p",  pNames, qNames, 1, ElectricityNaP, na.rm = TRUE )
  
  ## quantity indices for Electricity data with some NA prices and na.rm=TRUE
  printIndices( "q",  pNames, qNames, 1, ElectricityNaP, na.rm = TRUE )
  
  ## price indices for Electricity data with some NA prices and base=mean
  printIndices( "p",  pNames, qNames, 16, ElectricityNaP )
  
  ## quantity indices for Electricity data with some NA prices and base=mean
  printIndices( "q",  pNames, qNames, allObs, ElectricityNaP )
  
  ## price indices for Electricity data with some NA prices and na.rm=TRUE and base=mean
  printIndices( "p",  pNames, qNames, allObs, ElectricityNaP, na.rm = TRUE )
  
  ## quantity indices for Electricity data with some NA prices and na.rm=TRUE and base=mean
  printIndices( "q",  pNames, qNames, allObs, ElectricityNaP, na.rm = TRUE )
}

## Electricity data with some NA prices and quantities
## manipulating data of Electricity data
if( requireNamespace( 'micEcon', quietly = TRUE ) ) {
  ElectricityNaPQ <- ElectricityNaP
  for( i in 1:3 ) {
    ElectricityNaPQ[[ qNames[ i ] ]][ c( 2, ( i + 1 ) * 4, i * 8 ) ] <- NA
  }
  
  ## price indices for Electricity data with some NA prices and quantities
  printIndices( "p",  pNames, qNames, 1, ElectricityNaPQ )
  
  ## quantity indices for Electricity data with some NA prices and quantities
  printIndices( "q",  pNames, qNames, 1, ElectricityNaPQ )
  
  ## price indices for Electricity data with some NA prices and quantities and na.rm=TRUE
  printIndices( "p",  pNames, qNames, 1, ElectricityNaPQ, na.rm = TRUE )
  
  ## quantity indices for Electricity data with some NA prices and quantities and na.rm=TRUE
  printIndices( "q",  pNames, qNames, 1, ElectricityNaPQ, na.rm = TRUE )
  
  ## price indices for Electricity data with some NA prices and quantities and base=mean
  printIndices( "p",  pNames, qNames, 16, ElectricityNaPQ )
  
  ## quantity indices for Electricity data with some NA prices and quantities and base=mean
  printIndices( "q",  pNames, qNames, 16, ElectricityNaPQ )
  
  ## price indices for Electricity data with some NA prices and quantities and na.rm=TRUE and base=mean
  printIndices( "p",  pNames, qNames, allObs, ElectricityNaPQ, na.rm = TRUE )
  
  ## quantity indices for Electricity data with some NA prices and quantities and na.rm=TRUE and base=mean
  printIndices( "q",  pNames, qNames, allObs, ElectricityNaPQ, na.rm = TRUE )
}

## Electricity data with some NA prices, where quantities are partly zero
## manipulating data of Electricity data
if( requireNamespace( 'micEcon', quietly = TRUE ) ) {
  ElectricityNaP0Q <- Electricity
  for( i in 1:3 ) {
    ElectricityNaP0Q[[ pNames[ i ] ]][ c( 2, i * 4, i * 8 ) ] <- NA
    ElectricityNaP0Q[[ qNames[ i ] ]][ c( 2, i * 8 ) ] <- 0
  }
  
  ## price indices, base = 1
  printIndices( "p",  pNames, qNames, 1, ElectricityNaP0Q )
  
  ## quantity indices, base = 1
  printIndices( "q",  pNames, qNames, 1, ElectricityNaP0Q )
  
  ## price indices, base = mean, na.rm = TRUE
  printIndices( "p",  pNames, qNames, allObs, ElectricityNaP0Q, na.rm = TRUE )
  
  ## quantity indices, base = mean, na.rm = TRUE
  printIndices( "q",  pNames, qNames, allObs, ElectricityNaP0Q, na.rm = TRUE )
}