単回帰分析について。 【目次】 計算式等 ノート プログラムコード 回帰係数の詳細 参考 計算式等 単回帰分析は、２変量 {x,y} の関係を直線（一次式）で求める分析方法。モデル、式は次のとおり。 ■ 単回帰モデル ■ 単回帰式 は切片（intercept）、 は回帰…

クロス表の連関係数について記載。 【目次】 計算式等 計算例 プログラムコード 参考 計算式等 次のクロス表で考える。 Outcome Total 0 1 Factor 1 a b AB 2 c d CD Total AC BD N ファイ係数（Phi coefficient） -1～１の値を取り、絶対値が１に近いほど関…

相関係数について。 【目次】 計算式等 計算例 プログラムコード ピアソンの相関係数と順位相関係数 参考 計算式等 2変量を x, y とし、x の順位を z、y の順位を w として考える。タイ（同一順位）は平均順位を用いる。 ピアソンの積率相関係数（Pearson pr…

統計関係の記事いくつか書いているので、サイトマップ作りました。 差の評価 傾向性の評価 関係性の評価 多変量解析 生存時間解析 あとがき 差の評価 対応有無 目的変数 群数 検定 記事 対応なし 連続値 ２群 t検定（Student、Welch） 記事あり ２群以上 一…

クラスカル・ウォリス検定（Kruskal-Wallis test）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード 2標本でもKW検定 参考 クラスカル・ウォリス検定（KW検定）は３標本以上を評価するノンパラメトリックな手法。Wilcoxonの順位和…

Wilcoxonの順位和検定（Wilcoxon rank sum test）とマンホイットニーのU検定（Mann–Whitney U test）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード W統計量とU統計量の関係 参考 両方とも２標本の分布の位置のズレを評価するノ…

Cochran-Armitage test（コクラン-アーミテージ検定またはCA検定）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード その他 参考 CA検定は、割合の増加/減少傾向を評価する傾向性の検定。 検定の中身は、割合を目的変数とした回帰…

マンテル・ヘンツェルの検定（Mantel-Haenszel test）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード CochranとMantel and Haenszelの違い 参考 マンテル・ヘンツェルの検定は、 検定やFisher's Exact testと同じく、カテゴリー…

Fisher's exact testについて。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード 参考 Fisher's Exact testは 検定と同じく、カテゴリーデータの独立性を評価する手法。 cochineal19.hatenablog.com 帰無仮説、対立仮設 帰無仮説 対立仮設 要…

独立性の 検定について。 【目次】 帰無仮説、対立仮設 計算式等 ノート プログラムコード Cochran's rule 参考 検定はカテゴリデータのクロス表（分割表、クロステーブル）で用いられる検定。 帰無仮説、対立仮設 帰無仮説 対立仮設 要因と結果の間に関係が…

傾向性検定（ヨンキー検定）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード ヨンキー検定とケンドールの順位相関 参考 Jonckheere-Terpstra trend test（ヨンクヒール-タプストラ検定、またはヨンキー検定）は目的変数が連続値の…

LaTeXという組み込みソフトで数式を書けるようです。 [tex: 数式] で表します。 数式の作成には、以下のサイトさんがとても使いやすいです。手書きした数式をコードに変換してくれて、とても便利です。 webdemo.myscript.com 試しに適当な数式を書いてみまし…

一元配置分散分析（One-way ANOVA）について。 【目次】 帰無仮説、対立仮設 計算式等 計算例 プログラムコード 残差平方和と水準間平方和の違い 参考 t検定 は２群の平均の差を評価するもの。３群以上の場合は、分散分析（ANOVA：ANalysis Of VAriance）が…

2標本の平均の差の検定（t検定）について。 【目次】 帰無仮説、対立仮設 計算式等 プログラムコード 参考 帰無仮説、対立仮設 帰無仮説 対立仮設 群間の平均に差があるかどうかを評価する。帰無仮設の「群aと群bの平均に差がない」は「平均の差がゼロである…

タイトル通り、基本統計量の算出方法です。SQL Serverには、中央値を算出するmedian関数はありません ... （同じMicrosoftのExcelにはあるのに ...）（SASのproc sqlではmedian関数が使えるのに ...）とがっかりしましたが、どうやらpercentile_cont関数で対…

Rに限らず、プログラムファイルを見ると、いきなりコードを書き始める人もいますが、私は環境設定・データ取得・加工・出力はなるべくまとまりにすべきだと思っています。（そのほうがメンテナンス性が良いと思います。） ということで、R用テンプレートのメ…

Rscript.exeを使ったバッチ実行ツールの記事の続き。 記事その１はこちら。 【R言語＆Javascript】Rスクリプトをバッチ実行する（その１） - こちにぃるの日記 Rscript.exeで設定した引数の受け取り方法を解説。引数の設定は、次のとおり。 Rscript.exe [Rス…

Rscript.exeを使ったバッチ実行ツールの記事。 コマンドプロンプト上で、次のコマンドで実行できるが、Rスクリプトのファイルパスや引数情報を毎回手入力するのは面倒なので、javascriptを橋渡し役として使っている。（このコマンドだけではログ出力されない…

バックアップを作成するバッチをjavascriptで書いてみました。 使い方は簡単で、本記事のコードをJS形式で保存し、任意のフォルダに配置すればOKです。ダブルクリックで実行し、同一階層にあるファイルのバックアップをとります。 【処理内容（実行されると…

VBAからSQL Serverに接続してクエリ結果を取得します。 Microsoft ActiveX Data Objects (ADO)を使用したDBへのアクセス方法です。 本記事はWindows認証による方法です。 以下の図のように、C2～4セルにサーバ名、データベース名、SQL文を入れて実行すると、…

さて、R関係の初記事です。 もともとSASを使っていてSASマクロを使ったループが大変便利だったわけです。 Rはそこまで柔軟ではないものの、次の方法でループ処理できます。 loop1 <- c("A","B","C") A <- c(1,1,1) B <- c(2,2,2) C <- c(3,3,3) for (i in 1:…

Hatena Blogにコードを書いたとき、 背景色をつけられるかのテスト。 <pre style="background:色コード; color:色コード;"><code>data A2; set A1;run;</code></pre> 背景：青、文字：白 style="background: #0000ff; color: #ffffff; font-size: 80%;" data A2; set A1; run; 背景：薄緑、文字：濃い緑 style="background: #CCFFCC…

PDFのパスワード、アクロバットリーダーでは解除できなくて困ったのですが、 Pythonのpikepdfで簡単に解除できるようなのでメモ。 from pikepdf import Pdf mypdf = Pdf.open('パス/ファイル名.pdf', password='パスワード') outPDF = Pdf.new() outPDF.page…

Hatena BlogはHTML編集モードで直接HTMLをいじくれるようです。 ということでテストです（個人用）。 テーブル <table> <tr> <th>タイトル1</th> <th>タイトル2</th> </tr> <tr> <td>あ</td> <td>い</td> </tr> <tr> <td>う</td> <td>え</td> </tr> </table> タイトル1 タイトル2 あ い う え hrタグ <hr/> codeタグ&preタグ <pre><code> proc sql noprint; count(*) into :cnt </code></pre></hr/>…

Excelの文字切れ対策。 以下、簡易な方法。万能ではないが。 ①一度MSゴシックにしてフォントサイズでかくする 。②オートフィットしてから行幅を1.2倍（任意）する。（+15など行幅を足しても良い）③任意のフォントに戻す。 Sub 文字切れ対策() Dim baseFontNa…

成果物はA1セルをアクティブにするご作法。作業中ではないですよ。の無言の表現。 forループを使ったシンプルな対処法です。シート数 to 1（-1ずつ）で１シートずつA1セルをアクティブにしましょう。 ここで注意なのがwindowの固定。A1セルをselectしただけ…

アクティブなシートにシート名一覧を作成。 セル(2,2)から順々にシート名入れていく。 Sub シート名一覧化() Dim wb1 As Workbook, ws1 As Worksheet, i, cnt As Long Set wb1 = ThisWorkbook Set ws1 = wb1.Worksheets(ActiveSheet.Name) cnt = 2 For i = 1…