Example #1
0
File: phi.cpp Project: 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;
}
Example #2
0
File: phi.cpp Project: Horta/BASiCS
Phij Phi::phij(size_t j) const
{
  double phi0 = kappa.size() + 1;
  return Phij(kappa(j-1), phi0, sum(exp(kappa)));
}