Skip to content

Commit 5db7610

Browse files
authored
Update README.md
1 parent 9361181 commit 5db7610

File tree

1 file changed

+3
-8
lines changed

1 file changed

+3
-8
lines changed

README.md

Lines changed: 3 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -106,7 +106,7 @@ plot!(x, [target envelop], width = 2, label = ["target density" "envelop"])
106106

107107
![](img/example3.png)
108108

109-
### Elastic Net Eistribution
109+
### Elastic Net Distribution
110110

111111
The following example arises from elastic net regression and smoothing problems. In these cases, the integration constants are not available analytically.
112112

@@ -135,7 +135,7 @@ plot!(x, [target envelop], width = 2, label = ["target density" "envelop"])
135135

136136
![](img/example4.png)
137137

138-
### Tips for numerical stability
138+
### Tips for numerical stability, use of logdensity
139139

140140
Here are some tips:
141141

@@ -162,12 +162,7 @@ theta = 1.0
162162
# a complicated logdensity
163163
logf(v) = n * v - (n - k * alpha) * logsumexp([v, log(tau)]) - theta / alpha * ( (tau + exp(v) )^alpha )
164164

165-
# make two plots of logf and f
166-
p1 = plot(logf, -20, 20, label = "logf")
167-
p2 = plot(x -> exp(logf(x)), -20, 20, label = "f")
168-
plot(p1, p2, layout = (1, 2))
169-
170-
# runs sampler
165+
# run sampler
171166
δ = 0.1
172167
support = (-Inf, Inf)
173168
search = (0.0, 10.0)

0 commit comments

Comments
 (0)