Dear All,
I would like to ask for help on how to read different files automatically and do
analysis using scripts.
1. Description of the data
1.1. there are 5 text files, each of which contains cleaned data for the same
100 SNPs. Observations (e.g., position on gnome, alelle type, ...) for SNPs
are rows ordered by the SNP numbers,
1.2. there are 1 text file, containing the expression level of mRNAs 9 (and
other information), which are rows ordered by mRNA numbers,
So for each SNP the sample size is 5.
2. Description of aim
Take SNP 1 and mRNA 1 for example. Write a scrip that can
2.1 extract row 1 from each of the 5 text files for SNP data
2.2 extract row 1 from the mRNA text file
2.3 regress mRNA 1 counts on SNP 1
2.4 store regression results, SNP number, mRNA number and their positions on
gnome
3. key thing I need help
automatic extraction of observations for one particular SNP from SNP data files,
since there are actually 400 SNP data files, or more generally extract
information automatically from different files and store it properly for
operation.
Thanks a lot!
Chee
[[alternative HTML version deleted]]
