AssetPricing (version 0.0-10)

vsolve: Solve for expected value of assets.

Description

Solves a system of coupled differential equations for the expected value of a number q of (perishable) assets, with q running from 1 to qmax, given a pricing policy. Treats the system in a vectorized form and uses the method of Runge-Kutta.

Usage

vsolve(S, lambda, gprob, tmax=NULL, x, nstep,
                   alpha=NULL, salval=0, nverb=0)

Arguments

S

An expression, or list of expressions, or a function or list of functions, specifying the price sensitivity functions S_j(x,t). See Details.

lambda

A function (of residual time t --- see tmax) or a positive constant specifying the intensity of the (generally inhomogeneous) Poisson process of arrival times of groups of potential customers.

gprob

A function (to calculate probabilities) or a numeric vector of probabilities determining the distribution of the size of an arriving group of customers. Must be compatible with certain characteristics of s (see below). See detai

tmax

The maximum residual time; think of this as being the initial time at which the assets go on sale (with time decreasing to zero, at which point the value of each asset drops to the salvage value (salval), usual

x

An object of class flap (see xsolve()) specifying the (given) pricing policy. It has the form of a list of functions x_i(t), with i running from 1 to qmax<

nstep

The number of (equal) sub-intervals into which the interval [0,tmax] is divided in the process of solving the system of differential equations numerically.

alpha

A numeric scalar between 0 and 1 specifying the probability that an arriving group of size j > q (where q is the number of assets remaining for sale) will consider purchasing (all of) these remaining assets. It is

salval

A (non-negative) numeric scalar specifying the salvage value of an asset --- i.e. the quantity to which the value of an asset drops at residual time t=0. Usually salval is equal to 0.

nverb

Non-negative integer; if nverb > 0 then a progress report (actually consisting solely of the step number) is printed out every nverb steps of the Runge-Kutta procedure. (If nverb==0 the procedu

Value

• A list with components
• x
  The argument x which was passed to vsolve, possibly with its tlim attribute modified. It is an object of class flap.
• v
  An object of class flap whose entries are (spline) functions v_q(t) specifying the expected value of q assets at time t as determined by numerically solving the coupled system of differential equations.
• vdot
  An object of class flap whose entries are the derivatives (with respect to t) of the functions v_q(t) described above. The values of these derivatives are determined as the left hand side of the differential equations being solved.

Details

The components of the argument S may be provided either as expressions or functions. If the former, these expressions should be amenble to differentiation with respect to x and t via the function deriv3(). This is essentially a matter of convenience; the derivatives are not actually used by vsolve. The expressions are turned into functions by deriv3() in the same manner as is used by xsolve(). See the help for xsolve() for further information about the required nature of S.

The argument tmax (if specified) must be less than or equal to the tmax attribute of argument S if S is a piecewise linear price sensitivity function, and must also be less than or equal to the tlim attribute of argument x.

If tmax is not specified it will be set equal to the tmax attribute of argument S if S is a piecewise linear price sensitivity function, in which case this attribute must be less than or equal to the tlim attribute of argument x. (If this is not so then S and x are incompatible.) Otherwise tmax will be set equal to the tlim attribute of argument x.

The argument gprob determines the range of possible values of the size of an arriving group of customers. The maximum value of this group size is in effect that value of j for which the corresponding probability value is numerically distinguishable from zero. If the argument x is a doubly indexed list of functions (was created with type="dip") then the maximum value of group size as determined by gprob must be compatible with the indexing scheme of x. That is to say, it must be less than or equal to the jmax attribute of x, otherwise an error is given. Note that if single indexing is in effect (i.e. x was created with type="sip") then this attribute is equal to 1, but for single indexing x does not depend on group size and so no restriction is imposed.

References

Baneree, P. K. and Turner, T. R. A flexible model for the pricing of perishable assets. Omega, vol. 40, number 5, 2012, pages 533 -- 540, doi: 10.1016/j.omega.2011.10.001.

See Also

xsolve(), plot.flap()

Examples

#
# In these examples "nstep" has been set equal to 100
# to reduce the time for package checking on CRAN.
# A more "realistic" value is nstep=1000.  Likewise "qmax"
# has been set equal to 5 which is an unrealistically low
# value for the total number of assets.
#
S <- expression(exp(-kappa*x/(1+gamma*exp(-beta*t))))
attr(S,"parvec") <- c(kappa=10/1.5,gamma=9,beta=1)
lambda1 <- function(tt){
#	eps <- sqrt(.Machine$double.eps)
	84*(1-tt)
}

# Optimal pricing policy assuming customers arrive singly:
X <- xsolve(S=S,lambda=lambda1,gprob=1,tmax=1,qmax=5,nstep=100,nverb=50)
lambda2 <- function(tt){
#	eps <- sqrt(.Machine$double.eps)
	36*(1-tt)
}
# Expected values if the customers actually arrive in groups, using the
# (sub-optimal) pricing policy based on the (erroneous) assumption that
# they arrive singly.  Note that the two scenarios are ``comparable'' in
# that the expected total number of customers is 42 in each case.
V <- vsolve(S=S,lambda=lambda2,gprob=(5:1)/15,x=X$x,nstep=100,
            alpha=0.5,nverb=50)