rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] = 1]
rle30 = rle(temp_70year$TX >= 33) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 1]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 2]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 4]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 5]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 6]
rle30 = rle(temp_70year$TX >= 35) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 1]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 2]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 4]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 5]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 6]
rle30 = rle(temp_70year$TX >= 32) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 1]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 2]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 4]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 5]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 6]
temp_70year = read.csv('dailyMax_70y.csv', header = T)
rle30 = rle(temp_70year$TX >= 32) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 1]
rle30 = rle(temp_70year$T >= 32) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 1]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 2]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 4]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] == 5]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 1]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 2]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 4]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 5]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 6]
rle30 = rle(temp_70year$T >= 32) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 1 & rle30$lengths[rle30$values == T] <= 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 4 & rle30$lengths[rle30$values == T] <= 6]
rle30 = rle(temp_70year$T >= 30) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 1 & rle30$lengths[rle30$values == T] <= 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 4 & rle30$lengths[rle30$values == T] <= 6]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] > 6]
rle30 = rle(temp_70year$T >= 33) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 1 & rle30$lengths[rle30$values == T] <= 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 4 & rle30$lengths[rle30$values == T] <= 6]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] > 6]
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year = read.csv('rtd_70year.csv', header = T)
temp_70year= temp_70year[complete.cases(temp_70year), ]
plot(temp_70year$T[2600:3600], type = 'l')
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year= temp_70year[complete.cases(temp_70year), ]
plot(temp_70year$T[2600:3600], type = 'l')
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year= temp_70year[complete.cases(temp_70year), ]
plot(temp_70year$T[2600:3600], type = 'l')
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year= temp_70year[complete.cases(temp_70year), ]
plot(temp_70year$T[2600:3600], type = 'l')
d.duration = data.frame('duration' = NA, 'temp' = NA)[-1,]
t.min = min(temp_70year$T)
t.max = max(temp_70year$T)
t.seq = seq(t.min, t.max, 0.1)
n = length(t.seq)
mean = rep(NA, times = n)
median = rep(NA, times = n)
for (i in t.seq) {
tf = rle(temp_70year$T>= i)$values
count = rle(temp_70year$T>= i)$lengths
duration = count[tf == 'TRUE']
temp.duration = rep(i, times = length(duration))
neod.duration = data.frame(duration, temp.duration)
d.duration = rbind(d.duration, neod.duration)
}
for (j in 1:n) {
x = d.duration$duration[round(d.duration$temp.duration, 1) == t.seq[j]]
mean[j] = mean(x)
median[j] = median(x)
}
y = tapply(d.duration$duration, d.duration$temp.duration, mean, na.rm = T)
n = tapply(d.duration$duration, d.duration$temp.duration, length)
x = t.seq
sumry.drt = data.frame('temp' = t.seq, 'mean' = mean, 'median' = median)
sumry.drt$median[is.nan(sumry.drt$median)]<-NA
sumry.drt$median[is.nan(sumry.drt$median)]<-NA
sumry.drt = sumry.drt[complete.cases(sumry.drt), ]
rle30 = rle(temp_70year$T >= 33) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 1 & rle30$lengths[rle30$values == T] <= 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 4 & rle30$lengths[rle30$values == T] <= 6]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] > 6]
plot(sumry.drt$temp, sumry.drt$mean, log = 'y')
rle30 = rle(temp_70year$T >= 30) #
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 1 & rle30$lengths[rle30$values == T] <= 3]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] >= 4 & rle30$lengths[rle30$values == T] <= 6]
rle30$lengths[rle30$values == T][rle30$lengths[rle30$values == T] > 6]
plot(sumry.drt$temp, sumry.drt$mean, log = 'y')
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year= temp_70year[complete.cases(temp_70year), ]
rel_temp = rle(temp_70year$T >= 30)
rle_temp = rle(temp_70year$T >= 30)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle_T[rle_T >= 1 & rle_T <= 3]
rle_T[rle_T >= 4 & rle_T <= 6]
rle_T[rle_T > 6]
rle_T
View(temp_70year)
temp_70year = read.csv('dailyMax_70y.csv', header = T)
rle_temp = rle(temp_70year$T >= 30)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle_T[rle_T >= 1 & rle_T <= 3]
rle_T[rle_T >= 4 & rle_T <= 6]
rle_T[rle_T > 6]
sumry.drt = data.frame('temperature' = temp.drt, '1-3 d' = drt1_3, '4-6 d' = drt4_6, '> 6 d' = drt_m6)
sumry.drt = data.frame('temperature' = NA, '1-3 d' = NA, '4-6 d' = NA, '> 6 d' = NA)
View(sumry.drt)
sumry.drt = data.frame('temperature' = NA, '1-3 d' = NA, '4-6 d' = NA, '> 6 d' = NA)
sumry.drt = data.frame('temperature' = NA, '1-3d' = NA, '4-6d' = NA, '>6d' = NA)
View(sumry.drt)
sumry.drt = data.frame('temperature' = NA, '1_3d' = NA, '4_6d' = NA, 'm6d' = NA)
sumry.drt = data.frame('temperature' = NA, 'd1_3d' = NA, 'd4_6d' = NA, 'dm6d' = NA)
sumry.drt = data.frame('temperature' = NA, 'd1-3d' = NA, 'd4_6d' = NA, 'dm6d' = NA)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 1 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 6]
rlem6 = rle_T[rle_T > 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
rle_temp = rle(temp_70year$T >= 30)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 1 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 6]
rlem6 = rle_T[rle_T > 6]
rle_T
rle13
rle46
rlem6
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 1 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 6]
rlem6 = rle_T[rle_T > 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T == 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 6]
rlem6 = rle_T[rle_T > 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T == 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(temp_70year)
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 1 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 6]
rlem6 = rle_T[rle_T > 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 1 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 6]
rlem6 = rle_T[rle_T > 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_70year.csv', header = T)
View(schp_70y)
# loading data
schp_70y = read.csv('schp_70year.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
View(temp_70year)
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_70year.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
# save data
write.csv(temp_70year, 'schp_70year.csv', row.names = F, col.names = T)
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_70year.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_70year.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
View(schp_70y)
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_dailyMax70.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
# save data
write.csv(temp_70year, 'schp_70year_sorted.csv', row.names = F, col.names = T)
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_dailyMax70.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
# save data
write.csv(temp_70year, 'schp_70year_sorted.csv', row.names = F, col.names = T)
# summary heatwave duration
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_dailyMax70.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
# save data
write.csv(temp_70year, 'schp_70year_sorted.csv', row.names = F, col.names = T)
# summary heatwave duration
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$TX >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_dailyMax70.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
# save data
write.csv(temp_70year, 'schp_70year_sorted.csv', row.names = F, col.names = T)
# summary heatwave duration
sumry.drt = data.frame('temperature' = NA, 'd2_3' = NA, 'd4_5' = NA, 'dm6' = NA)[-1,]
for (i in 30:35) {
rle_temp = rle(temp_70year$TX >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle23 = rle_T[rle_T >= 2 & rle_T <= 3]
rle45 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d2_3 = length(rle23)
d4_5 = length(rle45)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d2_3, d4_5, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
write.csv(sumry.drt, 'sumry70y.csv', row.names = F, col.names = T)
temp_70year = read.csv('dailyMax_70y.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 28:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
View(sumry.drt)
library(tidyverse)
library(lubridate)
# loading data
schp_70y = read.csv('schp_dailyMax70.csv', header = T)
# subset and calibrate
temp_70year = schp_70y[, c(2, 12:16)]
temp_70year[,c(2:6)] = temp_70year[,c(2:6)]/10
# add date
temp_70year$date<-as.Date(as.character(temp_70year$YYYYMMDD),format="%Y%m%d")
temp_70year$year = format(temp_70year$date, format = '%Y')
temp_70year$vicennial = temp_70year$year
# save data
write.csv(temp_70year, 'schp_70year_sorted.csv', row.names = F, col.names = T)
# summary heatwave duration
sumry.drt = data.frame('temperature' = NA, 'd2_3' = NA, 'd4_5' = NA, 'dm6' = NA)[-1,]
for (i in 28:35) {
rle_temp = rle(temp_70year$TX >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle23 = rle_T[rle_T >= 2 & rle_T <= 3]
rle45 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d2_3 = length(rle23)
d4_5 = length(rle45)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d2_3, d4_5, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
for (i in 28:38) {
rle_temp = rle(temp_70year$TX >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle23 = rle_T[rle_T >= 2 & rle_T <= 3]
rle45 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d2_3 = length(rle23)
d4_5 = length(rle45)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d2_3, d4_5, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
temp_70year = read.c
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year = read.csv('schp_dailyMax70.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 28:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
View(temp_70year)
temp_70year = read.csv('dailyMax_70y.csv', header = T)
temp_70year = read.csv('schp_dailyMax70.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 28:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d1_3 = length(rle13)
d4_6 = length(rle46)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
View(temp_70year)
temp_70year = read.csv('schp_dailyMax70.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 28:35) {
rle_temp = rle(temp_70year$T >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d2_3 = length(rle23)
d4_5 = length(rle45)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
temp_70year = read.csv('schp_dailyMax70.csv', header = T)
View(temp_70year)
temp_70year = read.csv('schp_dailyMax70.csv', header = T)
sumry.drt = data.frame('temperature' = NA, 'd1_3' = NA, 'd4_6' = NA, 'dm6' = NA)[-1,]
for (i in 28:35) {
rle_temp = rle(temp_70year$TX >= i)
rle_T = rle_temp$lengths[rle_temp$values == T]
rle13 = rle_T[rle_T >= 2 & rle_T <= 3]
rle46 = rle_T[rle_T >= 4 & rle_T <= 5]
rlem6 = rle_T[rle_T >= 6]
temperature = i
d2_3 = length(rle23)
d4_5 = length(rle45)
dm6 = length(rlem6)
sumry.neo = data.frame(temperature, d1_3, d4_6, dm6)
sumry.drt = rbind(sumry.drt, sumry.neo)
}
View(sumry.drt)
