以下のようなグラフ 忘れないように記載しておきます データサンプル Group <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3) p…

data2 <- data.frame(c(rep("あり",15),rep("なし",15)), c(rep("有効",12),rep("無効",3),rep("有効",4),"無効","有効","有効",rep("無効",8))) colnames(data2) <- c("筋トレ","効果") xtabs(~ 効果 + 筋トレ, data = data2) 筋トレ 効果 あり なし 無効 3…

使用するパッケージ install.packages("ggplot2") library(ggplot2) データセットを作成します #乱数を使用のため、実施毎にグラフは異なります dat <- c(c(rnorm(15)),c(rnorm(15,3,1.2)),c(rnorm(15,10,2)),c(rnorm(15,17,1.5))) pre_post <- c(rep("前", …

データの準備 dat <- c(c(rnorm(15)),c(rnorm(15,3,1.2)),c(rnorm(15,10,2)),c(rnorm(15,17,1.5))) pre_post <- c(rep("前", 30),rep("後", 30)) treat <- c(rep("A",15),rep("B",15),rep("A",15),rep("B",15)) dataf <- data.frame(dat, pre_post, treat) h…

2018-07-19投稿、最終更新日2019-07-07 下のようなカテゴリカルデータをダミー変数に変更します 治療 効果 NO NO YES YES YES NO NO YES YES YES YES YES NO YES YES NO NO NO NO YES NO NO NO NO YES YES YES YES YES YES YES NO NO NO NO NO YES YES NO NO…

いつも勉強させていただいおります、下記サイトの共分散分析より グラフの描き方の備忘録です 我楽多頓陳館--雑学と統計学の館 以下のようなグラフになります データのサンプルは上記サイトから引用してdataに入れました # まずはこれを描いて考える・・・！…

サンプルはRのmtcarsを使用します．車に詳しくないので回帰分析に相応しいサンプルか分かりませんが、学習のため数値のみ使用させていただきます． サンプルから直接取り込み dat <- mtcars[,c(8,2,5,6)] #ここがポイント8列目を応答変数として、2列目･5列目…

投稿日:2018.2.13 最終更新日:2018.5.17 ロジット関数とロジスティック関数 準備として関数の特徴を押さえておきます ロジット関数 (標準)ロジスティック関数ロジット関数の逆関数= mathwords.net サンプルirisより 品種"virginica＝1"、"別の品種＝0"という…

決定係数 回帰係数がどの程度よく当てはまっているか、つまりXがYをうまく説明できているか、ということを明らかにすることは重要になります．決定係数は、その当てはまりを示す基準としてして使用されています． 観測値の平方和（全変動） 回帰による平方和…

投稿日2017.2.2 更新日2017.11.13 線形回帰分析 説明変数と目的変数との関係を直線でモデル化する回帰分析。 Xi <- c(35, 45, 55, 65, 75) Yi <- c(114, 124, 143, 158, 166) plot(Xi, Yi, pch = 16) abline(lm(Yi~Xi), col = 2) どのようにして赤のラインを…

標準偏回帰係数の求め方 投稿日2017.8.12 更新日2017.8.13 忘れないように載せておきます 例）予測変数が二つの場合の重回帰分析 母回帰方程式 Y = b0 + b1*X1i + b2*X2i + Ei 予測方程式 y = β0 + β1*x1i + β2*x2i +ei 標準偏回帰係数 sβ1、sβ2 sβ1 = ( yと…

トービット回帰直線 投稿日2017.8.4 打ち切りデータの場合、つまり天井効果や床効果が生じているデータの場合には、トービット回帰直線で分析します． 参考書はもちろん豊田秀樹 (著, 編集);回帰分析入門 (Rで学ぶ最新データ解析) ,東京図書 ,2012 Rのサンプ…

残差分析 Rのサンプル ChickWeight を使用y<-ChickWeight$weight[1:10]x<-ChickWeight$Time[1:10]summary(lm(y~x))plot(x,y,xlab = "生後日数",ylab = "体重")abline(30.327 ,7.030) #回帰直線の挿入 残差 ei をプロットすることで、回帰モデルからのズレを…

投稿日2017.2.3更新日2017.7.31 回帰係数の区間推定 標本回帰係数の不偏推定量から母回帰係数の信頼区間を求めてみます． 次式をRにペーストしてx<-read.table("clipboard",header=T) 次のデータをコピーしてRで実行します東京 福岡1019.4 1018.41005.7 1007…

Rのデータセットtrees を使用して説明します Girth=y , Height=x1 , Volume=x2 次のようにエクセルに準備します x<-read.table("clipboard",header=T) をコピーしてR Consoleにペーストします 分析したい一覧表をコピー（上記の表など） x<-read.table("clip…