我有下列排序基因數據的12700 X 307數據幀:在列的個別ID號碼,軌跡標識行(注意兩行每個基因座)改變交錯行中的數據幀
alist<-c("loci",185,186,187,188,189,190,191,"A549",1,1,1,1,1,1,1,"A549",0,0,1,1,1,0,1,"A588",1,1,1,1,1,1,1,"A588",0,0,0,0,0,0,1,"A794",1,1,1,1,1,1,1,"A794",1,0,1,0,1,1,0,"A081",1,1,1,1,1,1,0,"A081",1,1,1,1,1,1,1)
df <- data.frame(matrix(unlist(alist), nrow=9, byrow=T),stringsAsFactors=FALSE)
colnames(df) = df[1, ]
df<-df[-1, ]
我需要將其更改爲每個單獨兩行的數據框,每個軌跡一列。個體的第一行應該有第一個等位基因的存在/不存在條目,第二行存在該位置的第二個等位基因。
因此,它應該是這樣的:
blist<-c("individual","A549","A588","A794","A081","185",1,1,1,1,"185",0,0,1,1,"186",1,1,1,1,"186",0,0,0,1,"187",1,1,1,1,"187",1,0,1,1,"188",1,1,1,1,"188",1,0,0,1,"189",1,1,1,1,"189",1,0,1,1,"190",1,1,1,1,"190",0,0,1,1,"191",1,1,1,0,"191",1,1,0,1)
dfb <- data.frame(matrix(unlist(blist), nrow=15, byrow=T),stringsAsFactors=FALSE)
colnames(dfb) = dfb[1, ]
dfb<-dfb[-1, ]
它必須是相當做,能,但我沒有看到它。我會很感激任何想法。