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

おじさんは昔からJuliaばっかり書いてるんだよ

タグ：

posted at 23:59:09

ceptree @ceptree

Juliaに対する思いを吐露するTL

タグ：

posted at 23:58:22

計 @kei_p1

「でもそんなのやってる実験系の論文なんてないですよね？本当に必要だとは思えませんけど？？」みたいな返しをされた。はい、そうですか、お好きになさると良いと思いますという感じ。

タグ：

posted at 23:49:34

計 @kei_p1

先日は「ベイズ統計を使ってデータを纏める際の統計処理の正しさはR-hat値を使えばいいのか？」と訊かれたから、R-hat値はあくまでMCMCの収束判定の目安だからあなたが仮定した統計モデルの妥当性とはまた別で、そこは何かしらの手段で汎化誤差に関係する評価を行わないとわからないと答えたら、

タグ：

posted at 23:47:50

Massimo @Rainmaker1973

There's a waterfall in Romania whose aspect and colors make it look like it came out from a fairly tales book [read more, with more photos by Teea Cristea: buff.ly/2w5HGxg] pic.twitter.com/Giwg7sIOpV

タグ：

posted at 23:30:35

非公開

タグ：

posted at xx:xx:xx

2021年はJulia元年だから

タグ：

posted at 22:52:08

父さんな、Juliaコンサルとして生きていこうと思うんだ

タグ：

posted at 22:37:59

Makie @MakiePlots

Beautiful Makie examples by @LazarusAlon
lazarusa.github.io/Webpage/codeJu... pic.twitter.com/5YICaqReFm

タグ：

posted at 22:29:18

#超算数
いつぞやの、ネトウヨ公共教師の掛け算順序論者と対消滅してほしい。

タグ：超算数

posted at 22:07:28

Massimo @Rainmaker1973

This is how Justin and Jack, two employees working in the editorial department of the creative production company Wondros, made a giant commercial airline style paper airplane using the paper wrappers from company lunches on “Falafel Fridays" [read more: buff.ly/35e3yWM] pic.twitter.com/cmhBjeVTJg

タグ：

posted at 22:00:34

#超算数
トンデモ算数批判を差別に利用するクソ野郎ども
twitter.com/Hile285Banzai/... pic.twitter.com/YugkgZ0mG7

タグ：超算数

posted at 21:59:11

#超算数
トンデモ算数批判を差別に利用するクソ野郎どもtwitter.com/miyuabc123/sta... pic.twitter.com/cdagfj9yxc

タグ：超算数

posted at 21:57:41

#超算数
トンデモ算数批判を差別に利用するクソ野郎ども
twitter.com/AMATlls/status... pic.twitter.com/ikOIDWcpzx

タグ：超算数

posted at 21:56:27

非公開

タグ：

posted at xx:xx:xx

Minori @minotonefinland

ミリデコ@7.30Vショー&東京ゲームダ @Milideco_

そうなんです。しかもオシャレフード扱いです。パリ、ロンドン、コペンハーゲンのヒップスターたちがオープンするカフェやレストランの人気メニューとなっております。今これ作ったり食べたりしてると特定の界隈ではオシャレ人認定です。 twitter.com/aiamsterdam_nl...

タグ：

posted at 21:44:10

非公開

タグ：

posted at xx:xx:xx

至急求む！
なんでバツなんでしょう？
教えられない！w pic.twitter.com/ynGRMkq0pv

タグ：

posted at 21:08:11

#超算数なぜ彼らのことを私が知っているか、なぜ彼らが日本の状況に関与的だと思うのか。彼らを含めドイツ教育者について戦前の日本はかなり詳細に調べていたようなんです。

タグ：超算数

posted at 19:38:58

今井翔太 / Shota Imai@える @ImAI_Eruel

#超算数いちおう私の知っているところでは、19世紀に教会やラテン語を基準とする教育が、現代語による国民教育に転換した頃からすでにそういう傾向があります。言語表現から掛け算の式をつくるという順序を提唱したディースターヴェーク、作問を提唱したハルニシュの二人をドイツから紹介しています。

タグ：超算数

posted at 19:37:07

スタンフォード大学で使われている線形代数の応用に関する教科書の翻訳版が発売されるようです
既に海外では非常に評価が高い教科書で，AI・機械学習等，普通の線形代数の本では応用が見えにくい部分を学習できます
『スタンフォードベクトル・行列からはじめる最適化数学』
www.amazon.co.jp/dp/4065161967/... pic.twitter.com/uT5SSdPlWY

タグ：

posted at 18:50:45

#超算数場面、ことばの式、公式など、算数の理解の基盤として言語による理解を前提にしているとしか思えない概念が多い。この傾向は歴史的にどこまで遡るのか分からないが、戦後の算数教育界では普遍的と見ている。

タグ：超算数

posted at 18:28:00

#超算数算数の解答で一つの考え方しか認めない先生が多いという問題、指導要領の文言にも瑕疵がありそう。式が場面を表すというやつ。場面から式を抽出するとは読めなそうなのだ。

タグ：超算数

posted at 18:21:45

私が #超算数と呼称する算数教育界の不可解な指導法たちは、かならずしも日本固有の問題ではない。コモンコアのように日本が輸出した可能性が高い事例もあるが、日本とは無関係に掛順指導の国もある。ただし、諸外国には、日本の算数教育界が目を顰めるような*真っ当な*指導法も行われている。

タグ：超算数

posted at 17:43:32

#超算数前屋氏の記事は、芳沢氏の最近の記事への補足および近著の宣伝なのだろう。芳沢氏の言葉を引用して【「…ところが日本では、証明問題までを穴埋め式で答えさせている。暗記を試験しているにすぎない。『日本の常識は世界の非常識』でしかない」】とする。だが、くもわは米国由来かもしれない。

タグ：超算数

posted at 17:43:32

#超算数
前屋毅「子どもたちを『数学ぎらい』にしてしまった学校での教え方は、日本の教育を象徴している」Yahoo!ニュース、2021年1月10日。news.yahoo.co.jp/byline/maeyats...
【芳沢教授は最近、『AI時代に生きる数学力の鍛え方－思考力を高める学びとは－』（東洋経済新報社）を上述している。】上梓だろう。

タグ：超算数

posted at 17:43:30

just *wrote of course, was trying to squeeze into the tweet and accidentally a word. Also I'm tired.

タグ：

posted at 17:12:23

(PR: github.com/rr-debugger/rr...)

タグ：

posted at 17:09:59

Oh also, microarchitectural behavior of Intel's performance counters, but that's pretty much par for the course when talking about rr, so I forgot ;).

タグ：

posted at 17:09:06

This rr test I just to repro a #julialang CI issue may not be the most cursed I've ever written for rr but it's definitely up there: github.com/rr-debugger/rr.... Has a little bit of everything - buggy legacy Linux APIs, Spectre-like BTB poisoning, very specific rr internals...

タグ： julialang

posted at 17:07:44

ToruNiina @ToruNiina

@ponkamepon TOMLは設定ファイル用の言語なので、触れる機会が多いのはパッケージマネージャとかですね。有名所だと、RustのCargoやJuliaのPkg、pipやpoetryも使ってます。他だと、静的サイトジェネレータのHugoとかも使っています。

タグ：

posted at 16:45:29

非公開

タグ：

posted at xx:xx:xx

Oちゃん @osuke0208

息子のテストなんだけど、下の式が違うらしい

上の式に赤で修正させられてんだよなぁ。

意図はわかるけど、やるなら筆算も指導しろよな。

掛け算なら入れ替えると意味は変わるけど、足し算はどっちでもよくない？

#算数　 #小学校　 #テスト pic.twitter.com/PICZxLtB7R

タグ：テスト小学校算数

posted at 14:44:34

かとたく @goldfish0900101

ごまふあざらし(GomahuAzaras @MathSorcerer

プランクトン＝海中の小さな生き物ではありません。プランクトン＝海流に逆らって泳げない生き物です。なので遊泳力の乏しいサルパやクラゲはプランクトンという括りになります。

タグ：

posted at 14:18:43

1.5 は時代遅れ，時代は 1.6
#Julia言語

タグ： Julia言語

posted at 13:54:11

かとたく @goldfish0900101

オオサルパは尾索動物亜門に属するホヤの仲間で植物プランクトンを餌とする動物プランクトンです。植物プランクトンは二酸化炭素を体内に取り込むためそれを食べるサルパは糞や死骸となることで炭素を海底へと運搬し、生物ポンプの役割を担っています。サルパは地球温暖化の救世主とも言われています。

タグ：

posted at 13:50:28

❄雪下❄ @Mopepe51

@physics303 Juliaはあまり詳しくないんですが、
sort!(list, by=x->(x=="foo" ? "" : x))
という感じでソート時の比較関数に"foo"だと空文字列で比較するような関数を渡してやれば行けそうだと思います

タグ：

posted at 12:50:24

@sandcat_neco 愚直にそれでやってみました～ pic.twitter.com/Qwvuj4Evz5

タグ：

posted at 12:48:06

非公開

タグ：

posted at xx:xx:xx

こんな感じ？ pic.twitter.com/0iTwdL1piJ

タグ：

posted at 12:41:45

juliaで重複のない

list = ["x","y","z","foo","w","a"]

みたいなリストがあったときに，特定の要素"foo"を先頭に持ってくるソート

sorted_list = ["foo", "x","y","z","w","a"]

ってどうやって得ればいい？

（listの中身は文字列）

タグ：

posted at 12:12:20

Juliaの実装で相談に乗ってほしいことがあるんだけど，研究のコアな部分なので迂闊にツイートできなくて悩ましい...．

タグ：

posted at 11:55:28

非公開

タグ：

posted at xx:xx:xx

今更ですが、今日は寒い。何この積雪！

タグ：

posted at 09:42:33

倉 - @kura

オブジェクト指向でない言語をいじると、戸惑う。OOPとは何なのか？などと考えてしまう。これまで触れてきたいくつかのプログラミング言語を通じて学んできたのは、言語ではなくOOPなのだと知る。そして理解が中途半端だということも。という訳で、Julia苦戦中。

タグ：

posted at 09:38:16

長尾年恭・東海大教授の地震予測「地下天気図」は、もう可愛そうなくらい当たらない。
blog.goo.ne.jp/geophysics_lab...

タグ：

posted at 09:34:32

#Julia言語実1変数函数の数値積分は QuadGK.jl パッケージが定番だと思う。

自分で下手に数値積分のコードを書くより、これを利用した方が計算が速いことが多い。

許容誤差を大きめにして計算時間を減らしたければrtol=1e-3とかに設定して使えばよい。

github.com/JuliaMath/Quad... pic.twitter.com/KrfAIPGZJQ

タグ： Julia言語

posted at 09:34:24

高橋学・立命館大特任教授は、過去に雑誌でデタラメ地震予測を連発。

「2016年の2～3月頃に東京湾直下大地震で津波」（週刊女性2016年）
「2017年のGW明けに福岡で大地震」（週刊女性、日刊ゲンダイ2017年5月）
「2019年のGWに南海トラフ巨大地震」（日刊ゲンダイ2019年4月）
blog.goo.ne.jp/geophysics_lab...

タグ：

posted at 09:31:53

#Julia言語分布Gamma(α, θ)におけるlog xの期待値は

1/(θ^α Γ(α)) ∫_0^∞ e^{-x/θ} x^{α-1} log x dx
= (∂/∂α)log(θ^α Γ(α))
= log θ + ψ(α)

と書ける。ここでψ(s)=(log Γ(s))'はディガンマ函数と呼ばれる基本特殊函数で、適当なライブラリを使えばコンピュータで気楽に計算できる。 pic.twitter.com/XJnLieIlSw

タグ： Julia言語

posted at 09:31:11

「昨年（2020年）の震度5弱以上の地震は計9回とほぼ例年通り。非常に静穏な年でした」（村井俊治・東大名誉教授）

…この村井さんはちょうど1年前の週刊ポストで「2020年は東日本大震災と同規模の地震が起こる可能性がある」と言ってたんですよ。酷くないですか？
www.news-postseven.com/archives/20200...

タグ：

posted at 09:27:53

長年懸念されている「南海トラフ」地震は起こるのか？　2021年の巨大地震を専門家が予測（週刊SPA！）
hbol.jp/236110

…専門家としてこの記事が紹介しているのが、村井俊治氏、高橋学氏、長尾年恭氏。もう笑うしかありませんね。

タグ：

posted at 09:24:28

#Julia言語微訂正

❌E[f(X)|X～p]
⭕️E[f(X)|X～dist]

モンテカルロ法による近似で十分なら

E[f(X)|X～dist]
↓
X = rand(dist, L)
mean(f, X)

と翻訳できる。対応が見易い。

確率変数Xはモンテカルロ法では乱数列Xに対応。逆にこのイメージで確率変数を(測度論を経由せずに)理解してもよい。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 09:22:00

#Julia言語サンプルX上でのf(x)の平均が

mean(f, X)

とシンプルに書けるのはありがたい。

このスレッドのテーマはこれの一般化がJuliaで広範に使用可能になっているという話。これを知っていると、数学語の直訳でJuliaのコードを書ける場合が増える。しかも多くの場合に計算速度的にも効率的。

タグ： Julia言語

posted at 09:16:29

Juliaの標準ライブラリのDictよりメモリ消費が少なく、検索が（概ね）速い実装ができたかもしれん。 pic.twitter.com/hPblEGOcaJ

タグ：

posted at 09:15:38

#Julia言語確率分布 dist における確率変数(函数) f(X) の期待値 E[f(X)|X～p] のモンテカルロ法による計算は

X = rand(dist, L)
mean(f, X)

例えば、ガンマ分布でのlog Xの期待値は

α, θ = 10, 0.1
using Distributions
dist = Gamma(α, θ)
L = 10^7
X = rand(dist, L)
ElogX_mc = mean(log, X) pic.twitter.com/VKXvBWx6Ji

タグ： Julia言語

posted at 09:11:32

非公開

タグ：

posted at xx:xx:xx

#Julia言語

open(f -> ファイル f を読み込んで行う操作, ファイル名)

と等価な

open(ファイル名) do f
ファイル f を読み込んで行う操作
end

の形式のコードはよく書く。

1つ前のツイートの内容はこれの応用例。

タグ： Julia言語

posted at 08:10:46

Dr. KID @Dr_KID_

発熱はないものの、泣き止まないため受診した乳児を分析した研究。

5%ほどで重篤な疾患が隠れていたようです。

また、多くの場合は問診や身体所見から発見につながっており、小児科医の丁寧な診察が必要と考えられます。

pubmed.ncbi.nlm.nih.gov/19255012/

タグ：

posted at 08:09:00

やねうら王 @yaneuraou

やねうら王の新しい詰将棋エンジンで将棋図巧１番が解けたようなので記事にしました。

あと、ミクロコスモス(1525手詰)を詰ます脊尾詰が不詰を証明できない局面について、不詰を証明できたらしいのでその件についても書きました。

yaneuraou.yaneu.com/2021/01/12/yan...

タグ：

posted at 08:08:25

#Julia言語

using Base64
showimg(mime, fn; tag="img") = open(fn) do f
base64 = base64encode(f)
display("text/html", """<$tag src="data:$mime;base64,$base64" />""")
end

showimg("image/jpeg", "foo.jpeg"; tag="img width=80%")

でJupyterノートブック上に画像を取り込める。

タグ： Julia言語

posted at 08:08:16

#Julia言語

サンプル A の経験累積確率分布函数

ecdf(A, x) = count(≤(x), A)/length(A)

using Distributions してあれば

ecdf(A, x) = mean(≤(x), A)

でもよい。ブロードキャストで使う場合には配列AをRefで保護する。

x = range(0, 6; length=600)
y = ecdf.(Ref(A), x) pic.twitter.com/JpcoVZRsyg

タグ： Julia言語

posted at 08:04:15

#Julia言語 #数楽 #統計

〇〇達について互いに独立に△△を実行するプログラムのコードを書くことは易しい場合が多く、並列化も易しい。

しかし、△△の実行が各〇〇ごとに独立になっていない場合は面倒になり、効率的な並列化も難しくなる。

確率論でもi.i.d.だと話が易しくなる。

タグ： Julia言語数楽統計

posted at 07:46:31

#Julia言語 foo(f::函数, X::配列など繰り返し可能な何か) のXや for ループのインデックスの動く範囲に使えるものが非常に多彩であることを知るには

methods(iterate)

の結果を眺めるとよいです。

〇〇について△△を繰り返すことをやらせたい場合には、このスレッドの内容が役に立つ。 pic.twitter.com/wgQAoDqIbH

タグ： Julia言語

posted at 07:43:14

#Julia言語 m×n行列 A のすべての列(縦⃗ベクトル)に函数 f を作用させた結果を並べた1×n行列は

mapslices(f, A; dims=1)

これをただの1次元配列にしたければ

vec(mapslices(f, A; dims=1))

Aが縦ベクトルの配列なら f.(A) で似たようなことをできますが、同じことを行列でやりたいことは多い。 pic.twitter.com/vOi6fUYtF9

タグ： Julia言語

posted at 07:33:21

#Julia言語 open(args...; kwargs...) の結果に函数 io -> f(io) を作用させることも

open(f, args...; kwargs...)

の形式で可能です。io -> f(io) の定義が複雑で一回限りしか使わない場合には

open(args...; kwargs...) do io
f(io)の定義
end

がよく使われている。

タグ： Julia言語

posted at 07:19:59

#Julia言語

Xのすべての要素に手続き f を施すには

foreach(f, X)

返り値はforループと同じnothingになります。

Xの要素に函数 f を作用させた結果を集めたものは f.(X) だけではなく

map(f, X)

で作れる。

これら以外にもmaximum, minimum, count, sum, mean, …も似た使用法が可能。

タグ： Julia言語

posted at 07:15:12

#Julia言語二乗和は

X = randn(10^6)
sum(Base.Fix2(^, 2), X)
sum(abs2, X)

のようにも書ける。Fix2(f, a)は本質的に x->f(x, a) です。読み易さは x->f(x,a) の方が上のことが多い。

1.96より大きい要素の割合は

count(>(1.96), X)/length(X)

using Statistics
mean(>(1.96), X)

タグ： Julia言語

posted at 07:08:22

#Julia言語

Juliaでは

foo(f::函数, X::配列など)

の形式で、配列Xなど(generatorやiteratorを含む)のすべての要素に函数fを施した結果に "foo" の操作を施せる場合が多数あります。

例えばXの要素の絶対値の最大値と和はそれぞれ

maximum(abs, X)
sum(abs, X)

二乗和は

sum(x->x^2, X) twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 07:01:45

添付画像中の太文字部分をツイッターで検索してみたら、放流分は氷山の一角でしかなかった。😱

検索で見つけた発言は「報告」しておきました。😊 twitter.com/totb1984/statu...

タグ：

posted at 06:19:58

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

ihnorton @ihnorton

Here's a #JuliaLang version, which could even fit on one long line:

jl> signpost(x) = ccall(:kdebug_signpost,
Nothing, (UInt32, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}),
x, C_NULL, C_NULL, C_NULL, C_NULL)

jl> signpost(1000) pic.twitter.com/kUeTevUFJt

タグ： JuliaLang

posted at 02:22:10

oka ఒక @nowohyeah

@hinahypersonica ROOT について詳しいポッドキャストのエピソードと言えばご存知かもしれませんがこちらです。TCFM らしく低レイヤーの話をしているので、 ROOT を5年以上使っていても聞いたことのない話ばかりでした
turingcomplete.fm/15

タグ：

posted at 01:11:31

oka ఒక @nowohyeah