会社の新しいmac book pro に Julia をインストールする
自宅 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の間である。 この失得点差は、非常に僅差であると言えるだろう。ある意味人為的な操作があったとも考えられても不思議ではない。