phase space embeddings

tbep-tech · Sep 6, 2024 · 1acb87a · 1acb87a
1 parent ed543db
commit 1acb87a
Show file tree

Hide file tree

Showing 16 changed files with 16 additions and 6 deletions.
diff --git a/R/cumu_chl-max_TN_otb-subsegments.R b/R/cumu_chl-max_TN_otb-subsegments.R
@@ -35,7 +35,7 @@ epcwq3.sub <- epcwq3[ which( epcwq3$param=="Chla"
                              & epcwq3$subseg==subseg
                        ), ]
 epcwq3.sub$month <- floor_date( epcwq3.sub$date, unit = 'month' ) 
-chldat <- epcwq3.sub |> group_by(month) |> summarise( chl = max(value) ) |> as.data.frame()
+chldat <- epcwq3.sub |> group_by(month) |> dplyr::summarise( chl = max(value) ) |> as.data.frame()
 chldat$chl <- chldat$chl |> log10()
 wqdat <- inner_join( chldat, loaddat, by = 'month' )
 

diff --git a/R/cumu_chl-max_TN_otb.R b/R/cumu_chl-max_TN_otb.R
@@ -34,7 +34,7 @@ epcwq3.sub <- epcwq3[ which( epcwq3$param=="Chla"
                              # & epcwq3$subseg==subseg
                        ), ]
 epcwq3.sub$month <- floor_date( epcwq3.sub$date, unit = 'month' ) 
-chldat <- epcwq3.sub |> group_by(month) |> summarise( chl = max(value) ) |> as.data.frame()
+chldat <- epcwq3.sub |> group_by(month) |> dplyr::summarise( chl = max(value) ) |> as.data.frame()
 chldat$chl <- chldat$chl |> log10()
 wqdat <- inner_join( chldat, loaddat, by = 'month' )
 
@@ -75,7 +75,7 @@ plot( mean ~ TN, data = chl_TN, las = 1,
       ylim = c(0,2.5), yaxt = 'n',
       ylab = "Chlorophyll a (ug/L)", xlab = "", xlim = c(10,200)
 )
-mtext( "TN load (tons/month)", side = 1, line = 2, cex = 0.7 )
+mtext( "TN load (tons/month)", side = 1, line = 2, cex = 1 )
 axis( 2, at = 0:3, las = 1,
       labels = c( 10^0, 10^1, 10^2, 10^3 ) )
 # abline( h = axTicks(2), col = rgb(0,0,0,0.1) )

diff --git a/R/cumu_chl-mean_TN_otb-subsegments.R b/R/cumu_chl-mean_TN_otb-subsegments.R
@@ -35,7 +35,7 @@ epcwq3.sub <- epcwq3[ which( epcwq3$param=="Chla"
                              & epcwq3$subseg==subseg
                            ), ]
 epcwq3.sub$month <- floor_date( epcwq3.sub$date, unit = 'month' ) 
-chldat <- epcwq3.sub |> group_by(month) |> summarise( chl = mean(value) ) |> as.data.frame()
+chldat <- epcwq3.sub |> group_by(month) |> dplyr::summarise( chl = mean(value) ) |> as.data.frame()
 chldat$chl <- chldat$chl |> log10()
 wqdat <- inner_join( chldat, loaddat, by = 'month' )
 

diff --git a/R/cumu_chl-mean_TN_otb.R b/R/cumu_chl-mean_TN_otb.R
@@ -34,7 +34,7 @@ epcwq3.sub <- epcwq3[ which( epcwq3$param=="Chla"
                              # & epcwq3$subseg==subseg
                            ), ]
 epcwq3.sub$month <- floor_date( epcwq3.sub$date, unit = 'month' ) 
-chldat <- epcwq3.sub |> group_by(month) |> summarise( chl = mean(value) ) |> as.data.frame()
+chldat <- epcwq3.sub |> group_by(month) |> dplyr::summarise( chl = mean(value) ) |> as.data.frame()
 chldat$chl <- chldat$chl |> log10()
 wqdat <- inner_join( chldat, loaddat, by = 'month' )
 
@@ -75,7 +75,7 @@ plot( mean ~ TN, data = chl_TN, las = 1,
       ylim = c(0,2.5), yaxt = 'n',
       ylab = "Chlorophyll a (ug/L)", xlab = "", xlim = c(10,200)
 )
-mtext( "TN load (tons/month)", side = 1, line = 2, cex = 0.7 )
+mtext( "TN load (tons/month)", side = 1, line = 2, cex = 1 )
 axis( 2, at = 0:3, las = 1,
       labels = c( 10^0, 10^1, 10^2, 10^3 ) )
 # abline( h = axTicks(2), col = rgb(0,0,0,0.1) )

diff --git a/data/embed_pyro_otb.RData b/data/embed_pyro_otb.RData
diff --git a/data/embed_pyro_otb_ce.RData b/data/embed_pyro_otb_ce.RData
diff --git a/data/embed_pyro_otb_cw.RData b/data/embed_pyro_otb_cw.RData
diff --git a/data/embed_pyro_otb_ne.RData b/data/embed_pyro_otb_ne.RData
diff --git a/data/embed_pyro_otb_nw.RData b/data/embed_pyro_otb_nw.RData
diff --git a/data/embed_pyro_otb_se.RData b/data/embed_pyro_otb_se.RData
diff --git a/data/embed_pyro_otb_sw.RData b/data/embed_pyro_otb_sw.RData
diff --git a/data/pyro_routine.csv b/data/pyro_routine.csv
@@ -0,0 +1,10 @@
+name,lat,lon
+TB Pyrodinium Station 01,28.02,-82.675
+TB Pyrodinium Station 02,27.972,-82.673
+TB Pyrodinium Station 03,27.98,-82.62
+TB Pyrodinium Station 04,27.949,-82.561
+TB Pyrodinium Station 05,27.957,-82.622
+TB Pyrodinium Station 06,27.942,-82.675
+TB Pyrodinium Station 07,27.9055,-82.6046
+TB Pyrodinium Station 08,27.9244,-82.5496
+TB Pyrodinium Station 09,27.86111,-82.58911
diff --git a/data/ssa_pyro_otb_ne.RData b/data/ssa_pyro_otb_ne.RData
diff --git a/figs/cumu_chl-max_TN_otb.png b/figs/cumu_chl-max_TN_otb.png
diff --git a/figs/cumu_chl-mean_TN_otb.png b/figs/cumu_chl-mean_TN_otb.png
diff --git a/figs/ssa_pyro_plots.png b/figs/ssa_pyro_plots.png