library(lattice)
data <- read.table("cmi2005.txt", header=TRUE)

##Using FOR loop while to calculate avg temp
j <- 1
avg_temp <- array(0, dim=c(1,365))
sum_temp <-0
for(i in 1:max(data[,2])){
      while(data[j,2] == i) {
		sum_temp <- sum_temp + data[j,5]
		j <- j+1
	}
	avg_temp[1,i] <- sum_temp/24
	sum_temp <-0
}
plot(avg_temp[1:365])



##Using FOR loop while to calculate avg rh
k <- 1
avg_rh <- array(0, dim=c(1,365))
sum_rh <-0
for(i in 1:max(data[,2])){
      while(data[k,2] == i) {
		sum_rh <- sum_rh + data[k,6]
		k <- k+1
	}
	avg_rh[1,i] <- sum_rh/24
	sum_rh <-0
}
plot(avg_rh[1:365])

##Using FOR loop while to calculate avg SolarR
m <- 1
avg_solarR <- array(0, dim=c(1,365))
sum_solarR <-0
for(i in 1:max(data[,2])){
      while(data[m,2] == i) {
		sum_solarR <- sum_solarR + data[m,4]
		m <- m+1
	}
	avg_solarR[1,i] <- sum_solarR/24
	sum_rh <-0
}
plot(avg_solarR[1:365])


##Using FOR loop while to calculate avg wind
n <- 1
avg_wind <- array(0, dim=c(1,365))
sum_wind <-0
for(i in 1:max(data[,2])){
      while(data[n,2] == i) {
		sum_wind <- sum_wind + data[n,7]
		n <- n+1
	}
	avg_wind[1,i] <- sum_wind/24
	sum_wind <-0
}
plot(avg_wind[1:365])


##Using FOR loop while to calculate sum precip
o <- 1
avg_precip <- array(0, dim=c(1,365))
sum_precip <-0
for(i in 1:max(data[,2])){
      while(data[o,2] == i) {
		sum_precip <- sum_precip + data[o,8]
		o <- o+1
	}
	avg_precip[1,i] <- sum_precip/24
	sum_precip <-0
}

plot(avg_precip[1:365])

##Using FOR loop while to calculate sum precip and cumulative
cum_precip <- sum(avg_precip)


##Using FOR loop If to calculate frost days at 0C
Frost <- (data$DailyTemp.C)
fd <- numeric(length(Frost))
for(i in 1:length(Frost)){
ffd <- Frost[i]
if(ffd >= 0){
fd[i] <- 0
}else{
fd[i] <- ffd
}
}
plot(fd)

##Using FOR loop If to calculate frost days at -2C
Frost2 <- (data$DailyTemp.C)
fd2 <- numeric(length(Frost2))
for(i in 1:length(Frost2)){
ffd2 <- Frost2[i]
if(ffd2 >= -2){
fd2[i] <- 0
}else{
fd2[i] <- ffd2
}
}
plot(fd2)