theta_pv: Sampling plausible values(PVs) based on rejection sampling.
In takuizum/irtfun2: Some useful functions for IRT

Sampling plausible values(PVs) based on rejection sampling.

theta_pv(
  x,
  nofrands,
  eap_apply,
  const_apply,
  map_apply,
  n,
  maxtheta,
  mintheta,
  a,
  b,
  c,
  D,
  mu,
  sigma
)

`x`	response data
`nofrands`	the number of PVs
`eap_apply`	a numeric vector of estimated EAP.
`const_apply`	a numeric vector of const of a posterior distribution.
`map_apply`	a numeric vector of estimated MAP.
`n`	the number of subjects.
`maxtheta`	a max value of theta.
`mintheta`	a mininum value of theta.
`a`	an item slope parameter vector.
`b`	an item location parameter vector.
`c`	an item asymptote parameter vector.
`D`	a factor constant.
`mu`	p hyper parameter of a posterior distribution.
`sigma`	same as above.