I don't know of a package that calculates "lunar day" from a date. In theory you could do it from your dataset, by identifying the maximum and minimum values for phase, then converting phases to percentages and rounding up as a proportion of 29.53.
However, the lunar
package also calculates illumination (as a fraction of visible surface). I think this is a good proxy for lunar day and also gives you a physical value, rather than something more arbitrary.
Using your data, it's clear that new moon occurs near the start of the month:
library(tidyverse)
library(lunar)
sample_data <- read_csv("sample_data.csv")
sample_data %>%
mutate(Date = as.Date(fdate, "%d/%m/%Y"),
illum = lunar.illumination.mean(Date)) %>%
ggplot(aes(Date, illum)) + geom_point()

We could also fill in the missing dates, which makes the lunar cycle apparent:
all_dates <- data.frame(Date = seq.Date(min(as.Date(sample_data$fdate, "%d/%m/%Y")),
max(as.Date(sample_data$fdate, "%d/%m/%Y")),
by = "1 day")) %>%
mutate(illum = lunar.illumination.mean(Date))
all_dates %>%
ggplot(aes(Date, illum)) + geom_point()

Now, assuming your dataset has a column named catch
, we could begin analysis by joining the catch data with the full date range, then plotting catch and lunar illumination. This dataset may also be used for regression, correlation etc.
# simulated catch data
set.seed(123)
sample_data <- sample_data %>%
mutate(catch = rnorm(16, 100, 30))
all_dates %>%
left_join(mutate(sample_data, Date = as.Date(fdate, "%d/%m/%Y"))) %>%
select(Date, illum, catch) %>%
gather(variable, value, -Date) %>%
ggplot(aes(Date, value)) +
geom_point() +
facet_grid(variable~., scales = "free_y")
