달력

10

« 2024/10 »

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
2019. 1. 14. 13:33

AR(2) stationary condition 공부/통계학2019. 1. 14. 13:33

phi1 <- seq(-3,3,0.1)

phi2<-seq(-3,3,0.1)

phi <- merge(phi1,phi2)


flag<-rep(0,3721)


for(i in 1:3721){

z=polyroot(c(1, -phi$x[i],-phi$y[i]))

az=abs(z)

flag[i]=ifelse((az[1]>1 & az[2]>1),1,0)


c<-cbind(phi,flag)

plot(c$x,c$y,col=c$flag)


:
Posted by 시스메