plink2學(xué)習(xí)(一)
最近學(xué)東西很雜料滥,基本上是做啥學(xué)啥然眼,嘛,又懶得做葵腹。高每。屿岂。。就醬紫鲸匿。爷怀。。带欢。
- 下載plink2运授,下載解壓之后就是個(gè)可執(zhí)行文件了
- 把plink2配置到全局變量,參考我同事的筆記乔煞,懶得指路了
- 研究研究plink2的基本用法
plink2下載下來吁朦,包里是沒有example的,那我就只能看著別人的筆記看看這些命令都是啥意思了渡贾。
#以plink的命令為例逗宜,畢竟我么得example,
plink --file toy --freq --out toy_analysis
#plink以執(zhí)行時(shí)空骚,命令時(shí)用plink開頭的嘛纺讲,但是plink2就是以plink2開頭的
#--file 代表的輸入文件 只要所有輸入文件的前綴就可以
#--freq 這一層是表示要進(jìn)行運(yùn)算,后續(xù)有很多其他的運(yùn)算囤屹,用的時(shí)候就是換名字換成其他運(yùn)算熬甚,當(dāng)然還有各種參數(shù),這里freq指的是統(tǒng)計(jì)等位基因的頻率
#--out 指定輸出文件肋坚,默認(rèn)輸出‘plink.xyz','.xyz'就是其中一個(gè)擴(kuò)展名乡括,只要指定前綴就可以,當(dāng)然--out之后前面的'plink'就會(huì)被替換成你提供的前綴冲簿。建議前綴不要重復(fù)粟判,因?yàn)闀?huì)抹掉前面同名的文件。如果你不設(shè)置的話峦剔,應(yīng)該是會(huì)報(bào)錯(cuò)的
看起來就很簡單的操作,等我用起來之后再吐槽吧角钩,感覺不像看起來這么簡單吝沫,科研bug我本人。递礼。惨险。。