plot.cesEst: Plot RSSs of a CES Function Estimated by Grid Search

In micEconCES: Analysis with the Constant Elasticity of Substitution (CES) Function

Description

Plot a scatter plot, where the values of rho are on the x axis and the corresponding sums of the squared residuals obtained by a grid search for rho are on the y axis. Estimations that did not converge are marked with red.

Note that this method can be applied only if the model was estimated by a grid search for rho, i.e.\ cesEst was called with argument rho set to a vector of more than one values for rho.

Usage

## S3 method for class 'cesEst'
plot( x, negRss = TRUE, bw = FALSE, ... )

Arguments

x	object returned by cesEst if it was called with argument rho set a vector containing more than one value for rho so that a grid search was performed.
negRss	logical. Indicates whether the negative sum of squared residuals should be plotted in 3D plots (ignored in 2D plots).
bw	logical. Indicates whether 3D plots should be in black-and-white or colored.
...	All further arguments are passed to plot.default or persp.

Author(s)

Arne Henningsen and Geraldine Henningsen

See Also

cesEst.

Examples

   data( germanFarms, package = "micEcon" )
   # output quantity:
   germanFarms$qOutput <- germanFarms$vOutput / germanFarms$pOutput
   # quantity of intermediate inputs
   germanFarms$qVarInput <- germanFarms$vVarInput / germanFarms$pVarInput

   ## CES: Land & Intermediate Inputs
   cesLandInt <- cesEst( yName = "qOutput",
      xNames = c( "land", "qVarInput" ), data = germanFarms,
      rho = seq( from = -0.6, to = 0.9, by = 0.3 ) )

   # plot the rhos against the sum of squared residuals
   plot( cesLandInt ) 

Example output

Loading required package: minpack.lm
Loading required package: DEoptim
Loading required package: parallel

DEoptim package
Differential Evolution algorithm in R
Authors: D. Ardia, K. Mullen, B. Peterson and J. Ulrich


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/
Warning messages:
1: In nls.lm(par = start, fn = residFun, data = data, jac = jac, yName = yName,  :
  lmder: info = -1. Number of iterations has reached `maxiter' == 50.

2: In nls.lm(par = start, fn = residFun, data = data, jac = jac, yName = yName,  :
  lmder: info = -1. Number of iterations has reached `maxiter' == 50.

micEconCES documentation built on Jan. 7, 2021, 3:01 p.m.