Как загрузить результат арифметической операции с коэффициентами регрессии двух моделей - PullRequest
Новая
       87

Как загрузить результат арифметической операции с коэффициентами регрессии двух моделей

0 голосов
/ 08 марта 2019

Я хочу сделать начальную загрузку деления выходов двух регрессионных моделей, чтобы получить доверительные интервалы среднего значения результата операции.

# создание примера данных 

ldose <- rep(0:5, 2)
numdead <- c(1, 4, 9, 13, 18, 20, 0, 2, 6, 10, 12, 16)
sex <- factor(rep(c("M", "F"), c(6, 6)))
SF <- cbind(numdead, numalive = 20-numdead)
dat<-data.frame(ldose, numdead, sex, SF)
tibble::rowid_to_column(dat, "indices")

# создание функции для начальной загрузки 

out<-function(dat) {
     d<-data[indices, ] #allows boot to select sample
     fit1<- glm(SF ~ sex*ldose, family = binomial (link = log), start=c(-1,0,0,0))
     fit2<- glm(SF ~ sex*ldose, family = binomial (link = log), start=c(-1,0,0,0))
     coef1<-coef(fit1)
     numer<-exp(coef1[2])
     coef2<-coef(fit2)
     denom<-exp(coef2[2])
     resultX<-numer/denom
     return(mean(resultX))
}

#doing bootstrap 

results <- boot(dat, out, 1000)

# сообщение об ошибке

Ошибка в статистике (данные, оригинал, ...): неиспользованный аргумент (оригинал)

Заранее спасибо за любую помощь.

...