R - смоделировать результат модели - PullRequest
2 голосов
/ 24 апреля 2019

У меня есть модель glm в R, она вычисляет некоторую оценку.У меня есть исходные данные, коэффициенты, модель, но я не могу рассчитать результат самостоятельно.

Исходные данные

df <- read.table(text="
                 exe   velocity  days sequence term  amt
    1              1   1.666667   10         9  20   507.44 
", header=TRUE)

Модель выглядит следующим образом:

model <- glm(
    formula =
      exe ~
      velocity +
      days +
      sequence +
      term +
      amt
    data = df,
    family = binomial
  )

Коэффициенты этой модели:

(Intercept)                 velocity                     days                sequence                    term                     amt 
-0.8519663862            0.4389537214           -0.0011613631           -0.0433850503           -0.0050228834            0.0005687273

Если я не правильно понялЯ должен быть в состоянии повторить эту модель по этой

((-0.8519663862)+(0.4389537214*1.666667)+(-0.0011613631*10)+(-0.0433850503*9)+(-0.0050228834*20)+(0.0005687273*507.44))

Но результат моих расчетов равен -0,3343185, но правильный показатель, рассчитанный моделью, равен 0,4171902.Что я делаю не так?

...