plotgen {capn} | R Documentation |
Plot Generator for Shadow Price or Value Function
Description
The function draws shadowp or vfun-w plot from the simulation results of vsim
, psim
, or pdotsim
.
Usage
plotgen(simres, ftype = NULL, whichs = NULL, tvar = NULL, xlabel = NULL,
ylabel = NULL)
Arguments
simres |
A simulation results from |
ftype |
Plot type (ftype=NULL (default) or ftype="p" for shadow price; ftype="vw" for vfun-w plot) |
whichs |
A pisitive integer for indicating a specific stock for multi-sotck cases (ftype=NULL (default) or 1<= whichs <= the number of stocks) |
tvar |
An array of time variable if simulation result is a time-base simulation |
xlabel |
A character for x-label of a plot (xlabel=NULL (default); "Stock" or "Time") |
ylabel |
An array of characters for y-label of a plot (ylabel=NULL (default); "Shadow Price", "Value Function" or "W-value") |
Details
This function provides an one-dimensional plot for "shadow price-stock", "shadow price-time", "Value function-stock", "Value function-time", "Value function-stock-W value", or "Value function-time-W value" depending on input arguments.
Value
A plot of approximation resuts: shadow (accounting) prices, inclusive wealth, and Value function
See Also
Examples
## 1-D Reef-fish example: see Fenichel and Abbott (2014)
data("GOM")
nodes <- chebnodegen(param$nodes,param$lowerK,param$upperK)
simuDataP <- cbind(nodes,sdot(nodes,param),
dsdotds(nodes,param),dwds(nodes,param))
Aspace <- aproxdef(param$order,param$lowerK,param$upperK,param$delta)
# p-approximation
pC <- paprox(Aspace,simuDataP[,1],simuDataP[,2],
simuDataP[,3],simuDataP[,4])
# Without prividing W-value
GOMSimP <- psim(pC,simuDataP[,1])
# With W-value
GOMSimP2 <- psim(pC,simuDataP[,1],profit(nodes,param),simuDataP[,2])
# Shadow price-Stock plot
plotgen(GOMSimP)
plotgen(GOMSimP,ftype="p")
plotgen(GOMSimP,xlabel="Stock Size, S", ylabel="Shadow Price (USD/Kg)")
# Value-Stock-W plot
plotgen(GOMSimP2,ftype="vw")
plotgen(GOMSimP2,ftype="vw",xlabel="Stock Size, S", ylabel="Value Function")
plotgen(GOMSimP2,ftype="vw",xlabel="Stock Size, S", ylabel="Value Function")
## 2-D Prey-Predator example
data("lvdata")
aproxdeg <- c(20,20)
lower <- c(0.1,0.1)
upper <- c(1.5,1.5)
delta <- 0.03
lvspace <- aproxdef(aproxdeg,lower,upper,delta)
lvaproxc <- vaprox(lvspace,lvaproxdata)
lvsim <- vsim(lvaproxc,lvsimdata.time[,2:3])
# Shadow price-Stock plot
plotgen(lvsim)
plotgen(lvsim,ftype="p")
plotgen(lvsim,whichs=2,xlabel="Stock Size, S",ylabel="Shadow Price (USD/Kg)")
# Shadow price-time plot
plotgen(lvsim,whichs=2,tvar=lvsimdata.time[,1])
# Value Function-Stock plot
plotgen(lvsim,ftype="vw")
plotgen(lvsim,ftype="vw",whichs=2,
xlabel="Stock Size, S",ylabel="Shadow Price (USD/Kg)")
# Value Function-time plot
plotgen(lvsim,ftype="vw",tvar=lvsimdata.time[,1])
plotgen(lvsim,ftype="vw",whichs=2,tvar=lvsimdata.time[,1],
xlabel="Stock Size, S",ylabel="Shadow Price (USD/Kg)")