add tsv export line

script/03_sequencing_effort.R

@@ -78,8 +78,9 @@ write_tsv(samples, "finalresult/sequencing_effort/sample_seq_astan.tsv", na = ""
 sampling_effort <- read_tsv("finalresult/sampling_effort/sampling_events.tsv", col_select = c("sample_id", "dna_extracted")) |> 
   rename("sampling_event" = "sample_id")
 
-samples_wide <- samples |> 
+samples_wide <- samples |>
   select(marker, sampling_event) |> 
+  filter(!is.na(sampling_event)) |> 
   group_by(marker, sampling_event) |> 
   summarise(n = n()) |> 
   pivot_wider(names_from = "marker", values_from = "n")
@@ -94,3 +95,5 @@ samples_wide <- full_join(sampling_effort, samples_wide, by = "sampling_event")
     )
   ) |> 
   relocate(date,size_fraction,.after = sampling_event)
+
+write_tsv(samples_wide, "finalresult/sequencing_effort/sample_seq_astan_wide.tsv", na = "")