使用R除去NA之外，把其他元素换成0和1，怎么处理？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在下是R新手最近处理基因文件的时候要求将缺失值NA保留，其他碱基值置换为0或1，怎么处理啊？
文件格式大约就是这样:
C C G G NA NA T T G G NA…..
大约需要些什么函数以及思路如何，希望大神救济一下。
这是从师兄那里模仿来的代码，但是无法好好运行，结果出来只有行列数，为什么？
markers <- read.table(“~/Documents/176GATK.INDELsnpsFiltered0.05M0.25_outfile.txt”,header=TRUE,sep=”\t”,as.is=T,check.names=FALSE,row.names=NULL)

convert.snp <- function(x) {
alleles <- na.omit(unique(x))
y<- rep(NA,length(x))
y[which(x==alleles[1])] <- 0
y[which(x==alleles[2])] <- 1
return(y)}

M <- apply(markers[,-c(1:2)],1,convert.snp)
dim(M)
library(rrBLUP)
N <- M
A <- A.mat(N)
dim(A)
write.table(dim(A),”/Users/huli/Documents/176GATK.INDELsnpsFiltered0.05M0.25_outputfile.txt”,sep=”\t”,quote=F,col.names=TRUE,row.names=TRUE)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：怎么处理 marker marke Mark 文件格式 元素 如何

碱基值不是有四个么，你这用0 1怎么替换？

jiangbeilu 发表于 2015-6-10 10:47
碱基值不是有四个么，你这用0 1怎么替换？

除了NA以外随机替换，现在想想是有点不靠谱，但是姑且问一下，可能吗？

这个是可以做到的呀

用ifelse试试