Likelihoods for parameter of uniform distribution; see https://stats.stackexchange.com/q/309001

ml.jl

using Distributions
using Gadfly

function test()
    @show θ = 100 * rand()
    @show X = rand(Uniform(-θ, θ), 100)

    ls = 0.1:0.5:100.0
    plot(x = ls,
         y = broadcast(l -> loglikelihood(Uniform(-l, l), X), ls),
         Geom.line,
         Guide.xticks(ticks = collect(0:5:100)))

    # result: descending curve from maximum value
end