私のカウントデータにPoisson一般化線形モデルを使用し、データと適合モデルをggplotで実行しています。ggplot2とgeom_ribbon:eval(expr、envir、enclos)のエラー:オブジェクト 'Freq'が見つかりません。
マイデータ:
mod <- glm(Freq~YR, data = sub9, family = "poisson")
pred.df <- data.frame(YR = seq(min(sub9$YR), max(sub9$YR), length.out = 100))
pred <- predict(mod, newdata = pred.df, se.fit = TRUE)
pred.df$count <- exp(pred$fit)
pred.df$countmin <- exp(pred$fit - 2 * pred$se.fit)
pred.df$countmax <- exp(pred$fit + 2 * pred$se.fit)
ggplot(sub9,aes(x=YR,y=Freq)) +
scale_y_continuous(limits=c(0,75),breaks=c(10,20,30,40,50,60,70),expand=c(0,0)) +
scale_x_continuous(limits=c(1960,2018),breaks=c(1960,1965,1970,
1975,1980,1985,1990,1995,
2000,2005,2010,2015)) +
geom_point() +
geom_ribbon(data = pred.df, aes(ymin = countmin, ymax = countmax), alpha=0.3) +
geom_line(data = pred.df) +
xlab(" ") + ylab("Count")
私は、最初の1つのデータセットでこの手順を実行して成功した後:
structure(list(YR = c(1960, 1961, 1962, 1963, 1964, 1965, 1966,
1967, 1968, 1969, 1970, 1971, 1972, 1973, 1974, 1975, 1976, 1977,
1978, 1979, 1980, 1981, 1982, 1983, 1984, 1985, 1986, 1987, 1988,
1989, 1990, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999,
2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010,
2011, 2012, 2013, 2014, 2015, 2016), Freq = c(0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L,
0L, 0L, 2L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 4L, 1L, 2L, 4L, 5L,
3L, 2L, 5L, 2L, 14L, 6L, 5L, 5L, 10L, 13L, 10L, 5L, 8L, 7L, 6L,
10L, 12L, 14L, 2L, 16L, 15L)), .Names = c("YR", "Freq"), row.names = 58:114, class = "data.frame")
ここでのモデル、パラメータやgeom_ribbonを使用してプロットはまず、スクリプトです新しい類似のデータセットで同じ手順を試すときにエラーメッセージが表示されます。エラーメッセージ:
Error in eval(expr, envir, enclos) : object 'Freq' not found
私が行っている唯一のものは、データフレームの名前を交換し、データフレームに基づいて、「のFreq」に「カウント」から従属変数の名前を変更することです。私は間違って何をしていますか?