version 1.1, 2014/03/31 06:20:06 |
version 1.4, 2015/03/22 00:27:51 |
|
|
% $OpenXM$ |
% $OpenXM: OpenXM/src/R/r-packages/hgm/man/hgm.mleFBByOptim.Rd,v 1.3 2015/03/21 23:40:34 takayama Exp $ |
\name{hgm.mleFBByOptim} |
\name{hgm.z.mleFBByOptim} |
%%Todo, write documents for hgm.mleDemo, hgm.ssFB |
\alias{hgm.z.mleFBByOptim} |
|
%%Todo, write documents for hgm.z.mleDemo, hgm.ssFB |
%\alias{hgm.ncso3} |
%\alias{hgm.ncso3} |
%- Also NEED an '\alias' for EACH other topic documented here. |
%- Also NEED an '\alias' for EACH other topic documented here. |
\title{ |
\title{ |
|
|
distribution on S^d. |
distribution on S^d. |
} |
} |
\usage{ |
\usage{ |
hgm.mleFBByOptim(d=0,ss=NULL,data=NULL,start=NULL,lb=NULL,up=NULL,bigValue=10000) |
hgm.z.mleFBByOptim(d=0,ss=NULL,data=NULL,start=NULL,lb=NULL,ub=NULL,bigValue=10000) |
} |
} |
%- maybe also 'usage' for other objects documented here. |
%- maybe also 'usage' for other objects documented here. |
\arguments{ |
\arguments{ |
|
|
\item{ss}{Sufficient statistics} |
\item{ss}{Sufficient statistics} |
\item{data}{ |
\item{data}{ |
The argument data is a set of data on the d-dimensional sphere. |
The argument data is a set of data on the d-dimensional sphere. |
It's format is an n by (d+1) matrix where n is the number of data. |
Its format is an n by (d+1) matrix where n is the number of data. |
The arguments ss and data are exclusive. When data is given, ss must be NULL |
When data is given, ss must be NULL |
and ss is calculated from data by hgm.ssFB(data). |
and ss is calculated from data by hgm.ssFB(data). |
} |
} |
\item{start}{ |
\item{start}{ |
|
|
is used for the optimization. |
is used for the optimization. |
The output is used as a starting point for the holonomic gradient method. |
The output is used as a starting point for the holonomic gradient method. |
See nk_fb_gen_c.rr of \url{http://www.math.kobe-u.ac.jp/Asir}. |
See nk_fb_gen_c.rr of \url{http://www.math.kobe-u.ac.jp/Asir}. |
|
This function generates temporary work files whose names start with tmp. |
|
\code{data <- read.table(fileName,sep=",")} can be used to read CSV data |
|
from a file. |
} |
} |
\value{ |
\value{ |
The output format is that of the function optim(). |
The output format is that of the function optim(). |
Line 70 T.Koyama, H.Nakayama, K.Nishiyama, N.Takayama. |
|
Line 74 T.Koyama, H.Nakayama, K.Nishiyama, N.Takayama. |
|
%% ~Make other sections like Warning with \section{Warning }{....} ~ |
%% ~Make other sections like Warning with \section{Warning }{....} ~ |
|
|
\seealso{ |
\seealso{ |
\code{\link{optim}}, |
\code{\link{optim}} |
} |
} |
\examples{ |
\examples{ |
## ===================================================== |
## ===================================================== |
Line 83 T.Koyama, H.Nakayama, K.Nishiyama, N.Takayama. |
|
Line 87 T.Koyama, H.Nakayama, K.Nishiyama, N.Takayama. |
|
0.3276, |
0.3276, |
0.0063,0.0054,0.0762) |
0.0063,0.0054,0.0762) |
start <- c(0.1,0.1,1,1,1,-1,-1,-1,-1) |
start <- c(0.1,0.1,1,1,1,-1,-1,-1,-1) |
hgm.mleFBByOptim(d=d,ss=ss,start=start) |
hgm.z.mleFBByOptim(d=d,ss=ss,start=start) |
} |
} |
} |
} |
% Add one or more standard keywords, see file 'KEYWORDS' in the |
% Add one or more standard keywords, see file 'KEYWORDS' in the |