我正在编写一个简单的函数，使用比率估计器和数据集"；虹膜"；。这是我的代码：
library(boot)
rat=function(data,index)
{
X=data[index,1]
Y=data[index,2]
return(mean(Y)/mean(X)*Xbar)
}
Xbar=mean(iris$Sepal.Length)
N=nrow(iris)
n=0.7*N
samp<-sample(1:N,n,replace = F)
iris_samp<-iris[samp,]
output<-boot(iris_samp,rat,R=1000)
mean(output$t)
output$t0
我得到的结果是：
> mean(output$t)
[1] 3.057975
> output$t0
[1] 3.061799
我的印象是t0只是t的平均值，即每个样本的估计值。然而，这两个值在这里显然不匹配。我的编码是错误的还是我的理解是错误的？
提前感谢！
$t0是在原始数据上运行lm时的系数$t是自举的值。
这里是的详细信息