Return to hgm.se.hgm.Bingham.Rd CVS log | Up to [local] / OpenXM / src / R / r-packages / hgm / man |
File: [local] / OpenXM / src / R / r-packages / hgm / man / hgm.se.hgm.Bingham.Rd (download)
Revision 1.2, Mon Mar 31 00:49:51 2014 UTC (10 years, 5 months ago) by takayama
hgm.ncBingham is an alias to hgm.se.hgm.Bingham. Example. > library(hgm) Loading required package: deSolve > hgm.ncBingham(c(1,3,5)); [1] 17.640021 2.371808 4.166681 9.192789 > hgm.ncBingham(c(1,3,5,6)); [1] 39.918359 3.948743 6.039318 10.854341 15.742235 > hgm.ncBingham(c(1,3,5,6,10)); [1] 302.18060 18.39045 24.12517 34.42539 43.00941 165.81628 > hgm.ncBingham(c(1,3,5,6,10,12)); [1] 1359.53150 65.15867 80.34557 104.24621 121.92847 304.56343 623.78715 > hgm.ncBingham(c(1,3,5,6,10,12,13)); [1] 3746.9919 157.2332 188.3426 234.1278 265.9428 539.4425 924.3876 [8] 1292.3312 |
\name{hgm.ncBingham} %\alias{hgm.se.hgm.Bingham} %- Also NEED an '\alias' for EACH other topic documented here. \title{ The function hgm.ncBingham performs the holonomic gradient method (HGM) for Bingham distributions. } \description{ The function hgm.ncBingham performs the holonomic gradient method (HGM) for Bingham distributions with the deSolve package in R. } \usage{ hgm.ncBingham(th, d=rep(1,length(th)+1), logarithm=FALSE, ini.method="power", times=NULL, withvol=FALSE, ...) } %- maybe also 'usage' for other objects documented here. \arguments{ \item{th}{ A (p-1)-dimensional vector which specifies the first (p-1) components of the parameter vector of the Bingham distribution on the (p-1)-dim sphere. The p-th parameter is assumed to be zero.} \item{d}{ A p-dimensional vector which specifies the multiplicity of the parameter. The default is all-one vector. } \item{logarithm}{ If 'logarithm' is TRUE, then the result is log of the normalizing constant. } \item{ini.method}{ The method for computing the initial value. Only "power" is implemented now. } \item{times}{ a vector; times in [0,1] at which explicit estimates for G are desired. If time = NULL, the set {0,1} is used, and only the final value is returned. } \item{withvol}{ If 'withvol' is TRUE, then the normalizing constant with volume of sphere is returned. Otherwise that without volume is returned. Therefore, if 'withvol' is FALSE and the parameter is zero, then the normalizing constant becomes 1. } \item{...}{ Additional parameters for computing initial values. Details are omitted. } } \details{ The function hgm.ncBingham computes the normalizing constant of the Bingham distribution and its derivatives at any specified point. The initial value is computed by the power series expansion. % \code{\link[RCurl]{postForm}}. } \value{ The output is p-dimensional vector G. The first element of G is the normalizing constant and the following (p-1)-elements are partial derivative of the normalizing constant with respect to the first (p-1) components of the parameter 'th'. } \references{ \url{http://www.math.kobe-u.ac.jp/OpenXM/Math/hgm/ref-hgm.html} } \author{ Tomonari Sei } \note{ %% ~~further notes~~ } %% ~Make other sections like Warning with \section{Warning }{....} ~ \seealso{ %%\code{\link{oxm.matrix_r2tfb}} } \examples{ # Example 1. hgm.ncBingham(c(1,3,5)) } % Add one or more standard keywords, see file 'KEYWORDS' in the % R documentation directory. \keyword{ Normalization constant } \keyword{ Holonomic gradient method } \keyword{ HGM }