Fitting beta distributions in Scala 😉

beta_distribution_fit.scala

import scala.sys.process._

object BetaDistributionFit {

  val distName: String = "beta"

  def fitCommand(samples: Seq[Double]): Seq[String] =
    Seq("python", "-c", 
        s"""
        |from scipy import stats
        |samples = [${samples.mkString(",")}]
        |print(','.join(map(str, stats.$distName.fit(samples))))
        """.stripMargin)

  def fit(samples: Seq[Double]): Seq[Double] =
    fitCommand(samples).!!.split(",").map(_.toDouble)

}

More monads, please.