2021-04-19 用于筛选近交系数的R代码

之前写的R脚本用于筛选fit小于0.8的值

[lyc@200server ~]$ vim fit.R

data<-read.table("plink.hwe")
A<-data.frame(data[,])
A$Fit<-1-as.numeric(A[,7])/as.numeric(A[,8])
clean2fit<-subset(A,Fit>0.8)
write.table(clean2fit,"/home/lyc/clean2fit.hwe")

scp报错:not a regular file
原因是 这样是相当于下载文件夹,而非文件。上传文件夹同样用-r

解决办法是 scp后加参数 -r

data<-read.table("/Users/lyc/Desktop/head50.hwe")
A<-data.frame(data[,])
A$Fit<-1-as.numeric(A[,7])/as.numeric(A[,8])
write.table(A,"/Users/lyc/Desktop/allfit")
Fit<- A[,10]
hist(Fit,freq=T)
axis(1,c(0,0.5,1,1.5,2))
clean2fit<-subset(A,Fit>=0.8)
Clean2fit<-clean2fit[,10]
hist(Clean2fit,freq=T)
write.table(clean2fit,"/Users/lyc/Desktop/clean2fit.hwe")

修改之后的脚本,可以输出筛之前的fit和筛之后的fit


是烟草基因型vcf数据

以及对应的map、ped文件
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • [toc] 前提准备 系统信息 cat /etc/redhat-release 结果: Red Hat Enter...
    feihonInk阅读 1,346评论 0 1
  • 简明 VIM 练级攻略搜索文件夹下是否包含某个字符串systemctl鸟哥的 Linux 私房菜学习总结(超赞!!...
    _fhs阅读 349评论 0 0
  • 1. Linux命令 1.1 常用命令 显示文件或目录ls(选项)(参数)-l 表示long,长格式列出-a 表示...
    nimw阅读 6,272评论 1 12
  • 主目录 / 根 每个文件和目录都从根目录开始。 只有root用户拥有这个目录下的写权限。 请注意/ root是ro...
    杨健kimyeung阅读 433评论 0 0
  • 今天感恩节哎,感谢一直在我身边的亲朋好友。感恩相遇!感恩不离不弃。 中午开了第一次的党会,身份的转变要...
    迷月闪星情阅读 10,620评论 0 11