The most complete description of the state of a system at any time is given by knowledge of the probability density function, which describes the locus of possible states conditioned on any available measurement information. When employing radar returns, an admissible region approach provides a physics-based region of the right-ascension rate/declination rate space of possible Earth-bound orbit solutions. This work develops a method that employs a probabilistic interpretation of the admissible region and approximates the admissible region by a Gaussian mixture to formulate an initial orbit determination solution