Hello, If I understand correctly, just use ?cbind.
Rui Barradas Em 02-02-2017 13:33, greg holly escreveu:
Hi Rui; Is there any way to insert the chr ids in numeric as 1,2......,22 in the final output. Here is output from str(temp). So I need also chr ids in a column. 1 rs58108140 10583 G A -0.070438 0.059903 2 rs189107123 10611 C G -0.044916 0.085853 Regards, Greg On Wed, Feb 1, 2017 at 1:32 PM, Rui Barradas <ruipbarra...@sapo.pt <mailto:ruipbarra...@sapo.pt>> wrote: Hello, If what you want is to combine the files into one data.frame then there are 2 things you should see: 1) You create a variable named 'temp' and don't ever use it. 2) You never combine the data.frames you read in. Try instead the following. temp <- data.frame() for(i in 1:22) { infile<-paste("chr",i,"/Z-score.imputed",sep="") psT<-read.table(infile,header=T,as.is <http://as.is>=T,sep="\t") temp <- rbind(temp, psT) } str(temp) # to see what you have Hope this helps, Rui Barradas Em 01-02-2017 17:25, greg holly escreveu: Hi all; I have 22 directories named chr1, chr2,....,chr22. Each directory has a file named "Z-score.imputed". I would like to combine Z-score.imputed files into one. I wrote the following loop but did not get any results. Your helps are highly appreciated. regards, Greg temp<-c() for(i in 1:22) { infile<-paste("chr",i,"/Z-score.imputed",sep="") psT<-read.table(as.character(infile),header=T,as.is <http://as.is>=T,sep="\t") ps<-psT[psT$Var>0.6,] ratio=nrow(ps)/nrow(psT) print(ratio) } [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org <mailto:R-help@r-project.org> mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help <https://stat.ethz.ch/mailman/listinfo/r-help> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html <http://www.R-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code.
______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
