黒木玄 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 を見て下さい。
2019年12月27日(金)
@MathSorcerer ごまふあざらしさんのツイートで気づいたんですが，Dataframeを直接受け取れるのが珍しいものなんでしょうか？
ggplot2ライクに描画できるのがGadflyの魅力と聞いてたので確かにdfを受け取っていろいろできるのはggplot2的なんですが，他の描画系にはないんです？
タグ：
posted at 00:07:17
ごまふあざらし(GomahuAzaras @MathSorcerer
@389jan 他の描画系でDataFramesのオブジェクトを受け取れるかは未調査だったのでGadfly良さげでは？と思ってしまったのです．
ggplot2 は詳しくないです．ごめんささい．．．
タグ：
posted at 00:11:08
OokuboTact 大久保中二病中年 @OokuboTact
読書日記
『小学校算数科の精選・集約』 (1976年)
編集：川口廷、片桐重男、杉岡司馬、窪田騰
執筆者は当時の教師陣。
興味深い内容
＃超算数
タグ： 超算数
posted at 00:16:16
OokuboTact 大久保中二病中年 @OokuboTact
昭和５１年（１９７６年）と「現代」がどう繋がるのかわからないけど、
算数の「ゆとり教育」の発端なのかもしれない。
＃超算数
タグ： 超算数
posted at 00:23:31
@MathSorcerer @389jan #Julia言語
Windows 8.1環境で、Pythonさん御一行は~/.julia/conda以下にあります。スクショなしで申し訳ないのですが、Jupyter notebook上で
using Gadfly
plot(x=rand(10), y=rand(10))
でグラフが普通に表示されました。マウスで動かせますね！
タグ： Julia言語
posted at 00:38:12
OokuboTact 大久保中二病中年 @OokuboTact
１９７０年代〜８０年代は、「ゆとり教育」のはじまりであるけど、トンデモ算数にとっても重要な育成時代という気がする。
塾VS学校の戦いが激化したことと関係があるのかもしれない。
タグ：
posted at 00:38:49
ごまふあざらし(GomahuAzaras @MathSorcerer
飼い主， Windows 8.1 でもちゃんと動くの素直にすごいと思ってる twitter.com/genkuroki/stat...
タグ：
posted at 00:43:06
@MathSorcerer @389jan #Julia言語 Daraframeのプロットには
github.com/JuliaPlots/Sta...
の @ df マクロがちょっと便利です。@ df マクロはプロット専用ではなくて、Dataframeの中身を取り出すこと一般に使用できます。
タグ： Julia言語
posted at 00:43:42
TeX & LaTeX Advent Calendar 2019
adventar.org/calendars/4087
SATySFi Advent Calendar 2019
adventar.org/calendars/3929
『数式組版』を組む技術 Advent Calendar 2019
adventar.org/calendars/4719
数値計算 Advent Calendar 2019
qiita.com/advent-calenda...
タグ：
posted at 00:48:22
ごまふあざらし(GomahuAzaras @MathSorcerer
@genkuroki @389jan dfマクロしゅごい． pic.twitter.com/wbYCHJ15d8
タグ：
posted at 00:55:40
@MathSorcerer @389jan #Julia言語
using Plots
plotlyjs()
plot(sin)
などをJupyter notebookで表示するためには
juliagizmos.github.io/WebIO.jl/lates...
の注意に従う必要がありました。plotlyjs()はマウスでグリグリできるプロットを表示してくれます。
Gadflyでも同様の問題がないのかしら？と思いました。
タグ： Julia言語
posted at 01:04:58
@MathSorcerer @389jan #Julia言語 リンク先を間違った。マウスでぐりぐりなPlotlyjs関連で必要だった設定は
juliagizmos.github.io/WebIO.jl/lates...
の方に書いてあります。
タグ： Julia言語
posted at 01:07:30
私は学生時代の教員が良かったのでしょう。幸いこの理解には至ったつもりですが…
数学モデルの上でも信頼区間の95という数字が意味する心は良く分からないし、しかも現実の母集団が一般に厳密には正規分布に従わないので、95という数字に何の意味があるのかまるで分からないんですよね。 twitter.com/genkuroki/stat...
タグ：
posted at 01:57:33
何か尺度を定義しておく程度でしかなく、「経験的に95だと結構良い推定ができる」みたいな経験値が業界毎にあるだけなんでしょうか。（実際は業界内でも扱う母集団がしばしば異なるはずですけど）
タグ：
posted at 01:58:29
Celeste.jlとか66万5000CPUコア使ってるって言うし，普通にスパコンでJulia使ってるプロジェクトだと思うんだけど何か違う？
github.com/jeff-regier/Ce...
タグ：
posted at 02:10:44
> On a supercomputer, the VI procedure efficiently uses 665,000 CPU cores to construct an astronomical catalog from 50 terabytes of images in 14.6 minutes, demonstrating the scaling characteristics necessary to construct catalogs for upcoming astronomical surveys.
タグ：
posted at 02:13:09
Approximate inference for constructing astronomical catalogs from images
projecteuclid.org/euclid.aoas/15...
タグ：
posted at 02:13:44
非公開
タグ：
posted at xx:xx:xx
ODE Solver Multi-Language Wrapper Package Work-Precision Benchmarks (MATLAB, SciPy, #julialang, deSolve (R)) by @ChrisRackauckas
benchmarks.juliadiffeq.org/html/MultiLang... pic.twitter.com/E16nf7GSrY
タグ： julialang
posted at 02:31:19
非公開
タグ：
posted at xx:xx:xx
New Julia Benchmarks: @ChrisRackauckas published new reproducible benchmarks demonstrating Julia’s superior performance for ordinary differential equations (ODEs) compared with Python, Matlab and R.
benchmarks.juliadiffeq.org/html/MultiLang...
#julialang
タグ： julialang
posted at 07:49:28
食べ物を残したときに罪悪感を抱くのは健全ですが、無理して食べると高血圧、糖尿病、痛風、高脂血症になるので、躊躇なく残しましょう。反省と後悔を、次から作りすぎないように、と前向きな意識に転換しましょう。これが正しいのに、反省が面倒だから食べちゃう。困ったもんだ。 twitter.com/minesoh/status...
タグ：
posted at 08:24:39
#Julia言語
github.com/chakravala/Red...
Reduce.jl
を試してみました。Jupyter notebook上ならばMathJaxで数式をきれいに成形してくれるようになっていました。ただし、
Base.convert(::Type{Any}, x::RExpr) = x
としておく必要がありました。バグかも。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/xjJ5wiEO9i
タグ： Julia言語
posted at 08:30:05
Fredrik Bagge Carlso @baggepinnen
@JuliaComputing @ChrisRackauckas "SciPy Omitted due to failures at higher tolerances and because it's too slow to finish in a day!" hard to believe scipy still can't solve diffeqs 2019 🤣
タグ：
posted at 08:57:18
Want to learn how to build a network graph in Tableau? 😎 @ladataviz shows you how in three steps: tabsoft.co/2Z3WicT pic.twitter.com/ktNWF9bGPA
タグ：
posted at 09:00:31
#統計 補足
現実の未知の母集団分布については、
繰り返し多くの異なった標本について95%信頼区間を何回も計算した場合、真の値を区間内に含む95%信頼区間の割合は95%になる
という主張はナンセンスです。しかし、フィクションである数学的モデル内部では概ね成立しています。 pic.twitter.com/teqaaunTAo
タグ： 統計
posted at 09:06:15
#統計 東大出版会の『統計学入門』を見ると、「割合」の語が太字になっており、まるで、「確率」ではなく「割合」と言わなければいけないかのような印象を与える説明の仕方になっています。
これは統計学のイロハを理解していればありえない説明の仕方だと思います。 pic.twitter.com/1sLurmocMr
タグ： 統計
posted at 09:14:16
Dr. Chris Rackauckas @ChrisRackauckas
@baggepinnen @JuliaComputing Indeed, I was surprised when I found that, so I isolated out and found this code as a simple reproducer:
github.com/JuliaDiffEq/Sc...
Note that the problem it's failing on, ROBER, is the standard stiff ODE test problem since GEAR demonstrated you can solve stiff ODEs with BDFs.
タグ：
posted at 09:29:31
Re:RTs #Julia言語
JuliaのDifferentialEquations.jl(微分方程式の巨大パッケージ)は非常に便利。
「ボールの跳ね返りの瞬間」の記述もできる。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 10:02:56
#Julia言語 数値計算の背後にある数学は自分で実装してみた方が理解が進みやすいので、便利なパッケージがあっても、自分で実装してみることは大事なことだと思います。
Juliaでオイラー法を自分で書いて求めた近似解とDifferentialEquations.jlによる結果を比較してみたりすると勉強になるかも。 twitter.com/kado_judo0312/...
タグ： Julia言語
posted at 10:07:25
非公開
タグ：
posted at xx:xx:xx
関連 #Julia言語 の DifferentialEquations.jl は非常に優れたパッケージ。ソースコードを閲覧する価値もあると思う。
twitter.com/juliacomputing...
タグ： Julia言語
posted at 10:17:22
ぶっちゃけ、既存のパッケージを使って「〇〇できる」という段階に到達することは大したことがない話で、「それだけではない何か」について学ばなければ新しいものに手を出す価値は低いと思う。
タグ：
posted at 10:20:25
マイナスイオンの顛末を知らんのか、としか。権威：ノーベル賞学者のレナード、大手企業が商品化、日本国民の大多数が購入、でも結局マイナスイオンの化学種すら特定されず、濃度と効果の関係もまともなものが出なかった。 twitter.com/qqta5ydd/statu...
タグ：
posted at 10:23:11
STAPが否定されたのは、追試でどこも再現できず、本人のやらかした内容から捏造が強く疑われたから。そういう意味では同様の状態ですね。設置後1年後検査のデータ操作が強く疑われる状態です。偶然成立したブラインドテストに対する社長の発言によって確信が持てました。 twitter.com/qqta5ydd/statu...
タグ：
posted at 10:28:16
Dr. Chris Rackauckas @ChrisRackauckas
@moorepants @baggepinnen @JuliaComputing Indeed odeint doesn't seem to fail. When using Numba to accelerate it, rough estimates show it 500x slower than Julia, putting odeint in last place out of those that finish, but only slightly behind deSolve. Details:
gist.github.com/ChrisRackaucka...
We can add it to the benchmarks.
タグ：
posted at 10:34:31
Formatted string interpolation - Usage / First steps - JuliaLang discourse.julialang.org/t/formatted-st...
結局Juliaに%2f的なのを@sprintf以外でやる方法ってないのかな
タグ：
posted at 10:50:56
ごまふあざらし(GomahuAzaras @MathSorcerer
Luxor.jl で生成した gif の描画をJunoでできた．
- discourse.julialang.org/t/base-show-an...
- cormullion.github.io/pub/2019-06-23... pic.twitter.com/0RPyctj50Y
タグ：
posted at 11:10:39
ごまふあざらし(GomahuAzaras @MathSorcerer
ごま「このスクリーンショットの解説最高にわかりやすいでキュー」 twitter.com/genkuroki/stat...
タグ：
posted at 11:18:04
ごまふあざらし(GomahuAzaras @MathSorcerer
tspan is 何だったけれど微分方程式を解いている部分と描画部分別ブロックだった模様．
ということで動いたできゅ！の結果です． pic.twitter.com/PpsbUHOBVw
タグ：
posted at 11:35:05
#統計 ベイズ統計では「データは定数になる」というようなことが書いてある「統計学ハウツー本」がありますが、ベイズ統計であろうがなかろうが、無作為抽出で得たデータは運悪く偏っている可能性があるので、データは確率変数扱いしないとダメです。
twitter.com/genkuroki/stat...
タグ： 統計
posted at 11:43:10
そうすると、e-ポートフォリオを運営する一般社団法人教育情報管理機構もベネッセの傀儡であることが強く推認されますね。 twitter.com/y__hiroyuki/st...
タグ：
posted at 11:43:26
「オカンがな、掛算順序強制の数学的な根拠を忘れたらしくて」
「特徴教えて」
「なんか、右から係数を掛ける代数系らしいねん」
「ほな加群や！右から掛けんのやったら右加群や！」
「でも左からかけたらアカンらしいねん」
「ほな加群とちゃうなあ。整数環上の右加群は自明に左加群にもなるしな」
タグ：
posted at 11:51:32
#Julia言語 での algebraic data types の実装例
discourse.julialang.org/t/mlstyle-jl-a...
github.com/thautwarm/MLSt...
タグ： Julia言語
posted at 13:23:17
#Julia言語
github.com/thautwarm/MLSt...
【Performance Gain
When dealing with complex conditional logics and visiting nested datatypes, the codes compiled via MLStyle.jl is usually as fast as handwritten code. You can check the benchmarks for details.】
作者のアイコンが！(笑)
タグ： Julia言語
posted at 13:30:31
非公開
タグ：
posted at xx:xx:xx
@chiyosuke1002 @kaijikaito 失礼します。ちよ助さんにお伺いしたいのですが、面積なら順序はどちらでもいいということなのでしょうか？
平行四辺形は底辺×高さ、高さ×底辺、どちらでもいいのでしょうか？
タグ：
posted at 13:40:23
@baygon27288632 私は中高で理科を教えていますが、モルや濃度の概念や複数のパラメーターが絡むと混乱する生徒の中にも、順序指導やみはじ指導をはじめとする「計算や数値の意味より公式への代入」の影響で混乱していると考えられる生徒がいます。
タグ：
posted at 14:08:03
Amicitia sal vitae. @boots_fleck
これ、体重を乗せて背中をバネにして前に飛ぶように放つという、ロングボウ射手独特の射ち方そのものでちょっと動きがぎこちないけど再現度に驚いた。
この射ち方を習得した射手が、59kgのロングボウで約45度の角度で放つと、約300ヤード先（最大射程）まで飛ぶということである。 twitter.com/heikihenken/st...
タグ：
posted at 14:18:33
Dr. Chris Rackauckas @ChrisRackauckas
@moorepants @baggepinnen @JuliaComputing #scipy odeint did better than the quick benchmarks suggested. It's about 50x slower than #julialang, which makes it around #matlab and faster than R's deSolve. Reproducer:
github.com/JuliaDiffEq/Sc... pic.twitter.com/BUN4g6olpJ
posted at 14:46:15
Dr. Chris Rackauckas @ChrisRackauckas
@moorepants @baggepinnen @JuliaComputing Note that MATLAB is not in these plots right now because I need to fix some licensing issues. But the code for that is just commented out, so you can run it locally, and it'll be in the same spot as it is in the other plots. We'll update the full benchmarks soon.
タグ：
posted at 14:48:04
Interesting GIF but the strangest analogy. Google translates the caption as follows: "Indeed, rotation stops when there is something like anal beads" twitter.com/jagarikin/stat...
タグ：
posted at 15:19:23
#統計
頻度論 ベイズ
パラメータ 定数 確率変数
データ 確率変数 定数
という図式でベイズ統計について解説する困った人達(ググると見付かる)の話です。
頻度論とベイズで扱うデータは同じものなの前者でのみ確率変数になるらしい。ひどすぎ。
twitter.com/ML_deep/status...
タグ： 統計
posted at 16:46:03
#掛算 出尽くしてるだろうと思ってたが、これはw
さすがに初出ではなかろうか＞不完全性定理 twitter.com/kamo_hiroyasu/...
タグ： 掛算
posted at 17:31:14
ごまふあざらし(GomahuAzaras @MathSorcerer
．jmd で作成したJuliaつき文書をHTMLで変換しそれを別のサイトに移植できるか？ってのを動かしてたら
github.io と はてなブログにHTMLをそのままコピーするという方法はできた．
タグ：
posted at 17:37:37
@phys_suomusu ここで、勘違いしないで欲しいことは、「温度の電子状態への影響」を考慮することが困難なのであって、「温度の原子の運動への影響」は、第一原理分子動力学法により、考慮することが可能であるということです（今の場合、Born-Oppenheimer近似により、原子核と電子の運動を分離できると仮定している）
タグ：
posted at 17:55:35
せめて最初ぐらい「これは出来レースではないかも」と思わせてほしかったが、やはり予想通り。英語については民間試験利用以外の方法を考えるつもりはない模様。そうでなければ、評価だけでなく、現場に詳しい英語教育の専門家を最低一人は入れるはずだけど...。 twitter.com/kiitakashi/sta...
タグ：
posted at 18:14:07
#統計 ベイズ統計の解説が悲惨なことになっていることをググって確認したければ
www.google.com/search?q=%E9%A...
からリンクをたどるといいかも。 pic.twitter.com/PIQcx5JmOs
タグ： 統計
posted at 18:55:31
へへ熊(固定にthreads貼ってます) @kuma_256bit
暗黙の型宣言さん、冬コミ情報を何も発信してないっぽいな…サークル名すら分からん…隣がjulialangって言ってたからそこから探すか…
タグ：
posted at 19:07:06
大学入学共通テスト国語・数学の記述式問題の採点をすることになっていた学力評価研究機構の封筒。
私に対する抗議の手紙が入っていました。
電話番号がありません。 pic.twitter.com/FQhNc6TZs5
タグ：
posted at 19:15:07
うわあ、これは書いてて恥ずかしくないんかなぁ。ゲーデル…分かってないどころか名前どこかで知っただけやろ、この書き方（苦笑 twitter.com/adahuouoh/stat...
タグ：
posted at 19:47:49
今回の冬コミ(4日目南リ38a)新刊は「プログラミング言語Julia」 です。今までのやつに加筆（主に並行・並列計算）・修正（主にv1.3 対応）しつつ、（長かったので）タイトルを変えました。
webcatalog.circle.ms/Circle/14830757
タグ：
posted at 19:52:19
以下のPDFファイル（35頁，862KB, Google Drive）にて，カイ2乗検定にて，自由度が一元表のときはr-1，二元表のときは(r-1)×(c-1)になる理由を考えてみました．
drive.google.com/open?id=1WKvcy...
まあ，よく知られたことだとは思いますが...．
きっと，もっと短く説明できるのでしょけれども...無駄に長い...
タグ：
posted at 20:20:09
藤井聡太七段は王位戦予選トーナメント決勝で斎藤慎太郎七段に１２１手で勝ち、初の挑戦者決定リーグ入りを決めました。リーグ戦は紅白２組に分かれ、各６人の棋士が総当たりで対戦し、優勝者同士が木村一基王位への挑戦権をかけて戦います。
タグ：
posted at 20:20:50
#Julia言語 統計学での2×2の分割表に関するシミュレーションの計算で Memoize.jl がめっちゃ役に立ちました。さらに最近のバージョンアップ、どうもありがとうございます。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/hymd_mathemati... pic.twitter.com/uYeLFqQkYT
タグ： Julia言語
posted at 20:22:02
これまで勝てなかった斎藤慎太郎七段を破ってのリーグ入りはすごいですね。王将戦リーグに続いて、また観る将の楽しみが増えます。
本田奎五段と佐々木大地五段も、これからリーグ入りをかけた予選決勝ですね。 twitter.com/mainichi_shogi...
タグ：
posted at 20:24:06
ごまふあざらし(GomahuAzaras @MathSorcerer
#Julia言語
リモート用としてSlackも用意しています。 twitter.com/mathsorcerer/s...
タグ： Julia言語
posted at 21:10:38
Juliaでカルマンフィルタを実装してみた
#Julia言語
#カルマンフィルタ
qiita.com/fred55/items/9... #Qiita
posted at 21:20:07
julia言語 の DifferentialEquations.jl は
callback で condition() 指定すると、condition()のゼロクロスを検出して追い込んでくれるのか。便利そうですね。アッセンディング/ディッセンディングの区別もしてくれるみたい。精度が必用なら、追い込むロジックは要確認と。
タグ：
posted at 21:31:00
たまには宇宙っぽい Tweet を
軌道設計/解析の場面で、遠日点や近日点まで伝搬したいとき、動径 R の時間微分のゼロクロスを検出する方式を採用すれば、condition() のような汎用的な仕組みを適用できそうですね。
タグ：
posted at 21:31:00
直近の３ステップの動径 R を使って極大/極小を検出する方式は、特別扱いが多くなって見た目も使い勝手もちょっと残念な感じになりがちです。汎用的な軌道伝搬ライブラリーを設計するとき、ちょっと検討してみると良いかもよ(無保証)。
タグ：
posted at 21:31:01
Rumpの例題をJuliaで計算してみた。
Float64の答えは-1.1805916207174113e21
BigFloatの答えは-0.8273960599...以下省略
Juliaすげー。
タグ：
posted at 21:55:44
#Julia言語
* めんどくさいのでシンプルなモンテカルロシミュレーションでどういう感じになっているか様子を見てみること
* 同じような計算結果を綺麗にプロットするための調性を試行錯誤で行うこと
などでも、メモ化
ja.wikipedia.org/wiki/%E3%83%A1...
は非常に便利でした。
タグ： Julia言語
posted at 22:02:56
コードはこんな感じ
f(a,b)=333.75*b^6+a^2*(11*a^2*b^2 - b^6 - 121*b^4 - 2) + 5.5*b^8 + a/(2*b)
a=77617.0; b=33096.0
f(a,b)
→-1.1805916207174113e21
a_bf = BigFloat(a); b_bf = BigFloat(b)
f(a_bf,b_bf)
→-0.8273960599...以下省略（おそらく正しい値）
タグ：
posted at 22:14:44
謎水トラブルを防ぐのに科学の知識は必要ありません
「配管を40年延命するという装置の効果確認が1年後だけ。機材がリースではなく買い切り。赤水が出ても保証されない契約内容」
という状況は、消費者委員の知見で十分に「怪しさ」が理解できるはずです
タグ：
posted at 22:15:36
特定の製品を否定できないと言うなら
「この点を業者に確認してください」
「契約書のここが不利に書かれている、メリットとデメリットをよく考えて判断を」
「この問題はこの有識者が詳しい」
など、正しい理解を広める手段はたくさんあるはず
タグ：
posted at 22:15:39
現状では、日本技術士会が駆け込み寺になって、会員が「ボランティアで相談に乗る」などと言う手段でしか被害が防げていません
ボランティアがプライベートを犠牲にしているのに、給料が出ている行政が見て見ぬフリという状況がなぜ発生するのか？
疑問を禁じ得ません
タグ：
posted at 22:15:40
@phys_suomusu 分子動力学（Molecular Dynamics, MD）法というのは、こんな感じのですね。MD法で、その分子の相図を書くことができます。また、その分子が、任意の圧力と温度で、理想気体からどれだけずれている（PV/(Nk_BT)からのずれ）かも、調べることができます。
twitter.com/dc1394/status/...
タグ：
posted at 22:20:45
今日小5にこれやってもらったら受けてたみたいでよかった（笑）。数人に試した結果だが、なんかこれだと関係がよく見えるみたいなんだよなー。なん倍ですか、より。名前つけると把握しやすくなるのか。毎回つけるわけにもいかないが（笑） pic.twitter.com/VcSSG0uE66
タグ：
posted at 22:51:12
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数 ＃歴史的資料
１９７８年の資料から
「数と計算の具体的事例「余りあるわり算」第３学年を中心にして」 pic.twitter.com/zCnccO5nGu
posted at 23:35:47
