# NOT RUN {
# upper record times observed in a 100 length time series
L <- c(1, 4, 14, 40, 45, 90)
X <- series_record(L_upper = L, Trows = 100)
# now you can apply plots and tests for upper records to the X series
#N.plot(X)
#N_normal.test(X)
# if you also have lower record times
L_lower <- c(1, 2, 12, 56, 57, 78, 91)
X <- series_record(L_upper = L, L_lower = L_lower, Trows = 100)
# now you can apply plots and tests to the X series with both types of record times
#foster.plot(X, statistic = 'd')
#foster.test(X, statistic = 'd')
# apply to the 200-meter Olympic records from 1900 to 2020
or200m <- series_record(L_lower = Olympic_records_200m$t,
R_lower = Olympic_records_200m$value,
Trows = 27)
# some plots and tests
N.plot(or200m, record = c(0,1,0,0))
N.test(or200m, record = "lower", distribution = "poisson-binomial")
# }
Run the code above in your browser using DataLab