<p>Got it! :)</p>
<p>I switched the <code>summary(inv)</code> to <code>summary(inv, alpha=1)</code> to get all possible results and all my groups showed up.</p>
<p>So, the answer to my own question is that the groups are missing from the output because they don't have any statistically significant OTUs correlated with those groups.</p>

CC BY-SA 4.0
2024-03-27T14:52:21.740Z
2024-03-27T14:52:21.740Z
Rainman
<p>Based on presence absence data, I'm attempting to find indicator species for different sample sites. I have 25 different groups or &quot;clusters&quot; as I've labelled them. Running the R code below, only results for 13 of the 25 groups was returned.</p>
<p>Is the output truncated, i.e. is there a max output parameter in R I need to fix? Or does this imply that there simply aren't any statistically relevant indicator species in the groups not included in the output?</p>
<p>Code used:</p>
<pre class="lang-r prettyprint-override"><code>library(indicspecies)

data = read.csv(&quot;ISAI.txt&quot;, header= TRUE, sep= &quot;\t&quot;) #Load up data from file

abund = data[,3:ncol(data)] #Collect only numerical data, starts in 3rd column
cluster = data$Cluster #Groups are defined in the Cluster column of the imported data

inv = multipatt(abund, cluster, func = &quot;r.g&quot;, control = how(nperm=1000),duleg = TRUE) #Run analysis.

summary(inv) #print out results
</code></pre>
<p><strong>Notes</strong>:</p>
<p>For the <code>func</code> parameter, &quot;r&quot; was used due to data being in presence/absence format. The &quot;g&quot; corrects for groups with different numbers of members. The <code>duleg = True</code> prevents group site combinations from being assessed as I'm not interested in those.</p>

<r><features>
multiplatt from indicspecies package not producing results for all groups
# multiplatt from indicspecies package not producing results for all groups

(Asked by: username_0 on 2024-03-27)

Based on presence absence data, I'm attempting to find indicator species for different sample sites. I have 25 different groups or "clusters" as I've labelled them. Running the R code below, only results for 13 of the 25 groups was returned.
Is the output truncated, i.e. is there a max output parameter in R I need to fix? Or does this imply that there simply aren't any statistically relevant indicator species in the groups not included in the output?
Code used:
library(indicspecies)

data = read.csv("ISAI.txt", header= TRUE, sep= "\t") #Load up data from file

abund = data[,3:ncol(data)] #Collect only numerical data, starts in 3rd column
cluster = data$Cluster #Groups are defined in the Cluster column of the imported data

inv = multipatt(abund, cluster, func = "r.g", control = how(nperm=1000),duleg = TRUE) #Run analysis.

summary(inv) #print out results

Notes:
For the func parameter, "r" was used due to data being in presence/absence format. The "g" corrects for groups with different numbers of members. The duleg = True prevents group site combinations from being assessed as I'm not interested in those.

**Tags:** r, features

---

## 1 answers

### Answer by username_0 on 2024-03-27. Score: 1

Got it! :)

I switched the `summary(inv)` to `summary(inv, alpha=1)` to get all possible results and all my groups showed up.

So, the answer to my own question is that the groups are missing from the output because they don't have any statistically significant OTUs correlated with those groups.


---

