黒木玄 Gen Kuroki(@genkuroki)/2021年01月22日

項目反応理論の連続型モデルの改良と推定方法に関する論文のオンライン版が公開されました。
実装はJulia言語でやってます。

www.jstage.jst.go.jp/article/jjep/6...

タグ：

posted at 01:31:14

しぶてぃ @takuizum

RからもJuliaCallパッケージを使えば動くと思いますが，試してません。

タグ：

posted at 01:31:35

中学校のテスト問題の例で「定義」と答えさせたい問題↓

【漢字で習ったものは漢字で書きなさい。

(1)ことばが指し示すものが何であるかを，はっきりと簡潔に述べたものを，そのことばの[　　　]という。】

これだと日常用語としての「定義」でしかないと思う。

定義は簡潔である必要はない。 twitter.com/genkuroki/stat...

タグ：

posted at 04:24:30

【(5)変数χとyの関係が、y=a/χで表されるとき、yはχに
[　⑧　]するという。このとき、aを[　⑤　]という。

(6)反比例のグラフは、[　⑨　]である。】

⑧が「反比例」だと即答できる人は多いと思いますが、ここに引用した部分だけを見て、⑤と⑨が何であるかを答えられますか？

私はできません。 twitter.com/genkuroki/stat... pic.twitter.com/GUiHyndZnH

タグ：

posted at 04:29:16

漢字で書けることが【知識理解】と関係あるらしい。

中学生の保護者は知っていると思いますが、こういうテスト問題が本当に中学校の数学のテストで出されているんです。

こういうのは即刻やめさせないと、大島さんが指摘しているような問題はさらに悪化するでしょう。 pic.twitter.com/SsPRTzJm2l

タグ：

posted at 04:32:14

この問題については高等教育機関側にも責任があって、大島さんが指摘している問題に気付いた人は、講義の中で(特に教員免許取得に必要な講義の中で)、現実にすでに行われている教え方について批判的な考え方をはっきり述べるべきだと思います。

逆に肯定するようだと非常にまずい。

タグ：

posted at 04:35:02

【(6)反比例のグラフは、[　⑨　]である。】という問題について、

* 漢字で書かなくてもよい
* 反比例のグラフの正しい性質なら何でもよい
* 正解は無数にある
* 本質的に異なる答えを3通り以上書かせる

などすれば良い問題になる可能性もありますが、現実には答えを暗記して漢字で書かないとバツ。

タグ：

posted at 04:42:23

数学なのだから、論理的に正しかったり、文脈的に数学的に正しいことを言っていると解釈可能な答案にはすべてマルを付けなければいけません。

正解が無数にあるならそれらすべてに丸を付けなければいけない。

小学生のときからずっとそういう教育をしないと、数学的に正常な理解力は育ち難いと思う。

タグ：

posted at 04:45:59

#Julia言語おもしろすぎ！

少し待っていれば太陽系全体まで拡張したバージョンが出て来るのだろうか😊 twitter.com/asteroid_maru/...

タグ： Julia言語

posted at 05:54:19

M. Kanamaru @AstroshaperJP

@genkuroki ありがとうございます。

惑星の重力による摂動を入れて、小惑星の軌道を計算したいと思っています！粒が少ないので、全粒子ペアに対して愚直に重力を計算してもよさそうです。積分器はもう少しマシなものを書く必要があります。

タグ：

posted at 06:03:46

jkrumbiegel.github.io/pages/2020-10-...
Julia Helps To Bridge The Gap Between User and Creator

↑これは私が知る限りにおいて最も優れた #Julia言語の紹介文。

Juliaを使えば今まで単なるユーザーに過ぎなかった人達が開発者側に変貌する可能性がある！ twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 06:04:33

ほむほむ@アレルギー専門医：noteメン @ped_allergy

#Julia言語

「速く計算できて欲しい」という要求がある場合は多い。

速さが必要な部分をC, C++, Fortranなどで書かなければいけないことにされてしまうと、一般ユーザー側が手を出すことは困難になる。

Juliaを使えば、高速計算が必要な部分についてもユーザー側が開発者になれる！

タグ： Julia言語

posted at 06:04:33

ほむほむ@アレルギー専門医：noteメン @ped_allergy

【定期】mRNAワクチンは理論上、その人のDNAを修飾したり変えたりはしません。
しかし喫煙は、DNAを修飾して(メチル化といいます)、世代をこえ喘息の発症リスクを伝える可能性が示唆されています。
pediatric-allergy.com/2017/10/27/gra...

タグ：

posted at 07:56:54

DNAを修飾というと分かりにくいかもしれません。
たとえば、ある人そのものが変わっていなくてもメガネをしたり帽子をかぶったりすると印象が変わったりしますよね。
そんなDNA飾りをつけるような変化を『エピジェネティックな変化』といいます。
喫煙は、そんな変化を起こすということです😌

タグ：

posted at 08:01:14

非公開

タグ：

posted at xx:xx:xx

金子洋一神奈川20区（相模原市南区、座間 @Y_Kaneko

本の紹介とは別に，
記事の著者は Julia を使っているらしい．いいぞぉ〜．

計算機統計学・ベイズ統計学周辺でのお勧めの教科書10冊 qiita.com/s-yonekura/ite... #Qiita

タグ： Qiita

posted at 09:12:36

物価は経済の体温計。デフレに逆戻りしてしまったわが国経済。数十兆円規模の新たな財政支出と金融緩和をフルアクセルで取り組む必要がある。⇨去年12月の消費者物価指数前年同月比マイナス1.0％と大幅下落　www3.nhk.or.jp/news/html/2021...

タグ：

posted at 09:13:17

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
QuadGK, Optim, Roots, Memoize, NLsolve の使用例

QuadGK.jl 1変数函数の数値積分
Optim.jl 函数の最小化
Roots.jl 1変数函数の零点
Memoize.jl メモ化
NLsolve.jl 非線形(連立)方程式の解法

この手の基本的なパッケージの使い方を知っていると楽をできます。

タグ： Julia言語

posted at 09:33:57

#Julia言語

∫_0^1 ζ(s, x+t) dt = 1/((s-1)xˢ⁻¹) の数値的確認.

数値積分に QuadGK.quadgk を使用。私の場合には、自前で数値積分のコードを書くより、楽で同精度なら速いことが多いです。 pic.twitter.com/UpnRKTppuj

タグ： Julia言語

posted at 09:37:38

非公開

タグ：

posted at xx:xx:xx

#Julia言語

Optim.optimize を使ったベータ分布の最尤法

Distributions.jl の fit(Beta, sample) は最尤法でない推定法を採用しており、グラフを見れば分かるように誤差が大きくなる場合があります。 pic.twitter.com/JNC95VLc8X

タグ： Julia言語

posted at 09:40:02

#Julia言語

2×2の分割表のFisher検定の実装例

Roots.find_zeros を使って、Fisher検定と整合的なオッズ比の信頼区間を求めています。

Rのfisher.testが表示する信頼区間はP値と整合性がありません。 pic.twitter.com/0bNzzq160i

タグ： Julia言語

posted at 09:43:27

非公開

タグ：

posted at xx:xx:xx

#Julia言語

①与えられた平均と分散を持つWeibull分布をNLsolve.nlsolveを使って計算

②Weibull分布の最尤法をOptim.optimizeで実装

連立方程式をNLsolve.jlで解いたり、函数の最小化をOptim.jlで行ったりできます。もちろん、慎重な注意が必要な場合があります。 pic.twitter.com/WswYGX4bPy

タグ： Julia言語

posted at 09:48:06

めちゃくちゃ怖い。

タグ：

posted at 10:00:14

Yuki Nagai @cometscome_phys

Juliaでパッケージのcommit ハッシュ表示させるにはどうやればいいんだろ。文字列として取得したい

タグ：

posted at 10:06:08

非公開

タグ：

posted at xx:xx:xx

Termux上でJupyterを走らせてJuliaを書こうとしています。JuliaとJupyterのインストールにそれぞれ成功しました。しかし、IJuliaのプリコンパイルに失敗するためJupyterがJuliaカーネルを起動できません。困っています。#julialang #termux pic.twitter.com/holzSTE4Zh

タグ： julialang termux

posted at 11:49:01

I'm trying to code Julia by running Jupyter on Termux. I have successfully installed Julia and Jupyter respectively. However, Jupyter cannot boot the Julia kernel because IJulia precompile fails. How can I solve this? #Julialang #termux

タグ： Julialang termux

posted at 11:51:32

Python on Jupyter notebook and Julia itself can be run but precompiling IJulia faile inside julia. pic.twitter.com/lQNuoh8x9z

タグ：

posted at 11:59:46

M. Kanamaru @AstroshaperJP

３体問題できた（かわいい）

#Julia言語
 #julialang twitter.com/asteroid_maru/... pic.twitter.com/DdnmqoCLZN

タグ： julialang Julia言語

posted at 12:04:06

Dr. Tad @tak53381102

西浦先生らから

”「Go To Travel」前と期間中の旅行関連症例と観光関連症例の発生率を比較、観光起因COVID-19発生率は、対照の6/22-7/21, 7/15-7/19と比較し、それぞれ約8倍、2～3倍であった。国内観光客増加が旅行関連COVID-19症例増加に寄与したと考えられる”

www.mdpi.com/2077-0383/10/3...

タグ：

posted at 12:07:11

#Julia言語私のJulia v1.5.3では

pkg> st -m IJulia Conda ZMQ ZeroMQ_jll

の結果が添付画像のようになっています。

バージョンが違っていたら、例えば、

pkg> add ZMQ@1.2.1

のようにすれば(他のパッケージについても同様)、問題が解決する可能性があると思いました。 twitter.com/OPPO89694572/s... pic.twitter.com/sqEVYHNYfi

タグ： Julia言語

posted at 12:18:02

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

@genkuroki 早速の反応ありがとうございます。添付いただいたバージョンに揃えたあとビルド、using IJuliaしてみましたが結果変わらず。ZMQの初期化に失敗しているようですが、単体のusing ZMQは成功します。 pic.twitter.com/QaaBjg1VYN

タグ：

posted at 12:39:37

非公開

タグ：

posted at xx:xx:xx

みんなJulia使っとるな

タグ：

posted at 15:12:09

ceptree @ceptree

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

Juliaでやってみようという話は色々聞くので2-3年後どうなってるか楽しみじゃな

タグ：

posted at 15:13:09

Juliaが流行ればワシの本も100万部くらい売れそうやな

タグ：

posted at 15:16:00

ceptree @ceptree

分野で一番使われてるライブラリがJuliaって話を聞くのはいつだろう(流石にまだなさそう

タグ：

posted at 15:22:34

#Julia言語

計算結果をFisherText型のオブジェクトに格納するようにしてある。

添付画像

①FisherText型のオブジェクトのデフォルトでのテキスト表示とグラフ作画の仕方の設定

②Fisher検定の結果resultが自動的に清書されて表示され、plot(result)だけでグラフも表示

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/VexCcroUoo

タグ： Julia言語

posted at 15:42:46

#Julia言語

自分で定義した型のデフォルトでのテキスト表示の仕方の設定については以下のリンク先を参照。これ、非常に便利です。

docs.julialang.org/en/v1/manual/t...

docs.julialang.org/en/v1/base/io-...

タグ： Julia言語

posted at 15:47:36

#Julia言語

自分で定義した型の Plots.jl におけるデフォルトでのプロットの仕方の設定については以下のリンク先を参照。

これも非常に便利です。この機能があるので、PyPlot.jl とは別に Plots.jl も使いたくなる。

docs.juliaplots.org/latest/recipes/

タグ： Julia言語

posted at 15:50:22

訂正

❌FisherText
⭕️FisherTest twitter.com/genkuroki/stat...

タグ：

posted at 16:00:49

@OPPO89694572 #Julia言語私だったら、こういう場合には、「すべてを入れ直す」を実行します。

① ~/.julia を ~/.julia-bak にリネーム
② 最新安定板を入れる(入れ直す)
③ パッケージを全部入れ直す

これ、私は、多分、何十回もやっています。

タグ： Julia言語

posted at 18:35:29

非公開

タグ：

posted at xx:xx:xx

@genkuroki ありがとうございます。再インストールも試したのですが変わりませんでした。代わりにTermuxの中にUbuntu in termuxというUbuntu環境を入れて試しているところです。こちらではTermuxでうまく行かなかったPlutoも実行できたので望みがあります。

タグ：

posted at 20:11:47

白石祐太 @iasadafa

この解析、Rとpythonとjuliaとどれを使うべきかめっちゃ悩むな...

タグ：

posted at 20:19:12

```julia
using Franklin; newsite("$x")
```

をして git add -A && git commit -am "initial files"
git push したら簡単に

<自分のWebページ>/github.io でデプロイできてわろた・・・．数年前の苦労は何だったんだって気持ち

タグ：

posted at 21:03:36

やり方はこちら

franklinjl.org/workflow/deplo...

タグ：

posted at 21:04:22

成果物はこちら．過去のものを全消しして
テンプレートを push しました．

terasakisatoshi.github.io/index.html

タグ：

posted at 21:04:22

かずおじ@pythonerになりたい @kazu_oji

マニュアルにも書いてありますが，添付の写真のようにして
gh-pages を Branch にセットしてください．
こうしないと意図しない結果が得られます．

#Julia言語 pic.twitter.com/pxmOXsAuOi

タグ： Julia言語

posted at 21:06:13

最近はJuliaという言語が気になってます。
ニューラルネットとか簡単に書けるのかな。
早いうちにトライします。

タグ：

posted at 21:13:05

じりおん @zillione

計算と制御はJuliaに移行したけど、データ解析と作図はPythonを使ってしまう

タグ：

posted at 22:00:36

#Julia言語

昨日わかったこと

```
julia> goma goma kyu kkyu
```

というサンプルコードをREPLで動かすときに丁寧に

`julia>` の部分を含めずにコピペしてたけれど

`julia>` も含めても REPL 側でよしなに無視してくれる.

し・ら・な・か・た・・・

タグ： Julia言語

posted at 22:01:00