2017-06-08 5 views
0

通常、2つのグループに3つの条件の変数があります。この2つのグループ間に統計的な違いはありますか?例えば短いCSVについて分子遺伝学における2つのグループの違い

:5グラム/ 5グラム、4G/5グラムおよび4G/4G:

genotype,group,gene 
5g/5g,1,PAI 
5g/5g,1,PAI 
4g/5g,1,PAI 
4g/5g,1,PAI 
4g/5g,1,PAI 
4g/5g,1,PAI 
5g/5g,1,PAI 
4g/4g,1,PAI 
4g/4g,1,PAI 
4g/4g,2,PAI 
4g/5g,2,PAI 
4g/4g,2,PAI 
4g/4g,2,PAI 
4g/4g,2,PAI 
4g/5g,2,PAI 
5g/5g,2,PAI 
4g/4g,2,PAI 
5g/5g,2,PAI 
5g/5g,2,PAI 
5g/5g,2,PAI 
4g/5g,2,PAI 

遺伝子PAIは、3つの遺伝子型を持っています。

表形式でこのデータ存在する場合: - 絵下記だから、

enter image description here

、どのように私は1つのp値like thereを見つけることができますか? picture below

+0

明確にするために、何が違いますか? 2つのグループ間の各遺伝子型の割合? –

+0

はい、そうです。カイ二乗検定を使った違いを分析するために私が見ることができる(最後の絵)。そしてp値= 0.56。したがって、2つのグループに統計的な違いはありません。 – autumnrustle

答えて

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 

これは実際のデータではありませんが、同じことを行うことができます。

+1

火格子!しかし、私は警告 "chisq.test(テーブル(グループ、遺伝子型))があるたびに:カイ二乗近似は間違っているかもしれません"。それはテストデータセットでいくつかの値が非常に小さいbecaueですか?だから、いくつかの値が非常に小さい(時には私は実際には0または1の値を持っている)私はカイ二乗検定を使わないといけないのですか?私はfisher.testを使う必要がありますか? – autumnrustle

+1

良い点。カイ二乗近似は、小さなサンプルサイズでは信頼性が低くなります。予想されるカウントがある場合は、期待されるカウントのテーブルを生成することができます(一般的に期待通りのカウントがあれば 'with(df、chisq.test(table(group、genotype))) 5未満の場合、カイ2乗検定は推奨されません。代わりに、フィッシャーの正確なテスト(あなたはすでに見つかっています!)またはランダム化に基づいたカイ二乗検定です(df、chisq.test(テーブル(グループ、遺伝子型)、simulate.p.value = T)) ' 。たぶん両方を試して、p値が似ているかどうかを確認してください。 –

+0

大きな感謝!それは私があなたに多くのポイントを与えることはできません悲しいです) – autumnrustle

関連する問題