SNP分析命令
更新时间:2024-03-19 10:27:01 阅读量: 综合文库 文档下载
E:\\ > cd e: E:\\
E:\\ > cd plink-1
E:\\plink-1>plink –file test 1. Map 更新
Plink --sheep --file data --update-map position.txt --recode --out data1
Chrnew.txt -- update-chr --recode --out data2 Position: SNP code and position Chrnew:SNP code and Chr. 2.SNP merge
Plink --file data1 --merge data2.ped data2.map --recode --out merge 3.提取SNP位点
Plink --file data --extract 50kSNP.txt --recode --out data1 50kSNP.txt: 50k中的SNP名 4. Quality control
Call rate >98%/99%
Plink --file sheep --geno 0.02 --recode --out sheepgeno Plink --file sheepgeno --mind 0.01 --recode --out sheepmind
MAF>0.05
Plink --file sheepmind --maf 0.05 --recode --out sheepmaf Hardy-Weinberg equilibrium <0.0001
Plink --file sheepmaf --hwe 0.0001 --recode --out sheephwe
Exclude the SNP markers with either chromosome or both unknown Plink --sheep --file sheephwe --extract 4newsnp.txt --recode --out sheep4
Note: 制作4newsnp.txt(包含 chromosome 和 base-pair position 都为0的SNP) To identify sample duplication or half-sibs or closer Plink –sheep –file sheep4 –genome –max 0.85 Note:Check the genome file 5. LD quality control
Plink –sheep --file sheep4 –indep-pairwise 100 25 0.2 –out sheepld0.2 Plink --sheep --file sheep4 --indep-pairwise 100 25 0.05 --out sheepld0.05 Plink --file sheep4 --ld-window-r2 0.2 --out sheepldr0.2 输出结果为 data prunein 和 data prune out (质控时,要去除X染色体)
将data prune in 转化为ped和map
Plink --sheep --file 114hwe --extract 114sheep0.05.prune.in --recode --out sheepforpca
6. PCA-
PCA的三个文件:
Plink --sheep --file data(生成LD的文件) --extract data (LD).prune.in --recode --out sheepforpca 1sheepforpca.ped 改为5.ped 2sheepforpca.map 改为5.pedsnp
3将sheepforpca 制作成二进制文件 输出5b plink --file hapmap1 --make-bed --out hapmap1
结果为5b.farm即为ped文件的前6列, 将5b.farm 改名为5.pedind
Note: 5.pedind 文件中要将第六列-9换成familyID. 参数文件 Genotypename: 5.ped Snp name: 5.pedsnp Indivame: 5.pedind
Evecoutname: 5.pca.evec Evaloutname: 5.eval Altnormstyle: NO Numoutevec: 3 Numoutlieriter: 5 Numoutlierevec: 10 Outlier sigmathresh: 6.0 Qt mode: NO
将上述文件拷贝到eigensoft/bin 文件夹内 打开命令
Cd EIG5.01/bin 作图命令
./smartpca –p 5.par
./ploteig –I 5.pca.evec –c 1:2 –p AL:BSB:…Tiberan –x 5-0 即可得到PCA
在5.pca.evec文件中可以看到主成分占的比例。 7 原始SNP数据转化成map和ped文件
>data=read.csv(\>Ta=t(data)
> write.table(Ta,file=\检查命令:--Compound –genotypes
8: 近交系数计算, 多态性含量, Ho He 哈温P值) plink --file filename --het --out filename1
plink --file filename --homozyg --out filename1
plink --file filename --hardy --out filename1(Plink --file filename –hardy,9: ADZE软件计算Ar 和 pAr Plink 转化成structure
1Plink --file filename --recode-structure --out filename1 2用PGDspider 转化ped 文件为structure结构。
将plink转化的位点信息粘贴到PGDspider转化的文件中
全基因组关联分析
plink --file data --remove mylist.txt --recode --out filename plink --sheep --file filename --out name –assoc
plink --sheep --file filename --out name --assoc --adjust R软件中绘制manhattan 图 先安装qqman软件 > library(qqman)
> results<- read.table(\
> manhattan(results, ylim = c(0, 10), col = c(\
结果为plink.hwe) R软件中绘制qq-plot图 > library(qqman)
> results<- read.table(\> qq(results$P)
用其他关联分析方法:
plink --sheep --file mar --out mar-model --model --model-trend --adjust LD 分析(haploview)
1 info 文件生成:plink --file hu-M2 --recodeHV --out hu-M2HV R 安装 GenABEL
Install packages (GenABEL) Install packages (MASS)
Install packages(Gen ABEL.data) 加载安装包 Library (MASS)
Library(Gen ABEL.data) Library(GenABEL)
在使用GenABEL前需要准备4个文件
Ped、map、phen(当ped中含有多个表型时用到)、praw 1生成tped、tfram 文件
Plink –file name –transpose –recode –out gwa-gabel
当多个表型时还还需要—pheo phenol.phen –pheno-name 2制作praw文件
格式 id sex phen(sex:female=0, male=1 phen case=1 control=0) ―Sss12‖ 1 0
―Sss18‖ 1 1 GLM test
Testb<- scan.glm(?phen~CRSNP‘, family=binomial(),data=b.dat) 2 score test
Testb.qt<- qtscore(phen, data=b.dat, trait=‖binomial‖) Test.qt@lambda
正在阅读:
SNP分析命令03-19
感恩节感谢父母的作文06-13
小微企业融资难融资贵深层次原因解析03-20
现代推销学--第一次平时作业2022秋华南理工大学网络教育答04-08
语文课程与教学论的课后习题 一12-07
七年级下册数学《相交线与平行线》相交线 - 知识点整理03-16
服装厂有关计件工资04-04
小学数学教学反思21篇05-28
广东省2015年高新技术企业培育库拟入库企业名单(第1.2.3批全)06-03
深山养芒鼠助民走富路05-15
- 多层物业服务方案
- (审判实务)习惯法与少数民族地区民间纠纷解决问题(孙 潋)
- 人教版新课标六年级下册语文全册教案
- 词语打卡
- photoshop实习报告
- 钢结构设计原理综合测试2
- 2014年期末练习题
- 高中数学中的逆向思维解题方法探讨
- 名师原创 全国通用2014-2015学年高二寒假作业 政治(一)Word版
- 北航《建筑结构检测鉴定与加固》在线作业三
- XX县卫生监督所工程建设项目可行性研究报告
- 小学四年级观察作文经典评语
- 浅谈110KV变电站电气一次设计-程泉焱(1)
- 安全员考试题库
- 国家电网公司变电运维管理规定(试行)
- 义务教育课程标准稿征求意见提纲
- 教学秘书面试技巧
- 钢结构工程施工组织设计
- 水利工程概论论文
- 09届九年级数学第四次模拟试卷
- 命令
- 分析
- SNP
- 中国传媒大学南广学院文化管理学院
- 学生社团公益创业实践与职业素养提升的生态体系构建 -
- 一位野路子项目经理的PMP考试5A级成功经验
- 论文材料10-论文模板
- 2015-2016年广东省深圳市龙华区东王实验学校九年级上学期
- 足球基本技术第三次课教案 - 图文
- 清原满族自治土地整治规划 - 图文
- 企业有机加工操作规程模板
- 感官评价
- 植物生理学重点
- 7111斜巷修护安全技术措施
- 通风设施工技术操作规程
- 2019运用多媒体辅助教学手段提高数学的课堂教学效率语文
- WQS软件介绍 - 图文
- 商业银行零售业务市场营销策略探析
- Z010关于制定2014年度综合计划的通知
- 2011河北职称计算机考试操作题答案
- 2019-2020年高三一轮复习测试(二)数学理试题 含答案
- 恢复原状请求权的性质
- 证券投资实务课后答案