分子遗传学中两组之间的差异

时间:2017-06-08 17:50:29

标签: r statistics

通常我在2组中的3个条件中有一些变量。我如何才能发现,这两组之间是否存在统计差异?

例如短csv:

MouseUp

基因PAI有3种基因型:5g / 5g,4g / 5g和4g / 4g。

如果以表格形式提供此数据:

enter image description here

那么,我怎样才能找到下面的一个p值like there - 图片? picture below

1 个答案:

答案 0 :(得分:1)

这可能是你正在寻找的吗?

## simulating some data
df <- data.frame(genotype=sample(c("5g/5g","4g/5g","4g/4g"),16,replace=T), group=c(rep(1,8), rep(2,8)))

df
##    genotype group
## 1     4g/5g     1
## 2     4g/5g     1
## 3     4g/5g     1
## 4     5g/5g     1
## 5     5g/5g     1
## 6     4g/4g     1
## 7     4g/4g     1
## 8     4g/5g     1
## 9     4g/4g     2
## 10    4g/4g     2
## 11    4g/4g     2
## 12    4g/5g     2
## 13    5g/5g     2
## 14    5g/5g     2
## 15    4g/5g     2
## 16    4g/4g     2

## first a table
with(df, table(group, genotype))
##      genotype
## group 4g/4g 4g/5g 5g/5g
##     1     2     4     2
##     2     4     2     2

## and maybe a chi-squared test is what you want?
with(df, chisq.test(table(group,genotype)))
##    Pearson's Chi-squared test

## data:  table(group, genotype)
## X-squared = 1.3333, df = 2, p-value = 0.5134

请注意,这不是您的真实数据,但您可以做同样的事情。