
No.3ベストアンサー
- 回答日時:
preタグが使えないので見にくいですが、、、
g.dat <- as.factor(rep(c(1, 2), c(5, 6))) # グループを表す変数
x.dat <- c(1.3, 1.5, 1.5, 1.8, 2.0, 1.5, 1.5, 1.8, 1.8, 2.0, 2.2) # 説明変数
y.dat <- c(17.3, 15.6, 14.1, 13.8, 10.2, 17.8, 15.8, 16.0, 13.6, 11.2, 10.2) # 目的変数
mydata <- data.frame(Y = y.dat, X = x.dat, group = g.dat) # データフレームの作成
# 関数two.slopeを定義する
two.slope <- function(dat){
attach(dat) # 引数に指定されたデータフレームにパスを通す
# データフレームより各群のデータを抽出
Y1 <- Y[group==1]
X1 <- X[group==1]
Y2 <- Y[group==2]
X2 <- X[group==2]
# 2本の回帰直線を求める
model1 <- lm(Y1 ~ X1)
model2 <- lm(Y2 ~ X2)
# それぞれの残差を求め、残差平方和を求める
m1.resid <- sum(resid(model1) ^ 2)
m2.resid <- sum(resid(model2) ^ 2)
SSr <- m1.resid + m2.resid
# 各群の平均値
y1.M <- mean(Y1)
y2.M <- mean(Y2)
x1.M <- mean(X1)
x2.M <- mean(X2)
# 2つの回帰直線が平行であるとしたときの推定値と残差平方和
sum1 <- sum((X1 - x1.M) * (Y1 - y1.M))
sum2 <- sum((X2 - x2.M) * (Y2 - y2.M))
sum3 <- sum((X1 - x1.M) ^ 2)
sum4 <- sum((X2 - x2.M) ^ 2)
b <- (sum1 + sum2) / (sum3 + sum4)
a1 <- y1.M - b * x1.M
a2 <- y2.M - b * x2.M
sig1 <- sum((Y1 - (a1 + b * X1)) ^ 2)
sig2 <- sum((Y2 - (a2 + b * X2)) ^ 2)
SSrb <- sig1 + sig2
m <- length(Y1) # グループ1のデータ数
n <- length(Y2) # グループ2のデータ数
DF2 <- m + n - 4 # 自由度2
Fb <- (SSrb - SSr) / (SSr / (DF2)) # F統計量の計算
pb <- 1 - pf(Fb, df1 = 1, df2 = DF2) # p値の計算
cat("F = ", Fb, "\n")
cat("p = ", pb, "\n")
detach(dat) # データフレームのパスを解除
}
# 使用例
> two.slope(mydata)
F = 0.1161944
p = 0.7432014
一応、補足しておきますと、帰無仮説は「2本の回帰直線は平行である(傾きが同じである)」というもので、今回の場合はp = 0.74 < 0.05なので「2本の回帰直線は平行である」と結論づけられます。すなわち、傾きが異なるとはいえないということですね。
詳しくは参考文献
岩崎学「統計的データ解析入門 単回帰分析」東京図書 p88-89
を参照してください。
もっとちゃんと探せばこういう関数があると思いますけど(^_^;)
No.2
- 回答日時:
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
-
ネット速度にお困りなら…Wi-Fi 6
高速・安定のWiFi 6で、通信速度の不満を解消!接続機器が多くても安心。だから今買うならWi-Fi 6がおすすめ
-
回帰曲線の有意差の検定
数学
-
回帰直線の傾きの比較
数学
-
エクセルで回帰直線の傾きを検定するには
Excel(エクセル)
-
4
近似直線の傾き同士の相関
経済学
-
5
研究所に応募するときは、「御社」?
その他(就職・転職・働き方)
-
6
エクセルで計算すると2.43E-19などと表示される。Eとは何ですか?
Windows Me・NT・2000
-
7
回帰分析の論文の書き方について質問です。
数学
-
8
敬語についての質問です。 「無理だったら、大丈夫です。」を敬語にするとどうなりますか?
日本語
-
9
相関係数についてくるP値とは何ですか?
統計学
-
10
英語の正式名称と略称を繋げる場合、間に入れるのはコロン?セミコロン?
英語
-
11
ペアワイズとは?
その他(自然科学)
-
12
エクセルで正規分布かどうかを確かめる。
統計学
おすすめ情報
このQ&Aを見た人がよく見るQ&A
人気Q&Aランキング
-
4
切片あり回帰と切片なし回帰
-
5
散布図からの近似曲線の書き方...
-
6
残差について
-
7
回帰式と近似式について
-
8
統計の検定と信頼区間について
-
9
微分方程式と回帰分析の予測方...
-
10
重回帰分析で偏回帰係数を全て...
-
11
数学の問題について質問です。...
-
12
【統計】「相関係数の有意性の...
-
13
サンプル数の違うものの比較
-
14
相関を求める時の、はずれ値の...
-
15
アンケート調査と重回帰分析に...
-
16
質的データと量的データの相関...
-
17
相関係数Rの2乗について
-
18
DataSetの内容がDataGridViewに...
-
19
SPIの問題について質問です。 v...
-
20
重回帰分析におけるR2乗値
おすすめ情報