2013-11-26 41 views
-3

的名單上有基因匹配的基因組位置的這個名單[R數出現在數字矢量

$HES2 
[1] 6472685 6472960 6473183 6473556 6473721 6473880 6474011 6474016 6474021 
[10] 6474026 6474031 6474036 6474595 6475487 6475525 6475625 6475843 6475847 
[19] 6475851 6476040 6476416 6477257 6477498 6477543 6478033 6478181 6480077 
[28] 6481025 6481370 6481426 6483450 6483762 6483892 6484023 6484664 

$HGFAC 
[1] 3444133 3444265 3445567 3448970 3449333 3450621 

我想算在段倉的出現,讓我們說100個鹼基對的,並報告基因在這個特殊的段數,讓我告訴你虛擬數據提前

+1

你能告訴我們你到目前爲止嘗試過的嗎? – Henrik

+1

你是如何獲得150和589? – zx8754

+0

另外,請提供確切的輸入和輸出格式。從描述中,您提供的最後一個對象似乎混合了輸入(bin)和輸出(頻率)。 – flodel

回答

3

不知道

HES2 6474026-6474126 150 
HES2 6475038-6475138 589 

感謝,如果我失去了一些東西,但是這是你所追求的?

HES2<-c("6472685" ,"6472960" ,"6473183" ,"6473556" ,"6473721" ,"6473880" ,"6474011", "6474016" ,"6474021","6474026" ,"6474031" ,"6474036" ,"6474595", "6475487", "6475525" ,"6475625", "6475843", "6475847","6475851", "6476040" ,"6476416", "6477257", "6477498", "6477543" ,"6478033" ,"6478181" ,"6480077","6481025" ,"6481370", "6481426" ,"6483450" ,"6483762" ,"6483892", "6484023" ,"6484664") 

HES2<-as.numeric(HES2) 
bin_w<-100 
h=hist(x=HES2,breaks=bin_w,plot=FALSE) 
counts<-data.frame(from=h$breaks,to=h$breaks+bin_w,count=c(h$count,0)) 
counts[which(counts$count>0),] 

     from  to count 
1 6472600 6472700  1 
4 6472900 6473000  1 
6 6473100 6473200  1 
10 6473500 6473600  1 
12 6473700 6473800  1 
13 6473800 6473900  1 
+0

感謝特洛伊,是一個初步的,但我認爲它可以適用於列表 – user976991

+0

@ user976991,採取特洛伊的答案,然後在'plyr'包中看到'?ldply'。 –