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

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

since2016 ときどきTEXのメモ

sample

解析の練習で使うデータサンプル集 随時更新

sample1

症例 <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
性別 <- c("男性","男性","女性","男性","男性","女性","男性","男性","男性","男性","女性","女性","女性","女性","男性","女性","女性","女性","女性","男性")
治療 <- c("B","A","B","A","A","B","B","A","A","B","A","B","B","B","A","B","B","B","A","A")
治療前BP <- c(160,135,177,141,142,155,175,145,149,155,135,156,170,150,138,160,150,180,155,140)
治療後BP <- c(148,133,160,138,139,145,157,142,142,145,132,145,157,142,134,150,140,160,149,136)
BP変化量 <- c(12,2,17,3,3,10,18,3,7,10,3,11,13,8,4,10,10,20,6,4)
sample1 <- data.frame(症例,性別,治療,治療前BP,治療後BP,BP変化量 )

sample2
対応のある2群

ID <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10)
pre <- c(18, 8, 14, 15, 16, 13, 23, 20, 14, 20)
post <- c(23, 27, 15, 22, 15, 20, 20, 20, 18, 25)
sample2 <- data.frame(ID, pre, post)

sample3
独立した2群

グループ <- c("A","A","B","B","A","B","A","B","B","A","A","A","B","A","B","B","B","B","A","A")
ROMT <- c(18, 8, 23, 27, 14, 15, 15, 22, 15, 16, 13, 23, 20, 20, 20, 20, 18, 25, 14, 20)
sample3 <- data.frame(グループ, ROMT)

sample4
クロス表

sample4
判定 <- c(rep("非該当",30), rep("要支援", 25))
歩行 <- c(rep("可能", 20), rep("不可能", 10), rep("可能", 10), rep("不可能", 15))
sample4 <- data.frame(判定, 歩行)

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

   症例 性別 治療 治療前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   10 男性    B      155      145       10
11   11 女性    A      135      132        3
12   12 女性    B      156      145       11
13   13 女性    B      170      157       13
14   14 女性    B      150      142        8
15   15 男性    A      138      134        4
16   16 女性    B      160      150       10
17   17 女性    B      150      140       10
18   18 女性    B      180      160       20
19   19 女性    A      155      149        6
20   20 男性    A      140      136        4

上記テーブルの性別を切り出し
テキスト形式のベクトルに変換

data <- read.table("clipboard", header = T)
性別2 <- factor(data$性別)

install.packages("stringr")
library(stringr)

性別3 <- str_replace_all(性別2, pattern="k", replacement="*")
#性別3をwordで置換 スペースを,に変換

性別 <-c("男性","男性","女性","男性","男性","女性","男性","男性","男性","男性",
       "女性","女性","女性","女性","男性","女性","女性","女性","女性","男性")

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

\LaTeX

PDFにページ数を挿入する場合
PDFのファイル名はNo1,No2,No3の三枚
同じフォルダに入れる

\documentclass[uplatex,11pt]{jsarticle}
\usepackage[top=15truemm,bottom=15truemm,left=15truemm,right=15truemm]{geometry}     %余白

%図
\usepackage[dvipdfmx]{graphicx}
%図の位置調整 改行後に[H]
\usepackage{here}  

% 章が進むごとに図番号をリセットする
\makeatletter% プリアンブルで定義開始
\renewcommand{\thefigure}{\arabic{figure}}
\@addtoreset{figure}{section}
\makeatother % プリアンブルで定義終了

%図とキャプションを狭くする
\setlength\abovecaptionskip{5pt}

%---------------------------------------------------%
\begin{document} 
%---------- 1 ----------%
\begin{figure}[H] %現在の位置に
  \centering
    \includegraphics[width=17cm]{No1.pdf} %図の挿入
\end{figure}
%---------- 2 ----------%
\newpage
\begin{figure}[H] %現在の位置に
  \centering
    \includegraphics[width=17cm]{No2.pdf} %図の挿入
\end{figure}
%---------- 3 ----------%
\newpage
\begin{figure}[H] %現在の位置に
  \centering
    \includegraphics[width=17cm]{No3.pdf} %図の挿入
\end{figure}

\end{document}