ploting a oneBaseline figure

[Claudio Quezada R.]

Dear Nadège

If you want to plot some output that doesn't have alpha, you have to explicitly indicate it with plotAlpha = FALSE within credibilityIntervals(). Other way, this function is expecting to plot a 2 panels figure, one for trophic position, and another for alpha.

Cheers

Claudio

2018-06-27 6:08 GMT-04:00 <nadege....@hotmail.fr>:

In fact I also have this issue when using credibilityIntervals(), as alpha is not computed with only one baseline, and replaced by "NA".

TP_models$df looks like:

  model   group      consumer    lower    upper   median     mode alpha.lower alpha.upper alpha.median  alpha.mode combined
1    1b      22          TMAR 2.511890 4.098319 3.269893 3.266149          NA          NA           NA         NA     --1b   

> credibilityIntervals(TP_models$df, x = "group", xlab ="Community")
Error: Discrete value supplied to continuous scale

Le mercredi 27 juin 2018 12:03:41 UTC+2, nadege....@hotmail.fr a écrit :

Hi,

I have an issue with the calculation of the parametric version of TP with multiple species.

I only have one baseline, so that here is my code:

TP_List <- extractIsotopeData(TP_model, b1 = "Baseline1", b2 = NULL,
                             baselineColumn = "FG", consumersColumn = "Species",
                             groupsColumn = "Pond",
                             d13C = "D13C", d15N = "D15N")

This works nicely with the calculation of Bayesian TP using multiSpeciesTP.

Now I would like to compare these results with parametric estimates using parametricTP(). But because I only have one baseline, the function actually crashes on the first population and stops:

[1] "Parametric version of trophic position"
[1] "For consumer:  TMAR"
[1] "One baseline TP:  3.34"
[1] "Two baselines TP:  NaN NaN"
[1] "Full model TP. At the beginning:  NaN NaN"
Error in while (TP_p1 != TP_p2) { : missing value where TRUE/FALSE needed

Is there a way to make the function continue on other populations even with one baseline only?

Thank you for your help!

Nadège

--
To download tRophicPosition package visit https://github.com/clquezada/tRophicPosition
---
You received this message because you are subscribed to the Google Groups "tRophicPosition" group.
To unsubscribe from this group and stop receiving emails from it, send an email to trophicposition-support+unsubscr...@googlegroups.com.
To post to this group, send email to trophicposition-support@googlegroups.com.
Visit this group at https://groups.google.com/group/trophicposition-support.

To view this discussion on the web, visit https://groups.google.com/d/msgid/trophicposition-support/877b3e6c-5c95-4ea5-9362-21d68dccdd0d%40googlegroups.com.

For more options, visit https://groups.google.com/d/optout.

[nadege....@gmail.com]

That worked perfectly! Indeed I missed that argument in the help file, oops.

Thank you

Nadège

To unsubscribe from this group and stop receiving emails from it, send an email to trophicposition-support+unsub...@googlegroups.com.
To post to this group, send email to trophicposi...@googlegroups.com.

Visit this group at https://groups.google.com/group/trophicposition-support.