##TABLE S1
  #point to watch: for ITRDB dendro series, the column start.rec is used: this is first year with 8 samples. (Doesn't apply to end)

#ITRDB matrix
	download.file("http://www.meteo.psu.edu/~mann/supplements/MultiproxyMeans07/data/itrdbmatrix","temp.dat")
	test=scan("temp.dat",n=-1) #2426050
	test=t( array(test,dim=c(1210,length(test)/1210) )) ;dim(test) #2005 1210
	
	proxy_info=test[1:3,]
	test=test[4:nrow(test),] # from - 5 to 1996
	range(test[,1])
	test=ts(test[,2:ncol(test)],start=test[1,1])
	proxy=test;rm(test)
	dim(proxy) #[1] 2002 1209
	tsp(proxy) #- 5 1996

	proxy_info=t(proxy_info)
	proxy_info=proxy_info[2:nrow(proxy_info),]
	names(proxy_info)=c("long","lat","code")

##rtable
	rtable=scan("http://www.meteo.psu.edu/~mann/supplements/MultiproxyMeans07/data/rtable1209")
	rtable=array(rtable,dim=c(1209,7))

##LOAD DETAILS
	download.file("http://www.climateaudit.org/data/mann.2008/details.tab","temp.dat",mode="wb");load("temp.dat")
   		#coral proxies (7000,7001) may need to be inverted 
	coral=(details$code==7000)|(details$code==7001);sum(coral) #15

 
###MAKE TABLE S1: NUMBER OF PROXIES IN NETWORKS
     dendro=(details$code==9000)|(details$code==7500);sum(dendro) #1032
     annual=! (substr(as.character(details$code),4,4)=="1");sum(annual) #1158	#this digit controls definition of decadal
     NH=(details$lat>0)	;sum(NH) #1036	
     passing=!is.na(details$r1850_1995)	;sum(passing) #484

     countf=function(temp) {
	 count=array(NA,dim=c(19,3));count=data.frame(count);names(count)=c("dendro","other","total")
	 row.names(count)=paste(seq(0,1800,100),seq(99,1899,100),sep="-")
	 start0=seq(0,1800,100)
	 for (i in 1:19) {
	   count$dendro[i]=sum(details$start.rec[temp&dendro]<=start0[i])
	   count$other[i]=sum(details$start.rec[temp&!dendro]<=start0[i])
	   count$total[i]=sum(details$start.rec[temp]<=start0[i])
	  }
	countf=count[nrow(count):1,]
	countf
	}

     #set up list. Named per Table S1 2 columns across, then by geog, then by screened or not
	name1= c(outer(c("annual","all"),c("NH","SH","GLB"), function (x,y) paste(x,y,sep="_")))
	name1=c(outer(name1,c("full","screen"), function (x,y) paste(x,y,sep="_")))
	name1 
	  # [1] "annual_NH_full"    "all_NH_full"       "annual_SH_full"    "all_SH_full"       "annual_GLB_full"   "all_GLB_full"      "annual_NH_screen" 
	  # [8] "all_NH_screen"     "annual_SH_screen"  "all_SH_screen"     "annual_GLB_screen" "all_GLB_screen"   
	count=rep(list(NA),12);
	names(count)=name1

	count[[1]]= countf(annual&NH)
	count[[2]]=countf(NH)
	count[[3]]=countf(annual&!NH)
	count[[4]]=countf(!NH)
	count[[5]]=countf(annual)
	count[[6]]=countf(rep(TRUE,1209))
	count[[7]]= countf(annual&NH&passing)
	count[[8]]=countf(NH&passing)
	count[[9]]=countf(annual&!NH&passing)
	count[[10]]=countf(!NH&passing)
	count[[11]]=countf(annual&passing)
	count[[12]]=countf(passing)

	count
	   #this is close reproduction of Table S1