{-# LANGUAGE RankNTypes #-}
module Mcmc.Proposal.Slide
( slide,
slideSymmetric,
slideUniformSymmetric,
slideContrarily,
)
where
import Mcmc.Proposal
import Mcmc.Proposal.Generic
import Mcmc.Statistics.Types
import Statistics.Distribution.Normal
import Statistics.Distribution.Uniform
slidePFunction :: Mean Double -> StandardDeviation Double -> TuningParameter -> PFunction Double
slidePFunction :: Double -> Double -> Double -> PFunction Double
slidePFunction Double
m Double
s Double
t =
forall d a.
(ContDistr d, ContGen d) =>
d
-> (a -> Double -> a)
-> Maybe (Double -> Double)
-> Maybe (a -> Double -> Jacobian)
-> PFunction a
genericContinuous (Double -> Double -> NormalDistribution
normalDistr Double
m (Double
s forall a. Num a => a -> a -> a
* Double
t)) forall a. Num a => a -> a -> a
(+) (forall a. a -> Maybe a
Just forall a. Num a => a -> a
negate) forall a. Maybe a
Nothing
slide ::
Mean Double ->
StandardDeviation Double ->
PName ->
PWeight ->
Tune ->
Proposal Double
slide :: Double -> Double -> PName -> PWeight -> Tune -> Proposal Double
slide Double
m Double
s = forall a.
PDescription
-> (Double -> PFunction a)
-> PSpeed
-> PDimension
-> PName
-> PWeight
-> Tune
-> Proposal a
createProposal PDescription
description (Double -> Double -> Double -> PFunction Double
slidePFunction Double
m Double
s) PSpeed
PFast (Int -> PDimension
PDimension Int
1)
where
description :: PDescription
description = String -> PDescription
PDescription forall a b. (a -> b) -> a -> b
$ String
"Slide; mean: " forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> String
show Double
m forall a. [a] -> [a] -> [a]
++ String
", sd: " forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> String
show Double
s
slideSymmetricPFunction :: StandardDeviation Double -> TuningParameter -> PFunction Double
slideSymmetricPFunction :: Double -> Double -> PFunction Double
slideSymmetricPFunction Double
s Double
t =
forall d a.
(ContDistr d, ContGen d) =>
d
-> (a -> Double -> a)
-> Maybe (Double -> Double)
-> Maybe (a -> Double -> Jacobian)
-> PFunction a
genericContinuous (Double -> Double -> NormalDistribution
normalDistr Double
0.0 (Double
s forall a. Num a => a -> a -> a
* Double
t)) forall a. Num a => a -> a -> a
(+) forall a. Maybe a
Nothing forall a. Maybe a
Nothing
slideSymmetric ::
StandardDeviation Double ->
PName ->
PWeight ->
Tune ->
Proposal Double
slideSymmetric :: Double -> PName -> PWeight -> Tune -> Proposal Double
slideSymmetric Double
s = forall a.
PDescription
-> (Double -> PFunction a)
-> PSpeed
-> PDimension
-> PName
-> PWeight
-> Tune
-> Proposal a
createProposal PDescription
description (Double -> Double -> PFunction Double
slideSymmetricPFunction Double
s) PSpeed
PFast (Int -> PDimension
PDimension Int
1)
where
description :: PDescription
description = String -> PDescription
PDescription forall a b. (a -> b) -> a -> b
$ String
"Slide symmetric; sd: " forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> String
show Double
s
slideUniformPFunction :: Size -> TuningParameter -> PFunction Double
slideUniformPFunction :: Double -> Double -> PFunction Double
slideUniformPFunction Double
d Double
t =
forall d a.
(ContDistr d, ContGen d) =>
d
-> (a -> Double -> a)
-> Maybe (Double -> Double)
-> Maybe (a -> Double -> Jacobian)
-> PFunction a
genericContinuous (Double -> Double -> UniformDistribution
uniformDistr (-Double
t forall a. Num a => a -> a -> a
* Double
d) (Double
t forall a. Num a => a -> a -> a
* Double
d)) forall a. Num a => a -> a -> a
(+) forall a. Maybe a
Nothing forall a. Maybe a
Nothing
slideUniformSymmetric ::
Size ->
PName ->
PWeight ->
Tune ->
Proposal Double
slideUniformSymmetric :: Double -> PName -> PWeight -> Tune -> Proposal Double
slideUniformSymmetric Double
d = forall a.
PDescription
-> (Double -> PFunction a)
-> PSpeed
-> PDimension
-> PName
-> PWeight
-> Tune
-> Proposal a
createProposal PDescription
description (Double -> Double -> PFunction Double
slideUniformPFunction Double
d) PSpeed
PFast (Int -> PDimension
PDimension Int
1)
where
description :: PDescription
description = String -> PDescription
PDescription forall a b. (a -> b) -> a -> b
$ String
"Slide uniform symmetric; delta: " forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> String
show Double
d
contra :: (Double, Double) -> Double -> (Double, Double)
contra :: (Double, Double) -> Double -> (Double, Double)
contra (Double
x, Double
y) Double
u = (Double
x forall a. Num a => a -> a -> a
+ Double
u, Double
y forall a. Num a => a -> a -> a
- Double
u)
slideContrarilyPFunction ::
Mean Double ->
StandardDeviation Double ->
TuningParameter ->
PFunction (Double, Double)
slideContrarilyPFunction :: Double -> Double -> Double -> PFunction (Double, Double)
slideContrarilyPFunction Double
m Double
s Double
t =
forall d a.
(ContDistr d, ContGen d) =>
d
-> (a -> Double -> a)
-> Maybe (Double -> Double)
-> Maybe (a -> Double -> Jacobian)
-> PFunction a
genericContinuous (Double -> Double -> NormalDistribution
normalDistr Double
m (Double
s forall a. Num a => a -> a -> a
* Double
t)) (Double, Double) -> Double -> (Double, Double)
contra (forall a. a -> Maybe a
Just forall a. Num a => a -> a
negate) forall a. Maybe a
Nothing
slideContrarily ::
Mean Double ->
StandardDeviation Double ->
PName ->
PWeight ->
Tune ->
Proposal (Double, Double)
slideContrarily :: Double
-> Double -> PName -> PWeight -> Tune -> Proposal (Double, Double)
slideContrarily Double
m Double
s = forall a.
PDescription
-> (Double -> PFunction a)
-> PSpeed
-> PDimension
-> PName
-> PWeight
-> Tune
-> Proposal a
createProposal PDescription
description (Double -> Double -> Double -> PFunction (Double, Double)
slideContrarilyPFunction Double
m Double
s) PSpeed
PFast (Int -> PDimension
PDimension Int
2)
where
description :: PDescription
description = String -> PDescription
PDescription forall a b. (a -> b) -> a -> b
$ String
"Slide contrarily; mean: " forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> String
show Double
m forall a. [a] -> [a] -> [a]
++ String
", sd: " forall a. [a] -> [a] -> [a]
++ forall a. Show a => a -> String
show Double
s