version 1.9, 2015/03/26 11:54:13 |
version 1.12, 2016/02/14 00:21:50 |
|
|
% $OpenXM: OpenXM/src/R/r-packages/hgm/man/hgm.cwishart.Rd,v 1.8 2015/03/22 00:27:51 takayama Exp $ |
% $OpenXM: OpenXM/src/R/r-packages/hgm/man/hgm.cwishart.Rd,v 1.11 2016/02/13 22:56:50 takayama Exp $ |
\name{hgm.pwishart} |
\name{hgm.pwishart} |
\alias{hgm.pwishart} |
\alias{hgm.pwishart} |
%- Also NEED an '\alias' for EACH other topic documented here. |
%- Also NEED an '\alias' for EACH other topic documented here. |
|
|
} |
} |
\usage{ |
\usage{ |
hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
err,automatic,assigned_series_error,verbose) |
err,automatic,assigned_series_error,verbose,autoplot) |
} |
} |
%- maybe also 'usage' for other objects documented here. |
%- maybe also 'usage' for other objects documented here. |
\arguments{ |
\arguments{ |
Line 33 hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
|
Line 33 hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
|
} |
} |
\item{dp}{ |
\item{dp}{ |
Sampling interval of solutions by the Runge-Kutta method. |
Sampling interval of solutions by the Runge-Kutta method. |
|
When autoplot=1, it is automatically set. |
} |
} |
\item{q}{ |
\item{q}{ |
The second value y[0] of this function is the Prob(L1 < q) |
The second value y[0] of this function is the Prob(L1 < q) |
Line 40 hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
|
Line 41 hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
|
} |
} |
\item{mode}{ |
\item{mode}{ |
When mode=c(1,0,0), it returns the evaluation |
When mode=c(1,0,0), it returns the evaluation |
of the matrix hypergeometric series and its derivatives at x0. |
of the matrix hypergeometric series and its derivatives at q0. |
When mode=c(1,1,(m^2+1)*p), intermediate values of P(L1 < x) with respect to |
When mode=c(1,1,(m^2+1)*p), intermediate values of P(L1 < x) with respect to |
p-steps of x are also returned. Sampling interval is controled by dp. |
p-steps of x are also returned. Sampling interval is controled by dp. |
|
When autoplot=1, it is automatically set. |
} |
} |
\item{method}{ |
\item{method}{ |
a-rk4 is the default value. |
a-rk4 is the default value. |
Line 71 hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
|
Line 73 hgm.pwishart(m,n,beta,q0,approxdeg,h,dp,q,mode,method, |
|
If it is 1, then steps of automatic degree updates and several parameters |
If it is 1, then steps of automatic degree updates and several parameters |
are output to stdout and stderr. |
are output to stdout and stderr. |
} |
} |
|
\item{autoplot}{ |
|
autoplot=0 is the default value. |
|
If it is 1, then it outputs an input for plot. |
|
When ans is the output, ans[1,] is c(q,prob at q,...), ans[2,] is c(q0,prob at q0,...), and ans[3,] is c(q0+q/100,prob at q/100,...), ... |
|
} |
} |
} |
\details{ |
\details{ |
It is evaluated by the Koev-Edelman algorithm when x is near the origin and |
It is evaluated by the Koev-Edelman algorithm when x is near the origin and |
by the HGM when x is far from the origin. |
by the HGM when x is far from the origin. |
We can obtain more accurate result when the variables h is smaller, |
We can obtain more accurate result when the variables h is smaller, |
x0 is relevant value (not very big, not very small), |
q0 is relevant value (not very big, not very small), |
and the approxdeg is more larger. |
and the approxdeg is more larger. |
A heuristic method to set parameters x0, h, approxdeg properly |
A heuristic method to set parameters q0, h, approxdeg properly |
is to make x larger and to check if the y[0] approaches to 1. |
is to make x larger and to check if the y[0] approaches to 1. |
% \code{\link[RCurl]{postForm}}. |
% \code{\link[RCurl]{postForm}}. |
} |
} |
Line 109 The obtained value is not validated automatically. |
|
Line 116 The obtained value is not validated automatically. |
|
|
|
%% ~Make other sections like Warning with \section{Warning }{....} ~ |
%% ~Make other sections like Warning with \section{Warning }{....} ~ |
|
|
\seealso{ |
%\seealso{ |
%%\code{\link{oxm.matrix_r2tfb}} |
%%%\code{\link{oxm.matrix_r2tfb}} |
} |
%} |
\examples{ |
\examples{ |
## ===================================================== |
## ===================================================== |
## Example 1. |
## Example 1. |
Line 122 hgm.pwishart(m=3,n=5,beta=c(1,2,3),q=10) |
|
Line 129 hgm.pwishart(m=3,n=5,beta=c(1,2,3),q=10) |
|
## ===================================================== |
## ===================================================== |
b<-hgm.pwishart(m=4,n=10,beta=c(1,2,3,4),q0=1,q=10,approxdeg=20,mode=c(1,1,(16+1)*100)); |
b<-hgm.pwishart(m=4,n=10,beta=c(1,2,3,4),q0=1,q=10,approxdeg=20,mode=c(1,1,(16+1)*100)); |
c<-matrix(b,ncol=16+1,byrow=1); |
c<-matrix(b,ncol=16+1,byrow=1); |
|
#plot(c) |
|
## ===================================================== |
|
## Example 3. |
|
## ===================================================== |
|
c<-hgm.pwishart(m=4,n=10,beta=c(1,2,3,4),q0=1,q=10,approxdeg=20,autoplot=1); |
#plot(c) |
#plot(c) |
} |
} |
% Add one or more standard keywords, see file 'KEYWORDS' in the |
% Add one or more standard keywords, see file 'KEYWORDS' in the |