我不确定正确的词是什么,所以很抱歉术语严重错误。基本上我有大约 1000 个数据点,我想对 100 个数据点随机子采样 999 次,并对每个子样本执行相同的函数(广义最小二乘模型),并查看相关性为正的频率。
您可以使用
sample
从数据点中随机抽取样本。然后你可以使用 replicate
运行你的函数 n 次。
随机抽取 n=100 样本并计算平均值 10 次的示例:
> set.seed(1)
> datapoints <- runif(1000, max = 10000)
> result <- replicate(10, mean(sample(datapoints, 100)))
5194.298 5063.320 5064.992 4681.281 5008.011 4849.998 5320.206 5012.931 4900.636 4776.135