Displaying 1 result from an estimated 1 matches for "wsdist".
Did you mean:
sdist
2009 Feb 24
0
Wavelet Bootstrap Size Simulation
...ulations
wspval<-rep(0,nosims) #initialize vector of p-values
#outer loop for recording calculated p-values
for (j in 1:nosims){
series<-arima.sim(T, model=list(ar=arcoeff)) #series would be the observed
AR(1) series
test<-mos(series,1:length(series)) #observed value of test statistic
wsdist<-rep(0,bsreps) #initialize vector for wavelet bootstrap null
distribution
#innner loop for finding null distribution via wavelet bootstrap
for (i in 1:bsreps){
wsrep<-wavestrap(series, wf = "d4", infl=1)
wsdist[i]<-mos(wsrep,1:length(series))
}
wspval[j]<-sum(wsdi...