Dear all,
I am really sorry for not giving the input file because in my mail,I did not
explain my problem in a best way.
I have a file that is summary.txt(I have attached it) .we can read this file
using-
dfa=read.table("summar.txt",fill=T,colClasses = "character",header=T)
In V10 column I have ASCII values which I converted into decimal numbers
using this code-
dfa$V10=lapply(dfa[,4], function(c) as.numeric(charToRaw(c)))
Now I have a dataframe dfa with these columns something like this-
V7 V8 V9 V10
0 1 G 82
0 1 CGT c(90, 92,
96)
0 1 GA c(78, 92)
0 1 GAG c(90, 92,
92)
0 1 G 88
0 1 A 96
0 1 ATT c(90, 96,
92)
0 1 T 94
0 1 C 97
the values in column V10 corresponds to A,C,G T in column V9.I want only those
whose score is more than 91.so output of above should be-
V7 V8 V9 V10
0 1 GT c(90, 92, 96)
0 1 A c(78, 92)
0 1 AG c(90, 92, 92)
0 1 A 96
0 1 TT c(90, 96, 92)
0 1 T 94
0 1 C 97
Can you please tell me the solution.
Thanking you,
Warm Regards
Vikas Bansal
Msc Bioinformatics
Kings College London V7 V8 V9 V10
0 1 G `
0 1 T a
0 1 C a
0 1 A a
0 1 G _
0 1 G Z
0 1 C ^
0 1 C \
0 1 A Z
0 1 T a
0 1 g ^
0 1 A \
0 1 C _
0 1 G a
0 1 C X
0 1 C `
0 1 G a
0 1 G _
0 1 G a
0 1 G a
0 1 A a
0 1 G a
0 1 G a
0 1 C a
0 1 A a
0 1 C _
0 1 A X
0 1 g `
0 2 GG \\
0 2 GG [^
0 2 AA Y^
0 2 GG `]
0 2 AA a^
0 2 GG ]^
0 2 AA a\
0 2 GG a]
0 2 GG Z]
0 2 GG ]^
0 2 CC W\
0 2 CC a]
0 2 TT ``
0 2 GG a\
0 2 GG ``
0 2 aa aa
0 2 AA a^
0 2 CC b`
0 2 AA _\
0 2 CC ]`
0 2 TT ^\
0 2 CC Z`
0 2 Ac `a
0 3 AAA b`a
0 3 GGG aa]
0 3 AAA `[_
0 3 CCC a`_
0 3 TTT _]^
0 3 CCC aaa
0 3 CCC ^a`
0 3 CCC _``
0 3 AAA Z`]
0 3 CCC \a]
0 2 GG `_
0 2 GG `Y
0 2 AA a]
0 1 G Z
0 1 G _
0 1 T ^
0 1 T ^
0 1 C Y
0 1 A \
0 1 G ]
0 1 T [
0 1 T ^
0 1 C [
0 1 C \
0 1 A ^
0 1 C ^
0 1 A ^
0 1 C Y
0 1 T ^
0 1 C Y
0 1 C \
0 1 C Q
0 1 C X
0 1 T ^
0 1 C Z
0 1 C X
0 1 c a
0 1 a a
0 1 c a
0 1 c ^
0 1 t a
0 1 t _
0 1 c a
0 1 A a
0 1 g a
0 1 a a
0 1 c a
0 1 A a
0 1 G a
0 1 C a
0 1 A a
0 1 C a
0 1 C _
0 1 C a
0 1 T `
0 1 C `
0 1 C a
0 1 C ^
0 1 C `
0 1 G _
______________________________________________
R-help@r-project.org mailing list
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.
