r-statistics-fanの日記

統計好きの現場の臨床医の覚書のようなもの

グラードンのレイド時のCPの表(PL20,25)

レイドにグラードンが出ているではないか。 表を作成しておく AT DF HP 固体値% CP: PL 20 HP: PL 20 CP: PL 25 HP: PL 25 最大 CP: PL40 最大CP 順位 レイド上位% 80% get必要数 同90% 15 15 15 100.0 2328 117 2910 131 4074 1 0.5 322 460 15 14 15 97.8 …

ゴローニャのレイド時のCPの表

ゴローニャのレイド時のCPの表がほしいというリクエストが有ったので、表を作成しました。 アフィ貼っているわけでもなく、ほとんどblogにコメントも書かれないので、こういった反応があると嬉しいですね。 種族値は 【ポケモンGO】ゴローニャの巣とレア度&…

ホウオウのレイド時のCPの表(PL20,25)

2017.12.9午後3時 ゲットしたホウオウのCPが下がってるー。種族値を再々変更 #2017.12.9天候によってPL20とPL25の場合が有る上に、種族値も変わったらしいので再計算しておく。今日やっと1体目のホウオウをゲットしたら種族値が変わっていて表が役に立たず…

Fisher's exact testが"正確”かどうか

https://oku.edu.mie-u.ac.jp/~okumura/stat/fisher-chisq.html https://twitter.com/genkuroki/status/910360137256325121このあたりのFisher's exact testが"正確”かどうかという議論が勉強になるそれぞれの立場でそれぞれ一理ある前提として周辺度数の固…

ライコウのレイド時のCPの表~ポケモンGO

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ヨーロッパ出張で希少な自由時間にエンテイをゲット。1/1ゲットなので運が良かった。 ついでに運良くバリヤードもゲットした。 ということで、ライコ…

エンテイのレイド時のCPの表~ポケモンGO

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ヨーロッパ出張で希少な自由時間にエンテイをゲット。1/1ゲットなので運が良かった。 ついでに運良くバリヤードもゲットした。 ということで、ライコ…

複数人(同窓会対応)の生存曲線の描画とN年後に全員生存している確率を計算するwebアプリ

shinyアプリへのリンク ##追記:無料サービスを利用しているので接続が月25時間を超えると翌月まで使えなくなってしまうようです。自腹で課金する気はないので、つながらない場合はごめんなさい。 前回の記事 をshinyでwebアプリ化(ここをクリックすると…

スイクンのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 https://pokemongo-get.com/pokego01148/ を参考に、種族値HP200/AT180/DF235を採用し、スイクンのレイド用の表を作っておく。 AT DF HP 固体値% CP:…

同窓会で5年後に全員生存している確率は~H27の日本人完全生命表で計算する

#追記 ここをクリックでブラウザで各自計算できます 昔の記事で日本人の生存曲線を扱ったが、H27の国勢調査のデータによる完全生命表が出ているのでアップデートすることにした。http://www.mhlw.go.jp/toukei/saikin/hw/life/22th/ ここのデータを引っ張っ…

カビゴンのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

ラプラスのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

カイリキーのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

バンギラスのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

ルギアのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

フリーザーのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

ファイヤーのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

サンダーのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

ミュウツーのレイド時のCPの表

2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最大強化時CPを更新しました。 ミュウツーのEXレイドに備えてレイド時のCPの表を作っておく。 ついでに、バンギラス、ファイヤー、フリーザー、サンダー、ルギア、カイリキー、ラプ…

統計推理。どう計算したのだろう。比率の信頼区間

どう計算したのだろう。比率の信頼区間日経メディカルで、コデインの小児への処方が制限されたと言う記事を見た。PDMA(医薬品医療機器総合機構)が運営するMID-NETというデータベースを利用した解析結果が参考にされたらしい。http://www.jpa-web.org/dcms_…

100人を部屋に集めてお金をランダムな相手に渡し続けるとだんだんと貧富の差が生まれる~条件により印象通りの結果に:第2弾

GIGAZINEでの興味深い記事の検証第2弾。 http://gigazine.net/news/20170711-random-people-give-money-to-random-other-people/ 100人を部屋に集めてお金をランダムな相手に渡し続ける」とだんだんと貧富の差が生まれる というもの。記事の中では「45ドルを…

GIZAZINEの記事をなぞってみる:「100人を部屋に集めてお金をランダムな相手に渡し続ける」とだんだんと貧富の差が生まれる

gigazine.net#7/12追記 普通にソートした順でboxplotをすれば分かりやすかったので最下部に図を追加GIGAZINEで興味深い記事が。 「100人を部屋に集めてお金をランダムな相手に渡し続ける」とだんだんと貧富の差が生まれる というもの。記事の中では「45ドル…

#OsakaStan5参加 献本争奪じゃんけんで2回で決着した事象は珍しいものだったのか?

はじめて統計の勉強会に参加した。#7/2一部バージョン違いのコードだったので修正いわゆるアヒル本の勉強会である。 www.amazon.co.jpなんと著者の松浦先生もいらっしゃると聞いて駆けつけた次第である。何故アヒル本と呼ばれるのか? それは、この表紙のブ…

レイド戦後捕獲時のバンギラスCP/固体値と最大強化時CP/レア度順位の表

#2017.11.28 アップデートで最大強化ポケモンレベルが39から40になったので、最新の表は r-statistics-fan.hatenablog.com こちらをご参照ください 統計ブログ()だが、わりとポケモン記事へのアクセスが多い。 というわけで追加記事。レイド戦が始まり、バン…

データの出力時にrownamesを出したくない/ポケモンGO,ラッキーCP777になる固体値全部

ふつうに行列をコンソールに出力すると> lucky777 AT DF HP Lv [1,] 0 2 14 24 [2,] 0 3 11 24 [3,] 0 3 12 24 [4,] 0 4 8 24 [5,] 0 4 9 24 のようにrownamesがNULLであっても、[1,] のように何行目か 分かるようにコンソールに出力される。 これを消したい…

超ムカつくtibbleクラス:バージョンアップでエラー頻発

Rのバージョンアップで、エラーが出ることが多くなった。その多くはどうもtibbleクラスが関わっている! 色んなパッケージがtibbleではエラーを吐きまくる。読み込んだ段階でdeta.frameがtibbleになっていることがしばしばである。たとえばこんな感じ > str(…

ポケモンGO バンギラスの厳選はどこで妥協するか。Rで計算してみる

ポケモンgoで岩祭り開催中。 ヨーギラスがよく出るので、厳選したいが、どの程度で妥協すればよいだろうか。 バンギラスの各攻撃、防御、HP(スタミナ)毎での、最大CPやレア度を知りたい。 Rで計算してみる。#順番をCP強い方が先頭に来るように改定。CP順位…

R3.4.0にアップデート。Rが高速に。

R3.4.0が出たと聞いた。 JITバイトコンパイラになり、loopが高速になるらしい。 適当にフィボナッチ数を計算する関数を作って試してみる。 fibonacci <- function(n){ if(n == 1) 1 else if(n == 2) 2 else (fibonacci(n-2) + fibonacci(n-1)) } fibonacci2 …

探偵ナイトスクープ 得点した選手の背番号の合計は偶数が多い?

ナイトスクープでサッカーの試合で得点した人の背番号の合計が 偶数有利か奇数有利かということを検証していた。偶数背番号の選手も奇数背番号の選手もおなじ得点率なら同じになるだろう。 しかし、実際は異なる。 偶数の選手が得点することがたまたま多かっ…

受験番号と合格率の関係

受験番号と合格率の関係 rs_fan_jp 2016年4月3日 2016年4月3日20:20念のため追記。言うまでもなく受験番号から合格への直接的因果関係はないからね。統計上有意ですが、単なる見かけ上の関係です。間違っても受験番号1ゲットなどという不毛なことを目指さな…

レアなタクシーに遭遇する確率

レアなタクシーに遭遇する確率京都で学会があり、レアなタクシーに遭遇する機会があった。ヤサカタクシーは通常三つ葉なのだが、四つ葉と双葉 が存在する。それぞれ1200台中4台と2台しか無いらしい。今回同じ日に四つ葉と双葉を見た。なんだか良い気分。多分…