Hi all, After running this code (attaches is the input file): dat <- read.table("gene_prob.txt", sep = "\t") n <- length(dat$V1) print(n) print(dat$V1) I get this print out. ...... [8541] LOC552889 GPR15 SLC2A11 GRIP2 SGEF [8546] PIK3IP1 RPS27 AQP7 8548 Levels: 3.8-1 A2M A4GALT A4GNT AAAS AAK1 AAMP AANAT AARSD1 AASS ... hCG_1730474 What's the meaning of the last line? Is it an error? -- Gundala Viswanath -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: gene_prob.txt URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080529/99090d93/attachment.txt>
Hi all, After running this code: __BEGIN__ dat <- read.table("gene_prob.txt", sep = "\t") n <- length(dat$V1) print(n) print(dat$V1) __END__ With this input in "gene_prob.txt" __INPUT__ HFE 0.00107517988586552 NF1 0.000744355305599206 PML 0.000661649160532628 TCF3 0.000661649160532628 NF2 0.000578943015466049 GNAS 0.000578943015466049 GGA2 0.000578943015466049 ..... I get this print out. ...... [8541] LOC552889 GPR15 SLC2A11 GRIP2 SGEF [8546] PIK3IP1 RPS27 AQP7 8548 Levels: 3.8-1 A2M A4GALT A4GNT AAAS AAK1 AAMP AANAT AARSD1 AASS ... hCG_1730474 What's the meaning of the last line? Is it an error? How can I fix it? -- Gundala Viswanath
Gundala Viswanath <gundalav <at> gmail.com> writes:> > Hi all, > > After running this code (attaches is the input file): > > dat <- read.table("gene_prob.txt", sep = "\t") > n <- length(dat$V1) > print(n) > print(dat$V1) > > I get this print out. > > ...... > [8541] LOC552889 GPR15 SLC2A11 GRIP2 SGEF > [8546] PIK3IP1 RPS27 AQP7 > 8548 Levels: 3.8-1 A2M A4GALT A4GNT AAAS AAK1 AAMP AANAT AARSD1 AASS > ... hCG_1730474 > > What's the meaning of the last line? Is it an error? >Why should this be an error? I found hCG_173047 in your data set. For problem like this, try to locate the culprit be reducing the data set: first take the first half, then, if no error, the second, recursively. And please, don't post these huge sets again before having reduced the problem to the minimum. Dieter