Table Tests

Table from Microsoft Word

1	2
3	4
5	6

Table from simple HTML

Column 1	Column2
Value of π	3.14159

Table from Excel

1	2
3	4
5	6

Station ID	Location	Latitude	Longitude
10160355000	SKIKDA	36.93	6.95
10160360000	ANNABA	36.83	7.82
10160400001	CAP CARBON	36.80	5.10
10160402000	BEJAIA	36.72	5.07
10160403000	GUELMA	36.47	7.47
10160419000	CONSTANTINE	36.28	6.62
10160445000	SETIF	36.18	5.42
10160468000	BATNA	35.55	6.18
10160475000	TEBESSA	35.48	8.13
15260714000	BIZERTE	37.25	9.80
15260725000	JENDOUBA	36.48	8.80
62316539000	CAPO FRASCA	39.75	8.47
62316560000	CAGLIARI/ELMA	39.25	9.07
62316564000	CAPO CARBONAR	39.10	9.52

Deconstructing the Steig AWS Reconstruction

(This is merely a test of a full post.. The full post can be found at

http://www.climateaudit.org/?p=5151 )

So how did Steig et al. calculate the AWS reconstruction? Since we don’t have either the satellite data or the exact path they followed using the RegEM Matlab files, we can’t say for sure how the results were obtained. However, using a little mathemagic, we can actually take the sequences apart and then we can also calculate reconstructed values for those times when the reconstructions have been replaced by the measured temperature values. I realize that there are problems with the AWS data sets, but for our purposes that is irrelevant. This analysis can easily be repeated for the updated results if and when Prof. Steig does the recalcs.

Test Post

This is for testing various techniques:

Upload of an equation from a pdf converted to a jpg

Seems to work reasonably well. Now a slightly smaller png file:

Not bad.

Try WordPress latex:

Now redo the equation above:

Yes!! It works.

Post an R script as a block quote:

download.file(“http://www.climateaudit.org/data/antarctic_mann/Data.tab”,”temp2.dat”,mode=”wb”);load(“temp2.dat”)

slopes = function(tsdat) {
star = min(time(tsdat))
slopes = rep(NA,12)
for (i in 1:12){
dat.win = window(tsdat,start=c(star,i),deltat=1)
yr = time(dat.win)
slopes[i] = ifelse(sum(!is.na(dat.win)) < 2, NA, coef(lm(dat.win~yr))[2])
}
slopes}

slope(Data$aws[,8])

all.slope = function(alldat) {
nc = ncol(alldat)
slop = matrix(NA, nrow = nc, ncol = 12)
colnames(slop) = month.abb
for (j in 1:nc) {
slop[j,]= slope(alldat[,j])}
slop}

test3 = all.slope(Data$aws[,6])

Quotes copy out wrong. This needs fixing.

New test:

Attach text document as .doc file:

collation_functions

The blog has begun!

Welcome to WordPress.com. Hopefully this blog will be useful as a vehicle for math and stat ideas and as a repository for graphs and documents.