начальная загрузка для пропорции. анализ мощности - PullRequest
Новая
       23

начальная загрузка для пропорции. анализ мощности

0 голосов
/ 01 апреля 2020

Я адаптировал следующий код из здесь . Я пытаюсь смоделировать мощность для теста пропорции. Я продолжаю получать следующее сообщение об ошибке: «Ошибка в prop.test (Treat, Ctrl):« x »и« n »должны иметь одинаковую длину». Я не совсем понимаю, почему, так как я знаю, что мои 'x' и 'n' имеют одинаковую длину. Могу ли я попросить чью-то помощь? Спасибо! 

library('paramtest')
library('pwr')

df <- data.frame(trt = rep(c("smokers","patients"), each = 4),
                 values= c(83, 90, 129, 70, 86, 93, 136, 82))

p_func_boot <- function(data, indices) {
  sample_data <- data[indices, ]
  treat <- sample_data[sample_data$trt == 'smokers', 'values']
  ctrl <- sample_data[sample_data$trt == 'patients', 'values']

  ptest <- prop.test(treat,ctrl)
  stat <- ptest$statistic
  p <- ptest$p.value

  return(c(X_squared =stat, p=p, sig=(p < .05)))
}

power_proptest_boot <- run_test(p_func_boot, n.iter= 100, output='data.frame', boot=TRUE,
                             bootParams=list(data=df))
...