R で confit() 関数を使用する方法


R のconfint()関数を使用すると、近似回帰モデル内の 1 つ以上のパラメーターの信頼区間を計算できます。

この関数は次の基本構文を使用します。

confint(オブジェクト、パラメータ、レベル=0.95)

金:

  • object : 近似された回帰モデルの名前
  • parm : 信頼区間を計算するパラメータ (デフォルトは all)
  • level : 使用する信頼レベル (デフォルト値は 0.95)

次の例は、この関数を実際に使用する方法を示しています。

例: R で confit() 関数を使用する方法

R に、クラスの 10 人の生徒の学習時間数、受験した模擬試験の数、および最終試験の得点を示す次のデータ フレームがあるとします。

 #create data frame
df <- data. frame (score=c(77, 79, 84, 85, 88, 99, 95, 90, 92, 94),
                 hours=c(1, 1, 2, 3, 2, 4, 4, 2, 3, 3),
                 prac_exams=c(2, 3, 3, 2, 4, 5, 4, 3, 5, 4))

#view data frame
df

   score hours prac_exams
1 77 1 2
2 79 1 3
3 84 2 3
4 85 3 2
5 88 2 4
6 99 4 5
7 95 4 4
8 90 2 3
9 92 3 5
10 94 3 4

ここで、次の重線形回帰モデルを R に当てはめたいとします。

試験のスコア = β 0 + β 1 (時間数) + β 2 (実技試験)

lm()関数を使用して、このモデルを適応させることができます。

 #fit multiple linear regression model
fit <- lm(score ~ hours + prac_exams, data=df)

#view summary of model
summary(fit)

Call:
lm(formula = score ~ hours + prac_exams, data = df)

Residuals:
    Min 1Q Median 3Q Max 
-2.4324 -1.2632 -0.8956 0.4316 5.1412 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 68.4029 2.8723 23.815 5.85e-08 ***
hours 4.1912 0.9961 4.207 0.0040 ** 
prac_exams 2.6912 0.9961 2.702 0.0306 *  
---
Significant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 2.535 on 7 degrees of freedom
Multiple R-squared: 0.9005, Adjusted R-squared: 0.8721 
F-statistic: 31.68 on 2 and 7 DF, p-value: 0.0003107

モデルの概要には、近似された回帰係数が表示されることに注意してください。

  • 切片 = 68.4029
  • 時間 = 4.1912
  • prac_exams = 2.6912

これらの各係数の 95% 信頼区間を取得するには、 confint()関数を使用できます。

 #calculate 95% confidence interval for each coefficient in model
confined(fit)

                 2.5% 97.5%
(Intercept) 61.6111102 75.194772
hours 1.8357237 6.546629
prac_exams 0.3357237 5.046629

各パラメーターの 95% 信頼区間は次のように示されます。

  • 切片の 95% CI = [61.61, 75.19]
  • 時間の 95% CI = [1.84, 6.55]
  • prac_exams の 95% CI = [0.34, 5.05]

99% 信頼区間を計算するには、 level引数の値を変更するだけです。

 #calculate 99% confidence interval for each coefficient in model
confint(fit, level= 0.99 )

                 0.5% 99.5%
(Intercept) 58.3514926 78.454390
hours 0.7052664 7.677087
prac_exams -0.7947336 6.177087

また、特定のパラメーターの信頼区間のみを計算するには、 parm引数を使用して係数を指定するだけです。

 #calculate 99% confidence interval for hours
confint(fit, parm=' hours ', level= 0.99 )

          0.5% 99.5%
hours 0.7052664 7.677087

99% 信頼区間は時間変数についてのみ表示されることに注意してください。

追加リソース

次のチュートリアルでは、R の線形回帰に関する追加情報を提供します。

R で回帰出力を解釈する方法
R で単純な線形回帰を実行する方法
R で重回帰を実行する方法
R でロジスティック回帰を実行する方法

コメントを追加する

メールアドレスが公開されることはありません。 が付いている欄は必須項目です