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

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

since2016 ときどきTEXのメモ

グラフ

RでBland Altman Plotを描く

データ国際医療福祉大学 下井研究室信頼性の検討方法:相対信頼性と絶対信頼性例1(http://shimoi.iuhw.ac.jp/reliability_3_systematicbias.html)を使用させていただきましたx125.1 23.3 33.2 20.9 22.8 34.9 26.5 26.1 29.4 29.8 29.2 26.5 14.3 14.6 17.…

正規分布のグラフ

グラフ初歩の初歩 手探り状態でやってます 赤文字のみRで実行 正規分布#1正規分布の乱数を100個生成x<-rnorm(100)#ヒストグラム,freq=Fで確率密度,ylimでy軸範囲hist(x,freq=F,ylim=c(0,0.6))#枠の処理、lty線種は実践、btyで左と下のみbox(lty=1,bty = "l")…

二項分布のグラフ

こんな説明でどうでしょうか? コイン投げ(いかさまなし)表の出る確率確率1/23回コイン投げを実施してみた場合・・・表が出る回数(x)は0回、1回、2回、3回xの確率は𝐵𝑖(3 , 0.5) に従う.その確率分布を二項分布といいます. 村上 正康,安田 正実; 統計…

グラフに色をつける

library(RColorBrewer)#RColorBrewerパッケージのサンプルdisplay.brewer.all() ヒストグラムを塗ってみるどの色セットを使用するかを指定するcols <- brewer.pal(8,"Pastel1") # brewer.pal(何色、パレット名) y<-c(1,2,3,4,5,6,7)p<-c(2,3,4,5,4,3,2)q<-…

棒グラフの中央に散布図をプロットする方法 (pos.x)

y<-c(1,2,3,4,5,6,7)p<-c(2,3,4,5,4,3,2)q<-c(2,3,4,5,4,3,2)par(mar = c(5, 4, 1, 4)) #余白 底辺、左、上、右の順 pos.x <- barplot(q,ylim=c(0,6))points(pos.x, p)

時系列データの可視化01

下のデータをコピーして時系列グラフを作ってみます 心拍変動を三次元加速度を同時測定した結果です 被験者は私です RRI TEM X Y Z HF LFHF LF activity HR793 21.5 0.16 -1.09 -0.12 75.531 4.247 80.941 0.11 76800 21.5 0.16 -1.09 -0.12 75.436 4.078 80…

中心極限定理をグラフで理解する

#母集団の分布がどのような分布であっても、無作為抽出した標本における和の分布は、標本の大きさnが大きいときに正規分布に収束する。 #母集団の分布がどのような分布であっても、無作為抽出した標本における標本平均の分布は、標本の大きさnが大きいとき…

幹葉図、箱ひげ

# 箱ひげ、パーセンタイル説明のために偏ったデータを作成 x<-c( 5, 8, 11, 23, 33, 35, 41, 44, 44, 44, 44, 46, 47, 49, 50, 51, 52, 52, 52, 53, 55, 55, 55, 56, 56, 56, 58, 58, 58, 59, 59, 59, 60, 60, 60, 61, 61, 62, 63, 63, 64, 65, 65, 65, 65, …