# $Id: tempchallenge.R,v 1.2 2007/09/24 12:03:11 kaip Exp $
#
#
# Copyright (c) 2007 Kai Puolamaki <Kai.Puolamaki@iki.fi>
#
# Permission to use, copy, modify, and distribute this software for any
# purpose with or without fee is hereby granted, provided that the above
# copyright notice and this permission notice appear in all copies.
#
# THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
# WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
# MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
# ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
# WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
# ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
# OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
#


DATA <- read.table("ftp://ftp.engr.udayton.edu/jkissock/gsod/FIHELSIN.txt",
                   header=FALSE,na.strings="-99")
colnames(DATA) <- c("Month","Day","Year","Temp")

DATA <- DATA[!is.na(DATA[,"Temp"]),]

mintemp <- min(DATA[,"Temp"])
maxtemp <- max(DATA[,"Temp"])

XDATA <- data.frame(list(X=2*((DATA[,"Month"]-1)/12+(DATA[,"Day"]-1)/(12*31))-1,Y=(DATA[,"Temp"]-mintemp)/(maxtemp-mintemp)))
XDATA <- XDATA[order(XDATA[,"X"]),]
len <- dim(XDATA)[1]
rownames(XDATA) <- 1:len
XDATA[,"X"] <- jitter(XDATA[,"X"])
XXDATA <- XDATA
XXDATA[sample(len)[21:len],"Y"] <- NA
write.table(XXDATA,file="t613050-2-2007.txt")
write.table(XDATA,file="t613050-2-2007-orig.txt")