<p>I agree you have the right idea because that data doesn't look at all linear, so regression analysis looks really complicated. Your approach a straight comparison keeps it simple and is within the biological question you are asking, </p>

<p>However, I would use a non-parametric such as Wilcoxon instead of a T-test. I don't fully know the null distribution you are comparing it against. Wilcoxan doesn't require an assumption of normality. T-tests invariably require a transformation to normality and that gets complicated, but a log or sqrt transformation might be okay. If Wilcoxan produces a no significant difference you report it and move onto a T-test.</p>

<p>It is difficult to use outlier analysis with 9(?) observations, thats used if you have a way-out result and can afford to exclude it (>3 standard deviations). Outlier analysis is to remove a single result that is messing up your analysis and doesn't affect the bulk of the data. You alot of observations to be considering outlier analysis.</p>

CC BY-SA 4.0
2020-06-10T06:35:14.443Z
The actual dataset is 130 substrates thankfully. But I still think you are correct that it would not be normally distributed. Thus would a glm work here with either a negative binomial or poisson as the distribution?
Lamma
CC BY-SA 4.0
2020-06-10T06:43:35.253Z
I think negative binomial would be best as I don't think the variance is equal to the mean, but I guess Icould just test that assumption.
Lamma
CC BY-SA 4.0
2020-06-10T10:01:18.663Z
That would need fitting and at thiss point you are into formal regression analyses
M__
CC BY-SA 4.0
2020-06-10T10:06:25.510Z
I have changed the analysis now to have the frequence of occurance of each enzyme and lables them via the enzyme substrate so my model is glm.nb(substrate ~ frequency). I think this is better as it generates a mean and variance for each substrate thus allowing a much better identification of if the degridation of any singular substrate is enriched?
Lamma
CC BY-SA 4.0
2020-06-10T10:07:20.277Z
Okay thats a different question
M__
CC BY-SA 4.0
2020-06-09T15:23:15.293Z
2020-06-09T15:23:15.293Z
M__
<p>I have the following dataframe that is comprised of enzyme substrates as the variable and then the unique number of times enzymes capiable of degrading this substrates appeared in my genomes and then,total count of times all of those enzymes appeard in my samples. For example, a 7 unique enzymes were found to be capable of degrading chitin, and they occured a total of 298 times across my genomes.</p>
<p>I wish to see if any particular substrate has significantly larger occurence in both unique and total, than the rest of the substrates. Would outlier detection be appropriate for this or would a simple T-test do?</p>
<p><a href="https://i.sstatic.net/DdLVZ.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/DdLVZ.png" alt="enter image description here" /></a></p>

<r>
Is outlier detection the right choice for determining if a value lies significantly away from the mean in a population in gene frequency?
# Is outlier detection the right choice for determining if a value lies significantly away from the mean in a population in gene frequency?

(Asked by: username_0 on 2020-06-09)

I have the following dataframe that is comprised of enzyme substrates as the variable and then the unique number of times enzymes capiable of degrading this substrates appeared in my genomes and then,total count of times all of those enzymes appeard in my samples. For example, a 7 unique enzymes were found to be capable of degrading chitin, and they occured a total of 298 times across my genomes.

I wish to see if any particular substrate has significantly larger occurence in both unique and total, than the rest of the substrates. Would outlier detection be appropriate for this or would a simple T-test do?

[![enter image description here](https://i.sstatic.net/DdLVZ.png)](https://i.sstatic.net/DdLVZ.png)

**Tags:** r

---

## 1 answers

### Answer by username_1 on 2020-06-09. Score: 0

I agree you have the right idea because that data doesn't look at all linear, so regression analysis looks really complicated. Your approach a straight comparison keeps it simple and is within the biological question you are asking,

However, I would use a non-parametric such as Wilcoxon instead of a T-test. I don't fully know the null distribution you are comparing it against. Wilcoxan doesn't require an assumption of normality. T-tests invariably require a transformation to normality and that gets complicated, but a log or sqrt transformation might be okay. If Wilcoxan produces a no significant difference you report it and move onto a T-test.

It is difficult to use outlier analysis with 9(?) observations, thats used if you have a way-out result and can afford to exclude it (>3 standard deviations). Outlier analysis is to remove a single result that is messing up your analysis and doesn't affect the bulk of the data. You alot of observations to be considering outlier analysis.

> **Comments:**
> **username_0**: The actual dataset is 130 substrates thankfully. But I still think you are correct that it would not be normally distributed. Thus would a glm work here with either a negative binomial or poisson as the distribution?  
> **username_0**: I think negative binomial would be best as I don't think the variance is equal to the mean, but I guess Icould just test that assumption.  
> **username_1**: That would need fitting and at thiss point you are into formal regression analyses  
> **username_0**: I have changed the analysis now to have the frequence of occurance of each enzyme and lables them via the enzyme substrate so my model is glm.nb(substrate ~ frequency). I think this is better as it generates a mean and variance for each substrate thus allowing a much better identification of if the degridation of any singular substrate is enriched?  
> **username_1**: Okay thats a different question  

---

