Does anyone know of R code to implement the Cayula-Cornillon SST front detection algorithm?

More Matthew A. Birk's questions See All

How to convert between buffering capacity quantities - physiologists and oceanographers?

Physiologists measure buffering capacity of a cell or fluid in terms of mM / pH (or mEq / pH). However, oceanographers measure buffering capacity of seawater in mmol / kg seawater (total...

09 October 2018 4,818 0 View

Does pH adjustment technique matter for m-cresol purple solution used in spectrophotometric pH measurement?

08 September 2014 6,699 4 View

How to normalise data on OPM package from R Studio ?

Good afternoon, I recently used OmniLog from BIOLOG for my experimentations : I tested the metabolism of different strains on 2 types of plates. I have 16 strains of 3 different groups...

02 March 2021 3,584 1 View

How to plot the probability paper of Gumbel distribution using R programing?

In R-studio, there are many commands of Gumbel package. Arguments are also different. I`m asking about the alpha parameter of the Copula which must be greater than 1. If this is the one used to...

25 February 2021 5,229 3 View

Bayesian decision trees and bayesian regression trees in Python and/or R?

I am looking for implementation of bayesian decision trees and bayesian regression trees in Python and/or R.. Do you know if the existing libraries support this

23 February 2021 6,689 3 View

Are there packages (R or Python) for simple models (bucket?) to assess soil water content?

Hi everyone. I am researching droguht effects on soil water drought. I am looking for already implemented functions in R or Python to estimate soil water content from precipitation,...

21 February 2021 2,084 1 View

ANCOVA controlling baseline in RCT - how to perform power analysis and sample size estimation?

Hi, I have a simple question. I am hoping to perform a power analysis/sample size estimation for an RCT. We will be controlling for baseline symptoms, and using post-treatment or change scores as...

19 February 2021 2,648 23 View

MVN, how to conciliate contradictory results from Multivariate normality and Multivariate outliers test?

Dear colleagues, I'm working with Path analyses in lavaan and MVN packages. There are some results that for me are confusing and features from the MVN that I do not know how to set. My dataset is...

18 February 2021 190 3 View

Would you please recommend a function in R software that could be used with FUNGuild output to merge different rows that have similar guilds?

Dear colleagues, I have an output of FUNGuild tool (ASVs abundance in different samples and assigned to guilds) and would like to merge the ASVs that had similar ecological guilds using R...

16 February 2021 4,292 7 View

How to remove the species data from rda data frame in R (only remain the sites/samples point)?

Dear everyone, When I using the rda function in r, I met this question. As you can see in the figure, there are a lot of deep blue points which refers to the species. But I would like remain...

15 February 2021 8,982 2 View

How to find out or construct(?) the right priors from data of previous experiment to use in a bayesian glm?

I have the following equation, Eating choice (Yes/No)~ Box type (A/B/C). So the question is whether a certain box type affects eating choice in groups of dogs. The eating choice of each member in...

10 February 2021 6,024 9 View

XRPD Quantification - Building a calibration curve of drug excipient?

I am aiming to utilise XRPD to quantify a drug excipient within a drug product and have been having a little trouble building the calibration curve to allow for this quantification up until this...

10 February 2021 6,145 3 View

Ivan Allaman

I can help only.

Matthew A. Birk

Ivan, how would you like to help? Do you already have access to code? Or are you interested in writing your own code?

If you are familiarity with programming in R, WE can write code and to publish it. What do you think?

At the moment, I don't have the time to commit to that project, unfortunately.

Jon Lopez

Hi there,

I would also be interested on having a R code to estimate frontal detection algorithms like the Cayula-Cornillon algorithm. If anyone have access to that, I would really apreciate it.

I know that the Marine Geospatial Ecology tools for arcGIS do that. Unfortunately, I have no access to that software.

Thanks in advance,

Jon

I just found this package in R. Its called Cyclones and is clim.pact and akima packages dependent. I didn't try yet but looks interesting.

http://finzi.psych.upenn.edu/library/cyclones/html/oceancurrents.html

Hope it helps,

Thanks Jon, I wasn't too impressed by the YouTube videos showing its performance, but it may be an option. the MGET you referred to before may be an option. It uses Python on the backend. I'm looking into running the Python code in R. I'll let you know if I find a workable solution.

Wencheng Lau-Medrano

Hello everyone. Well, some months ago, I have started to work in this topic and I have developed an (official) R package, called 'grec'. Although it is a CRAN package, I am still working on it so the method that it uses is not precisely the SIED algorithm, but produces a good result.

I suggest you to download the developer version (on Github)

devtools::install_github("LuisLauM/grec")

Then, run the demo:

demo(topic = "detectFronts", package = "grec")

Thank you so much!!

Should try that.

I have also recently found the boaR package, which uses the Belkin's algorithm for SST and CHL fronts detection.

https://rdrr.io/github/galuardi/boaR/man/boaR-package.html

Yeah. That's a very useful package, but it performs their loops in R and grec uses the package 'imagine', which uses Rcpp (C++) in order to speed up the application of filters.

Nevertheless, I have done some changes and now 'grec' already has the Contextual Median Filter (contextualMF function), but only on the Github version, because I am checking and correcting some bugs with NA (the typical problem when you work with C++).