Построение предсказанных значений GLM с категориальным и непрерывным предиктором - PullRequest
Новая
       5

Построение предсказанных значений GLM с категориальным и непрерывным предиктором

1 голос
/ 25 марта 2020

У меня есть обобщенная линейная модель, в которой один предиктор является категориальной переменной (двоичная 1/0), а другой - непрерывной переменной.

Вот мои данные и мой код: 

library(reprex)
library(tidyverse)
library(datapasta)
library(ggplot2)

# my data frame: rainfall_binary is a binary vector (factor) while temp is continuous (it's temperature)

df <- data.frame(
  VH_average = c(-18.4527033816948,
                 -16.2644305598873,-16.1847107297772,-16.1971205524948,
                 -16.5239874732068,-17.2211302093816,-18.2875256347705,
                 -13.7489056675713,-14.5000673290099,-15.4042266341501,-15.1817907735829,
                 -15.6276170790517,-18.3260089724533,-15.1530169022796,
                 -16.1034334250424,-15.8582055282567,-16.2250640523659,
                 -16.9422869158206,-17.5646332225128,-15.4220277527455,
                 -15.7061506787604,-16.0053241375835,-15.6587460003135,
                 -16.848757403998,-17.9766074787419,-15.2396630934534,
                 -15.7573442344174,-15.8994493522684,-15.5906833828475,
                 -16.8160028280318),
  rainfall_binary = c(1,1,1,1,1,1,1,1,1,1,1,
                      1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0),
  temp = c(4.05,4.05,4.05,4.05,4.05,
           4.05,7,7,7,7,7,7,4.9,4.9,4.9,4.9,4.9,4.9,5.7,
           5.7,5.7,5.7,5.7,5.7,2.25,2.25,2.25,2.25,2.25,2.25)
)


# build model

GLM.REPREX <- glm(data=df,VH_average~rainfall_binary+temp)

# create predictors


df_newtemp <- 
  data.frame(temp= seq(0,20, length=100))


df_predicts <- predict(GLM.REPREX, newdata=df_newtemp, int = 'c')
#> Error in eval(predvars, data, env): oggetto "rainfall_binary" non trovato

Создано в 2020-03-25 пакетом Представить (v0.3.0)

Он говорит, что не находит объект. Я пробовал другими способами, но я новичок в R, и я даже не уверен, правильно ли то, что я делаю. Любые предложения приветствуются! Спасибо

...