2017-07-29 20 views
-1

私自身の関数 "foo"を使ってARIMAモデルを推定しようとしました。エラーコードはRの関数プロットがどこにも見つかりませんでした。本当に私を悩ます。どのように起こったか考えてみませんか?私はstats:::plotまだ動作しませんしようとしました。エラー:関数プロットが見つからないR

```{r warning=FALSE, message=FALSE} 
require(forecast) 
foo <- function(N, y1=1, p=0.8, q=1-p, seed=NULL, stepwise=TRUE, verbose=0) { 
if (!is.null(seed)) set.seed(seed) 
y <- c(y1, rep(NA,N-1)) 
for (i in 2:N) y[i] <- p*y[i-1] + q*rnorm(1) 
if (verbose) plot(y, type='l', xlab='time') 
return (list(y=y, aa=auto.arima(y, stepwise=stepwise))) 
} 
v <- foo(N=100, seed=10, verbose=1); v$aa$coef 
``` 
+0

'plot'コールの直前に' search() 'を置き、' package:stats'があなたの呼び出しツリーにあることを確認してください – MichaelChirico

+0

@MichaelChiricoありがとう。もう一度試しました。 (プロット) – xxyyzz

+1

あなたは私が言ったことをしなかった:p – MichaelChirico

答えて

0
require(forecast) 
Loading required package: forecast 
> foo <- function(N, y1=1, p=0.8, q=1-p, seed=NULL, stepwise=TRUE, verbose=0) { 
+  if (!is.null(seed)) set.seed(seed) 
+  y <- c(y1, rep(NA,N-1)) 
+  for (i in 2:N) y[i] <- p*y[i-1] + q*rnorm(1) 
+  if (verbose) plot(y, type='l', xlab='time') 
+  return (list(y=y, aa=auto.arima(y, stepwise=stepwise))) 
+ } 
> v <- foo(N=100, seed=10, verbose=1); v$aa$coef 
    ar1 
0.8946974 

あなたのコードはうまく機能..それもプロットをプロットします。だから、おそらくprobleは、予測パッケージの読み込みです。または、3つのコロンではなく2つのプロットの前にstat :: :: plotを試してみてください。

関連する問題