版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

[求助] perl程序咨詢，有關(guān)pedigree信息分析的。

各位鄉(xiāng)親父老，小女又來啦~~~~

這次的input是：
#Fid Iid Pid Mid Sex Stat Geno
5    2908 0    0    1 3    A G T C G T
5    2108 6256 4532  2 1       G T C A T C
20 4556 2089 2935  1 2    T T G A C C
5    6256 2908 6266  1 1    A C C T T C
5    4561 1145 2108  2 1    C C T G A G
117 1123 0    2354  1 2    C C T G A C
基本上，從2908— 6256— 2108— 4561就是一個家族的四代人，這個就是我想找到的信息。從這個聯(lián)系來分析我的整個數(shù)據(jù)文件，有多少個家庭里有四代人的信息，有多少家族只有兩代人的信息，諸如此類。

現(xiàn)在我有的程序，已經(jīng)可以把每個相同家庭里的人，放在一個group里，接下來的分析，搞不懂了，特求助各位。

這里附上我的上半段程序：

CODE:

use warnings;

use strict;

open (FN, "file1") or die $!;

open(OUT,">file2")or die $!;

my %hash;

while (){

    chomp;

    if (/(\d+)\s([\d\w]+\s[\d\w]+\s[\d\w]+)/){

        push @{ $hash{$1} },$2;            ###This is to save multiple elements into a key.

       }

}

回復(fù)此樓

» 猜你喜歡

268求調(diào)劑已經(jīng)有7人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑已經(jīng)有4人回復(fù)
302求調(diào)劑已經(jīng)有11人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) 已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
313求調(diào)劑已經(jīng)有3人回復(fù)
南昌大學(xué)材料專碩311分求調(diào)劑已經(jīng)有6人回復(fù)
316求調(diào)劑已經(jīng)有6人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有7人回復(fù)

1樓 2012-09-10 16:40:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

繼續(xù)等

回復(fù)此樓

2樓2012-09-11 15:53:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

等不到有緣人啊~~~

贊一下

回復(fù)此樓

3樓2012-09-12 14:25:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

繼續(xù)咯~~~

回復(fù)此樓

4樓2012-09-13 14:58:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ch_sunshine

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 89.5
帖子: 31
在線: 10.9小時
蟲號: 2000904
注冊: 2012-09-14
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

若是單純想看看2、4代家系的個數(shù)信息，倒是有其他方法能得到。（比如：2代家系有N個）

贊一下

回復(fù)此樓

5樓2012-09-14 17:06:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

引用回帖:

5樓: Originally posted by ch_sunshine at 2012-09-14 17:06:43
若是單純想看看2、4代家系的個數(shù)信息，倒是有其他方法能得到。（比如：2代家系有N個）

比如呢？

回復(fù)此樓

6樓2012-09-14 19:46:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ch_sunshine

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 89.5
帖子: 31
在線: 10.9小時
蟲號: 2000904
注冊: 2012-09-14
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
silicare: 金幣+6, 應(yīng)助指數(shù)+1, 謝謝參與，鼓勵新蟲發(fā)帖交流 2012-09-17 17:28:23

引用回帖:

6樓: Originally posted by 淡藍(lán)幽香 at 2012-09-14 19:46:55
比如呢？...

用MERLIN軟件關(guān)于ibd和kinship評估功能，可以根據(jù)你的.ped文件內(nèi)容分析出家系結(jié)構(gòu)（包括家系大小，世代數(shù)等一些信息）。你可以看一下官網(wǎng)http://www.sph.umich.edu/csg/abecasis/Merlin/

贊一下

回復(fù)此樓

7樓2012-09-17 13:59:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

引用回帖:

7樓: Originally posted by ch_sunshine at 2012-09-17 13:59:48
用MERLIN軟件關(guān)于ibd和kinship評估功能，可以根據(jù)你的.ped文件內(nèi)容分析出家系結(jié)構(gòu)（包括家系大小，世代數(shù)等一些信息）。你可以看一下官網(wǎng)http://www.sph.umich.edu/csg/abecasis/Merlin/...

多謝，我去試一下
對了，請問怎么個你加分？

贊一下

回復(fù)此樓

8樓2012-09-18 16:10:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ch_sunshine

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 89.5
帖子: 31
在線: 10.9小時
蟲號: 2000904
注冊: 2012-09-14
專業(yè): 生物信息學(xué)

引用回帖:

8樓: Originally posted by 淡藍(lán)幽香 at 2012-09-18 16:10:03
多謝，我去試一下
對了，請問怎么個你加分？...

不知道我是剛申的這個。希望能幫到你

贊一下

回復(fù)此樓

9樓2012-09-18 17:33:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淡藍(lán)幽香

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 243.5
散金: 1
帖子: 24
在線: 16.8小時
蟲號: 644125
注冊: 2008-11-03
專業(yè): 生物信息學(xué)

引用回帖:

9樓: Originally posted by ch_sunshine at 2012-09-18 17:33:07

不知道我是剛申的這個。希望能幫到你...

引用回帖:
8樓: Originally posted by 淡藍(lán)幽香 at 2012-09-18 16:10:03
多謝，我去試一下
對了，請問怎么個你加分？...

不知道我是剛申的這個。希望能幫到你

評分顯示：
當(dāng)前回帖不是應(yīng)助回帖，您不能對此進(jìn)行金幣獎勵哦

請問還有其他方法加金幣么。。。。

贊一下

回復(fù)此樓

10樓2012-09-19 15:18:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主淡藍(lán)幽香的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +8	吃吃吃才有意義 2026-03-19	8/400	2026-03-21 00:49 by 劉國森
[考研] 22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4	一瓶苯甲酸 2026-03-14	4/200	2026-03-20 20:36 by fen_rao
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 086500 325 求調(diào)劑 +3	領(lǐng)帶小熊 2026-03-19	3/150	2026-03-20 18:38 by 盡舜堯1
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 生物學(xué)調(diào)劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015

24小時熱門版塊排行榜

[求助] perl程序咨詢，有關(guān)pedigree信息分析的。

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] perl程序咨詢，有關(guān)pedigree信息分析的。