##REPLICATION LOEHLE 

	Data=array(NA,dim=c(1995,18)); proxy_sm=array(NA,dim=c(1995,18))
	for (i in 1:18) {
		test=proxy[[i]]
		test=test[!is.na(test$t),]
		f=approxfun(test[,1],test[,2])
		Data[,i]=f(1:1995)
		proxy_sm[,i]=filter(Data[,i],rep(1/30,30) )	
	}
	proxy_sm=scale(proxy_sm,scale=FALSE); dimnames(proxy_sm)[[2]]=id

	annual=ts(apply(proxy_sm,1,mean,na.rm=TRUE),start=1)
	count=ts(apply(!is.na(Data),1,sum,na.rm=TRUE),start=1)
		
	
##LOEHLE VERSION
	url<-"http://www.ncasi.org/programs/areas/climate/LoehleE&E2007.csv"
	test=read.csv(url,skip=2)
	recon=ts(test[,2],start=test[1,1])  #17 to 1980
	mean(recon[(1902:1980)-16]) # 0.01586732

##PLOT
	par(mar=c(3,3,2,1))
	ts.plot(annual)
	lines(recon,col="red")



##12 CRU NH #potential target 1961-1990
	url<-"http://www.cru.uea.ac.uk/cru/data/temperature/hadcrut3nh.txt"  #1850 2007 
	h<-scan(url)
	start0<-h[1]
	N<-length(h)
	h<-array(h,dim=c(27,N/27))
	h<-t(h)
	CRU<-ts(h[,14],start=start0,end=2006) # 1850 2006 
	delta=mean(CRU[(1902:1980)-1849]);delta #  -0.1592278