Вязание документа приводит к ошибке в seq.int - PullRequest
0 голосов
/ 02 ноября 2019

При вязании моего документа следующий код останавливает вязание и приводит к ошибке. Следует отметить, что при нормальной работе этот код работает путем отображения необходимого графика. Эта проблема возникает только при вязании документа.

Я попытался вставить echo = false, надеясь просто обойти проблему, но в результате получило то же сообщение об ошибке.

Код:

qplot(wt, mpg, data=mtcars, colour=factor(am)) +
  geom_smooth(aes(ymin = lcl, ymax = ucl), data=grid, stat="identity") +
  ggtitle("MPG, Weight, & Transmission Relationship")

Ошибка:

Quitting from lines 131-134 (MotorTrend.Rmd) 
Error in seq.int(U[1L], U[2L], length.out = nx + 1) : 
  'length.out' must be a non-negative number
Calls: <Anonymous> ... lapply -> FUN -> <Anonymous> -> f -> <Anonymous> -> f

Execution halted

1 Ответ

0 голосов
/ 02 ноября 2019

Мне нужно было лучше определить мои термины, поэтому я добавил дополнительный код:

model <- lm(mpg ~ wt + factor(am) + wt*factor(am), data=mtcars)
grid <- with(mtcars, expand.grid(
  wt = seq(min(wt), max(wt), length = 20),
  am = levels(factor(am))))
grid$mpg <- stats::predict(model, newdata=grid)

err <- stats::predict(model, newdata=grid, se = TRUE)
grid$ucl <- err$fit + 1.96 * err$se.fit
grid$lcl <- err$fit - 1.96 * err$se.fit
qplot(wt, mpg, data=mtcars, colour=factor(am)) +
  geom_smooth(aes(ymin = lcl, ymax = ucl), data=grid, stat="identity") +
  ggtitle("MPG, Weight, & Transmission Relationship")
...