黒木玄 Gen Kuroki(@genkuroki)/2018年03月19日

#統計以上では触れませんでしたが、その他にも

* 真の分布をカンニングしなければ計算できない量と真の分布をカンニングせずにサンプルのみから計算できる量の違いを強調しているか

という規準も結構重要。現実の統計学は真の分布(＝母集団分布)が不明の場合を扱わなければいけません。

タグ：統計

posted at 23:55:40

メガネ⭐︎動作確認中 @yasu_ando

#統計まとめ：「スタイン推定量」「階層モデル」「WAIC」がまとめて出て来る文脈では、解説者が「サンプルサイズ1」の場合と「サンプルサイズn」の場合を明瞭に区別し、looパッケージで計算した階層モデルに関するサンプルサイズ1のWAICの値を参照していないかどうかに注意を払う必要がある。

タグ：統計

posted at 23:52:04

Jupyter上でRも動かせるらしく、次のチュートリアルに向けて移行しなければと誓った。Juliaも気になってる。

タグ：

posted at 23:49:51

#統計例：n次元の多変量正規分布に従う確率変数 (Y_1,…,Y_n) はたった1回の独立試行で生成された乱数ベクトルとみなされるので、n次元の多変量正規分布のサイズ1のサンプルとみなされる。

例：Y_1,…,Y_nのそれぞれが同一の正規分布に従う独立な確率変数であるとき、そのサンプルサイズはnである。

タグ：統計

posted at 23:48:15

#統計サンプルサイズ1の設定とサンプルサイズnの設定の違いを認識することは階層モデルについて理解を深めるときにかなり大事です。

サンプルサイズの定義：同じ確率分布による独立試行で生成されたとみなされるY_1,…,Y_nのサンプルサイズはnである。

タグ：統計

posted at 23:45:42

#統計以上の確率モデル達について最尤法を実行するには積分で書かれた確率モデル

p(y|σ,w) = ∫ p(y|μ,σ)p(μ|w)dμ

を直接扱う必要があります。これはp(y|σ,w)の形がわかっていない一般の場合には結構大変です。最尤法よりもベイズ推定法の方が階層モデルを扱い易い理由もこれでわかりました！

タグ：統計

posted at 23:43:02

2次元の波動方程式まとめ。

色々試したので１つにまとめました。
twitter.com/i/moments/9757...

タグ：

posted at 23:41:45

ceptree @ceptree

MATLAB、ユニコード使えないの地味に痛い pic.twitter.com/WFZeeTIWgI

タグ：

posted at 23:41:31

あ〜る菊池誠(反緊縮)公式 @kikumaco

#統計確率モデル

p(y_1,…,y_n|μ_1,…,μ_n,σ,w) = Π_{i=1}^n p(y_i|μ_i,σ)p(μ_i|w)

の事後分布から確率モデル

p(y|σ,w) = ∫ p(y|μ,σ)p(μ|w)dμ

の事後分布が得られるおかげで、前者のモデルに関するMCMCによって後者のモデルも解けるわけです。後者は積分が入っているので直接の扱いは難しい。

タグ：統計

posted at 23:40:37

あら、僕も何か狙われてるのか。知らなかった。怖い怖い。そりゃあ通報しないといかんわねえ twitter.com/don_jardine/st...

タグ：

posted at 23:38:38

#統計こんなことになっている理由は、確率モデル

p(y_1,…,y_n|μ_1,…,μ_n,σ,w) = Π_{i=1}^n p(y_i|μ_i,σ)p(μ_i|w), σは既知

の(μ_1,…,μ_n,w)の事後分布のwへの射影は確率モデル

p(y|σ,w) = ∫ p(y|μ,σ)p(μ|w)dμ, σは既知

のwの事後分布に一致するという事情になっているからです。

タグ：統計

posted at 23:37:32

みんな波動方程式好きすぎでは pic.twitter.com/ctptrjZOOZ

タグ：

posted at 23:34:12

#統計所謂「スタイン推定量」の話題の面白さに目を奪われて、サンプルサイズ1ではなく、サンプルサイズnの設定で正しくWAICを計算しなければいけないことを見逃すと、looパッケージなどで機械的に階層モデルのWAICを計算して失敗してしまうかもしれません。この辺は要注意。

タグ：統計

posted at 23:30:56

#統計以上の解説では最初に正規分布を仮定しました。それは読者にイメージを持ってもらいやすいようにそうしただけではなく、所謂「スタイン推定量」の設定でもあるからです。「スタイン推定量」の設定はサンプルサイズ1の設定です。

階層モデルではサンプルサイズnの設定に変わる。

タグ：統計

posted at 23:29:03

非公開

タグ：

posted at xx:xx:xx

#統計階層モデルのWAICの計算では積分で表わされた p(y|σ,w) が必要です。それの公式が知られていなければ数値積分で計算する必要があります。結構面倒だし、コンピューターでの計算時間も跳ね上がります。

タグ：統計

posted at 23:23:20

the_spoiler @don_jardine

今菊池さんや林さんのネットリンチを呼び掛けている人(たち)が一時期「俺たちほどまともにデマと戦ってきた人間はいない」という趣旨のことを言っていて危うく顎が外れそうになった。

タグ：

posted at 23:21:51

#統計 Rのlooパッケージは階層モデルについては、このサンプルサイズ1の推定とみなした場合のWAICなどを計算する仕様になっています。だから、階層モデルでWAICを計算したと言われても、looパッケージのWAICなら「それ意味ないです」と教えてあげる必要があるわけです。

タグ：統計

posted at 23:20:30

#統計注意が必要なのは、階層モデルの設定もやろうと思えばサンプルサイズ1の設定ともみなせることです。すなわち、確率モデルは

p(y_1,…,y_n|μ_1,…,μ_n,σ,w) = Π_{i=1}^n p(y_i|μ_i,σ)p(μ_i|w)

であるとみなすと、Y_1,…,Y_nを与えることはサイズ1のサンプルを与えることになってしまいます。

タグ：統計

posted at 23:18:26

#統計 Normal(μ,σ)の密度函数をp(y|μ,σ)と書くと、階層モデルの設定の確率モデルは

p(y|σ,w) = ∫ p(y|μ,σ)p(μ|w)dμ

と積分で表わされることになる。今回は私にしては珍しく、異なるpを全部同じpで書いているので読者は注意して下さい。独立変数名で区別しています。(やはり、やめた方がいいかな？)

タグ：統計

posted at 23:14:57

@koumori_2011

財務省の増税工作は陰謀じゃなくて業務で粛々とやってるのは常識だよね。(*'ω'*)

タグ：

posted at 23:14:50

ぎーち(ブレイク兄) @BREAK_BROTHER

#統計階層モデルの設定では、Y_i達は同一の分布に従う独立な確率変数とみなされるので、サンプルサイズはnとみなされます。サンプルサイズが1とnでは大違いなので、この辺の違いには敏感になるべきだと思います。

タグ：統計

posted at 23:12:02

pic.twitter.com/fDb6nax5RY

タグ：

posted at 23:10:49

ぎーち(ブレイク兄) @BREAK_BROTHER

衛星設計をする時によく見るJAXAの技術資料たち。

全部で5049ページあります。

これ何が凄いかって、共通文書で一般の方もダウンロード出来るんですよ。

衛星開発に興味のある方は是非読んでみてくださいね！(地獄への誘い) pic.twitter.com/GhYATVTQD4

タグ：

posted at 23:10:43

#統計個体差を考慮する階層モデルでは、μ_i達に次の制限を付けるところが、一つ前のツイートの設定と違っています。

階層モデルの設定：μ_i達は同一の分布 p(μ|w) (wはパラメーター) に従う独立な確率変数である。

この設定では、Y_i達は同一の分布に従う確率変数とみなされるようになります！

タグ：統計

posted at 23:09:58

#統計例えば、Y_1,…,Y_n達は独立で各Y_iは正規分布Normal(μ_i, σ)に従っているとします。未知のμ_1,…,μ_nを推定するという設定はY_i達が異なる分布に従うので「サンプルサイズ1」の推定になることに注意が必要です。数ベクトル(Y_1,…,Y_n)がサイズ1のサンプルになる。続く

タグ：統計

posted at 23:04:31

#統計詳しくは以下のリンク先で解説しておきました。

twitter.com/genkuroki/stat...

タグ：統計

posted at 23:01:24

#仲正昌樹氏の連載を読む
（続き）　むしろ、「俺はユークリッド幾何学と不完全性定理の関係についても述べている」と反論するならわかる。

あるいは私が「仲正のバカがユークリッド幾何学と不完全性定理の関係を誤解している」と言ったら、もっともな反論なんだけどなあ
＃ポストモダン

タグ：ポストモダン仲正昌樹氏の連載を読む

posted at 23:01:15

#統計「階層モデル」「WAIC」という用語を使う人を見たら、必ず確認しておくべきこと：

「WAICはどうやって計算しましたか？もしかしてRStanを使っていて、looパッケージで計算したりしていませんか？」

回答がYesなら「そのWAICはサンプルサイズ1と解釈した場合の意味のない数値ですよ」となる。

タグ：統計

posted at 22:59:40

天むす名古屋 Temmus @temmusu_n

＃仲正昌樹氏の連載を読む
仲正氏はユークリッド幾何学と不完全性定理を結び付ける話などをしたことがない、と私のツイートに反論している。
でも私は「仲正氏がユークリッド幾何学と不完全性定理の関係を無視している」と言っているんだから当然だと思うんだけどなあ。
＃ポストモダン pic.twitter.com/eCmE51quUY

タグ：ポストモダン仲正昌樹氏の連載を読む

posted at 22:56:36

亠(なべぶた)の第一画が第2画に直角でないと漢字テストで誤答とされる実例があるとか。活字を手書き文字の規範とみなす傾向の一部なのだろうが、こんな細かいところまで統制を強めるのは、看板屋さんとか一定の書体を実現しなければいけない人たちだけがやればいいこと。隷書でも左傾が多いよなぁ。

タグ：

posted at 22:48:16

非公開

タグ：

posted at xx:xx:xx

#Julia言語 #Python #Rstats

JUlia + PYThon + R = JUPYTeR notebook

私は Ruby も Jupyter notebook で使えるようにしてあります。

タグ： Julia言語 Python Rstats

posted at 20:53:08

非战闘員 @dontsentouin

これでPython,Julia,Rのトロイカ体制に…

タグ：

posted at 20:47:30

ゆーま @reisin0930

juliaを使ってみたいので会社のpcに入れていたら時がたっていた

タグ：

posted at 20:43:45

Yuki KUBOTA @aza_yuuki1021

山口貴士 aka無駄に感じが悪いヤマベン @otakulawyer

旧帝大がこぞって「人間力」という言葉で人間を評価しようとしてるの、世紀末みがある。 pic.twitter.com/pHS14saSKy

タグ：

posted at 20:38:25

山口貴士 aka無駄に感じが悪いヤマベン @otakulawyer

記者会見の全体がこちらで公開されています。私が通訳もしています。
www.youtube.com/watch?v=AJO4di...

タグ：

posted at 20:34:09

日本脱カルト協会は、３月１５日、法務大臣に対し「松本智津夫（麻原彰晃）死刑囚以外の12人の死刑を執行せず、いずれは無期懲役に減刑すべき」という要請を行い、本日、外国特派員協会で記者会見をしました。日本脱カルト協会事務局長としてのお仕事です。

thepage.jp/detail/2018031...

タグ：

posted at 20:21:25

Hiroyasu Kamo @kamo_hiroyasu

藤井光／Hikaru Fujii @fujiihikaru

十年前の自分が証明を知っているはずですが、連絡がとれなかったので現在の自分が証明しました。4時間かかりました。

タグ：

posted at 20:13:16

英語のアルファベットに生き残れなかった12文字。thやwに取って代わられたり、発音が変わって不要になったりと事情は様々ですが、面白いのは”&”が19世紀まで27番目のアルファベットに入ることがあったという話。この12文字が今でも英語の本に出てきたら楽しそうなんですが。
mentalfloss.com/article/31904/...

タグ：

posted at 19:56:22

Hakky(´･∀･) @St_Hakky

MINO@「作例で学ぶ Substanc @MINO_8601

「ベイズ推論による機械学習」を読んだ | 10001 ideas wp.me/p95uaT-5a via

タグ：

posted at 19:41:46

ツイートで知って早速入れたChromeの拡張機能「Read Ruby」めっちゃ便利だわ…。難しめの英単語にルビ振ってくれるやつ。平易な英文だと翻訳しなくてもそれなりには読めるけどたまにある難しめな単語をいちいち調べるの面倒だったからすごく助かる。インストールはChromeのウエブストアから。 pic.twitter.com/32EQuXezhz

タグ：

posted at 19:25:53

ニカイドウレンジ @R_Nikaido

ゲーム紹介したらストアとか公式へのリンクをスレッドで繋げてやらんとな。 www.sokobond.com

まだ序盤しかやってないけど、終盤鬼のように難しくなるであろうことが容易に想像できる。

タグ：

posted at 18:51:05

ニカイドウレンジ @R_Nikaido

倉庫番系パズルゲーム『Sokobond』、めっちゃ秀逸。
・「原子を動かして他の原子と結合させて実在の物質を作る」という目的が秀逸
・原子ごとに手の数が決まっているのがパズルに見事活かせているのが秀逸
・結合していくごとにプレイヤーの形が変わるのでそれで色々ジレンマが生まれてるのが秀逸 pic.twitter.com/rMwn0pYywE

タグ：

posted at 18:42:09

道草 @econdays

タイラー・コーエン「Twitterではフェイクニュースの方が拡散しやすい？」 econ101.jp/%e3%82%bf%e3%8... (経済学101)

タグ：

posted at 18:36:03

じゃん @jean_lose

Juliaって言語初めて聞いた

タグ：

posted at 18:18:00

#Julia言語そうそう、parse(Int, "-123") のようなコードを見ると、「え？パースしなきゃいけないの？遅いんじゃない？」と思う人がいるかもしれませんが、速いです。

tk3369.wordpress.com/2018/02/04/an-...

のParse Intの節を見て下さい。

タグ： Julia言語

posted at 17:35:23

NOBIE.jl @NOBIE41

@genkuroki C言語の最高の発明はprintfですよね。

タグ：

posted at 17:32:50

#Julia言語浮動小数点数の文字列変換で @ sprintf を使うより効率的な方法を知っている人がいれば教えて下さい。まだそういう基本的なことさえ理解していない。

タグ： Julia言語

posted at 17:32:50

#Julia言語もとツイートが消えちゃいましたが、整数の文字列への変換法をゆるぽしていました。

タグ： Julia言語

posted at 17:31:09

twitter.com/t_zepp/status/...

#Julia言語整数の文字列への変換なら dec(x) で逆向きは parse(Int, s)

s = dec(-123) # s = "-123"
a = parse(Int, s) # a = -123

docs.julialang.org/en/stable/stdl...

浮動小数点数については @ sprintf "%f" -123.456 よりも効率的な方法があるかどうか知らない。

タグ： Julia言語

posted at 17:30:14

パステルしこめ @passhiko

役所は教えてない、申請すればもらえるお金・戻ってくるお金 pic.twitter.com/d61qIykZNy

タグ：

posted at 17:13:02

パステルしこめ @passhiko

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

役所は教えてくれない、申請すればもらえるお金・戻ってくるお金 pic.twitter.com/6UQF5Azow3

タグ：

posted at 17:03:23

加藤一二三【公式】 @hifumikato

@zakozako291 処理はpandas側で行われるので，特に速くならないと思います。むしろJulia<=>Python間のデータ変換コストの分遅くなるかと。

タグ：

posted at 16:37:11

うひゃ！うひゃ！うひゃ！
めちゃめちゃかわいい子猫ちゃんたちだよー(o^^o)
寝るとき思い出しそう(≧∀≦) pic.twitter.com/scT3JLGHMx

タグ：

posted at 16:21:52

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

JuliaHub @JuliaHub_Inc

Engineers at @fugroroames use #MachineLearning in #JuliaLang to identify network failures and potential failures 100x faster, by analyzing high resolution, detailed 3D Maps of the network (see below)

More info : goo.gl/EQSHno pic.twitter.com/xwVy7aNiJQ

タグ： JuliaLang MachineLearning

posted at 15:37:44

仲正昌樹氏のネット連載は、エゴサーチをして批判を発見して条件反射的にメモした文章の集合という印象。
それならツイッターを始めればいいのに。
ツイートを集めて後でまとめればいいから簡単。仲正氏が高く評価する千葉雅也氏もツイートして纏めて本にしている

#ポストモダン　 #現代思想

タグ：ポストモダン現代思想

posted at 14:41:05

“2013年4月30日アベノミクスでデフレ脱却はなるか

▽2%の物価目標は世界標準

そもそもインフレ目標はインフレを防止するために考案された仕組みであり、ハイパーインフレや財政ファイナンスを防止するのに最適なのである”

👍‼️

www.iist.or.jp/jp-m/2013/0218...

タグ：

posted at 14:35:13

中村圭志　神話・無神論・宗教学 @seattlelubbock

“受動喫煙対策
自民が厚労省に歩み寄れ” 👏
共同通信・柳沼勇弥

www.saga-s.co.jp/articles/-/153...

タグ：

posted at 14:19:52

ホーキング博士が亡くなったとき、原理主義者が「無神論者は今ごろ業火に焼かれている」とツイートして、ドーキンズ教授に批判されています。ドーキンズは、なぜ宗教の形をとったときだけヘイトスピーチが許されるのか疑問に思っています。一般的にこれは今後どんどん大きな問題になるでしょう。 pic.twitter.com/qSz7aymA1Y

タグ：

posted at 14:17:15

うおッ！「2012年9月」の時点で消費増税に警鐘を‼️
もう熱烈なファンになるしかない。😍

www.iist.or.jp/jp-m/2012/0211...

タグ：

posted at 14:16:14

Grahamian | データアナリスト @grahamian2317

“2017.05.17 【論説】　東日本大震災６年・被災地差別／正しい知識の普及を
柳沼勇弥

この論説は共同通信社によって配信され、山形新聞、茨城新聞、日本海新聞、山陰中央新報、大分合同新聞（順不同）などの地方紙に掲載されました。”

不覚にもノーマークであった。

www.tokyo-kuwano.com/blog/blog027.h...

タグ：

posted at 14:11:42

一般化線形モデルとかベイズみたいな基本やらないでscikit-learn, DeepLearnigうぇーいな今のAI業界やばいと思う。

タグ：

posted at 13:42:21

TJO @TJO_datasci

金子洋一神奈川20区（相模原市南区、座間 @Y_Kaneko

マーケティングとか実社会・実世界系のモデリングやると一定の割合で必ず混合効果モデルとか階層ベイズモデルとかやらされる羽目になるというのが経験値から来る肌感覚なんだけど、機械学習ベースのデータサイエンス講座の類でその辺取り上げてるのってかなり稀なので、そこはいつも気になってる

タグ：

posted at 13:38:06

消費増税分のわずか1/5しか社会保障の増額にあてられていれないのだから、慶応大の土居氏の「財務省への不信から税財政改革を遅らせて、困るのは財務省ではなく国民である」との主張は奇奇怪怪そのものである。困るのは見返りの法人減税がなくなる財界くらいのものだろう。　toyokeizai.net/articles/-/212...

タグ：

posted at 13:27:45

（続き）　仲正氏は、尾上氏について「学者としてまともな仕事もせず」と判断しているけど、これも根拠が不明。
知り合いなのかな？

＃仲正昌樹　　＃ポストモダン　＃現代思想

タグ：ポストモダン仲正昌樹現代思想

posted at 13:24:06

今回の仲正昌樹氏の連載では、私以外に奈良大の社会学者の尾上正人氏のツイートを批判している。これは仲正氏についてのツイートではなく、千葉雅也氏についてのツイートみたいだ。仲正氏は「自称社会学者」と揶揄しているけど、どういう基準で「自称」と判断しているんだろう？ pic.twitter.com/9xDVOO9snO

タグ：

posted at 13:18:03

tₖᶠ @tkf

Julia は @ edit や @ less で REPL からソースまでジャンプ出来てとても便利なんだけど、引数の型を指定しないといけないのが面倒だったので peco を使って絞り込み検索できるようにしてみた
github.com/tkf/Interactiv... pic.twitter.com/ANZgUPn0NG

タグ：

posted at 13:08:02

「全裸のシュワルツネッガーをからかうチンピラ」って、仲正昌樹氏にちょっかい出す今の俺かも。
でもツイート検索しても、今は俺以外いなんだよな。そんなチンピラ

twitter.com/shinkai35/stat...

＃現代思想　＃ポストモダン

タグ：ポストモダン現代思想

posted at 13:05:41

コーヒーパパさん @0418291

金子洋一神奈川20区（相模原市南区、座間 @Y_Kaneko

スペイン人の知り合いに「スペイン人って、毎日パエリア食べて闘牛見てるの？」と聞いたら「そんな事ないよ。日本人だって毎日米食べてアニメ見てないでしょ？」と返ってきたので「日本人は毎日米食ってアニメ見てるよ」と返した。「ファッキンクレイジー！」って言われた。

タグ：

posted at 12:57:30

勘違いしてはいけない。潜在成長率は成長率の上限ではなく、いわば「過去の経済から計算した標準的成長率」だ。少しプラスになったからといって、政府が政策を引き締めてはデフレ体質から抜け出せない。⇒２９年の需給ギャップ、９年ぶりプラス　デフレ脱却判断に環境整う www.sankei.com/politics/news/...

タグ：

posted at 12:50:56

Miura Hideki @miura1729

あ〜る菊池誠(反緊縮)公式 @kikumaco

TOTOって流しやトイレの排水の設計のために流体計算のためのフレームワークを並列計算機(TUBAMEだったかな？)に実装して論文を出したりして何言にハイテク企業と言うイメージ。LIXILの営業にTOTOみたいに学術論文出さないんですか？って聞いたら嫌な顔された

タグ：

posted at 12:29:55

金融緩和の効果でデフレ脱却がようやく見えてきたということなんだけど、消費税を上げずに積極財政をやってりゃとっくに脱却できていた可能性が高いわけよね。これは失政だよね。だから、アベノミクスの上をいく経済政策は可能なのに、なぜかアベノミクスを完全否定したがる政治家が多いよね

タグ：

posted at 10:51:10

波瑠 @shizupon2004

@nycity1022 @takamatsu_yabu @JapaneseSchool @genkuroki 王様。1年の知的の子に…何回も書き直しさせてた。 pic.twitter.com/zuh1jxj0hG

タグ：

posted at 10:33:49

久保拓弥 @KuboBook

"On Sociology: 計量分析を用いた論文も社会的な構築物なのか？：経済学と政治学の再現性ポリシーの違い" on-sociology.blogspot.com/2018/03/blog-p... pic.twitter.com/9slGMBYD34

タグ：

posted at 09:39:49

名取宏（なとろむ） @NATROM

過剰診断がなくても、単に診断が前倒しされただけでも、5年生存率は「改善」する。リードタイムバイアスという。がん検診の有効性を5年生存率で評価してはいけない。[がん検診の考え方　なぜ「生存率」ではだめなのか canscreen.ncc.go.jp/kangae/kangae3... ]。

タグ：

posted at 09:29:35

名取宏（なとろむ） @NATROM

年間10万人中20人が診断され、うち10人が5年以内に亡くなる疾患があったとする。5年生存率は50％（10÷20）。ここで有効性ゼロだが新たに20人の過剰診断が生じる検診を行うと、年間に10万人中40人が診断され、うち10人が5年間中に亡くなることになる。5年生存率は75％（30÷40）に見かけ上「改善」する。 pic.twitter.com/SzWnh73NV0

タグ：

posted at 09:29:13

円形の境界があるとして計算。複雑なものもいいけど、こういうシンプルなものもいいですね pic.twitter.com/pZpDu6VBZX

タグ：

posted at 09:14:06

名取宏（なとろむ） @NATROM

過剰診断が増えれば5年生存率は改善します。5年生存率の改善は過剰診断に対する反論にはなりません。 htn.to/BA9Lnr

タグ：

posted at 08:40:54

2次元の波動方程式。
四方は壁に覆われているものとして計算。1次元よりも複雑な動きが見れて面白い。 pic.twitter.com/8CAl5gFl3O

タグ：

posted at 08:37:25

片瀬久美子 @kumikokatase

天むす名古屋 Temmus @temmusu_n

再度言明しておきますが、私は辛さんが福島デマを広めた行為を問題視していますが、一方で石井氏らによる彼女への中傷は全く擁護できません。 twitter.com/kumikokatase/s...

タグ：

posted at 08:29:54

Dr. Chris Rackauckas @ChrisRackauckas

#超算数ページは485の間違いだった。

タグ：超算数

posted at 07:49:05

Dr. Chris Rackauckas @ChrisRackauckas

@ceptree @genkuroki @termoshtt Discontinuity tracking is standard in most state-dependent DDE solvers for this reason (see RADAR5). On top of that, none of the MATLAB DDE solvers can handle stiff equations. Giving it a rating of Fair is quite generous when the facts are laid out. Only R's deSolve is worse.

タグ：

posted at 05:16:26

Dr. Chris Rackauckas @ChrisRackauckas

@ceptree @genkuroki @termoshtt 3rd order stepping is too low to get good efficiency out of an RK method, and residual control strategies fully rely on adaptivity error control which throws dt down too low. The good thing is these are flexible and can apply anywhere, but they are not efficient (try it yourself)

タグ：

posted at 05:13:53

@ceptree @genkuroki @termoshtt My explanation in the post still holds. MATLAB's dde23 is low order. MATLAB's ddesd is only residual control RK4 without discontinuity tracking. Both of these have analogues in DifferentialEquations.jl, but aren't recommended to be used because they benchmark poorly.

タグ：

posted at 05:12:46

ゴルゴ・サーディーン @golgo_sardine

.@any_k_tg さんの「差別の被害者が別の差別に加担しうる、というきくまこさんの意見と、不思議な反差別界隈の皆さんの反応。」togetter.com/li/1209920 をお気に入りにしました。

タグ：

posted at 02:33:29

@youkokko いきなりで失礼します。
【掛け算と足し算に順序ないの？】
何をごらんになったのかハッキリしないですが、おそらく、
「掛算と足算に順序ない」という論は「Ａ×ＢでもＢ×Ａでもいいだろ」という主張です。
加減乗除が混在した式では、乗除が優先、で正しいです。#掛算

タグ：掛算

posted at 02:19:47

Taku @Less_role

ゴルゴ・サーディーン @golgo_sardine

julia一回コンパイルしたらすごく速い。fortranに負けてない

タグ：

posted at 02:09:23

@hama_t_titech そもそも【どちらを先に書いても良いというのが「×」という記号の使い方なのです】
という事を順序アリ論者が認めないから
「トランプのように配るだとか」
と言わなければならなかったのですが。 #掛算

タグ：掛算

posted at 01:10:49

twitter.com/ito_haru/statu...

#統計私もそう思います。だから深刻な問題だと思ってこの話題を繰り返しています。

著者が「間違ってました。ごめんなさい」と言うべきだと思うのですが、私には説得不可能だと感じています。

周囲に「騙されないように注意！」と警告するしかない感じ。

タグ：統計

posted at 01:02:33