なげなわ法を使用した特徴選択とcovTest::covTest
、p.values を取得するための共分散検定を実施しました。私はそのようなものから例を借りますcovTest
:
require(lars)
require(covTest)
set.seed(1234)
x=matrix(rnorm(100*10),ncol=10)
x=scale(x,TRUE,TRUE)/sqrt(99)
beta=c(4,rep(0,9))
y=x%*%beta+.4*rnorm(100)
a=lars(x,y)
covTest(a,x,y)
$results
Predictor_Number Drop_in_covariance P-value
1 105.7307 0.0000
6 0.9377 0.3953
10 0.2270 0.7974
3 0.0689 0.9334
7 0.1144 0.8921
2 0.0509 0.9504
9 0.0508 0.9505
8 0.0006 0.9994
4 0.1190 0.8880
5 0.0013 0.9987
$sigma
[1] 0.3705
$null.dist
[1] "F(2,90)
covTest の結果は、上位ヒット機能の p 値を示しました。私の質問は、予測子 1 の係数やそのStd.err
およびなどのこれらの機能の係数を取得する方法95%CI
です。これらの見積もりを の対応する見積もりと比較したいと思いglm
ます。