{-# LANGUAGE DeriveDataTypeable #-}
{-# LANGUAGE DeriveGeneric #-}
module ELynx.Tree.Distribution.TimeOfOrigin
( TimeOfOriginDistribution (..),
cumulative,
density,
quantile,
)
where
import Data.Data
( Data,
Typeable,
)
import ELynx.Tree.Distribution.Types
import GHC.Generics (Generic)
import qualified Statistics.Distribution as D
data TimeOfOriginDistribution = TOD
{
TimeOfOriginDistribution -> Int
todTN :: Int,
TimeOfOriginDistribution -> Double
todLa :: Rate,
TimeOfOriginDistribution -> Double
todMu :: Rate
}
deriving (TimeOfOriginDistribution -> TimeOfOriginDistribution -> Bool
forall a. (a -> a -> Bool) -> (a -> a -> Bool) -> Eq a
/= :: TimeOfOriginDistribution -> TimeOfOriginDistribution -> Bool
$c/= :: TimeOfOriginDistribution -> TimeOfOriginDistribution -> Bool
== :: TimeOfOriginDistribution -> TimeOfOriginDistribution -> Bool
$c== :: TimeOfOriginDistribution -> TimeOfOriginDistribution -> Bool
Eq, Typeable, Typeable TimeOfOriginDistribution
TimeOfOriginDistribution -> DataType
TimeOfOriginDistribution -> Constr
(forall b. Data b => b -> b)
-> TimeOfOriginDistribution -> TimeOfOriginDistribution
forall a.
Typeable a
-> (forall (c :: * -> *).
(forall d b. Data d => c (d -> b) -> d -> c b)
-> (forall g. g -> c g) -> a -> c a)
-> (forall (c :: * -> *).
(forall b r. Data b => c (b -> r) -> c r)
-> (forall r. r -> c r) -> Constr -> c a)
-> (a -> Constr)
-> (a -> DataType)
-> (forall (t :: * -> *) (c :: * -> *).
Typeable t =>
(forall d. Data d => c (t d)) -> Maybe (c a))
-> (forall (t :: * -> * -> *) (c :: * -> *).
Typeable t =>
(forall d e. (Data d, Data e) => c (t d e)) -> Maybe (c a))
-> ((forall b. Data b => b -> b) -> a -> a)
-> (forall r r'.
(r -> r' -> r) -> r -> (forall d. Data d => d -> r') -> a -> r)
-> (forall r r'.
(r' -> r -> r) -> r -> (forall d. Data d => d -> r') -> a -> r)
-> (forall u. (forall d. Data d => d -> u) -> a -> [u])
-> (forall u. Int -> (forall d. Data d => d -> u) -> a -> u)
-> (forall (m :: * -> *).
Monad m =>
(forall d. Data d => d -> m d) -> a -> m a)
-> (forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d) -> a -> m a)
-> (forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d) -> a -> m a)
-> Data a
forall u.
Int
-> (forall d. Data d => d -> u) -> TimeOfOriginDistribution -> u
forall u.
(forall d. Data d => d -> u) -> TimeOfOriginDistribution -> [u]
forall r r'.
(r -> r' -> r)
-> r
-> (forall d. Data d => d -> r')
-> TimeOfOriginDistribution
-> r
forall r r'.
(r' -> r -> r)
-> r
-> (forall d. Data d => d -> r')
-> TimeOfOriginDistribution
-> r
forall (m :: * -> *).
Monad m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
forall (c :: * -> *).
(forall b r. Data b => c (b -> r) -> c r)
-> (forall r. r -> c r) -> Constr -> c TimeOfOriginDistribution
forall (c :: * -> *).
(forall d b. Data d => c (d -> b) -> d -> c b)
-> (forall g. g -> c g)
-> TimeOfOriginDistribution
-> c TimeOfOriginDistribution
forall (t :: * -> *) (c :: * -> *).
Typeable t =>
(forall d. Data d => c (t d)) -> Maybe (c TimeOfOriginDistribution)
forall (t :: * -> * -> *) (c :: * -> *).
Typeable t =>
(forall d e. (Data d, Data e) => c (t d e))
-> Maybe (c TimeOfOriginDistribution)
gmapMo :: forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
$cgmapMo :: forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
gmapMp :: forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
$cgmapMp :: forall (m :: * -> *).
MonadPlus m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
gmapM :: forall (m :: * -> *).
Monad m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
$cgmapM :: forall (m :: * -> *).
Monad m =>
(forall d. Data d => d -> m d)
-> TimeOfOriginDistribution -> m TimeOfOriginDistribution
gmapQi :: forall u.
Int
-> (forall d. Data d => d -> u) -> TimeOfOriginDistribution -> u
$cgmapQi :: forall u.
Int
-> (forall d. Data d => d -> u) -> TimeOfOriginDistribution -> u
gmapQ :: forall u.
(forall d. Data d => d -> u) -> TimeOfOriginDistribution -> [u]
$cgmapQ :: forall u.
(forall d. Data d => d -> u) -> TimeOfOriginDistribution -> [u]
gmapQr :: forall r r'.
(r' -> r -> r)
-> r
-> (forall d. Data d => d -> r')
-> TimeOfOriginDistribution
-> r
$cgmapQr :: forall r r'.
(r' -> r -> r)
-> r
-> (forall d. Data d => d -> r')
-> TimeOfOriginDistribution
-> r
gmapQl :: forall r r'.
(r -> r' -> r)
-> r
-> (forall d. Data d => d -> r')
-> TimeOfOriginDistribution
-> r
$cgmapQl :: forall r r'.
(r -> r' -> r)
-> r
-> (forall d. Data d => d -> r')
-> TimeOfOriginDistribution
-> r
gmapT :: (forall b. Data b => b -> b)
-> TimeOfOriginDistribution -> TimeOfOriginDistribution
$cgmapT :: (forall b. Data b => b -> b)
-> TimeOfOriginDistribution -> TimeOfOriginDistribution
dataCast2 :: forall (t :: * -> * -> *) (c :: * -> *).
Typeable t =>
(forall d e. (Data d, Data e) => c (t d e))
-> Maybe (c TimeOfOriginDistribution)
$cdataCast2 :: forall (t :: * -> * -> *) (c :: * -> *).
Typeable t =>
(forall d e. (Data d, Data e) => c (t d e))
-> Maybe (c TimeOfOriginDistribution)
dataCast1 :: forall (t :: * -> *) (c :: * -> *).
Typeable t =>
(forall d. Data d => c (t d)) -> Maybe (c TimeOfOriginDistribution)
$cdataCast1 :: forall (t :: * -> *) (c :: * -> *).
Typeable t =>
(forall d. Data d => c (t d)) -> Maybe (c TimeOfOriginDistribution)
dataTypeOf :: TimeOfOriginDistribution -> DataType
$cdataTypeOf :: TimeOfOriginDistribution -> DataType
toConstr :: TimeOfOriginDistribution -> Constr
$ctoConstr :: TimeOfOriginDistribution -> Constr
gunfold :: forall (c :: * -> *).
(forall b r. Data b => c (b -> r) -> c r)
-> (forall r. r -> c r) -> Constr -> c TimeOfOriginDistribution
$cgunfold :: forall (c :: * -> *).
(forall b r. Data b => c (b -> r) -> c r)
-> (forall r. r -> c r) -> Constr -> c TimeOfOriginDistribution
gfoldl :: forall (c :: * -> *).
(forall d b. Data d => c (d -> b) -> d -> c b)
-> (forall g. g -> c g)
-> TimeOfOriginDistribution
-> c TimeOfOriginDistribution
$cgfoldl :: forall (c :: * -> *).
(forall d b. Data d => c (d -> b) -> d -> c b)
-> (forall g. g -> c g)
-> TimeOfOriginDistribution
-> c TimeOfOriginDistribution
Data, forall x.
Rep TimeOfOriginDistribution x -> TimeOfOriginDistribution
forall x.
TimeOfOriginDistribution -> Rep TimeOfOriginDistribution x
forall a.
(forall x. a -> Rep a x) -> (forall x. Rep a x -> a) -> Generic a
$cto :: forall x.
Rep TimeOfOriginDistribution x -> TimeOfOriginDistribution
$cfrom :: forall x.
TimeOfOriginDistribution -> Rep TimeOfOriginDistribution x
Generic)
instance D.Distribution TimeOfOriginDistribution where
cumulative :: TimeOfOriginDistribution -> Double -> Double
cumulative = TimeOfOriginDistribution -> Double -> Double
cumulative
cumulative :: TimeOfOriginDistribution -> Time -> Double
cumulative :: TimeOfOriginDistribution -> Double -> Double
cumulative (TOD Int
n Double
l Double
m) Double
x
| Double
x forall a. Ord a => a -> a -> Bool
<= Double
0 = Double
0
| Bool
otherwise = Double
te forall a. Floating a => a -> a -> a
** forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
n
where
d :: Double
d = Double
l forall a. Num a => a -> a -> a
- Double
m
te :: Double
te = Double
l forall a. Num a => a -> a -> a
* (Double
1.0 forall a. Num a => a -> a -> a
- forall a. Floating a => a -> a
exp (-Double
d forall a. Num a => a -> a -> a
* Double
x)) forall a. Fractional a => a -> a -> a
/ (Double
l forall a. Num a => a -> a -> a
- Double
m forall a. Num a => a -> a -> a
* forall a. Floating a => a -> a
exp (-Double
d forall a. Num a => a -> a -> a
* Double
x))
instance D.ContDistr TimeOfOriginDistribution where
density :: TimeOfOriginDistribution -> Double -> Double
density = TimeOfOriginDistribution -> Double -> Double
density
quantile :: TimeOfOriginDistribution -> Double -> Double
quantile = TimeOfOriginDistribution -> Double -> Double
quantile
density :: TimeOfOriginDistribution -> Time -> Double
density :: TimeOfOriginDistribution -> Double -> Double
density (TOD Int
nn Double
l Double
m) Double
x
| Double
x forall a. Ord a => a -> a -> Bool
< Double
0 = Double
0
| Bool
otherwise = Double
n forall a. Num a => a -> a -> a
* Double
l forall a. Floating a => a -> a -> a
** Double
n forall a. Num a => a -> a -> a
* Double
d forall a. Floating a => a -> a -> a
** Double
2 forall a. Num a => a -> a -> a
* Double
t1 forall a. Floating a => a -> a -> a
** (Double
n forall a. Num a => a -> a -> a
- Double
1.0) forall a. Num a => a -> a -> a
* Double
ex forall a. Fractional a => a -> a -> a
/ (Double
t2 forall a. Floating a => a -> a -> a
** (Double
n forall a. Num a => a -> a -> a
+ Double
1.0))
where
d :: Double
d = Double
l forall a. Num a => a -> a -> a
- Double
m
n :: Double
n = forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
nn
ex :: Double
ex = forall a. Floating a => a -> a
exp (-Double
d forall a. Num a => a -> a -> a
* Double
x)
t1 :: Double
t1 = Double
1.0 forall a. Num a => a -> a -> a
- Double
ex
t2 :: Double
t2 = Double
l forall a. Num a => a -> a -> a
- Double
m forall a. Num a => a -> a -> a
* Double
ex
quantile :: TimeOfOriginDistribution -> Double -> Time
quantile :: TimeOfOriginDistribution -> Double -> Double
quantile (TOD Int
n' Double
l Double
m) Double
p
| Double
p forall a. Ord a => a -> a -> Bool
>= Double
0 Bool -> Bool -> Bool
&& Double
p forall a. Ord a => a -> a -> Bool
<= Double
1 =
-Double
1.0 forall a. Fractional a => a -> a -> a
/ Double
d forall a. Num a => a -> a -> a
* forall a. Floating a => a -> a
log (Double
t1 forall a. Fractional a => a -> a -> a
/ Double
t2)
| Bool
otherwise =
forall a. HasCallStack => [Char] -> a
error forall a b. (a -> b) -> a -> b
$
[Char]
"PointProcess.quantile: p must be in [0,1] range. Got: "
forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> [Char]
show Double
p
forall a. [a] -> [a] -> [a]
++ [Char]
"."
where
d :: Double
d = Double
l forall a. Num a => a -> a -> a
- Double
m
n :: Double
n = forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
n'
t1 :: Double
t1 = Double
l forall a. Num a => a -> a -> a
* (Double
1.0 forall a. Num a => a -> a -> a
- Double
p forall a. Floating a => a -> a -> a
** (Double
1.0 forall a. Fractional a => a -> a -> a
/ Double
n))
t2 :: Double
t2 = Double
l forall a. Num a => a -> a -> a
- Double
p forall a. Floating a => a -> a -> a
** (Double
1.0 forall a. Fractional a => a -> a -> a
/ Double
n) forall a. Num a => a -> a -> a
* Double
m
instance D.ContGen TimeOfOriginDistribution where
genContVar :: forall g (m :: * -> *).
StatefulGen g m =>
TimeOfOriginDistribution -> g -> m Double
genContVar = forall d g (m :: * -> *).
(ContDistr d, StatefulGen g m) =>
d -> g -> m Double
D.genContinuous