|
24 | 24 | integral, _ = quadgk(t -> pdf(dist, t), 0.130, 20.0) |
25 | 25 | @test integral ≈ 1 atol = 1e-4 |
26 | 26 |
|
27 | | - dist = Wald(; ν = 2.0, η = 0.0, α = 0.3, τ = .130) |
28 | | - integral, _ = quadgk(t -> pdf(dist, t), .130, 20.0) |
| 27 | + dist = Wald(; ν = 2.0, η = 0.0, α = 0.3, τ = 0.130) |
| 28 | + integral, _ = quadgk(t -> pdf(dist, t), 0.130, 20.0) |
29 | 29 | @test integral ≈ 1 atol = 1e-4 |
30 | 30 | end |
31 | 31 | end |
|
63 | 63 | using StableRNGs |
64 | 64 |
|
65 | 65 | rng = StableRNG(34) |
66 | | - |
| 66 | + |
67 | 67 | Θ = (; ν = 2.0, η = 0.3, α = 0.3, τ = 0.130) |
68 | 68 | dist = Wald(; Θ...) |
69 | 69 | data = rand(rng, dist, 10_000) |
70 | 70 |
|
71 | | - νs = range(Θ.ν * .8, Θ.ν * 1.2, 100) |
| 71 | + νs = range(Θ.ν * 0.8, Θ.ν * 1.2, 100) |
72 | 72 | LLs = map(ν -> loglikelihood(Wald(; Θ..., ν), data), νs) |
73 | 73 | _, max_idx = findmax(LLs) |
74 | | - @test νs[max_idx] ≈ Θ.ν rtol = .02 |
| 74 | + @test νs[max_idx] ≈ Θ.ν rtol = 0.02 |
75 | 75 |
|
76 | 76 | ηs = range(0, Θ.η * 1.2, 100) |
77 | 77 | LLs = map(η -> loglikelihood(Wald(; Θ..., η), data), ηs) |
78 | 78 | _, max_idx = findmax(LLs) |
79 | | - @test ηs[max_idx] ≈ Θ.η atol = .05 |
| 79 | + @test ηs[max_idx] ≈ Θ.η atol = 0.05 |
80 | 80 |
|
81 | | - αs = range(Θ.α * .8, Θ.α * 1.2, 100) |
| 81 | + αs = range(Θ.α * 0.8, Θ.α * 1.2, 100) |
82 | 82 | LLs = map(α -> loglikelihood(Wald(; Θ..., α), data), αs) |
83 | 83 | _, max_idx = findmax(LLs) |
84 | | - @test αs[max_idx] ≈ Θ.α rtol = .02 |
| 84 | + @test αs[max_idx] ≈ Θ.α rtol = 0.02 |
85 | 85 |
|
86 | | - τs = range(Θ.τ * .8, Θ.τ, 100) |
| 86 | + τs = range(Θ.τ * 0.8, Θ.τ, 100) |
87 | 87 | LLs = map(τ -> loglikelihood(Wald(; Θ..., τ), data), τs) |
88 | 88 | _, max_idx = findmax(LLs) |
89 | | - @test τs[max_idx] ≈ Θ.τ rtol = .02 |
| 89 | + @test τs[max_idx] ≈ Θ.τ rtol = 0.02 |
90 | 90 | end |
91 | 91 |
|
92 | 92 | @safetestset "simulate" begin |
|
0 commit comments