setwd("C:/Temp/Workstation")

#ALA

ALA = read.table("ALA.txt", header = T)
ALA05 = ALA[,3]
ALA10 = ALA[,4]
ALA20 = ALA[,5]
ALA35 = ALA[,6]
ALA50 = ALA[,7]

cpgram(ALA05, main = "ALA, teplota v 5 cm pod zemí")
cpgram(ALA10, main = "ALA, teplota v 10 cm pod zemí")
cpgram(ALA20, main = "ALA, teplota v 20 cm pod zemí")
cpgram(ALA35, main = "ALA, teplota v 35 cm pod zemí")
cpgram(ALA50, main = "ALA, teplota v 50 cm pod zemí")

spec.pgram(ALA05, spans = NULL, main = "ALA, teplota v 5 cm pod zemí")
spec.pgram(ALA10, spans = NULL, main = "ALA, teplota v 10 cm pod zemí")
spec.pgram(ALA20, spans = NULL, main = "ALA, teplota v 20 cm pod zemí")
spec.pgram(ALA35, spans = NULL, main = "ALA, teplota v 35 cm pod zemí")
spec.pgram(ALA50, spans = NULL, main = "ALA, teplota v 50 cm pod zemí")


#ALA1

ALA1 = read.table("ALA1.txt", header = T)
ALA105 = ALA1[,3]
ALA110 = ALA1[,4]
ALA120 = ALA1[,5]
ALA135 = ALA1[,6]
ALA150 = ALA1[,7]

cpgram(ALA105, main = "ALA1, teplota v 5 cm pod zemí")
cpgram(ALA110, main = "ALA1, teplota v 10 cm pod zemí")
cpgram(ALA120, main = "ALA1, teplota v 20 cm pod zemí")
cpgram(ALA135, main = "ALA1, teplota v 35 cm pod zemí")
cpgram(ALA150, main = "ALA1, teplota v 50 cm pod zemí")

spec.pgram(ALA105, spans = NULL, main = "ALA1, teplota v 5 cm pod zemí")
spec.pgram(ALA110, spans = NULL, main = "ALA1, teplota v 10 cm pod zemí")
spec.pgram(ALA120, spans = NULL, main = "ALA1, teplota v 20 cm pod zemí")
spec.pgram(ALA135, spans = NULL, main = "ALA1, teplota v 35 cm pod zemí")
spec.pgram(ALA150, spans = NULL, main = "ALA1, teplota v 50 cm pod zemí")


#Fiedler

Fiedler = read.table("Fiedler.txt", header = T)
Fiedler10 = Fiedler[,4]
Fiedler60 = Fiedler[,5]
Fiedler200 = Fiedler[,3]

cpgram(Fiedler10, main = "Fiedler, teplota v 10 cm pod zemí")
cpgram(Fiedler60, main = "Fiedler, teplota v 60 cm pod zemí")
cpgram(Fiedler200, main = "Fiedler, teplota vzduchu ve 2m")

spec.pgram(Fiedler10, spans = NULL, main = "Fiedler, teplota v 10 cm pod zemí")
spec.pgram(Fiedler60, spans = NULL, main = "Fiedler, teplota v 60 cm pod zemí")
spec.pgram(Fiedler200, spans = NULL, main = "Fiedler, teplota vzduchu ve 2m")


#Srazky_manual

srazky = read.table("srazky.txt", header = T)
rada = srazky[,1]
cpgram(rada, main = "Periodogram srážek za období od ledna 2001 do prosince 2007")
spec.pgram(rada, spans = NULL, main = "Specifický periodogram srážek za období od ledna 2001 do prosince 2007")