教育资源为主的文档平台

当前位置: 查字典文档网> 所有文档分类> 高等教育> 生物学> admixture-manual

admixture-manual

上传者:李小滨
|
上传时间:2015-04-24
|
次下载

admixture-manual

admixture软件使用手册

Admixture1.23SoftwareManualDavidH.AlexanderJohnNovembre

August22,2013KennethLange

admixture软件使用手册

Contents1Quickstart

2Reference

2.1HowdoIchoosethecorrectvalueforK?..........

2.1.1Background.......................

2.2HowdoIplottheQestimates?................

2.3DoIneedtothinthemarkersetforlinkagedisequilibrium?

2.4HowmanymarkersdoIneedtosupplytoADMIXTURE?.

2.5HowdoIchangetherandomseed?’.............

2.6OptimizationMethods.....................

2.7Terminationcriteria......................

2.8Accelerationmethods.....................

2.9Bootstrapping..........................

2.10Supervisedanalysis.......................

2.11Multithreadedmode......................

2.12Penalizedestimation......................

3Citationandfurtherinformation.....................................................................................................................1334556777888101010Notetousers:thematrixformerlyreferredtoasF,containingthepopulationallelefrequencies,isnowknown(asofv1.20)asP,andisoutputtoa.P?le.ThisistobringournotationinlinewiththatoftheSTRUCUTREpapers.

1Quickstart

ADMIXTUREisaprogramforestimatingancestryinamodel-basedmannerfromlargeautosomalSNPgenotypedatasets,wheretheindividualsareunrelated(forexample,theindividualsinacase-controlassociationstudy).

ADMIXTURE’sinputisbinaryPLINK(.bed),ordinaryPLINK(.ped),orEIGENSTRAT(.geno)formatted?lesanditsoutputissimplespace-delimited?lescontainingtheparam-eterestimates.

TouseADMIXTURE,youneedaninput?leandanideaofK,yourbeliefofthenumberofancestralpopulations.Youshouldalsohavetheassociatedsupport?lesalongsideyourmaininput?le,inthesamedirectory.Forexample,ifyourprimaryinput?leisa.bed?le,youshouldhavetheassociated.bim(binarymarkerinformation?le)and.fam(pedigreestub?le)?lesinthesamedirectory.Ifyourprimaryinput?leisa.pedor.geno?le,acorrespondingPLINKstyle.map?leshouldbeinthesamedirectory.

1

admixture软件使用手册

IfyouhaveanbinaryPED(.bed)formatted?leinyourcurrentdirectory

%ls

hapmap3.bedhapmap3.bimhapmap3.fam

andyoubelievethattheindividualsinthesamplederivetheirancestryfromthreeancestralpopulationsthenrunadmixturelikethis:

%admixturehapmap3.bed3

ADMIXTUREwillstartrunning.Hopefullyitwill?nishsoon,anditwillthenoutputsomeestimates:

%ls

hapmap3.bedhapmap3.bimhapmap3.famhapmap3.3.Qhapmap3.3.P

Thereisanoutput?leforeachparameterset:Q(theancestryfractions),andP(theallelefrequenciesoftheinferredancestralpopulations).Notethattheoutput?lenameshave‘3’inthem.Thisindicatesthenumberofpopulations(K)thatwasassumedfortheanalysis.This?lenameconventionmakesiteasytorunanalysesusingdi?erentvaluesofKinthesamedirectory.

IfyouhaveaPLINK.ped“12”coded?le—generatedbyacommandlike

%plink--filehapmap--recode12--outhapmap

—thenyoufollowbasicallythesameinstructions:

%admixturehapmap3.ped3

NotethatADMIXTUREinferstheformatofyourinput?lebasedonthe?leextension(.bed,.ped,or.geno).

YoucanrunADMIXTUREonainput?lelocatedinanotherdirectory,butitsoutputwillalwaysbeputinthecurrentworkingdirectory.Forexample:

%pwd

/home/dalexander/AdmixtureRuns

%ls

(nooutput—currentdirectoryisempty)

%ls~/Data

hapmap3.pedhapmap3.map

%admixture~/Data/hapmap3.ped3

(waitforittorun)

2

admixture软件使用手册

%ls

hapmap3.3.Qhapmap3.3.P

%ls~/Data

hapmap3.pedhapmap3.map

Ifyoualsowantedstandarderrors,theninsteadofthelastcommandyoushouldhaveused

%admixture-B~/Data/hapmap3.ped3

Thiswillperformpointestimationandthenwillalsouseabootstrappingproceduretocalculatethestandarderrors.Notethat(point-estimation&bootstrapping)takesconsid-erablylongerthanpoint-estimationalone,soyouwillhavetobepatient.Eventuallyitwill?nish,yieldingpointestimatesandstandarderrors:

%ls

hapmap3.3.Qhapmap3.3.Phapmap3.3.Q_se

The“se”?leisinthesameunadorned?leformatasthepointestimates.

Ifyouranalysesaretakingalongtime,?rstconsiderwhy.Isyourdatasethuge?Doyoureallyneedtoanalyzeallthemarkers,orcanyouthinthemarkerset(asdescribedinSection2.3)?Ifyoureallydoneedtoanalyzeallthemarkers,orifyouranalysisisslowishforotherreasons(verylargeK,extensivecross-validationorbootstrapping)youmightconsiderrunningADMIXTUREinmultithreadedmode.Ifyourcomputerhasfourprocessors,forexample,youcoulduseacommandlike

%admixture~Data/huge_dataset.bed3-j4

tosplitADMIXTURE’sworkamongfourthreads—whichcanmakeADMIXTURErunalmostfourtimesasfast.

2

2.1ReferenceHowdoIchoosethecorrectvalueforK?

UseADMIXTURE’scross-validationprocedure.AgoodvalueofKwillexhibitalowcross-validationerrorcomparedtootherKvalues.Cross-validationisenabledbysimplyaddingthe--cv?agtotheADMIXTUREcommandline.Inthisdefaultsetting,thecross-validationprocedurewillperform5-foldCV—youcanget10-foldCV,forexample,using--cv=10.Thecross-validationerrorisreportedintheoutput.Forexample,ifinourbashshellweran

3

admixture软件使用手册

%forKin12345;\

doadmixture--cvhapmap3.bed$K|teelog${K}.out;done

(i.e.,ranADMIXTUREwithcross-validationforKvalues1,2,3,4and5),thenwecouldquicklyviewtheCVerrors:

%grep-hCVlog*.out

CVerror(K=1):0.55248

CVerror(K=2):0.48190

CVerror(K=3):0.47835

CVerror(K=4):0.48236

CVerror(K=5):0.48985

wherethenumberinparenthesesisthestandarderrorofthecross-validationerrorestimate.Wecaneasilyplotthesevaluesforcomparison,asinFigure

内容需要下载文档才能查看

1,whichmakesitfairlyclearthatK=3isasensiblemodelingchoice.

Figure1:Cross-validationplotforthehapmap3dataset

2.1.1Background

ThosefamiliarwithSTRUCTUREknowthatitprovidesameansofidentifyingthe“best”valueforK,thenumberofpopulations,basedoncomputingthemodelevidenceforeachpossibleKvalue.ThemodelevidenceforKisde?nedas

Pr(G|K)=f(G|Q,P,K)π(Q,P|K)dQdP

4

版权声明:此文档由查字典文档网用户提供,如用于商业用途请与作者联系,查字典文档网保持最终解释权!

下载文档

热门试卷

2016年四川省内江市中考化学试卷
广西钦州市高新区2017届高三11月月考政治试卷
浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
广西钦州市钦州港区2017届高三11月月考政治试卷
广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
广西钦州市高新区2016-2017学年高二11月月考政治试卷
广西钦州市高新区2016-2017学年高一11月月考政治试卷
山东省滨州市三校2017届第一学期阶段测试初三英语试题
四川省成都七中2017届高三一诊模拟考试文科综合试卷
2017届普通高等学校招生全国统一考试模拟试题(附答案)
重庆市永川中学高2017级上期12月月考语文试题
江西宜春三中2017届高三第一学期第二次月考文科综合试题
内蒙古赤峰二中2017届高三上学期第三次月考英语试题
2017年六年级(上)数学期末考试卷
2017人教版小学英语三年级上期末笔试题
江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
四川省简阳市阳安中学2016年12月高二月考英语试卷
四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
安徽省滁州中学2016—2017学年度第一学期12月月考​高三英语试卷
山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷

网友关注

2017国考国家物资储备局面试真题及解析参考答案 (2月24日上午)
2017国考国税面试真题及解析参考答案 (3月2日)
2017国考黄河水务局面试真题及解析参考答案 (2月26日)
2017国家公务员面试国税系统面试真题(考生回忆篇)
面试题库:2017国考面试每日一练结构化面试模拟题3.2
2017国考海事局面试真题及解析参考答案 (2月25日)
2017国考中央办公厅面试真题及解析参考答案 (2月24日上午)
面试题库:2017国考面试每日一练结构化面试模拟题3.6
2017国考甘肃能源局面试真题及解析参考答案 (2月24日上午)
2017国考国安局面试真题及解析参考答案 (2月25日上午)
2017国考工信部、统计局等面试真题及解析参考答案 (2月24日)
2017国考国税系统面试真题及解析参考答案 (3月3日)
2017国考铁路公安面试真题及解析参考答案 (2月24日)
2017国考铁路公安面试真题及解析参考答案 (2月27日)
2017国考工信部面试真题及解析参考答案 (2月25日上午)
2017国考铁路公安面试真题及解析参考答案 (2月26日上午)
2017国家公务员面试国税系统面试真题(3月1日)
面试题库:2017国考面试每日一练结构化面试模拟题3.7
2017国考水利部面试真题及解析参考答案 (2月24日下午)
2017国考广西保监局面试真题及解析参考答案 (2月24日上午)
2017国考统计局面试真题及解析参考答案 (2月24日上午)
2017国考工信部面试真题及解析参考答案 (2月24日上午)
2017国考统考面试真题及解析参考答案 (2月23日)
面试题库:2017国考面试每日一练结构化面试模拟题答案3.6
2017国考国家卫计委面试真题及解析参考答案 (2月24日上午)
2017国考国家知识产权局面试真题及解析参考答案 (2月25日上午)
面试题库:2017国考面试每日一练结构化面试模拟题答案3.2
2017国考黄河水务局面试真题及解析参考答案 (2月25日上午)
2017国考海事局面试真题及解析参考答案 (2月27日)
2017国考司法部(面试真题及解析参考答案 (2月24日上午)

网友关注视频

《空中课堂》二年级下册 数学第一单元第1课时
【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,天津市
七年级英语下册 上海牛津版 Unit3
冀教版小学数学二年级下册第二单元《有余数除法的整理与复习》
三年级英语单词记忆下册(沪教版)第一二单元复习
沪教版八年级下册数学练习册21.3(3)分式方程P17
【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,湖北省
沪教版牛津小学英语(深圳用) 四年级下册 Unit 4
【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,广东省
沪教版牛津小学英语(深圳用) 五年级下册 Unit 10
冀教版小学英语四年级下册Lesson2授课视频
二年级下册数学第二课
七年级英语下册 上海牛津版 Unit5
沪教版八年级下次数学练习册21.4(2)无理方程P19
人教版二年级下册数学
【部编】人教版语文七年级下册《老山界》优质课教学视频+PPT课件+教案,安徽省
冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣
19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T502436
3月2日小学二年级数学下册(数一数)
30.3 由不共线三点的坐标确定二次函数_第一课时(市一等奖)(冀教版九年级下册)_T144342
飞翔英语—冀教版(三起)英语三年级下册Lesson 2 Cats and Dogs
七年级英语下册 上海牛津版 Unit9
青岛版教材五年级下册第四单元(走进军营——方向与位置)用数对确定位置(一等奖)
精品·同步课程 历史 八年级 上册 第15集 近代科学技术与思想文化
19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T3763925
《小学数学二年级下册》第二单元测试题讲解
【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,辽宁省
3.2 数学二年级下册第二单元 表内除法(一)整理和复习 李菲菲
外研版英语七年级下册module3 unit2第一课时