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

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

since2016 ときどきTEXのメモ

グラフ

ヒストグラムのビン調整

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

相関係数のイメージ

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

正規分布の色塗り (polygon関数 )

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

共分散分析を考えるためのグラフ

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

独立している2群を比較するグラフ ggplot2

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

stripchartの横にエラーバー

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

正規分布の重ね描き

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

級内相関係数 Case1

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

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

移行しました y2pt.com

回帰直線と交互作用

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

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

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

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

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

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

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

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

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

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

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

2変量の正規分布をグラフでイメージ(persp)

また、ここで勉強させていただきました. http://cse.naro.affrc.go.jp/minaka/R/R-binormal.html 忘れないように要点のみ転記させていただます.まさに備忘録. 今回はRの関数perspを使用して、密度関数の数式から3Dのグラフを描いてみます 確率変数x1…

2変量の正規分布をグラフでイメージ(scatterplot3d)

ここで勉強させていただきました. http://cse.naro.affrc.go.jp/minaka/R/R-binormal.html 忘れないように要点のみ転記させていただます. 必要なパッケージをインストールします install.packages("mvtnorm") library(mvtnorm) install.packages("scatterp…

信頼区間のプロット

同じサイズのデータサンプルからt分布を利用した信頼区間の作図 まずは3×4の場合(サンプルサイズ3を4回実施する) x <- matrix(NA,nrow=3,ncol=4) #3×4の空セル for (i in 1:4){ #列数分乱数を代入 x[,i] <- rnorm(3) #標準正規分布の乱数を行数分繰り返…

クラスター分析

Rのサンプルattitudeを使用して、クラスタ分析(階層的方法)を勉強します. attitudeは、管理者態度のデータです.無作為に選ばれた35名の雇用者よりアンケート.好意的な割合が数値化されています. rating全般的評価、 complaints雇用者からの苦情処理、 …

指数分布

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

残差分析

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

F分布の確率密度関数

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

t分布からの信頼区間

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

x軸操作いろいろ

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

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")…

グラフに色をつける

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…

幹葉図、箱ひげ

# 箱ひげ、パーセンタイル説明のために偏ったデータを作成 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, …