黒木玄 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年02月25日(日)
サイズnの1つのサンプルからの推定の計算に数秒かかってしまうと1000回の繰り返しに1~2時間もかかってしまいます。数秒が0.1秒になってくれると100秒ですみます。この違いは大きいです。
各々は短時間ですむ計算であっても、統計がらみの確率がからむ場合にはたくさんやる必要が出て来る。
タグ:
posted at 23:57:02
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
須山敦志 Suyama Atsushi @sammy_suyama
上の「モデル構築→推論導出→実験→評価」のサイクルに若手を投入すると信じられない速さで機械学習を習得してくれます.やる気や基礎能力に応じて1サイクル当たりのハードルの高さ(モデルの拡張度合いとか)を調整するのがポイントだと思っています.
タグ:
posted at 23:07:45
#Julia言語 を使っていて「極めて遅い」場合に出会ったら、
nbviewer.jupyter.org/gist/genkuroki...
を見るといいかも。
タグ: Julia言語
posted at 23:03:45
須山敦志 Suyama Atsushi @sammy_suyama
本当はStanを始めとする確率的プログラミング言語を使って「アルゴリズムを導出する」部分をすっ飛ばすのが良いのですが,当時は推論がめちゃくちゃ難しいモデルを扱わなきゃならないことが多かったので,自力で効率の良い推論を探索する必要がありました.
タグ:
posted at 22:58:46
須山敦志 Suyama Atsushi @sammy_suyama
モデルを書く→アルゴリズムを導出する→実験する→失敗する→モデルを修正する→アルゴリズムを再導出する→・・・
というサイクルを毎日永遠と繰り返していたので,計算速度を損なわずにカジュアルに内部が変更できる必要があったんですね.
タグ:
posted at 22:55:23
twitter.com/inert_sh/statu...
#Julia言語 的には「Pythonの3倍」は「極めて遅い」ので、ソースコードを公開すれば誰かが改善案を示してくれる可能性が高いです。
遅くなる原因でありがちなのは、
1. 計算の主要部分を函数の中に入れていない。
2. @ code_warntype での赤字警告を放置している。
タグ: Julia言語
posted at 22:42:28
非公開
タグ:
posted at xx:xx:xx
OokuboTact 大久保中二病中年 @OokuboTact
筑波付属小算数の本はたくさん出ている。山本良和と田中博史の本が代表的なのかなあ?
代表的というのは、「特徴がよく出ている」という意味合いなんだけど
タグ:
posted at 22:38:02
須山敦志 Suyama Atsushi @sammy_suyama
既存の機械学習ライブラリやAPIをガチャガチャ組み合わせてアプリを作る場合はPythonが便利ですが,モデルや推論アルゴリズムを直接いじってオーダーメイドする場合は今のところJulia一択です.
タグ:
posted at 22:34:30
DLAmodelで、コロイド粒子(粒子数N=1200)が、原点にいるコロイドにくっついたときにできる塊(クラスター)(フラクタル性を持つ)の一例を作った
Juliaでも102sかかる
写真はN=300,600,900,1200でクラスターが成長する様子が分かる pic.twitter.com/zn2gbFon6o
タグ:
posted at 22:34:12
OokuboTact 大久保中二病中年 @OokuboTact
小学校の算数は、草の根運動のTOSS VS エリート育成の筑波付属小の構図で見るとわかりやすい。
どちらも酷いのだけど、筑波付属小教育を受けた生徒の多くは有名大学に合格している人が多いと思う。彼(彼女)達は筑波付属小の算数教育を素晴らしいと思っているんだろうか?
タグ:
posted at 22:33:56
@mamorumatsuo @vecchio_ciao [{()}]は日本の中学高校にローカルな慣習です。日本の中学高校の数学の教科書と中学高校の先生が書いた文書の他には、見たことがありません。そして、大学入試の採点をするのは大学の先生であって中学高校の先生ではありません。「大学の先生がそう言っている」が説得に役立つと良いのですが。
タグ:
posted at 22:33:36
Juliaのplotsは初めて使ったけど、出力されるグラフの設定を細かく指定したくても、日本語のサイトがあんまりなくて英語を読まなきゃいけなかったのが大変だったなぁ
タグ:
posted at 22:18:17
OokuboTact 大久保中二病中年 @OokuboTact
超算数タグで何度か話題になった山本良和(筑波付属小)の本を買ってみた。
記号の書き順とかいろいろ書いてある
タグ:
posted at 22:16:24
OokuboTact 大久保中二病中年 @OokuboTact
駅ビルの5Fにあるかなり大きな本屋に行ったけど、ガラガラ。
本屋は委託商売だけど、人件費は赤字な気がする
タグ:
posted at 22:13:22
3man chessまとめ
・真ん中の円部分は通過できる
・ビショップは鬼神の如き働き
・ルークは中心3列に入ると円を描いて殺戮を始める
・気を許すとすぐ駒が落ちる
・許さなくても落ちる
・開始1時間程度でチェス盤酔いが始まる
・チェックと言われても誰がどうやってチェックされたかわからない
・最高 pic.twitter.com/LZaLwxy9Cz
タグ:
posted at 22:04:02
なされるのであれば、っていうか、既になされているわけですが。教員は勝手気儘だから拒否しやすいんだけど、事務職員にも寄付を求めていて、そちらのほうが同調圧力が強いのではないかと懸念しますね twitter.com/f_nisihara/sta...
タグ:
posted at 21:34:26
@hankagosa 乱数の種を変えれば上手くいくのですが、やや気になるのは、実際にpoisson_log_rngに渡っているパラメータをチェックする方法はないのかな?って事です。
lambda[n]がそのまま渡ってる訳ではなさそうなので。
タグ:
posted at 20:49:51
@Akimasa_K お手数おかけしてすみません。途中でそのエラーが消えるなら無視でOKかと。各パラメータの初期値は-2から2の乱数で初期化されるのですが、その組み合わせによってはどこかのexpの内側が大きくなりすぎるようです。気になるなら小さめの初期値を指定します。
タグ:
posted at 20:20:47
rstanのバージョンは2.17.3 です。
エラー自体は、poisson_log_rng が内部的にexpを呼ぶときに、結果が32bit整数の範囲を超えないようにするためのエラーチェックに引っかかってるぽいのですが、
参考: discourse.mc-stan.org/t/help-with-po...
そこまで大きな値を渡してはいないはずなのです。
タグ:
posted at 20:03:44
ごまふあざらし(GomahuAzaras @MathSorcerer
実質DellXPSを購入して自分でいろいろいじってUbuntuを入れるという解決策しか知らないこの頃ですし.
タグ:
posted at 20:03:24
@hankagosa RStan本の演習問題5-(6)をやっているんですが、chap05/exercise/ex6.R を実行すると
poisson_log_rng: Log rate parameter is 21.2, but must be less than 20.7944
ってエラーが出ます。
でも、rstan::extract(fit) してlambdaの最大値を調べてみても、そこまで大きく無いのです。
タグ:
posted at 20:01:12
むぅ、Stanから謎のエラーが返ってくる。
「poisson_log_rngのパラメータは20.7944以下であるべきなのに、21.2が入ってきたよ」ってのは、どうやら32bit整数の制約みたいだが、
MCMCサンプルでは最大2.518しか突っ込まれてないっぽい...
タグ:
posted at 19:52:31
Fumiaki Nishihara(西原 @f_nisihara
昨晩、阪大の入試ミス対応について色々ツイートしたのですが、それの簡単なまとめです。 - 阪大は入試ミス対応のために教職員に寄付をつのるべきでない id.fnshr.info/2018/02/25/han...
タグ:
posted at 19:46:04
ウィキペディアを信用するならば「日本などの一部の国では、数式における括弧の入れ子は [{()}] の順で用いられてきた。しかし、世界的には {[()]} の順で用いられる方式が多数派である。(中略)現在は「丸括弧」・「角括弧」・「波括弧」と表記されている」──誰が大中小を勝手に決めたんだ? twitter.com/mamorumatsuo/s...
タグ:
posted at 19:03:04
#Julia言語
github.com/JuliaLang/juli...
を見るまで知らなかったのは resize! の存在。
途中で配列のサイズを増やすようなコードを書きたい場合には、push!で要素を追加するのではなく、resize!で一挙に配列のサイズを増やすと効率的なようです。
タグ: Julia言語
posted at 18:50:36
@bicycle1885 同じ掛け算でも × は結合則なくて、 AST レベルで区別つけてるの面白い気がします twitter.com/tkf/status/967...
タグ:
posted at 18:26:23
やっとJuliaで DLA model (ブラウン運動する粒子が凝集して大きなクラスターを形成するモデル)を実装できた
きちんとした画像を出力できるようにしたらツイートしよ
タグ:
posted at 18:15:43
@tkf 結合則が無ければa * b * cとは書かずに (a * b) * cなどと書くので,それはそれで良いのではないでしょうか。括弧をつければちゃんと構文木がとれます。 pic.twitter.com/yWgyOeaaYY
タグ:
posted at 17:53:11
Julia の和や積って AST のレベルで結合則があるのがな、って思ったけど lisp でも (* a b c) って書くしな。 pic.twitter.com/3V3znwts2m
タグ:
posted at 17:44:04
数学クラスタの皆様。
2016年ぐらいからポケモンのカビゴンのパロディキャラクタの「カンビゴン」というのを考えているのですが、名前だけでいいキャラ付けが思いつきません。どうかお知恵をください🙏 pic.twitter.com/LkC4s35Jxp
タグ:
posted at 17:15:46
#Julia言語 HCubature.jlという数値積分パッケージではDataStructures.jlというパッケージを使っていた。
github.com/stevengj/HCuba...
github.com/JuliaCollectio...
Julia言語では、適切に定義された型は速度劣化の原因にならないので、適切な型を定義して使うようにするのがよいのだと思う。
タグ: Julia言語
posted at 17:01:10
斉藤ひでみ・現職教師(西村祐二) @kimamanigo0815
【読売新聞 2018.2.25】
給特法改正に向けての議論がいよいよ本格化する見込みだ。 pic.twitter.com/nj9iX3Tm5W
タグ:
posted at 16:59:24
#Julia言語 配列の要素に代入する操作と比較すると、配列にpushする操作はずっと複雑なことをやっているので当然遅くなる。sizehintを使ってもその点は変わらない。
github.com/JuliaLang/juli...
Q. sizehintを使ってもどうしてあまり速くならないの?
A. push!が遅いから。
タグ: Julia言語
posted at 16:57:24
@Gelada Pseudosphere in Dynamic Geometry
Altitudes : huit.re/PShaut
Medians : huit.re/PSmedianes
Circle circ : huit.re/PScCirc
Move uA, uB, uC and A, B , C and te pseudosphere by clic-mouse pic.twitter.com/d2ZGP5nB9i
タグ:
posted at 16:34:21
接吻数問題 と 24 次元リーチ格子 - tsujimotterのノートブック tsujimotter.hatenablog.com/entry/kissing-...
タグ:
posted at 16:07:49
TikZ で図を描いていらっしゃるからこういう言明になったのだと思うのだけど,別に TikZ でなくても「図を描く」と執筆ペースは落ちるんじゃないかという気がする.
タグ:
posted at 16:05:45
それから、自主避難支援で最も重要なのはカウンセリングなんですが、これもなかなかはっきり言われない。そして、「自称支援者」による「囲い込み」から救い出す必要があるんだけど、そんなことは怖くて誰もできない。でも、必要なのはそういうことです
タグ:
posted at 15:43:02
福島からの自主避難とそれ以外の地域からの自主避難とは明確に分けて考えるべきだよね。これを言う勇気のない人が多いですけど、これは誰かが言わなくてはならないんだよ
タグ:
posted at 15:40:34
numirが充実してきたのでnumpyやjuliaで面白いと思ったコードを元にexampleを追加しています。今はchar RNNの言語モデルと、NMFの音源分離があります。しかし、まだまだ小規模なのでnumpyやjuliaにあるこの機能が欲しい人はぜひissue立ててください github.com/libmir/numir/t...
タグ:
posted at 15:24:43
tsujimotter 日曜数学者 @tsujimotter
また、先ほどのRTで言及されていたリーチ格子についてはこちらなんかいかがでしょう(キャノンボール問題とリーチ格子が関係するというのは初耳でした)
接吻数問題 と 24 次元リーチ格子 - tsujimotterのノートブック tsujimotter.hatenablog.com/entry/kissing-...
タグ:
posted at 15:09:08
tsujimotter 日曜数学者 @tsujimotter
仮面ライダービルドの24話の数式はリュカのキャノンボール問題だったみたいですね。関連記事はこちら。
リュカのキャノンボール問題 - tsujimotterのノートブック tsujimotter.hatenablog.com/entry/cannon-b...
タグ:
posted at 15:07:41
CPU pipelines: when more is less www.juliabloggers.com/cpu-pipelines-... #julialang (Re-post)
タグ: julialang
posted at 14:50:36
沢山の投票ありがとうございました!私が学生のときは主に Mathematica (学内無償)とサブとして C 言語を使っていました。その時と状況はだいぶ変わっているようですね。他の言語として Haskell、C++、Fortran、Julia などご意見を頂きました。ありがとうございました。
タグ:
posted at 13:42:30
Julialangでlibrary使わずにLogistic Regression (L2) やってみた。150行くらいで書けた。でも意図している結果と結構違って、機会があったら識者の意見を仰ぐためqiitaあたりに何かかくかもしれません。。。
タグ:
posted at 12:43:13
非公開
タグ:
posted at xx:xx:xx
JuliaでNLsolve使いたいのに一番簡単な例( www.cas.cmc.osaka-u.ac.jp/~paoon/misc/ju... )でやってもNaNが出てしまうんですが… pic.twitter.com/SpDycm6P0D
タグ:
posted at 11:55:02
バイクで転倒して頭を打ったけれどどこも痛くないですでも念のため来ました、と言う若者がレントゲン撮ったら頸椎破裂骨折しててすごく驚いた事はあるから、高エネルギーの事故は要注意ですね。
タグ:
posted at 11:43:47
きむりん(小六女子:Explicit) @ura_kimrin
Julialangでライブラリ使わずにLogistic Regression (L2) やってみた。150行くらいで書けた。でも意図している結果と結構違って、機会があったら識者の意見を仰ぐためqiitaあたりに何かかくかもしれません。。。
タグ:
posted at 11:38:56
非公開
タグ:
posted at xx:xx:xx
まー、阪大を含む国立大学職員は、東北の大震災のときも給料10%くらいさっ引かれたからなぁ。あれは誰のミスでもない天災だったのに。理不尽なことこの上ない。 twitter.com/f_nisihara/sta...
タグ:
posted at 11:11:23
非公開
タグ:
posted at xx:xx:xx
入試でミスを無くすことなど絶対にできない。再発防止策なんか、いくら議論したって無駄。ヒューマンエラーは必ず起きる。ミスが嫌なら入試を廃止するしかない。ミスが判明した時のリカバリーを議論しておくほうが建設的
タグ:
posted at 10:43:03
この理論はあまりうまく成功しなかったのですが、この方向性は10次元の超ひも理論として復活し活発に研究されています。10次元なのは、8次元のE8格子という性質の良い格子があり、これに落とした2次元を足して10次元という背景もあります。8次元の充填問題はE8格子、これに対応する符号もあります。
タグ:
posted at 10:22:13
リーチ格子の構成法の中で一番明快なのは、26次元ミンコフスキー空間に埋め込む方法です。最初の関係式により(0,1,2,…23,24,70)をワイルベクトルにとって、適切な射影で2次元落とすことで構成できます。初期のひも理論は26次元空間の理論でしたが、その背景の一つにはこの関係があります。
タグ:
posted at 10:20:50
惑星探査機ボイジャーの通信にも使われたゴレイ符号という符号があります。これは3ヶ所までエラーを修復できる完全符号なのですが、これはリーチ格子と同じ方法で構成できます。そしてこの符号の置換は、有限単純群の分類でも現れる散在単純群のマシュー群で与えられます。
タグ:
posted at 10:20:04
2次元円盤を平面に出来るだけ多く充填する問題を考えると、これは三角格子で置くのが最良です。この3次元版はケプラー問題で、90年代にやっと証明されました。この高次元版は実は8次元と24次元だけ(!)が解決しており、24次元の最密充填はリーチ格子の配置です(証明は2016年)。
タグ:
posted at 10:19:07
#仮面ライダービルド 今週のオープニングの話数の式は、1+4+9+…+N^2が平方数になるのはN=24の場合(和は70^2)のみ、という結果(Lucas予想)です。これはリーチ格子という24次元空間(!)の格子の背景にある関係の一つで、充填問題、符号理論、代数、さらにはひも理論まで幅広い分野と関係してます。
タグ: 仮面ライダービルド
posted at 10:17:42
GitHub に AWS キーペアを上げると抜かれるってほんと???試してみよー! on @Qiita qiita.com/saitotak/items...
タグ:
posted at 10:02:31
非公開
タグ:
posted at xx:xx:xx
阪大の例の寄付は本部から直接の依頼ではなく「学部長一同」からの依頼という形式なんですが、じゃあ本部としては経費をどうやって捻出する気なんだよ、と問い詰めたくなるわけです。この取り繕いかたがまたセコい
タグ:
posted at 09:45:03
もう10年以上前だけど、学生のときに某大学から入試の採点を依頼されたことがある。学生に入試の採点させてる。そうしないと回らないほど、文科省は我々の業務を圧迫している。 twitter.com/tori_esa0110/s...
タグ:
posted at 08:57:38
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
周波数依存重み付きLQRのJuliaによるシミュレーション
概要 クリスマスが終わっても寄稿できるようなので枠を頂きます。タイトルの通り周波数依存重み付きLQR(FWLQR: Frequency Weighted LQR)の紹介とJulia Contro… ift.tt/2GHDfKS
タグ:
posted at 05:18:25
#数楽 以上に話題はクラスター代数の話の特殊な場合。ランク2のクラスター代数に関する情報私による詳しい解説が
genkuroki.github.io/documents/clus...
にあります。高校生~大学1年生向けのフリーズパターンの解説が次の場所にある。
genkuroki.github.io/documents/2012...
タグ: 数楽
posted at 04:45:11
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
呪いの画像に対する反応として「『改行』で覚えればいける」「columnは柱の意味だから」「縦列駐車横行」「A列車で1行こう」「『行高』『列幅』っていうから大丈夫」って役に立つ覚え方がたくさん寄せられてきてるからあれはむしろ幸福の画像だったのでは???
タグ:
posted at 02:25:24
Morphs III (renaissance): alteredqualia.com/xg/examples/mo... #webgl #deepwarp pic.twitter.com/8l8KxTprtp
posted at 02:20:16