示例#1
0
文件: phi.cpp 项目: Horta/BASiCS
vec Phi::get(void) const
{
  vec phi(kappa.size() + 1);
  phi(0) = phi.size();
  phi.tail(kappa.size()) = phi(0) * exp(kappa) / sum(exp(kappa));
  return phi;
}
示例#2
0
文件: phi.cpp 项目: Horta/BASiCS
Phij Phi::phij(size_t j) const
{
  double phi0 = kappa.size() + 1;
  return Phij(kappa(j-1), phi0, sum(exp(kappa)));
}