#############################################

Prob1 <- read.table (file = "C:\\Users\\Owner\\Documents\\Agron 590\\R\\cmi2005.txt", header = TRUE) # reading data for a .txt file saved in my pc. I created a data.frame with this data


####################### Problem 1

## 1


m1 <- Prob1[1:4400,]
m2 <- Prob1[4401:8757,]


LFrost <- NULL

for(i in nrow(m2):1){
if (m2$DailyTemp.C[i]<0) {LFrost<-m2[i,]}
};LFrost

FFrost <- NULL

for(i in 1:nrow(m1)){
if (m1$DailyTemp.C[i]<0) {FFrost<-m1[i,]}
};FFrost




####################### Problem 2 

## 1

MeanTemp <- numeric (length (365))

for (i in 1:365 ) {
	MeanTemp [i] <- mean (Prob1 [Prob1$doy == i, "DailyTemp.C"])
	}
print (MeanTemp)

MeanSolar <- numeric (length (365))

for (i in 1:365 ) {
	MeanSolar [i] <- mean (Prob1 [Prob1$doy == i, "solarR"])
	}
print (MeanSolar)

MeanRH <- numeric (length (365))

for (i in 1:365 ) {
	MeanRH [i] <- mean (Prob1 [Prob1$doy == i, "RH"])
	}
print (MeanRH)

MeanWS <- numeric (length (365))

for (i in 1:365 ) {
	MeanWS [i] <- mean (Prob1 [Prob1$doy == i, "WindSpeed"])
	}
print (MeanWS)

MeanP <- numeric (length (365))

for (i in 1:365 ) {
	MeanP [i] <- mean (Prob1 [Prob1$doy == i, "precip"])
	}
print (MeanP)

MeanMV <- data.frame(MeanTemp, MeanRH, MeanSolar, MeanWS, MeanP)
MeanMV


## 2

SumTemp <- numeric (length (365))

for (i in 1:365 ) {
	SumTemp [i] <- sum (Prob1 [Prob1$doy == i, "DailyTemp.C"])
	}
print (SumTemp)

SumSolar <- numeric (length (365))

for (i in 1:365 ) {
	SumSolar [i] <- sum (Prob1 [Prob1$doy == i, "solarR"])
	}
print (SumSolar)

SumRH <- numeric (length (365))

for (i in 1:365 ) {
	SumRH [i] <- sum (Prob1 [Prob1$doy == i, "RH"])
	}
print (SumRH)

SumWS <- numeric (length (365))

for (i in 1:365 ) {
	SumWS [i] <- sum (Prob1 [Prob1$doy == i, "WindSpeed"])
	}
print (SumWS)

SumP <- numeric (length (365))

for (i in 1:365 ) {
	SumP [i] <- sum (Prob1 [Prob1$doy == i, "precip"])
	}
print (SumP)

SumMV <- data.frame(SumTemp, SumRH, SumSolar, SumWS, SumP)
SumMV

####################### Problem 3

CSTemp <- sum (Prob1 [1:365,5])
CSTemp

CSSolar <- sum (Prob1 [1:365,4])

CSRh <- sum (Prob1 [1:365,6])

CSWs <- sum (Prob1 [1:365,7])

CSP <- sum (Prob1 [1:365,8])

CumulativeSum <- data.frame (CSTemp, CSSolar, CSRh, CSWs, CSP)
CumulativeSum