# Issues with Unifrac calculations

**URL:** <https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183>\
**Category:** MetaPhlAn\
**Created:** [April 21, 2023, 7:13pm UTC](https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183 "2023-04-21T19:13:55Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![anitalijia](https://avatars.discourse-cdn.com/v4/letter/a/5f8ce5/32.png) [@anitalijia](https://forum.biobakery.org/u/anitalijia)\
**Post date:** [April 21, 2023, 7:13pm UTC](https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183/1 "2023-04-21T19:13:55Z")

</div>

I am trying to calculate unifrac beta diversity. I used the SGB tree in the Newick format, but I received the error messages. PS: I had no issues calculating other alpha and beta diversities.  
Thanks!

Error in FUN(X[[i]], …) : attempt to set index 0/0 in SET\_STRING\_ELT  
Calls: → lapply → FUN  
In addition: Warning message:  
In FUN(X[[i]], …) : NAs introduced by coercion  
Execution halted

---

<div class="post-metadata">

**Author:** ![PeterCx](https://avatars.discourse-cdn.com/v4/letter/p/f04885/32.png) [@PeterCx](https://forum.biobakery.org/u/PeterCx)\
**Post date:** [April 26, 2023, 10:24am UTC](https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183/2 "2023-04-26T10:24:05Z")

</div>

I had this problem myself a few days ago.

Are you using the most recent tree “mpa\_vOct22\_CHOCOPhlAnSGB\_202212.nwk”?

In my case I was using the tree from Jan21.

---

<div class="post-metadata">

**Author:** ![anitalijia](https://avatars.discourse-cdn.com/v4/letter/a/5f8ce5/32.png) [@anitalijia](https://forum.biobakery.org/u/anitalijia)\
**Post date:** [April 26, 2023, 12:59pm UTC](https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183/3 "2023-04-26T12:59:53Z")

</div>

I tried both nwk version, and got the same error messages. : (

---

<div class="post-metadata">

**Author:** ![aitor.blancomiguez](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/aitor.blancomiguez/32/86_2.png) [@aitor.blancomiguez](https://forum.biobakery.org/u/aitor.blancomiguez)\
**Post date:** [April 27, 2023, 8:18am UTC](https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183/4 "2023-04-27T08:18:18Z")

</div>

Hi @anitalijia which metaphlan database version are you using when running metaphlan?, you can check it in the first lines of your profiles

---

<div class="post-metadata">

**Author:** ![anitalijia](https://avatars.discourse-cdn.com/v4/letter/a/5f8ce5/32.png) [@anitalijia](https://forum.biobakery.org/u/anitalijia)\
**Post date:** [February 9, 2024, 3:32pm UTC](https://forum.biobakery.org/t/issues-with-unifrac-calculations/5183/5 "2024-02-09T15:32:00Z")

</div>

Hi @aitor.blancomiguez : I used mpa\_vJan21\_CHOCOPhlAnSGB\_202103.nwk, and my metaphlan used the same version. I am still experiencing issues with it.
