Вычисление lm () в цикле - PullRequest
       3

Вычисление lm () в цикле

1 голос
/ 09 февраля 2020

Цель: общая цель проблемы состоит в том, чтобы вычислить доверительный интервал (CI) для различных размеров выборки (n = 2,4..1024) rnorm, 10000 раз, а затем подсчитать, сколько раз каждый из них терпит неудачу ( это, вероятно, требует счетчика и оператора if / else). Наконец, результаты должны быть нанесены на график

Я пытаюсь вычислить CI средних для нескольких моделей размеров выборки, однако я сначала пытаюсь разбить код для одной конкретной c размера выборки a = 8.

У меня проблема в том, что я не знаю, как создать линейную модель для каждой строки. Кто-нибудь знает, как я могу это сделать? Вот что у меня пока так:

a <- 8
n.sim.3 <- 10000

for ( i in a) {
  r.mat <- matrix(rnorm(i*n.sim.3), nrow=n.sim.3, ncol = a)
  lm.tmp <- apply(three.mat,1,lm(n.sim.3~1)  # The lm command is where I'm stuck I don't think this is correct)
  confint.tmp <- confint(lm.tmp)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...