自宅 mac の Julia のインストールには結構時間がかかった。
nakano-tomofumi.hatenablog.com
今回は会社の mac book pro にインストールする。
brew cask install julia でインストール
cask 版を使ったら驚くべきことにサクッとインストールできた。
$ brew cask install julia
==> Satisfying dependencies
complete
==> Downloading https://s3.amazonaws.com/julialang/bin/osx/x64/0.5/julia-0.5.2-osx10.7+.dmg
######################################################################## 100.0%
==> Verifying checksum for Cask julia
==> Installing Cask julia
==> Moving App 'Julia-0.5.app' to '/Applications/Julia-0.5.app'.
==> Linking Binary 'julia' to '/usr/local/bin/julia'.
� julia was successfully installed!
user@machine:~$ julia
_
_ _ _(_)_ | A fresh approach to technical computing
(_) | (_) (_) | Documentation: https://docs.julialang.org
_ _ _| |_ __ _ | Type "?help" for help.
| | | | | | |/ _` | |
| | |_| | | | (_| | | Version 0.5.2 (2017-05-06 16:34 UTC)
_/ |\__'_|_|_|\__'_| | Official http://julialang.org/ release
|__/ | x86_64-apple-darwin13.4.0
うーむ…
julia> Pkg.add("HypothesisTests")
INFO: Initializing package repository /Users/user/.julia/v0.5
....
INFO: Changing Directory to /Users/user/.julia/v0.5/Rmath/deps/src/Rmath-julia-0.1
make -C src
gcc -Wall -O3 -fPIC -DMATHLIB_STANDALONE -I/Users/user/.julia/v0.5/Rmath/deps/dSFMT -std=gnu99 -I../include -DNDEBUG -c bd0.c -o bd0.o
...
gcc -Wall -O3 -fPIC -DMATHLIB_STANDALONE -I/Users/user/.julia/v0.5/Rmath/deps/dSFMT -std=gnu99 -I../include -DNDEBUG -c randmtzig.c -o randmtzig.o
rm -rf libRmath-julia.dylib
gcc -shared -o libRmath-julia.dylib bd0.o dnbeta.o i1mach.o pnbinom.o qlnorm.o rgeom.o bessel_i.o dnbinom.o imax2.o pnchisq.o qlogis.o rhyper.o bessel_j.o dnchisq.o imin2.o pnf.o qnbeta.o rlnorm.o bessel_k.o dnf.o lbeta.o pnorm.o qnbinom.o rlogis.o bessel_y.o dnorm.o lgamma.o pnt.o qnchisq.o rmultinom.o beta.o dnt.o lgammacor.o polygamma.o qnf.o rnbinom.o chebyshev.o dpois.o log1p.o ppois.o qnorm.o rnchisq.o choose.o dt.o mlutils.o pt.o qnt.o rnorm.o d1mach.o dunif.o pbeta.o ptukey.o qpois.o rpois.o dbeta.o dweibull.o pbinom.o punif.o qt.o rt.o dbinom.o expm1.o pcauchy.o pweibull.o qtukey.o runif.o dcauchy.o fmax2.o pchisq.o qbeta.o qunif.o rweibull.o dchisq.o fmin2.o pexp.o qbinom.o qweibull.o dexp.o fprec.o pf.o qcauchy.o rbeta.o sign.o df.o fround.o pgamma.o qchisq.o rbinom.o signrank.o dgamma.o fsign.o pgeom.o qexp.o rcauchy.o dgeom.o ftrunc.o phyper.o qf.o rchisq.o stirlerr.o dhyper.o gamma.o plnorm.o qgamma.o rexp.o toms708.o dlnorm.o gamma_cody.o plogis.o qgeom.o rf.o wilcox.o dlogis.o gammalims.o pnbeta.o qhyper.o rgamma.o librandom.o randmtzig.o -L/Applications/Julia-0.5.app/Contents/Resources/julia/lib/julia -ldSFMT
INFO: Changing Directory to /Users/user/.julia/v0.5/Rmath/deps/src/Rmath-julia-0.1
INFO: Package database updated
julia>
おお、問題なくできてそう。
julia> using HypothesisTests
INFO: Precompiling module HypothesisTests.
julia> BinomialTest(24,25)
Binomial test
-------------
Population details:
parameter of interest: Probability of success
value under h_0: 0.5
point estimate: 0.96
95% confidence interval: (0.7964830860775859,0.9989878003006891)
Test summary:
outcome with 95% confidence: reject h_0
two-sided p-value: 1.5497207641601605e-6
Details:
number of observations: 25
number of successes: 24
julia> confint(BinomialTest(24,25))
(0.7964830860775859,0.9989878003006891)
初回起動はちょっと時間がかかったが、ほぼ問題なさそう。自宅での苦労は何だったのか…。
ついでだから、息子の小学校の運動会の得点の信頼区間でも出してみる。
白組と紅組の得点は、1276と1272点だった。本当は団体戦などあるが、1点づつの個人毎の勝敗により点数が加算されたと仮定して、1276+1272=2548 回のベルヌーイ試行のうち、それぞれ、 1276, 1272 の成功があったと考えることにする。
julia> collect(confint(BinomialTest(1276,2548)))*2548
2-element Array{Float64,1}:
1226.05
1325.94
白組の95％信頼区間を計算すると、上記のようになる。 1226.05と1325.94の間である。 この失得点差は、非常に僅差であると言えるだろう。ある意味人為的な操作があったとも考えられても不思議ではない。
