# -*- coding: utf-8 -*-
"""
Created on Wed May 18 11:23:36 2011
Test SPIRIT and Frahst on Contol Time Series: Cyclic
@author: - Musselle
"""
from ControlCharts import Tseries
from CUSUM import cusum
import numpy as np
import matplotlib.pyplot as plt
from Frahst_v3_1 import FRAHST_V3_1
from SPIRIT import SPIRIT
from utils import analysis, QRsolveA, pltSummary
#===============================================================================
# Runscript
#===============================================================================
#============
# Initialise
#============
'Create time series string'
series_1 = Tseries(0)
series_2 = Tseries(0)
series_3 = Tseries(0)
series_4 = Tseries(0)
series_5 = Tseries(0)
# S 1
# CYCLIC - Normal - Cyclic
series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, \
noise_type = 'gauss')
series_1.normalEt(100, base=6, noise=0.00000000000000001, noise_type='gauss')
# ============ # Initialise # ============ "Create time series string" series_1 = Tseries(0) series_2 = Tseries(0) series_3 = Tseries(0) series_4 = Tseries(0) series_5 = Tseries(0) # S 1 # CYCLIC - Normal - Cyclic series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, noise_type="gauss") series_1.normalEt(100, base=6, noise=0.00000000000000001, noise_type="gauss") series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, noise_type="gauss") # S 2 # Cyclic series_2.cyclicEt(300, base=6, noise=0.1, amp=2, period=5, phase=1.5, noise_type="gauss") # S 3 # Cyclic series_3.cyclicEt(300, base=6, noise=0.1, amp=2, period=5, phase=0.5, noise_type="gauss") # S 4 # Cyclic series_4.cyclicEt(300, base=1, noise=0.1, amp=2, period=5, phase=0.2, noise_type="gauss") # S 5 # Cyclic series_5.cyclicEt(300, base=5, noise=0.1, amp=2, period=5, phase=0.3, noise_type="gauss")
import matplotlib.pyplot as plt #=============================================================================== # Runscript #=============================================================================== #============ # Initialise #============ 'Create time series string' series_1 = Tseries(0) # NORMAL - normalEt(self, size, base=0, noise, noise_type = 'gauss') series_1.normalEt(50, 0, 1, 'gauss') series_1.normalEt(50, 3, 1, 'gauss') # CYCLIC - cyclicEt(self, size, base=0, noise=2, amp=10, period=25, \ # phase = 0, noise_type = 'gauss'): series_1.cyclicEt(100, 3, 1, 5, 50) # NORMAL - series_1.normalEt(50, 3, 1, 'gauss') series_1.normalEt(50, -3, 1, 'gauss') # UP - upEt(self,size,base=0,noise=1,gradient=0.2, noise_type = 'gauss') series_1.upEt(100, -3, 1, 0.2, 'gauss') # Normal series_1.normalEt(50, 17, 1, 'gauss')
# -*- coding: utf-8 -*-
"""
Created on Wed May 18 11:23:36 2011
Test SPIRIT and Frahst on Contol Time Series: Cyclic
@author: - Musselle
"""
from ControlCharts import Tseries
from CUSUM import cusum
import numpy as np
import matplotlib.pyplot as plt
from Frahst_v3_1 import FRAHST_V3_1
from SPIRIT import SPIRIT
from utils import analysis, QRsolveA, pltSummary
#===============================================================================
# Runscript
#===============================================================================
#============
# Initialise
#============
'Create time series string'
series_1 = Tseries(0)
series_2 = Tseries(0)
series_3 = Tseries(0)
series_4 = Tseries(0)
series_5 = Tseries(0)
# S 1
# CYCLIC - Normal - Cyclic
series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, \
noise_type = 'gauss')
series_1.normalEt(100, base=6, noise = 0.00000000000000001, noise_type ='gauss')
import matplotlib.pyplot as plt #=============================================================================== # Runscript #=============================================================================== #============ # Initialise #============ 'Create time series string' series_1 = Tseries(0) # NORMAL - normalEt(self, size, base=0, noise, noise_type = 'gauss') series_1.normalEt(50, 0, 1,'gauss') series_1.normalEt(50, 3, 1,'gauss') # CYCLIC - cyclicEt(self, size, base=0, noise=2, amp=10, period=25, \ # noise_type = 'gauss'): series_1.cyclicEt(100, 3, 1, 5, 50, 'gauss') # NORMAL - series_1.normalEt(50, 3, 1,'gauss') series_1.normalEt(50, -3, 1,'gauss') # UP - upEt(self,size,base=0,noise=1,gradient=0.2, noise_type = 'gauss') series_1.upEt(100, -3, 1, 0.2, 'gauss') # Normal series_1.normalEt(50, 17, 1,'gauss')