{
"cells": [
{
"cell_type": "code",
"execution_count": 75,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"g (generic function with 1 method)"
]
},
"execution_count": 75,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# コーシー分布\n",
"f(y) = 1/(π*(1+y^2))\n",
"# 正規分布\n",
"\n",
"g(σ,x,μ) = 1/(√(2π*σ^2))*exp((x-μ)^2/2σ^2)"
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {},
"outputs": [],
"source": [
"using Plots"
]
},
{
"cell_type": "code",
"execution_count": 105,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
" \n",
" \n"
]
},
"execution_count": 105,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"lables = Array{String}(1,3)\n",
"\n",
"# コーシー分布\n",
"X = Float64[]\n",
"for i in -5:.1:5\n",
" push!(X,f(i))\n",
"end\n",
"#plot(X)\n",
"x = linspace(-5, 5, 100)\n",
"d = Normal(0, 1)\n",
"y = pdf(d, x)\n",
"labels[1] = string(\"コーシー分布\")\n",
"labels[2] = string(\"正規分布\")\n",
"\n",
"plot(x, [X,y], linewidth=2, alpha=0.6, label=labels)"
]
},
{
"cell_type": "code",
"execution_count": 109,
"metadata": {},
"outputs": [
{
"ename": "LoadError",
"evalue": "\u001b[91mMethodError: no method matching g(::Array{Float64,1})\u001b[0m\nClosest candidates are:\n g(::Any, \u001b[91m::Any\u001b[39m, \u001b[91m::Any\u001b[39m) at In[75]:5\u001b[39m",
"output_type": "error",
"traceback": [
"\u001b[91mMethodError: no method matching g(::Array{Float64,1})\u001b[0m\nClosest candidates are:\n g(::Any, \u001b[91m::Any\u001b[39m, \u001b[91m::Any\u001b[39m) at In[75]:5\u001b[39m",
"",
"Stacktrace:",
" [1] \u001b[1minclude_string\u001b[22m\u001b[22m\u001b[1m(\u001b[22m\u001b[22m::String, ::String\u001b[1m)\u001b[22m\u001b[22m at \u001b[1m./loading.jl:522\u001b[22m\u001b[22m"
]
}
],
"source": [
"log(g(y)/f(y))g(y)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 0.6.3",
"language": "julia",
"name": "julia-0.6"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "0.6.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}