黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2018年03月22日(木)
2枚目の画像の全文はこちらになります。Python(numpy)の方についてnp.array([[1, 2,3]])とすればいけるというご指摘がありましたが、MATLAB的には、Pythonで１次元配列、2次元配列の違いを意識しないといけないのは、難しいだろと言いたいようです。 pic.twitter.com/K2znbNsmmE
タグ：
posted at 00:01:54
JuliaはPythonやRubyなんかと比較するべき言語であって、型や最適化がそこそこ強いので処理内容によっては追いつける場合はあるけど、CやC++やRustと同レベルの最適化を目指した設計には初めからなっていないし、そういう観点で語っても意味はない気がする
タグ：
posted at 00:23:02
twitter.com/kazatsuyu/stat...
#Julia言語 原理的にはJulia言語は速さでCには勝てないと私も信じています。
しかし、現実にはやってみないとわからない。
数学的に必要なEi函数がJuliaに実装されていなかったので、CのコードをそのままJuiaに翻訳して試してみたら、Cと同じ速さでした。証拠に続く
タグ： Julia言語
posted at 00:35:49
#Julia言語 証拠は
nbviewer.jupyter.org/gist/genkuroki...
にあります。Cで書かれたEi函数を極めて忠実にJuliaに翻訳しています。C版とJulia版で計算速度を比較したら、気にするような違いはありませんでした。もとのCのコードは2種類で片方ではJuliaの方が微小に速かったです。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 00:38:51
#Julia言語 で特殊函数のライブラリを書くことは、実用的に十分に考えられる選択肢です。
E1函数をJuliaで実装して最適化するという内容の極めて面白い話を
nbviewer.jupyter.org/github/steveng...
で読めます。Python scipyで使われているFortranで書かれたライブラリより、Julia版は5～6倍速くなっています。
タグ： Julia言語
posted at 00:42:11
いくらJITがあってもCには勝てないだろうという考え方は多くの場合妥当ですが、事例によってはあっさり抜くので油断できない。古くはCommon Lispなどで blog.practical-scheme.net/shiro/20100620... twitter.com/genkuroki/stat...
タグ：
posted at 00:54:14
まあ、実際に試してもいない人が何を言っても説得力ゼロ。
私がやっている程度に実証的なきちんとコードを全公開して話を続けないと意味がない。
みんなが欲している情報は、馬鹿でもわかる一般論ではなく、実際にうまく行く場合(うまく行かない場合)に関する具体的な情報だと思います。
タグ：
posted at 00:55:25
すでに普及している特殊函数のライブラリの中にも最適化が不十分なものがあるという事実は、若い人達にとっては名を残すチャンスかも。(←Julia言語とは無関係な話)
「2018年の〇〇氏による最適化によって△△函数のライブラリ群が3倍高速化した」という感じに名を遺すチャンス。
タグ：
posted at 00:58:13
①公立教員を残業代の支払いの対象外とする #給特法 （公立の義務教育諸学校等の教育職員の給与等に関する特別措置法）や
②私立学校や民間企業で、残業してもいいよという協定 #３６協定 （時間外労働・休日労働に関する労使協定）
について、わかりやすく書かれていますね。
news.yahoo.co.jp/byline/konnoha...
posted at 00:59:33
Ei函数をわざわざ #Julia言語 で実装しようとしたモチベーションはJuliaでお手軽に素数定理を数値的に確認しようとしたら、li(x)がJuliaに見当たらなかったからです。li(x)=Ei(log(x))なのでEi函数を実装すればよいわけです。
結果的に添付画像をうまく描けました。楽しい！
pic.twitter.com/AchyIxQxBx
タグ： Julia言語
posted at 01:10:03
幾つか特殊函数の数値計算のライブラリのソースコードを見てみたのですが、ものすごく面倒そうに見えたのは、オーバーフローとアンダーフローへの対策です。函数の性質と実際に使われる場面について非常によく理解していないとものすごく難しい。
タグ：
posted at 01:14:33
twitter.com/kazatsuyu/stat...
単なる思い込みは怖い。
「 #Julia言語 はC,C++,Rustのような速さを出せるような設計になっていないことは明らかである」なんてJulia言語開発陣に言ったらフルボッコにされることは確実(笑)。
論よりコード
タグ： Julia言語
posted at 01:23:30
引き続き、「プログラミング言語 Julia の紹介」を読むもしかして読んだことあるかも)
www.slideshare.net/kentaroiizuka/...
タグ：
posted at 01:25:41
なんと言われようとJuliaが動的型付けであることに起因する関数初回呼び出し時のコストの存在は自明ですし、言語設計者がそれをわかっていないとは思えませんが twitter.com/genkuroki/stat...
タグ：
posted at 01:34:31
斉藤 淳 『アメリカの大学生が学んでいる @junsaito0529
帰国するときに、大学に就職しないで起業した理由の一つを思い出した。「誰でも分かるタイポ」を学会や研究会で指摘するオーディエンスの存在。「漢字が間違ってますね」レベルの指摘を大勢の前でして溜飲を下げる、書評でどうでもいい史実の誤りを指摘する、そういうのが嫌でたまらんのだった。
タグ：
posted at 01:43:52
twitter.com/kazatsuyu/stat...
#Julia言語 JITコンパイルのオーバーヘッドの問題解決は将来の問題。
まだ完全実用には程遠いのですが、すでにAtcoder方面向けにコメントしたことですが、Juliaもコンパイルしてバイナリを作れます。
github.com/JuliaLang/Pack...
タグ： Julia言語
posted at 01:44:46
#統計 #Julia言語
MCMCパッケージでよくテスト用に実行されている所謂 Eight Schools を #Julia言語 で扱ってみました。
nbviewer.jupyter.org/gist/genkuroki...
「八学校」への最尤法とベイズ推定法の適用
ベイズ推定はMambaを使い、最尤法は自前のコードで行いました。続く
posted at 02:11:14
#統計 モデル(A)で最尤法を機械的に試してみたら、ρ=0と推定されました。そこでモデル(A)の尤度函数を描いてみたら添付画像のようになりました。(尤度函数は最大値が1になるように正規化してあります。)
確かに尤度函数はρ=0で最大になっています。 pic.twitter.com/8jY6nXdlej
タグ： 統計
posted at 02:31:00
#統計 もしかしたら、
nbviewer.jupyter.org/gist/genkuroki...
「八学校」への最尤法とベイズ推定法の適用
は #Julia言語 でのMCMCパッケージであるMamba.jlのシンプルな使用例として貴重かも。今回はスライスサンプラーを使いました。
#Julia言語 でもStan.jl経由でStanを使えます。
posted at 02:37:43
非公開
タグ：
posted at xx:xx:xx
斉藤 淳 『アメリカの大学生が学んでいる @junsaito0529
掛け算や足し算の順番、漢字のトメ、ハネなどのどうでも良い間違いを指摘され続けて育った子どもが大人になると、他人のどうでも良い間違いを指摘し、揚げ足をとってマウンティングするようになる。日本社会で一番嫌いなのが、そういう大人の存在だなぁ。
タグ：
posted at 06:33:52
ネット上でいろんな変種をたくさん見てきた。
www.boj.or.jp/mopo/mpmsche_m... pic.twitter.com/JuvNisc5an
タグ：
posted at 08:07:47
時系列の統計モデルをAICでモデル選択することが自分のなかで納得できていなくて、その路線の本やブログ記事を安易にRTできないでいる。AICは未来の予測の良さじゃないんだよなぁ。何を測ってるか自覚ありますかね。
タグ：
posted at 10:22:00
不景気になって低賃金労働者のクビが切られると「実質賃金は“上がる”」もの、民主党政権下でも、昭和恐慌下でも、実質賃金が上がっているのはそういう事、と言ったら、「ならば自分はアベノミクスよりも恐慌を選ぶ」と言ってのけた「アベ嫌い」の人が現れた時はさすがに驚いた。 アタマオダイジニー。
タグ：
posted at 11:33:52
あしたこれ行く予定なので今日は Julia 言語に入門しよう（１㍉も知らない・・・
data-refinement.connpass.com/event/81109/
タグ：
posted at 11:36:17
今日次男が持って帰ってきた算数のテストでは「長方形をすべてえらべ」という設問の模範解答が正方形を含めないものになっていた。授業でもやはりそのように習ったらしい。昔（僕が小学生の頃）は違ったと思う。
何故こんなになってしまったんだ。日本の初等教育・中等教育はもうダメかも知れない。
タグ：
posted at 11:39:59
非公開
タグ：
posted at xx:xx:xx
@benthed00peness コメントありがとうございます。僕の主張はこれです。
twitter.com/hankagosa/stat...
もし未来の時点における予測分布への評価を議論している論文などがありましたら教えていただきたいです…。
タグ：
posted at 11:52:02
@y__hiroyuki 高校教員の話しでは、比率計算の意味が分かっていない生徒が目立つそうです。小学生には「きはじ」で距離=速さ×時間を覚えさせるそうで、苦労して意味を分からせるより手軽な公式暗記型教育の蔓延でしょう。先生もその教育で育っていますから、ますます物理は苦手って事になりますね。暗澹たる思い・・
タグ：
posted at 11:53:09
結果出してこその独立性 中原伸之・元日銀審議委員：日経www.nikkei.com/article/DGXMZO...
今の日銀では政策委員は外部との接触が制限されていると聞く…事務方が作った経済見通しや政策に賛成せざるを得なくなっている…事務方が審議委員の行動を制限するなど言語道断の越権行為で日銀法の趣旨にそぐわない。
タグ：
posted at 12:28:19
「＜東北大雇い止め＞東北18国公立大の無期雇用転換ルール 東北大のみ実施せず」
“現在、５年超の有期職員がいない山形大と福島大も「要件を満たせば５年を超えられ、無期転換を申し込める」などとしている。”
headlines.yahoo.co.jp/hl?a=20180322-...
タグ：
posted at 12:39:46
matlabとjuliaとpythonでベクトルと行列の実装が違うらしい。matlabには行列としての列ベクトル、行ベクトルはあるが、ベクトルがない。juliaはベクトルを転置すると、行列でない行ベクトルになるが、pythonはベクトルを転置しても変わらない、であってる？
タグ：
posted at 12:48:31
@y__hiroyuki @genkuroki ＃小学校 の ＃算数 では、執拗に「長方形と正方形は別のもの」という刷り込みをされていますね。
twitter.com/vecchio_ciao/s...
posted at 12:53:17
@vecchio_ciao @y__hiroyuki うちでも、算数プリントの「長方形を選べ」という問題の回答欄が、正方形2個を合わせて長方形4個なのに、2つしかなくて、結局、4つに増やして回答していたことがありました。
後で「丸になっていた」と教えてもらった。たぶん採点基準を統一していない。
教科書準拠の算数プリントが色々酷すぎ。
タグ：
posted at 13:00:26
仮に安倍政権がコケて金融引き締め&緊縮財政になったとしたら、当然再びデフレで大不況の就職難時代が来るわけだけど、さすがにそうなれば何が問題だったのか今よりも多くの人が気付くわけで、いわゆる“リベラル”さん達の主張の結果が大不況だと理解されたら… 今度こそ真の極右政権生まれるかもよ。
タグ：
posted at 13:16:15
＜東北大雇い止め＞東北１８国公立大の無期雇用転換ルール 東北大のみ実施せず | 河北新報オンラインニュース www.kahoku.co.jp/tohokunews/201... @kahoku_shimpoさんから 諸悪の根源は文科省天下りの大槻理事。里見総長は大槻理事を引責辞任させ、違法な無期労働契約移行妨害を速やかに止めよ!!
タグ：
posted at 13:21:40
＜東北大雇い止め＞東北１８国公立大の無期雇用転換ルール 東北大のみ実施せず | 河北新報オンラインニュース / ONLINE NEWS sp.kahoku.co.jp/tohokunews/201...
タグ：
posted at 13:22:55
文科省から天下りの大槻理事が諸悪の根源。里見総長は、労契法を故意に脱法することに固執している大槻理事を引責辞任させ、他の国立大学法人並にきちんと無期雇用転換ルールの適用を! このままだと東北大学は再びブラック企業大賞にノミネートされ、指定国立大からも twitter.com/nichidai_union...
タグ：
posted at 13:27:34
大学入試数学に関する迷信は、多くの大学の先生がtwitterとかで否定してくれていて非常にありがたい。同様に、小学校の数学に関する迷信（掛け算の順序とか）に関して、中学入試を採点する中学校の先生たちからメッセージが欲しいなぁ。
タグ：
posted at 14:14:35
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
すごい発明をしたかも…💡英語で技術ブログ書くとき、HackmdとGoogle翻訳を組み合わせると、Markdown+英語で書きながら、内容をGoogleのAIにチェックさせられる。英語として流暢かはさておき、少なくともGoogleが理解できる英語にはなる。この手法はGitHubにIssueを書く際も応用可 pic.twitter.com/u5sA9prv91
タグ：
posted at 15:18:01
これにGrammarlyのChrome拡張を組み合わせれば文法の添削もできて最強の英文ライティングIDEになる気がしてきた。が、試してはいない😁
www.grammarly.com
タグ：
posted at 15:27:10
twitter.com/hankagosa/stat...
#統計 リンク先のリンク先のリンク先がリンク切れなのですが、
web.archive.org/web/2017042704...
で読めます。
タグ： 統計
posted at 15:57:29
#統計 リンク切れの
web.archive.org/web/2017042704...
は
qiita.com/advent-calenda...
からもリンク切れ。
もしも取り下げられたのだとしたら、もったいないと思いました。このレベルで何か議論になるならそれだけで社会貢献になると思う。
タグ： 統計
posted at 16:09:30
"RやPythonでデータ分析をしている人がパフォーマンスの問題に悩まされた時に, 救いを求めて手を出すことがあると言われているプログラミング言語, Julia"
www.mwsoft.jp/programming/ju...
タグ：
posted at 16:10:44
ゆえに、S(q)=-∫q(y)log q(y)dyとおくと、p(y)のサンプルY_iに関する対数尤度
Σ_{i=1}^n log p(Y_i)
の-1/n倍はn→∞で
D(q||p) + S(q)
という
p(y)によるq(y)の予測誤差 + p(y)によらない量
に収束します。これを
pの予測誤差+定数
と略記しましょう。
タグ：
posted at 17:13:20
#統計 Kullback-Leibler情報量によって、「モデルの確率分布によるシミュレーションにおける真の分布の予測誤差」と呼ぶべき量を表せること(Sanovの定理)に関する易しめの解説が次のリンク先にあります。
genkuroki.github.io/documents/2016...
雑談終わり
タグ： 統計
posted at 17:13:28
ねぇ、このモーメント、悪意あると思うんだけどw
まぁ、ニュースの見出しだけ見て、福島県の広さを知らない方のご意見かな。
もうずっと放射線量変化ないし、誰も見てないモニタリングポストをいつまで置いておくのか疑問だったから、撤去賛成だけどな。
twitter.com/i/moments/9766...
タグ：
posted at 17:15:21
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
僕もずっとそう思って（いたがって？）ました。色々と此れまでも片鱗はあって、まぁそれでも所詮学校の先生だし、教育学系やと高校生に毛が生えた程度で凄いみたいやからまぁ致し方ないと想うてもいましたが、実際に娘が掛算で被害に遭い問い質して実態を知りここ数年の新規採用な先生達の言動を観察し twitter.com/y__hiroyuki/st...
タグ：
posted at 18:16:14
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
てると、なんだろう、文章が読めないとかAIガァとかいうよりは何キロも手前の悲惨な状態です。いやもあほんまに教えたらあかんやろという...本人達に悪気はなさそうなのですが、2年ほど一緒に働いた若者などはまるで良い子ロボットで己の頭で考えるというのはパターンを探すことのような為体で...辛い
タグ：
posted at 18:18:42
娘が漢字苦手で、しかも先取り学習をしたら「まだ習ってない！」っで怒鳴りながらおこるんだよね。原因は担任が習ってない漢字を書いた子を立たせて怒鳴るからだったけど、クソ担任はどれだけクレームが入っても先行学習した子に罰を与え続けたんだよね。
タグ：
posted at 18:24:13
んで、娘と違って精神的に成熟した子はくもんと学校を使い分けして6年までの漢字を先行してるのを、やっと娘が気づいて今さら漢字をやる気になった。先生が他のクラスの担任に恋をして担任をやる気を失ってくれたせいで呪縛が溶けてよかった
タグ：
posted at 18:26:28
マスコミが批判すれば「下衆」、逆に肯定的に擁護したら「人格者」。つまりは、テレビ中毒、マスコミ中毒なんだろうな…。 twitter.com/xWAypYyLj09DhK...
タグ：
posted at 19:32:05
生徒が持ってきた情報をもとに現在 Trilinear Coordinates について教材研究中。位置ベクトルの話……？ ベクトルの題材でどこまで使えるだろうか。
タグ：
posted at 20:03:02
@ysmemoirs ベクトルとの相性は barycentric coordinates のほうが良いです。trilinear coordinate が l:m:n のとき b.coord. は al:bm:cn （a,b,cは辺長）と変換は簡単ですので、t.coord. からいちいち b.coord. に変換して計算しても手間は大きくありません。
タグ：
posted at 20:15:04
@ysmemoirs ベクトルとの相性の良さは、barycentric coordinates で X=l:m:n と書けるとき、Xが無限遠点でなければ位置ベクトルで (l+m+n)X=lA+mB+nC と書けることです。
タグ：
posted at 20:17:23
掛け算の順番とか、習っていないことは使って解くと×をつけられるとか、厳しい出来事が多く見受けられます。全国の教育学部の教員養成での教育はどうなっているんだろうか。大学の教員養成の教科の教育は教育学部から切り離すべきかもしれません。 twitter.com/y__hiroyuki/st...
タグ：
posted at 22:07:21
雪国の旅のお供におひとついかが？
#CuBag
#Cuboard
#私の作品もっと沢山の人に広がれ祭り pic.twitter.com/PT3DnWQMGC
タグ： CuBag Cuboard 私の作品もっと沢山の人に広がれ祭り
posted at 22:45:24
非公開
タグ：
posted at xx:xx:xx
twitter.com/less_role/stat...
Anaconda3のJupyterで #Julia言語 を使っています。私も最初はdead kernelの永久ループに悩まされました。
ENV["JUPYTER"]=raw"jupyterバイナリのフルパス"
の状態でPkg.add("IJulia")したらうまく行きました。
RをJupyterで使うことにはトラブルの経験無し。
続く
タグ： Julia言語
posted at 23:06:42
#Julia言語 私によるインストールの詳細な記録が
nbviewer.jupyter.org/gist/genkuroki...
にあります。
JuliaカーネルがJupyterで死ななくなる方法は第1.4節にあります。(結構みんなはまっている。)
RをJupyterで使えるようにしたときの記録は第1.14節にあります。
タグ： Julia言語
posted at 23:06:43
#統計 #Julia言語
nbviewer.jupyter.org/gist/genkuroki...
最小二乗法による回帰の場合のAICと予測誤差の分布(NEW!)
nbviewer.jupyter.org/gist/genkuroki...
対数尤度の比較によるモデル選択の簡単な例(2017-11-01)
nbviewer.jupyter.org/gist/genkuroki...
AICと汎化損失の簡単な計算例(2017-11-06)
公式がある場合も数値計算ですませた。
posted at 23:25:37
読書猿『独学大全』14刷26万部、『文章 @kurubushi_rm
告白しますが、実は動画を見るのが苦手で、特に学習系の動画は、読めば１０秒で終わる情報を得るのに、数分間相手のペースに合わさなくてはならないのが耐え難く感じます。
#マシュマロを投げ合おう
marshmallow-qa.com/kurubushi_rm?u... pic.twitter.com/7ZcinoFO7B
タグ： マシュマロを投げ合おう
posted at 23:40:00
