timeserie<-runif(100)+1:100-(seq(1,10,length=100)^2)
model<-NULL
for(i in 1:100){
model<-cbind(model,runif(100)+1:100-(seq(1,10,length=100)^2.01))
}
pomic.timeseries(timeserie,model,eps=10^-20,fullmsd=TRUE,
plotting=TRUE,half_window_size=10,
check_whole=TRUE,check_diffs=TRUE)
Run the code above in your browser using DataLab