#TO EXTEND SERIES BY PERSISTENCE
extend.persist<-function(tree) {
	extend.persist<-tree
	for (j in 1:ncol(tree) )  {
		test<-is.na(tree[,j])
		end1<-max ( c(1:nrow(tree)) [!test])
		test2<-( c(1:nrow(tree))>end1) & test
		extend.persist[test2,j]<-tree[end1,j]
		}
	extend.persist
	}	

if(method =="online") {
	proxy=read.table("http://www.climateaudit.org/data/mbh99/proxy.txt",sep="\t",header=TRUE)
	proxy=ts(proxy[,2:ncol(proxy)],start=1000) #1000 1980
	proxy=extend.persist(proxy)
	m0=apply(proxy[903:981,],2,mean);sd0=apply(proxy[903:981,],2,sd)
	if(scale.method==FALSE) sd0=FALSE
	proxy=ts(scale(proxy,center=m0,scale=sd0),start=1000) } else {

	load("d:/climate/data/mbh99/proxy.tab");dim(proxy)
	proxy=ts(proxy[(1000:1980)-tsp(proxy)[1]+1,],start=1000) #1000 1980
	proxy=extend.persist(proxy)
	m0=apply(proxy[903:981,],2,mean);sd0=apply(proxy[903:981,],2,sd)
	if(scale.method==FALSE) sd0=FALSE
	proxy=ts(scale(proxy,center=m0,scale=sd0),start=1000)
    }