# Lefse without any subclass- still valid?

**URL:** https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233
**Category:** LEfSe
**Created:** [February 26, 2020, 4:53am UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233 "2020-02-26T04:53:29Z")
**Posts on this page:** 6
**Page:** 1

<div class="post-metadata">

### Author: ![Anushka\_Khasnobish](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/anushka_khasnobish/32/131_2.png) [@Anushka\_Khasnobish](https://forum.biobakery.org/u/Anushka_Khasnobish)
#### Post date: [February 26, 2020, 4:53am UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233/1 "2020-02-26T04:53:29Z")

</div>

A similar but not identical query already exists in the forum however it wasnt explained in detail so I am posting my query here.

I am trying to use the lefse for finding group of taxa which can successfully classify two disease groups from each other. As far as I understand, Lefse uses factorial KW rank sum test to select features that are differentially abundant between the two classes (say disease 1 and disease 2) and it used Wilcoxon test to further screen features which are consistently differentially abundant between the subclasses.

However in my dataset, there are no subclasses. In that case does Lefse skips Wilcoxon test ? And if it does, then are the selected biomarkers at the end of Lefse considered successful biomarkers or not?  
Thank you for the help.

---

<div class="post-metadata">

### Author: ![darmecian](https://avatars.discourse-cdn.com/v4/letter/d/47e85d/32.png) [@darmecian](https://forum.biobakery.org/u/darmecian)
#### Post date: [February 26, 2020, 6:13am UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233/2 "2020-02-26T06:13:18Z")

</div>

> [@LefSe: what happens if no subclass](https://forum.biobakery.org/t/lefse-what-happens-if-no-subclass/86):
>
> Hello, I’ve seen this question asked a few times on the previous lefse forum but never answered succinctly. If one does not provide a subclass, does LefSe run the wilcoxon? If so, what is it testing and will the result differ from only the KW? If one does provide a subclass, will the results from the wilcoxon show up in the “.res” file? Thanks in advance! John

I think you have an answer here.

---

<div class="post-metadata">

### Author: ![Anushka\_Khasnobish](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/anushka_khasnobish/32/131_2.png) [@Anushka\_Khasnobish](https://forum.biobakery.org/u/Anushka_Khasnobish)
#### Post date: [February 26, 2020, 6:32am UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233/3 "2020-02-26T06:32:03Z")

</div>

Thank you for the response. The answer that you linked only partly answers my question so if you can please kindly explain the following:

1. As I understand, in absence of subclass, the LefSe will perform per-feature KW, skip wilcoxon and perform LDA scoring and ranking. Is that right?
2. If so, then are the selected biomarkers at the end of Lefse considered successful biomarkers or not?

Thank you for the help.

---

<div class="post-metadata">

### Author: ![darmecian](https://avatars.discourse-cdn.com/v4/letter/d/47e85d/32.png) [@darmecian](https://forum.biobakery.org/u/darmecian)
#### Post date: [March 6, 2020, 3:57am UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233/4 "2020-03-06T03:57:46Z")

</div>

Sorry for my late reply. I understand what you mean, and I could not find answer in the thread I specified, as you stated. From my current understanding,

1. Absence of subclass, LEfSe will perform per-feature KW, and **NOT** skip wilcoxon unless you specified --wilc option to 0. Instead, LEfSe seem to perform both KW and Wilcoxon between class. I was examining format\_input.py and noticed this. I made another thread regarding this (if it is valid for lefse to perform KW and Wilcoxon both between class).  
[Question about LEfSe input\_format.py when specifing no subclass](https://forum.biobakery.org/t/question-about-lefse-input-format-py-when-specifing-no-subclass/253)

2. I think yes, considered they went through KW and LDA, however, in this case (–wilc 0), the p-value from KW, which is considered to be reported whether wilcoxon is performed or not, is not reported. I additionally asked this problem in the thread above.

Sincerely,

---

<div class="post-metadata">

### Author: ![sma](https://avatars.discourse-cdn.com/v4/letter/s/dec6dc/32.png) [@sma](https://forum.biobakery.org/u/sma)
#### Post date: [April 10, 2020, 7:18pm UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233/5 "2020-04-10T19:18:23Z")

</div>

Hi there,  
I’d agree with darmecian here. For 1, LEfSe will still perform Wilcoxon unless otherwise specified. For 2, I’d say so. Even in the absence of subclass tests LEfSe still filters for features that a) pass the KW test and b) has strong LDA score support. The two together should provide enough evidence for biomarkers.  
Thanks,  
Siyuan

---

<div class="post-metadata">

### Author: ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)
#### Post date: [November 18, 2020, 2:04pm UTC](https://forum.biobakery.org/t/lefse-without-any-subclass-still-valid/233/6 "2020-11-18T14:04:05Z")

</div>

Hi @sma @darmecian- what LDA score we should consider as “strong”? Can I use 1.5 instead of 2.0?

Thanks,  
DC7
