从t.test获得1000置信区间

Question

(case (car '(c d))
  ((a e i o u) 'vowel)
  ((w y) 'semivowel)
  (else 'consonant)) ; ==> consonant

我想运行1000次并创建一个获得1000个间隔的矩阵。我尝试了一些循环，但每次我得到相同的1000个间隔。但是，每次它应该给我一个不同的间隔，因为我在执行t.test之前每次采样。

Answer 1

您可以使用replicate：

执行此操作

x <- c(1:100)
y <- c(89:300)

myCI = function(x,y) {
  s1 <- sample(x, 30)
  s2 <- sample(y, 30)
  mytest <- t.test(s1, s2)
  mytest$conf.int
}

CIs = t(replicate(1000, myCI(x,y)))