はてなブックマーク

Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学 http://abrahamcow.hatenablog.com/entry/2014/09/11/024924 「時系列データの相関係数はあてにならない……のか？　教えて下さい」なんだけど... 私のコメントが気に触ることが多いようなのですが（特に悪意はないつもりなんですけど，すみませんね） 私は，経済学とか時系列についてはよく知らないのですが，「これは「見せかけの相関（擬似相関；spurious correlation）」の例だ」ということならば，偏相関係数を考えればよいのではないでしょうかね？？社会学などでは当たり前のように使われていると思うのですが。 > set.seed(1) > x = cumsum(rnorm(100)) > set.seed(2) > y = cumsu

Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学 筆者が恥ずかしかったのか，今では非公開扱いになってしまって，ほかの人が閲覧できなくなっている。 で， 「ダイナマイト・プロットとはこのページにある最初のグラフのようなもの。」 というのが意味不明（参照不能）になっているので，「ダイナマイトプロットって何？」に答えておこう。 まあ，以下のようなもの。昔の漫画（若い人には分かるまい）のダイナマイトの点火装置（上のバー）を押すとダイナマイトが爆発するというような図である。こういうと，今の若い人にも「だせ〜〜〜〜」とわかるだろう。 というのが捕捉で，以下が当時の記事。 ============================================================================= http://abra

Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学 「A/Bテストの数理」への批判 http://abrahamcow.hatenablog.com/entry/2013/12/28/210035 だけどねえ。ちょっといただけない。 ではこの場合、有意水準はいくつくらいにしたらいいんでしょうか？ とてもよい質問です。 たぶん20％くらいでいいんじゃないかと思います。 その根拠は天気予報です。 降水確率0％～20％くらいならだいたいみんな傘持ってかないじゃないですか。 40％～100％ならだいたいみんな傘持っていくんじゃないでしょうか。 30％だと人によってはちょっと迷うと思います。 だから20％で切っちゃいましょう。 そのくらいのざっくりした感じで十分だと思います。 有意水準は，リスク・ベネフィットの観点を考慮するというのはあたりまえな

Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学 相変わらず，この記事が参照されるようで。 以下のプログラムを見て分かる人はそもそも，こんな記事を読みに来るとも思えないので，一応，式を書いておきましょうかねぇ。 二点 (x1, y1) と (x2, y2) を通る直線と，もう二組みの点 (x3,y3) と (x4,y4) を通る直線の交点の座標 （x，y） は，まず， a1 = (y2-y1)/(x2-x1) a3 = (y4-y3)/(x4-x3) を計算して，その後， x = (a1*x1-y1-a3*x3+y3)/(a1-a3) y = (y2-y1)/(x2-x1)*(x-x1)+y1 ということです。下のプログラムそのまんまですね。。。下の方に追記 (x1,y1)と(x2,y2)を通る直線と，(x3,y3)と(x4,y4)を

算額（その1692） 新潟県三条市 三条本成寺 北越蒲原郡三十番神 寛政12年(1800) 神谷定令門人 東都本郷 松下清六郎興昌 藤田嘉言(1807)：続神壁算法 http://www.wasan.jp/jinpeki/jinpeki.html 涌田和芳，外川一仁：三条本成寺の紛失算額，長岡工業高等専門学校研究紀要 第49巻（2013）. https://kinpoku.nagaoka-ct.ac.jp/lib/kiyo/vol_46-50/vol_49/49_07wakuta2.pdf キーワード：菱形，等円 #Julia, #SymPy, #算額, #和算 菱形の中に 2 本の斜線を引き，等円を数個容れる（全部で偶数個）。菱形の対角線の長い方が 4 寸，短い方が 3 寸，等円の個数が全部で 10 個のとき，等円の直径はいかほどか。 解法は涌田，外川（『続神壁算法解義』）によるが，Sy