Arrays 将数据重塑为CSV

Arrays 将数据重塑为CSV,arrays,r,matrix,Arrays,R,Matrix,我正在编写一个脚本来创建csv文件进行分析。运行脚本时,OPA 5701和6561分别提供1个csv文件和1个csv文件。这是脚本的两部分之间唯一的区别 ##Samplesheet for GS0005701 rows<-unique(samples$Sample_Name) samplesheet<-rows opa.panels<-sort(unique(samples$Pool_ID)) for ( i in 1:length(opa.panels)){ samps&

我正在编写一个脚本来创建csv文件进行分析。运行脚本时,OPA 5701和6561分别提供1个csv文件和1个csv文件。这是脚本的两部分之间唯一的区别

##Samplesheet for GS0005701
rows<-unique(samples$Sample_Name)
samplesheet<-rows
opa.panels<-sort(unique(samples$Pool_ID))
for ( i in 1:length(opa.panels)){
  samps<-samples[samples$Pool_ID == opa.panels[i],]
  idx<-match(samps$Sample_Name,rows)
  samplesheet<-cbind(samplesheet,samps$Sentrix_ID[idx],samps$Sentrix_Position[idx])
}
colnames(samplesheet)[2:(length(opa.panels)*2+1)]<-c("SentrixBarcode_A","SentrixPosition_A","SentrixBarcode_B","SentrixPosition_B","SentrixBarcode_C","SentrixPosition_C","SentrixBarcode_D","SentrixPosition_D")[1:(length(opa.panels)*2)]
colnames(samplesheet)[1]<-"Sample_Name"
idx<-match(rows,samples$Sample_Name)
samplesheet<-cbind(samplesheet,samples[idx,c("Sample_Group","NorTum","Sample")])
ss_header<-c("[Header]","Investigator Name,Sander","Project Name,HNPCC_NA_MYH","Experiment Name,OPA1+2+3+4","Date,5062012","[Manifests]")
for (i in 1:length(opa.panels)) ss_header<-c(ss_header,paste(LETTERS[i],opa.panels[i],sep=","))
ss_header<-c(ss_header,"[Data]")
writeLines(ss_header,"Samplesheet5701.csv")
write.table(samplesheet,file="Samplesheet5701.csv",sep=",",row.names=FALSE,quote=FALSE,append=TRUE,na="")

##Samplesheet for GS0006561-OPA
rows2<-unique(samples2$Sample_Name)
samplesheet2<-rows2
opa.panels2<-sort(unique(samples2$Pool_ID))
for ( j in 1:length(opa.panels2)){
  samps2<-samples2[samples2$Pool_ID == opa.panels2[j],]
  idx2<-match(samps2$Sample_Name,rows2)
  samplesheet2<-cbind(samplesheet2,samps2$Sentrix_ID[idx2],samps2$Sentrix_Position[idx2])
}
colnames(samplesheet2)[2:(length(opa.panels)*2+1)]<-c("SentrixBarcode_A","SentrixPosition_A","SentrixBarcode_B","SentrixPosition_B","SentrixBarcode_C","SentrixPosition_C","SentrixBarcode_D","SentrixPosition_D")[1:(length(opa.panels)*2)]
colnames(samplesheet2)[1]<-"Sample_Name"
idx2<-match(rows2,samples2$Sample_Name)
samplesheet2<-cbind(samplesheet2,samples2[idx2,c("Sample_Group","NorTum","Sample")])
ss_header<-c("[Header]","Investigator Name,Sander","Project Name,HNPCC_NA_MYH","Experiment Name,OPA1+2+3+4","Date,5062012","[Manifests]")
for (j in 1:length(opa.panels2)) ss_header<-c(ss_header,paste(LETTERS[j],opa.panels2[j],sep=","))
ss_header<-c(ss_header,"[Data]")
writeLines(ss_header,"samplesheet6561.csv")
write.table(samplesheet2,file="Samplesheet6561.csv",sep=",",row.names=FALSE,quote=FALSE,append=TRUE,na="")
我知道这是一个很难回答的问题,但我希望有人能给我一个提示,说明一个代码如何创建data.frame,另一个代码创建矩阵


非常感谢

这个问题的正确答案是在索引部分进行交换

idx<-match(samps$Sample_Name,rows)

idxCan您可以用dput(head(samples))的输出替换数据的图片(以及损坏的复制和粘贴部分)这个问题太混乱了。你应该针对一个小的可复制的例子。但要回答data.frame vs matrix的问题,很可能是因为一个输出包含数字和字符值的混合,而另一个只包含数字或字符(后者更可能是给定的输入)。感谢您的回答,我已经找到了导致此问题的原因。这与
match`语句和交换
idx有关。您是否可以回答自己的问题,这样就可以找到问题的好答案,并且人们可以看到问题已经得到了回答。
idx<-match(samps$Sample_Name,rows)
idx<-match(rows,samps$Sample_Code)