理学療法士がまとめたノート

統計学備忘録(R言語のメモ)

since2016 ときどきTEXのメモ

対応のあるデータのグラフ

移行しました y2pt.com

回帰直線と交互作用

下記のサイトに移行しました y2pt.com

基準値の変更

下記のサイトに移転いたしました y2pt.com

Tex 図を並べる

TeX

図を縦に二つ並べる \begin{figure}[H] \centering \includegraphics[width=40mm]{ファイル01名.jpg} \caption{図の説明01} \end{figure} \begin{figure}[H] \centering \includegraphics[width=40mm]{ファイル02名.jpg} \caption{図の説明02} \end{figure} …

カッパ係数

投稿日2018.3.12 最終更新日2020.10.13 個人的によく利用させていただいております以下のHPをもとに、今回はカッパ係数について少し勉強してみます 統計学入門−第5章 まずはHPに掲載してある次のサンプルデータを使用して、Rを使って処理してみます 分類数が…

ANOVA君 で分散分析

投稿日2017.8.25, 更新日2020.8.17 ANOVA君はフリーの統計ソフトウェア「R」で動作する分散分析関数です。ダウンロード、使い方、リリース情報などは「井関龍太のページ」を必ずご確認ください. ANOVA君 - 井関龍太のページ このブログではanovakun version …

分散分析(繰り返し vs 反復)

登校日:2017.11.06, 更新日:2020.08.12 「繰り返し」と「反復」は区別なく混同して使用される場合もあります.要点のみ記載しておきます。 繰り返し 実験内容としては「水準の設定を最初から施して既定のサンプル数を繰り返してとる」ことになります. 「繰り…

2標本のグラフ、記述統計

下のサイトに移転しました y2pt.com

単回帰分析を一括で実行する方法

tokei.net 全人類が分かる統計学より to-kei.net 上記サイトに自分なりの注釈をつけました。 sampleもまったく同じです。 番号 <- c(1:30) 年齢 <- c(22,23,24,25,27,28,28,29,30,31,32,32,33,33,34,36,37,37,38,39,40,42,46,49,50,53,56,58,64,65) 血圧 <- …

エラーバー付きのグラフ ggplot2

下記のサイトに移転しました エラーバー付きのグラフ ggplot2 - 統計学備忘録

簡単なエラーバーの描き方

下のサイトに移転しました y2pt.com

比較した箱ひげにアスタリスクを入れる

下記のサイトに移転しました y2pt.com

sample

解析の練習で使うデータサンプル集 随時更新 sample1 症例 <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20) 性別 <- c("男性","男性","女性","男性","男性","女性","男性","男性","男性","男性","女性","女性","女性","女性","男性","女性","女性…

テーブルをベクトルに変換

症例 性別 治療 治療前BP 治療後BP BP変化量 1 1 男性 B 160 148 12 2 2 男性 A 135 133 2 3 3 女性 B 177 160 17 4 4 男性 A 141 138 3 5 5 男性 A 142 139 3 6 6 女性 B 155 145 10 7 7 男性 B 175 157 18 8 8 男性 A 145 142 3 9 9 男性 A 149 142 7 10 1…

ブログ名の変更

旧)統計学備忘録since2016 新)統計学備忘録(R言語のメモ) ブログ名だけで趣旨が分かるように変更しました

PDFにページ数を挿入する場合

TeX

PDFにページ数を挿入する場合 PDFのファイル名はNo1,No2,No3の三枚 同じフォルダに入れる \documentclass[uplatex,11pt]{jsarticle} \usepackage[top=15truemm,bottom=15truemm,left=15truemm,right=15truemm]{geometry} %余白 %図 \usepackage[dvipdfmx]{gr…

層別した散布図と回帰直線

以下のサイトに移転しました y2pt.com

開設3周年

はてなブログさんから素敵なメッセージが届きました

表1の作り方 "tableone"

下記のサイトに移転いたしました y2pt.com

if関数 層別カテゴリー化した変数を追加

下記のサイトに移転いたしました よろしくお願いいたします y2pt.com

並べ替え検定と正規近似

2017-07-11投稿, 2019.7.4更新 正確なp値 特定の確率分布をもとに推定を行うのではなく、母集団の未知のパラメータやサンプリング誤差が入らないため計算上も正しいp値が得られる. 並べ替え検定 例)x群とy群を比較します x <- c(5, 9) #平均 = 7 y <- c…

Rで簡単(二元配置分散分析:データには対応なし、繰り返し数は5)

投稿2016.11.17 最終更新日2019.6.5 検定は帰無仮説が大切です! 商品のの主効果:帰無仮説=カップ麺でもインスタント麺でも評価の母平均は等しい スープの主効果:帰無仮説=スープが違っても評価の母平均は等しい 商品とスープの主効果:帰無仮説=商品と…

相関行列と共分散行列

相関行列 x1 <- c(93, 89, 115, 90, 75) x2 <- c(121, 136, 121, 161, 125) x3 <- c(101, 115, 118, 122, 106) #データフレームにすること dat <- data.frame(x1, x2, x3) cor(dat) x1 x2 x3 x1 1.0000000 -0.2220449 0.3981155 x2 -0.2220449 1.0000000 0.6…

対応のある一元配置分散分析(繰り返しのない二元配置分散分析)、多重比較

以下のサイトに移転いたしました y2pt.com

Fisherの直接法

以下のサイトに移転しました y2pt.com

確率分布の計算関数

下記のサイトに移転しました y2pt.com

生存時間解析の基本

パッケージのインストール install.packages("survival") library(survival) survivalのデータセットheartを使用します head(heart) start, stop, event: Entry and exit time and status for this interval of time age: age-48 years year: year of accept…

subset関数で層別化したデータセット

以下のようなデータセットを男女で層別したセットに変更します id 年齢 性 身長 1 39 女 162.6 2 39 女 161.4 3 41 男 174.1 4 41 女 154.7 5 42 男 184.7 6 43 男 169.3 7 43 男 163.2 8 43 男 164.4 9 44 男 167.8 10 44 男 168.5 11 44 女 154.8 12 45 男…

CSVファイルの取り込み

下記のサイトに移転いたしました y2pt.com

データフレームの列を入れ替えて、欠損値のある行を除外する

忘れないうちに書いときます 欠損値のあるデータセットを用意します x <- LETTERS[1:10] y <- c(1,2,NA,2,4,NA,5,1,2,5) z <- c(1,2,2,1,1,2,2,1,1,1) xyz <- data.frame(x,y,z) x y z 1 A 1 1 2 B 2 2 3 C NA 2 4 D 2 1 5 E 4 1 6 F NA 2 7 G 5 2 8 H 1 1 9…