我得到t的所有值都等于0.999999999999995 \ n无法计算置信区间
我的数据是:
$data
y x2 x3
1 2 1 1
2 0 2 1
3 -2 3 1
4 -4 4 1
5 -6 5 1
6 5 1 2
7 3 2 2
8 1 3 2
9 -1 4 2
10 -3 5 2
11 8 1 3
12 6 2 3
13 4 3 3
14 2 4 3
15 0 5 3
16 11 1 4
17 9 2 4
18 7 3 4
19 5 4 4
20 3 5 4
21 14 1 5
22 12 2 5
23 10 3 5
24 8 4 5
25 6 5 5
我这样做了:
b<-boot(data=d,statistic=lmcoefs,R=5)
> b
ORDINARY NONPARAMETRIC BOOTSTRAP
Call:
boot(data = d, statistic = lmcoefs, R = 5)
Bootstrap Statistics :
original bias std. error
t1* 1 -9.436896e-15 4.831384e-15
t2* -2 1.110223e-15 1.169693e-15
t3* 3 1.776357e-15 8.881784e-16
然后我用boot.ci计算了置信区间:给出了
c<-boot.ci(b,conf=0.95,type="all")
[1] "All values of t are equal to 0.999999999999995 \n Cannot calculate confidence intervals"
我得到这个错误我无法摆脱。我也不明白为什么。
答案 0 :(得分:0)
问题是我没有做足够的bootstrap复制,因此我的数据没有足够的变化,这就是我遇到这个错误的原因。
谢谢@joran!