------------------------------------------------------------------------------- help forcircdpvm-------------------------------------------------------------------------------

Density probability plot for von Mises distribution fitted to circular data

circdpvmvarname[ifexp] [inrange] [,a(#)param(numlist)generate(newvar1 newvar2)line(line_options)graph_options]

Description

circdpvmgives a density probability plot for the fit of a von Mises (a.k.a. circular normal) distribution to a circular variable on a scale between 0 and 360 degrees. By default circvm is used to fit the distribution, estimating the two parameters vector mean mu and concentration parameter kappa. Both observed and expected densities are rotated so that each set is centred on the vector mean. Each is nevertheless labelled in terms ofvarname.

RemarksTo establish notation, and to fix ideas with a concrete example: consider an observed variable theta, whose distribution we wish to compare with a von Mises distributed variable phi. That variable has density function

f(phi), distribution functionP = F(phi) and quantile functionQ(P). (The distribution function and the quantile function are inverses of each other.) Clearly, this notation is fairly general and also covers other distributions, at least for continuous variables.The particular density function

f(theta | parameters) most pertinent to comparison with data for phi can be computed given values for its parameters, either estimates from data on theta, or parameter values chosen for some other good reason. In the case of a von Mises distribution, these parameters would usually be the vector mean and the concentration parameter kappa.The density function can also be computed indirectly via the quantile function as

f(Q(P)). In the case of von Mises distributions, it is convenient to centre variables at the vector mean, so thatP= 0.5 corresponds to the vector mean. In practicePis calculated as so-called plotting positionsp_iattached to values of a sample of sizenwhich have ranki. One simple rule usesp_i= (i- 0.5) /n. Most other rules follow one of a family (i-a) / (n- 2a+ 1) indexed bya.Plotting both

f(theta | parameters) andf(Q(P=p_i)), calculated using plotting positions, versus observed theta gives two curves. In our example, the first is von Mises by construction and the second would be a good estimate of a von Mises density if theta were truly von Mises with the same parameters. The match or mismatch between the curves allows graphical assessment of goodness or badness of fit. What is more, we can use experience from comparing frequency distributions, as shown on histograms, dot plots or other similar displays, in comparing or identifying location and scale differences, skewness, tail weight, tied values, gaps, outliers and so forth.Such

density probability plotswere suggested by Jones and Daly (1995). See also Jones (2004). For further discussion see dpplot.

Options

a()specifies a family of plotting positions, as explained above. The default is 0.5. Choice ofais rarely material unless the sample size is very small, and then the exercise is moot whatever is done.

param()specifies two parameter values which give an alternative reference distribution. The first is the vector mean and the second is the concentration parameter kappa. (By default these parameters are estimated from the data.)

generate()specifies two new variable names to hold the results of densities estimated from the data directly (asf() given parameters) and indirectly (asf(Q(P)) given parameters).

line(line_options)are options of twoway mspline and twoway line, which may be used to control the rendition of the density function curve.

show()specifies a numlist of axis labels to be shown, overriding the default.

graph_optionsare options of twoway.

Examples

. circdpvm wallasp. circdpvm wallasp, show(0(45)315)

ReferencesJones, M.C. 2004. Hazelton, M.L. (2003), "A graphical tool for assessing normality,"

The American Statistician57: 285-288: Comment.TheAmerican Statistician58: 176-177.Jones, M.C. and F. Daly. 1995. Density probability plots.

Communicationsin Statistics, Simulation and Computation24: 911-927.

AuthorNicholas J. Cox, University of Durham, U.K. n.j.cox@durham.ac.uk

Also seeOn-line: help for circvm, circpvm, dpplot (if installed)