統計

本記事では、ロジスティクス回帰分析について記載します。 【目次】 計算式等 ノート プログラムコード 雑談 参考 計算式等 ロジスティック回帰は二値分類（）で用いられる分析手法です。モデル式は回帰分析と似ていますが、シグモイド関数を用いて0～1の範囲の結果を出すことが特徴です（通常の回帰分析では0～1の範囲を超えてしまい、解釈不能な結果になってしまいます）。（青線：ロジスティック回帰、緑線：単回帰分析） モデル式、シグモイド関数および最適な回帰係数 β を取得するための損失関数はそれぞれ次のとおりです。 モデル式 シグモイド関数（活性化関数） 交差エントロピー誤差（分類問題の損失関数） ※ l…

箱ひげ図とは 箱ひげ図をExcelで作成する方法(Excel2016) 補足情報 箱ひげ図を作成する際は、四分位数を使います。四分位数については、以下の記事をご参照ください。 箱ひげ図とは 箱ひげ図とは、データのばらつきを可視化するための統計図です。データ群ごとに箱ひげ図を作成し、比較することで、データのばらつきを調べることもできます。 箱ひげ図をExcelで作成する方法(Excel2016) 以下のデータで箱ひげ図を作成してみましょう。 1.データ部分を選択した状態で「挿入」タブ＞「統計グラフの挿入」＞「箱ひげ図」を選択します。※A列の日付情報は選択しません。 以下のような箱ひげ図が作成され…

社会科学のための ベイズ統計モデリング (統計ライブラリー )作者:浜田 宏,石田 淳,清水 裕士発売日: 2019/12/01メディア: 単行本（ソフトカバー） 「社会科学のための」と謳われているように具体的な政策分析事例を取り上げ、加えて「文系人間」でも数式展開を追えるよう丁寧に説明される。統計モデリングに関しては、いまも一般化線型モデル（GLM）を最尤法で推定するのが一般的かつ説明も容易であると思われるが、本書をひと通り読むと、（GLMでは得られない）ベイズ統計モデリングの「威力」を徐々に垣間見ることができるようになる。（各種の記述統計から回帰分析、Ｐ値へと進む）一般的な頻度主義統計学に…

本記事では、重回帰分析について記載します。 【目次】 計算式等 ノート プログラムコード 雑談 参考 単回帰分析については次の記事を参照ください。 cochineal19.hatenablog.com 計算式等 重回帰分析は、多変量 { } の関係を直線で求める分析方法です。モデル式は次のとおりです。 は切片（intercept）、 は回帰係数（regression coefficient）であり、回帰係数は x が１増えた時に y がどのくらい増えるかを表します。回帰係数の式は、「回帰係数」＝「説明変数間の分散共分散行列の逆行列」×「説明変数と目的変数の共分散行列」です。 なお、正規方程式か…

本記事では、単回帰分析について記載します。 【目次】 計算式等 ノート プログラムコード 雑談 参考 計算式等 単回帰分析は、２変量 {x,y} の関係を直線（一次式）で求める分析方法です。モデル式は次のとおりです。 は切片（intercept）、 は回帰係数（regression coefficient）であり、回帰係数は x が１増えた時に y がどのくらい増えるかを表します。また、説明変数で表現できないものを残差 と言い、この残差は独立に正規分布 に従います。 単回帰分析は次の分散分析表で表現することができます。この分散分析表は一元配置分散分析で示したものと同じです。 ■ 分散分析表 変…

この記事で、RESASとは別にV-RESASというウェブサイトも公開されていることを初めて知りました。RESASは内閣府のまち・ひと・しごと創生本部が設けているビッグデータが公開されているサイトですが、PCでの閲覧を前提としたサイトなので、スマートフォンやタブレットでは閲覧し辛いというところが難点でした。 ですが、このV-RESASはスマートフォンでも閲覧可能ですし、国全体や都道府県ごとのデータをビジュアルで即座に視認出来るという点では優れたシステムだと思います。後はどれだけデータの充実が図られるか次第ですね。 統計とそれに基づく分析は、国や地方の政策立案、また、企業の経営戦略にも大きな影響を…

本記事では、クロス表の連関係数について記載します。 【目次】 計算式等 ノート プログラムコード 雑談 参考 計算式等 次のクロス表で考えます。 Outcome Total 0 1 Factor 1 a b AB 2 c d CD Total AC BD N ファイ係数（Phi coefficient） -1～１の値を取り、絶対値が１に近いほど関連性が大きい。 ２×２のクロス表で用いられ、本係数の絶対値はクラメールの連関係数と一致する。 ピアソンの連関係数 （Pearson's contingency coefficient）または一致係数 ※ は連続修正なしの値を用いる。 ０～１の値を取り、…

本記事では、相関係数について記載します。 【目次】 計算式等 ノート プログラムコード 雑談 参考 計算式等 2変量を x,y とし、x の順位を z、y の順位を w として考えます。タイ（同一順位）は平均順位を用います。 ピアソンの積率相関係数（Pearson product-moment correlation coefficient） スピアマンの順位相関係数（Spearman's rank correlation coefficient） ケンドールの順位相関係数（Kendall rank correlation coefficient） ※ただし、i を基準として、i＜ j の組合…

本記事では、クラスカル・ウォリス検定（Kruskal-Wallis test）について記載します。 【目次】 帰無仮説、対立仮設 計算式等 ノート プログラムコード 雑談 参考 Wilcoxonの順位和検定とマンホイットニーのU検定（以降、Wicoxonの順位和検定等と言います）に続く記事です。 cochineal19.hatenablog.com Wilcoxonの順位和検定等は２標本の分布の位置のズレを評価する検定でしたが、クラスカル・ウォリス検定（KW検定と略称します）は３標本以上の分布の位置のズレを評価するノンパラメトリックな手法です。 帰無仮説、対立仮設 帰無仮説 全ての標本の確率分…

本記事では、Wilcoxonの順位和検定（Wilcoxon rank sum test）とマンホイットニーのU検定（Mann–Whitney U test）について記載します。 【目次】 帰無仮説、対立仮設 計算式等 ノート プログラムコード 雑談 参考 両方とも２標本の分布の位置のズレを評価するノンパラメトリックな手法です。これらは等価な手法であり、同じ結果が得られます。 一方、パラメトリックな手法としては t検定が代表的です。 cochineal19.hatenablog.com 帰無仮説、対立仮設 帰無仮説 ２標本の確率分布に差がない（ズレていない）。 対立仮設 ２標本の確率分布に差があ…

統計の基礎についてまとめてみようと思います. どの分野でも必要ですし知っていると役に立つことばかりです。確率というとサイコロの話とかがまず思い浮かぶでしょうが,大学では連続的な関数を扱うこともできます。連続型確率変数と離散型確率変数、この二つを扱いますが,離散型の∑は連続型の積分と対応しています。離散型の扱いには慣れていると思うので主に連続型で説明をします. 密度関数 密度関数というものを考えます. 密度関数の性質は特に連続型の場合、(離散型の場合は積分ではなくシグマになります.)まあつまりどういうことかっていうと, たとえばｘが1以上3以下の確率は という形で表されます。 密度関数はとなるの…

本記事では、Cochran-Armitage test（コクラン-アーミテージ検定またはCA検定）について記載します。 【目次】 帰無仮説、対立仮設 計算式等 ノート プログラムコード 雑談 参考 CA検定は、割合の増加/減少傾向を評価する傾向性の検定です。 検定の中身は、割合を目的変数とした回帰分析であり、回帰直線の傾きが０であるかどうかを見ています。 なお、目的変数が連続値の場合は、Jonckheere-Terpstra trend test（ヨンクヒール-タプストラ検定、またはヨンキー検定）等が用いられます。 【統計】ヨンキー検定（傾向性検定） - こちにぃるの日記 帰無仮説、対立仮設 …

健康食品サイト男性の美容行為について、乳酸菌の中でもサプリメントや健康食品を徹底解剖。1日の必要量国内製造のプレミアムペット新着サプリメントなど情報満載です。利用者の口コミ・評価の統計全体重に対する割合はサプリメントの活用法3。医師のアドバイスのもと脂肪の吸収を抑制するとされるケルセチン。このサイトの使い方ベース・サプリメントをお選びください。飲んでクセがなく人間の免疫の観点から、アイアンマンマガジンでは、病気予防サプリなどがあります。けんこうしょくひん、マイアカウントの情報がありません。100%愛知産のふきエキスが、妊娠12週なのですが600円を頂戴致します。アカツメクサ花部中国で製造されて…

コンピ指数から見た京成杯のレース傾向です。 今回の京成杯のコンピ指数と同程度の指数配置であった過去のレース結果を照らし合わせて、出走馬の三連複と馬連の馬券絡み率を表にしました。 表はあくまでも、今回の京成杯と同じような指数のレース結果がどういう傾向であったかという統計データであって、 京成杯自体の過去実績ではないのですが傾向に近い結果になることが多いので、読者様のレース予想の後押しになればと思います。 （表の見方としては三連複フォーメーションと馬連フォーメーションをイメージしたものであり、各欄のどこかで10％を超えている部分のある馬が体感的には馬券に絡みやすいです） コンピ 指数 馬番 3連複…

どうも たついち（@4YIns）です。 2020年は、在宅勤務という今までは考えられなかったテレワークなど、働き方が大きく変わった1年でした。2021年に入っても再び在宅勤務に移行するなどと今なお状況は大きくは変わっていません。 そうなってくると、持ち運びに便利なPCやタブレットを求めるようになりました。僕の中では小型でロングバッテリーなものに注目していたのです。 何に注目していたかって？ 昨年頃から大々的に日本国内でも普及している感のあるChromebookです。 Chromebookは、コスパが非常に優れており、動作性が快適という特徴を持っています。 そんな中でも5万円を切る価格帯でありな…

誰が見ても、「先手先手」ではなく「後手後手」である。菅首相が追い詰められているということである。だが、首相を非難してもしょうがない。それが現実であるから、現実の中で最善手を模索するしかない。 統計学の教える所では、コロナについて天気予報のように先を読む材料がいくつかある。その一つが、感染者数と死亡者数の関係である。現在のところ、東京の死亡者数は、昨年に比べて大して増えていない。これで安心してよいのか。 統計学者は、第２波の例を引用して、感染してから重症化して死亡に至る期間にずれがあるから、そのずれを考慮して統計の数値を見なければいけないと教える。 そのずれとは、第2波の例から40日であるとする…

領域範囲階層レベル体系 第１位＝第４位＝第７位＝第１０位 第２位＝第５位＝第８位＝第１１位 第３位＝第６位＝第９位＝第１２位 業界体系 日本大学＝日本大学高等部＝日本大学中等部＝日本大学初等部 日本大学＝法律事務所＝サイマル＝グローバルキッズ（未定） 日本大学＝マルハニチロ＝日本水産＝グローバルキッズ（未定） 職業体系 大学教授＝高等学校教諭＝中学校教諭＝小学校教諭 音楽教師＝法律家＝通訳＝幼稚園教諭 美術教師＝食品衛生管理者＝食品衛生責任者＝保育士 企業体系 営業職＝一般事務職＝経理職＝企画職＝広報職＝マーケティング職＝人事職 順位体系 第１位＝人間力あり＝社会力あり＝学際力あり＝教育力あり…

いきなりですがクイズです。 干しシイタケを日本一消費するのはどこでしょう？ 正解は 消費量 １位大分市、 ２位鹿児島市、 ３位名古屋市。 出汁をとったり煮物で炊いたり、栄養価たっぷりで和食に欠かせない食材。「干ししいたけ」です！ １位の大分市は、いわずと知れた、干ししいたけの一大生産地です。大分を一度訪れたことのある人はご存じかもしれませんが、お土産屋さんには干ししいたけの商品がずらりと並び、車で山道を通るとしいたけの原木を見かけることがあります。 【クイズ】「干し椎茸」消費量 3位名古屋、2位鹿児島、1位は…? | sotokoto online（ソトコトオンライン） では、干しシイタケ生産…

1週間ぶりの更新となりました。何かとバタバタ続きでしたがエントリーは継続しておりました。東京都では緊急事態宣言が再度発令されましたが依然として高止まりしています。お願いベースだとそこまでの強制力は無くなってきていますね。 今週1週間の結果というとだいぶ損を出してしまいました。こちらの思惑と逆方向に進んでいくので難しい限りです。米国の影響は強く発言一つで反転することもしばしば。 www.satokibi2.work www.satokibi2.work 含み損の結果 米国の雇用統計調査で相当な含む損を抱えていたので損切りを実行することにしました。 たしか前回時点では▲4万くらいでしたがズルズルと…

モンゴル国北部フブスグル県で2021年1月12日にマグニチュード6.5～6.8の地震が発生、震源地周辺で建造物への亀裂や地割れ等が生じました。また現地報道からは、今後より大規模な地震が発生した際の懸念も露わになっています。

犬を座敷に上げて飼っては行け無い理由が有る、犬の山水画師 家を抜け出した時に、他家の座敷に上がり込み、見事な山水画を描いて仕舞う。山水画は桂林が有名。犬の山水画師で有る。 最近、毛筆の様なボールペンが販売されて居る。習字の練習がボールペンで出来る。漢字の練習は呆け防止にも成る。朝鮮は反日の為、日帝残滓で漢字を廃止し、論理的な思考が出来無く成り、理性的な行動が取れ無く成り。常識の欠落、記憶力の減退、学力の低下、知能の低下、言語能力の退化が始まっって、認知症が増えて居るものと思い込んで居たら、何と日本人が世界一で有った。日本は自由平等博愛をＧＨＱから教わり、平和で世界一安全で衛生的なくにで有る事が…

2021年1月15日金曜日 ポンドドル・ポン円・ユーロドル・ユーロ円 どうも最近チキン利食いが酷くなってきて、益々値幅が伸びないんじゃないか？もしかして逆に行かれるんじゃないか？と不安が増してきています。原因の心当たりはひとつ、それはロットが上がってきているからです。これも慣れなんですが、心配になりだすとどんどん意識がそのように増していきますので、この週末でもう一回考えを整理します。つまりやるべき事、決める事、やらないように注意する事、逆に行かれたらする事、最悪の事。だからこのトレードで統計的に確率がまあまあ高いと思われるからそうそること。メンタルの要素とはこういう事なんです。 それで金曜日の…

本記事では、ロジスティクス回帰分析について記載します。 【目次】 計算式等 ノート プログラムコード 雑談 参考 計算式等 ロジスティック回帰は二値分類（）で用いられる分析手法です。モデル式は回帰分析と似ていますが、シグモイド関数を用いて0～1の範囲の結果を出すことが特徴です（通常の回帰分析では0～1の範囲を超えてしまい、解釈不能な結果になってしまいます）。（青線：ロジスティック回帰、緑線：単回帰分析） モデル式、シグモイド関数および最適な回帰係数 β を取得するための損失関数はそれぞれ次のとおりです。 モデル式 シグモイド関数（活性化関数） 交差エントロピー誤差（分類問題の損失関数） ※ l…