0
我有一个奇怪的情况,当我在控制台中运行它时,一大块代码正在为我工作,但是当我试图编织它时R Markdown文件,我收到错误。下面是代码:R中的错误:找不到对象,代码不工作R Markdown
{r,echo=FALSE, message=FALSE, fig.width=3, fig.height=3}
library(cdcfluview)
library(mosaic)
library(ggplot2)
library(dplyr)
neflu <- get_flu_data("census", "1", "ilinet", years=1997:2016)
neflu[neflu=="X"]<-NA
wk_new <- rep(NA, nrow(neflu))
season_new <- rep('0', nrow(neflu))
for(i in 1:nrow(neflu)){
if(neflu$WEEK[i] == 40){
wk_num <- 1
}
wk_new[i] <- wk_num
wk_num <- wk_num+1
if(neflu$WEEK[i] < 40) {
season_new[i] <- paste(neflu$YEAR[i]-1, '-', neflu$YEAR[i])
} else {
season_new[i] <- paste(neflu$YEAR[i], '-', neflu$YEAR[i]+1)
}
}
neflu_new <- cbind(neflu, wk_new, season_new)
neflu_agg <- neflu_new %>%
mutate(weighted_ili = as.numeric(X..WEIGHTED.ILI))%>%
group_by(REGION.TYPE, season_new) %>%
summarize(peak_ili = max(weighted_ili, na.rm=TRUE),
cum_ili = sum(weighted_ili,na.rm=TRUE),
peak_wk = wk_new[which.max(weighted_ili)])
neflu_agg2 <- subset(neflu_agg, peak_wk != 5, select=c(REGION.TYPE,season_new,peak_ili,cum_ili, peak_wk))
plot(peak_wk, peak_ili, main= "Peak ILI vs Peak Week", xlab= "Peak Week", ylab= "Peak ILI")
nelm1 <- lm(peak_ili~1+peak_wk)
abline(nelm1, col='red')
plot(peak_ili, cum_ili, main= "Cumulative ILI vs Peak ILI", xlab= "Peak ILI", ylab= "Cumulative ILI")
nelm2 <- lm(cum_ili~1+peak_ili)
abline(nelm2, col='blue')
,这里是错误消息我得到:
The object peak_wk not found Calls:<Anonymous> . . .
WithCallingHandlers -> withVisible -> eval -> eval -> plot execution halted".
我曾试图改变这个代码这一点,但我仍然得到错误:
{r,echo=FALSE, message=FALSE, fig.width=3, fig.height=3}
library(cdcfluview)
library(mosaic)
library(ggplot2)
library(dplyr)
neflu <- get_flu_data("census", "1", "ilinet", years=1997:2016)
neflu[neflu=="X"]<-NA
wk_new <- rep(NA, nrow(neflu))
season_new <- rep('0', nrow(neflu))
for(i in 1:nrow(neflu)){
if(neflu$WEEK[i] == 40){
wk_num <- 1
}
wk_new[i] <- wk_num
wk_num <- wk_num+1
if(neflu$WEEK[i] < 40) {
season_new[i] <- paste(neflu$YEAR[i]-1, '-', neflu$YEAR[i])
} else {
season_new[i] <- paste(neflu$YEAR[i], '-', neflu$YEAR[i]+1)
}
}
neflu_new <- cbind(neflu, wk_new, season_new)
weighted_ili <- mutate(neflu_new, weighted_ili = as.numeric(X..WEIGHTED.ILI))
peak_ili <- max(weighted_ili, na.rm=TRUE)
cum_ili <- sum(weighted_ili,na.rm=TRUE)
peak_wk <- wk_new[which.max(weighted_ili)]
neflu_agg <- neflu_new %>%
mutate(weighted_ili = as.numeric(X..WEIGHTED.ILI))%>%
group_by(REGION.TYPE, season_new) %>%
summarize(peak_ili = max(weighted_ili, na.rm=TRUE),
cum_ili = sum(weighted_ili,na.rm=TRUE),
peak_wk = wk_new[which.max(weighted_ili)])
neflu_agg2 <- subset(neflu_agg, peak_wk != 5, select=c(REGION.TYPE,season_new,peak_ili,cum_ili, peak_wk))
plot(peak_wk, peak_ili, main= "Peak ILI vs Peak Week", xlab= "Peak Week", ylab= "Peak ILI")
nelm1 <- lm(peak_ili~1+peak_wk)
abline(nelm1, col='red')
plot(peak_ili, cum_ili, main= "Cumulative ILI vs Peak ILI", xlab= "Peak ILI", ylab= "Cumulative ILI")
nelm2 <- lm(cum_ili~1+peak_ili)
abline(nelm2, col='blue')
这里是我现在得到的错误:
Error in FUN(X[[i]],...): only defined on data frame with all
numeric variables Calls:
任何帮助只是让这个代码编织将不胜感激。我相信我只是错过了一些非常明显的东西,但我无法弄清楚。
在您的第一个代码块中,您可能需要引用peak_wk作为neflu_agg $ peak_wk。我的理解是'dplyr'返回一个data.frame。您引用peak_wk的代码似乎正在提供一个名为neflu_agg的数据框。 – lmo
你的意思是我应该在代码的总结部分这样做?所以neflu_agg $ peak_wk = wk_new [which.max(weighted_ili)])? 如果是这样,我收到一条新的错误消息,说cum_ili = sum(weighted_ili,na.rm = TRUE)中有一个意外的'=', – Nathan
不,我指的是非'dplyr'代码,如'plot'功能。注意错误信息:'plot execution halted'。 – lmo