# NOT RUN {
df <- read_df_raw(df_wbw,0.02,7)
df_slice <- dplyr::filter(df, df$n_period == 1)
# Check that this has enough data! needs more than 10 rows to be meaningful
nrow(df_slice)
# needs median temperature
temp <- median_temp(df_wbw)
res <- p_iv.week(df_slice, temp = temp, N_c = 60)
# }
Run the code above in your browser using DataLab