黒木玄 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 を見て下さい。
2021年01月24日(日)
非公開
タグ:
posted at xx:xx:xx
1.5は上手く行きました。で1.6betaもと続けてやったらば add IJulia で滞って進まず... #知らんけど な感じ。何でしょね。らまぁ1.5は使えるのかなぁ。 twitter.com/genkuroki/stat...
タグ: 知らんけど
posted at 22:44:28
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
マクロビなど節制こそが大事であると説く食養が子どもの成長を損なうのは、財政赤字を回避するための緊縮財政が日本の成長を損なう様子と似ていると思った。
成長には必要な器官に十分な栄養を供給することこそ大事であると思う。
タグ:
posted at 20:11:32
非公開
タグ:
posted at xx:xx:xx
ごまふあざらし(GomahuAzaras @MathSorcerer
@_cormullion 👍 It works ! Great!! pic.twitter.com/QX8UK3vKI9
タグ:
posted at 19:39:40
非公開
タグ:
posted at xx:xx:xx
#渡辺明 王将が勝利
【 第70期 #王将戦 七番勝負第2局 】
○渡辺明 王将 VS ●永瀬拓矢 王座
第70期王将戦七番勝負の第2局が1月23、24日(土、日)に大阪府高槻市の「山水館」で行われ、渡辺王将が120手で永瀬王座を破り、七番勝負の成績を2勝0敗としました。
buff.ly/3sPcNYF
posted at 18:56:17
@physics303 @MathSorcerer #Julia言語
TensorOperations.jl の使い方はこんな感じ。 pic.twitter.com/PsUAxkzwRd
タグ: Julia言語
posted at 18:20:17
@physics303 @MathSorcerer そうだ。言おうと思って忘れていたこと。
Einstein's notationでテンソル計算するには
github.com/Jutho/TensorOp...
が便利でかつ速いです。
@ ncon マクロを使えば可変ndimsの場合も楽に扱えると思います。 #Julia言語
タグ: Julia言語
posted at 18:05:17
aarch64上のJuliaでPlotsをインストールするときにハマった話とその解決法って需要ある?
#julialang pic.twitter.com/VelQnFcBrF
タグ: julialang
posted at 17:48:54
@physics303 @MathSorcerer 「素直な方法の直訳」を書くための部品を公式ドキュメントで探して使うとシンプルなコードになりやすいです。意外なほど部品が揃っている。
小ループの内側程度なら効率の良いコードを書く必要はない。
実質的な大ループの内側のコードだけは公式ドキュメントのPerformace Tipsに従う。
タグ:
posted at 17:17:23
@misaki70533995 少しモヤモヤしたから #LuaLaTeX に #luamplib 突っ込んで #MetaPost で同じ事やってみる... 此れなら #Julialang は要らんのだよなぁ... pic.twitter.com/5M0acSCzyr
タグ: Julialang LuaLaTeX luamplib MetaPost
posted at 16:54:50
こいつのリツイートで、シュムちゃんがオーストリアの財務大臣や銀行頭取やってたのを知らんのかバカめ、みたいなのがきたけど、我が国の論外クラスの連中みても、そんな肩書きでは何もわかんないのは自明だと思う。ちなみに、シュムちゃんはいずれも無能ですぐクビになってる。
タグ:
posted at 16:40:55
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
そしてこんな感じで表示が止まってるのやけどimportは終了なのやろか... とかたまに使う身としては毎度判らんようになって困るの... twitter.com/tsatie/status/... pic.twitter.com/6NmcP0lZKe
タグ:
posted at 15:22:37
@physics303 @MathSorcerer #Julia言語 ととのいました♫
B[circshift(~)...] = v の行がポイント
nd = ndims(A)
for i in 1:nd
v = vec(sum(A, dims=[1:i-1; i+1:nd]))
B[circshift([:; ones(Int, nd-1)], i-1)...] = v
end
B
gist.github.com/genkuroki/57ac... pic.twitter.com/NKMwhU9C6B
タグ: Julia言語
posted at 15:01:11
非公開
タグ:
posted at xx:xx:xx
#統計 警告:ど素人なので、jackknife法やbootstrap法の実装が標準的なものと違う可能性があります。1つ前のツイートの結果は私の実装の場合の結果です。そこからもリンクをはったように(添付画像にも全コードが示されている)、以下の場所でソースコードを読めます。
gist.github.com/genkuroki/efb0... pic.twitter.com/CLi4Q2aWOR
タグ: 統計
posted at 14:13:14
#統計 コンピュータで計算してみた。
⑤中央値の推定にjackknife法を使ってはいけない。
添付画像1のようにjackknife法の95%信頼区間に中央値の真の値が含まれる確率(標本分布に関する確率)が70%程度になっていたりする。bootstrap法ならそう悪くない。
gist.github.com/genkuroki/efb0... pic.twitter.com/xGh1PQJZpA
タグ: 統計
posted at 14:06:46
addコマンドのレジストリ更新がとろいのでctrl-cでキャンセルしたらGeneralレジストリがぶっ壊れたようだ
#julialang
タグ: julialang
posted at 13:23:49
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
julia言語
for文使って,綺麗にしたい.どうしたらええ?
Not(i) は using InvertedIndices の機能.Tとetaは大きさが同じ配列. pic.twitter.com/PSC4s1zkcH
タグ:
posted at 11:47:48
富谷(助教);監修 シン仮面ライダー @TomiyaAkio
統計自信あるニキ/ネキに質問したいんだけど、誤差の推定に使うJackknife法とbootstrap法、それぞれの長所と短所って何なんでしょうか
(Jackknifeが便利というのはおいておいて)
タグ:
posted at 11:13:42
ごまふあざらし(GomahuAzaras @MathSorcerer
fontface("Zapinfo") を使うと Luxor.jl での文字出力がめっちゃカッコよくなる.
#JuliaLang
#Julia言語 pic.twitter.com/Urf4mATw4T
posted at 11:00:26
#統計
* Stein推定では、調節された事前分布+MAP法の平均二乗誤差が最尤法よりも小さくなるという話。推測の誤差を小さくするために事前分布を利用できることは現代では常識の1つになっている。
nbviewer.jupyter.org/github/genkuro...
タグ: 統計
posted at 10:57:10
#統計
* 各標本分布の型ごとに、最大期待損失を最小化するために事前分布+ベイズ法を使える(このスレッドの話題)。
* ベルヌイ分布モデルにおいては、片側検定のP値と事後分布での確率がぴったり誤差無しで一致するように事前分布を設定できる。
nbviewer.jupyter.org/gist/genkuroki...
タグ: 統計
posted at 10:54:21
うちのおじいちゃん(97歳)、ボケ防止で数学解いたり法律の勉強したりしてて、たまに質問してくるのだけど、いろんな意味で難しいんだよね。。 pic.twitter.com/QSwqBytinM
タグ:
posted at 10:42:48
#統計 しかし、それだと、
どんなにダメな考え方の持ち主であっても、その人の主観内で期待損失を最小化していれば合理的だとみなす
という意味での合理性しか扱えない。
GelmanさんとShaliziさんは伝統的な主観主義ベイズ哲学は統計学にとって有害だったと言っている。
twilog.org/genkuroki/sear...
タグ: 統計
posted at 10:38:26
#統計 理学部数学科出身でも騙され易くなる理由の1つは、ベイズ統計の数学が結構複雑で難しいせいで、歴史的理由で不幸にも教科書に書かれるようになってしまった方針以外の道が見え難くなってしまうことです。
コンピュータの助けを借りながら、自分の頭で考えないとまともな理解は不可能。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 10:38:23
俺はシミュレーションをPythonなりJuliaなりでやった人が、計算結果をわざわざCSVで吐き出してOriginに読み込ませてグラフを修正するという作業は無駄かと思うので、直でmatplotlibとかでやりたい派ですわね。
タグ:
posted at 10:09:17
juliaだとこれがよさそうだな
1から始める Juliaプログラミング 進藤 裕之 www.amazon.co.jp/dp/433902905X/... @amazonJPより
タグ:
posted at 10:03:49
Do #JuliaLang users document and (continuously) test their packages? I dived into the General registry. The answer is in giordano.github.io/blog/2021-01-2... (spoiler alert: they do!)
タグ: JuliaLang
posted at 10:02:34
このスレッドでは、事前分布とベイズ法を最大期待損失を下げるための道具として使っている。しかも期待値は標本分布に関する期待値。
もちろん、これ以外にも事前分布とベイズ法の使い方がある。
タグ:
posted at 09:44:10
「事前分布と事後分布は主観確率の分布であり、信念の度合いを表す」という主張は、「50×100における50は単価を表し、100は数量を表す」という掛算順序指導に賛成の困った人達の主張に似ている。
事前分布と事後分布を信念の度合いだと解釈する必要はない。
ものすごく当たり前の話。
タグ:
posted at 09:39:55
「事前分布と事後分布は主観確率の分布であり、信念の度合いを表す」だとか、「標本分布を使うのは頻度主義であり、ベイズ統計ではない」だとか、「ベイズ統計ならばP値ハッキングの問題がなくなる」というようなおかしなことを言っている主観主義ベイジアン達は私には単なるアホにしか見えない。
タグ:
posted at 09:36:42
数の使い方や解釈が多彩であるのと同様に、事前分布の使い方も多彩でよい、と言えばわかり易いかな?
「Aという考え方もあるし、Bという考え方もあるし、Cという考え方もある。目的に応じて使い分けよう」という穏健な立場が数学に応用では普通だと思います。
普通から離れるとおかしくなる。
タグ:
posted at 09:22:47
#統計
標本分布がn=120の二項分布で、推定法が事前分布Beta(a, a)のベイズ法で、損失がKL情報量で定義されている場合の期待損失の動画が以下のリンク先にあります。
損失がKL情報量の場合の最尤法a=0では、真の成功確率qが0または1からちょっとだけ離れた値のときに期待損失が非常に大きくなる。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 09:08:46
キューブワールド2地点間の往復最短ルート
www.quantamagazine.org/the-crooked-ge...
立方体の測地線(展開図で2点間を直線で結ぶと作図できる)。一見シンプルに見えるけど奥が深くて「十二面体には他の頂点を経由せずに元の地点に戻れる経路あるか」問題がつい最近になって証明されたといった話
タグ:
posted at 09:03:11
#統計 損失を二乗誤差ではなく、KL情報量で定義した場合の同様の計算については
nbviewer.jupyter.org/gist/genkuroki...
および以下のリンク先を参照。
試行回数nを固定する場合にはa=0.5のJeffreys事前分布が出て来ます。小さな成功回数kを固定する場合にはa<0.5となり、k→∞でa→0.5となるっぽい。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 08:41:46
#統計 データの取得法によって、標本分布は変化するので、期待二乗誤差も変わります。だから、最適な事前分布も変わることに注意。
これは、事前分布を「成功確率に関するデータ取得前の信念」と解釈しなければいけないとする頭の固い主観主義ベイジアン達には受け入れられないことでしょう。 pic.twitter.com/k8FVs8ntJI
タグ: 統計
posted at 08:36:36
#統計 試行回数nを固定せずに、成功回数kを固定してちょうどk回成功するまでに行った試行回数nがデータの場合に、同様の方法で最大期待二乗誤差を最小化するaも求めてみました。
添付画像はk=8の場合です。a=1.63程度で最小化されている。 pic.twitter.com/0uae0viZgz
タグ: 統計
posted at 08:32:32
#統計 真の成功確率qが0または1に近い場合には、推定値pも0または1に近くなりやすく、二乗誤差(p-q)²の期待値も小さくなる。
大きな損失が起こり易いのは真の成功確率qが0と1から離れている場合で、そういう場合の損失を減らすにはa=√(n/4)のBeta(a, a)を事前分布に採用すればよいらしい。 pic.twitter.com/1Lrl2RBniB
タグ: 統計
posted at 08:16:13
#統計 事前分布Beta(a, a)のa=0のimproper事前分布の場合が最尤法になります。この手のゲームで最尤法は最適戦略になり難いですね。
添付画像を見ると、n=25のとき、最大期待二乗誤差を最小化するaは2.5になっています。
一般に√(n/4)がそのようなaになっているようです。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/avr1QEWIuE
タグ: 統計
posted at 08:11:52
Great story on @ClimateMachine and #julialang, featuring @AlanEdelmanMIT in @csmonitor!
"There was no way we could have done it with another language" - Prof. Ferrari at MIT. "Julia paid off for us better than they would have imagined" - Prof. Schneider at Caltech. twitter.com/csmonitor/stat...
タグ: julialang
posted at 08:10:53
#Julia言語 #統計
真の成功確率qのベルヌイ試行n回分のデータから、成功確率の推定値pを事前分布Beta(a, a)のベイズ法で求め、二乗誤差(p - q)²に比例した罰金を支払う
というゲームで期待損失の(qを動かす)最大値を最小化する計算をやってみた(パラメータaを最適化)。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/R2V7LiVSFR
posted at 08:06:48
見返してたらJuliaが多値を返せるので、返り値がエラーコードにとられちゃう問題が避けられる事に言及するのをすっかり忘れてました。近いうちに追記しておきます。
タグ:
posted at 07:59:04
#Julia言語
github.com/genkuroki/Isin...
では、2D Ising モデルでの周期境界条件を ifelse で処理した場合と「複数のforループのべた書き」で処理した場合の速度の比較のために、Holy trait パターンっぽい書き方をしています。
同じ問題を異なるアルゴリズムで解けるようにしておくときに便利。
タグ: Julia言語
posted at 07:33:57
#Julia言語 関連
Holy traitには、問題の型ごとにデフォルトの解法アルゴリズムを変えるというような使い方もできる。
Juliaのコードを読むときに、Holy traitについて知らないと困る場合が結構あります。
何も知らないと「中身が空の型が定義されまくっているけど、これなに?」となる。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 07:26:26
@michikusa_heya twitter.com/sekibunnteisuu...
こう書いたけど、中にはなかなか気づかないケースもあって、悩むこともありま
す。
なるべく教えないで自分で気づいてほしいけど、生徒によっては教えることもありますね。
タグ:
posted at 07:20:12
#Julia言語 Holy traitの解説
Juliaでは抽象型と具体型の階層が木構造になってしまうのですが、Holy traitを使えば「木構造しか不可」という問題を回避してコードを書けます。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 07:15:51
@genkuroki このパターンは是非取り入れたいと思ってます。
毒攻撃はボスと雑魚で効き目を変えたかったりするので、そういったあたりですかね
パラメトリック型もややこしいので当初触れない予定だったんですが、〇〇キラー攻撃が〇〇族に効く、みたいな処理を一元化するのに良いなあ、とか思うと終わらないです笑
タグ:
posted at 07:13:23
「これまでの知識重視から思考力がどーたら」というおざなりな記事じゃなくてちゃんと問題を分析していて珍しいな、と思ったら、清 史弘さんの記事だった。
twitter.com/f_sei/status/1...
タグ:
posted at 06:56:03
大学入学共通テスト、思考力は測れたのか?予備校講師が検証(HARBOR BUSINESS Online)
#Yahooニュース
news.yahoo.co.jp/articles/bb2b3...
タグ: Yahooニュース
posted at 06:54:13
#Julia言語 #Jupyter #ZeroTier
私はZeroTierで自宅サーバーに外部の端末から接続しています。
①ZeroTier www.google.com/search?q=ZeroT... を使えるようにする。
②サーバーにZeroTierをインストールして設定
③端末側にもZeroTierをインストールして設定
端末からZero TierのIPアドレスにアクセス twitter.com/oppo89694572/s...
posted at 06:45:37
科学計算が得意な julia っていう言語が生えてきたのかー
確かに計算量の多い分野だと実行速度速いほうが何かと良いですよね
project.nikkeibp.co.jp/idg/atcl/19/00...
タグ:
posted at 06:44:59
#Julia言語 情報
Juliaでは Tim Holy's trait trick もよく使われていて、他人が書いたパッケージのソースコードの解読でも必須の予備知識になっています。RPGとも相性が良さそう。
この辺の話も入れるとさらに大変な分量になりそう。
ahsmart.com/pub/holy-trait...
github.com/mauro3/SimpleT... twitter.com/muuumin20/stat...
タグ: Julia言語
posted at 06:19:32
道草学習の部屋@オンライン家庭教師&教育 @michikusa_heya
@sekibunnteisuu 生徒によってはなかなか自ずから気づいてもらえず、時々待てなくて、ヒントを出すか、教えてしまうことがあり、反省するところです。
タグ:
posted at 04:06:38
ちなみに今回の記事ではスキルが一つしか実装できてなくて、まだまだ先は長い、と言うか、終わりが見えなくて泣きそうになってます。
まだ書けてない事の方が多いんじゃないか?全5回くらいで終わらせるはずだったのに。
どうして、どうしてこんな事に・・・
タグ:
posted at 03:32:03
「Julia言語で入門するプログラミング」第6回を書きました。コンストラクタ、抽象型と具体型、型の階層構造、例外処理、といったあたりの話題に触れています。
はてなブログに投稿しました #はてなブログ #Julia言語
Julia言語で入門するプログラミング(その6) - SE教育パ…
muuuminsan.hatenablog.com/entry/2021/01/...
posted at 03:22:48
tau18analytics @tau18analytics
@OPPO89694572 great ! ...
there are also hosted environments
(which are also customizable), so that
your home PC is not required (if you
prefer) ...
i have had great experiences with cocalc.com ... @cocalc_com ...
タグ:
posted at 03:04:56
ごまふあざらし(GomahuAzaras @MathSorcerer
Drawing a Reuleaux triangle using Luxor.jl
#Julialang #Julia言語 pic.twitter.com/DQ0Rg1gRNl
posted at 01:19:33
りそな銀行が明細の表示期間を過去無限に引き伸ばしたらしい。これは大変素晴らしいニュースだと思う。僕はこれまで、過去の明細を遡れない銀行の仕様について大いに不満を抱いてて、ことあるごとにツイートしてきた。ゆうちょ・三菱・住友・みずほも早急に対応してほしい www.nikkei.com/article/DGXZQO...
タグ:
posted at 01:15:06
@fartelengelbert @nick_mayhall @Viral_B_Shah @NJIT Numba is great, but julia is far better. Numba supports only small part of functions. Overhead between ei. pandas and numba is very annoying. In #JuliaLang there is not such issue.
タグ: JuliaLang
posted at 00:33:54
ケインズやアダムスミスは、「うんそうだよな」「だいたいそうだよな」とほぼ納得しながら読めるが、シュムペーターって「いやそんなことないんじゃね?」「おまえ、技術ってなんだか知ってるの?」みたいな首傾げだらけで先に進む気がなかなか起きない
タグ:
posted at 00:32:43
本当は、「一部の長老とその取り巻きが中心で運営している」ではない、学会のようなものを作りたいとは考えています。良いものはよい、悪いものは悪いと言える学会を。
文科省が誤れば、それを提言できるような学会を。
タグ:
posted at 00:10:44
中学・高校の教員、塾・予備校の先生達で情報の共有をしていきます。閉じた世界ですが、その今あげた関係者であればwelcome です。閉じた世界なので、内々での資料をコピーして外部に出さないことを条件に公開していきます。
タグ:
posted at 00:06:45
最初のツイートで、教材、書籍、資料等の検討はFBでと書きましたが、2月に入ってから方向性を考えてみようと思います。参考書のここが変じゃない?のようなものも集めていきますが、例えば、参考書の一部分がおかしくても、それを全体に波及させたり、著者への攻撃にならないようにする世界を作ります。
タグ:
posted at 00:03:44