我的標(biāo)題可能信息量不大,但這是一個(gè)暴露我問(wèn)題的例子:
我有這個(gè)數(shù)據(jù)框:
df=data.frame(cond1=c(1,1,1,2,2,2,3,3,3,1,1,1,2,2,2,3,3,3),
group=c("F","V","M","F","V","M","F","V","M","F","V","M","F","V","M","F","V","M"),
gene=c("A","A","A","A","A","A","A","A","A","B","B","B","B","B","B","B","B","B"),
value=c(1,2,3,4,5,6,7,8,9,1,3,2,4,3,2,2,3,4))
df
cond1 group gene value
1 1 F A 1
2 1 V A 2
3 1 M A 3
4 2 F A 4
5 2 V A 5
6 2 M A 6
7 3 F A 7
8 3 V A 8
9 3 M A 9
10 1 F B 1
11 1 V B 3
12 1 M B 2
13 2 F B 4
14 2 V B 3
15 2 M B 2
16 3 F B 2
17 3 V B 3
18 3 M B 4
我想得到的是,對(duì)于每個(gè)基因,有多少不同的cond1的值對(duì)應(yīng)于F組,小于它們對(duì)應(yīng)于V的值,它們對(duì)應(yīng)于M的值。
在前三個(gè)品系中,我們的第二個(gè)基因是A。對(duì)應(yīng)于組F=1、V=2、M=3的值。因此,對(duì)于cond1=1組的A基因,F(xiàn)<V<M。
我對(duì)基因A的預(yù)期輸出是3,因?yàn)樗衏ond1組的值都滿足F<V<M。我對(duì)基因B的預(yù)期輸出是1,因?yàn)橹挥衏ond1=3組的值符合F<V<M。
理想情況下,我希望的輸出是一個(gè)數(shù)據(jù)幀,其中包含gene和cond1之和,不滿足我的標(biāo)準(zhǔn):
gene count
1 A 3
2 B 1
如果你能給我一些建議,告訴我該怎么做,我將不勝感激
檢查所有數(shù)據(jù)是否按遞增順序排列,并計(jì)算每個(gè)
gene
存在多少這樣的值。