[BACK]Return to hgm.cwishart.Rd CVS log [TXT][DIR] Up to [local] / OpenXM / src / R / r-packages / hgm / man

Diff for /OpenXM/src/R/r-packages/hgm/man/hgm.cwishart.Rd between version 1.10 and 1.12

version 1.10, 2015/03/27 02:36:30 version 1.12, 2016/02/14 00:21:50
Line 1 
Line 1 
 % $OpenXM: OpenXM/src/R/r-packages/hgm/man/hgm.cwishart.Rd,v 1.9 2015/03/26 11:54:13 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.
Line 12 
Line 12 
 }  }
 \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 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
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.12
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>