hi there :)
i want to use barplot with if else but i dont know how to do it ?
i tried this but it is not working with me
SNP <- read.table("my.txt")
>SNP[,2]
[1] 1175 483 240 170 99 79 76 45 38 35 21 16 14 19 16
[16] 3 3 3 10 2 1 6 8 6 8 2 0 5 1 1
[31] 1 0 6 2 0 13 0 5 0 5 0 1 0 0 0
[46] 0 5 5 0 1 0 0 0 0 0 0 0 0 0 0
[61] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 7
[76] 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0
[91] 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0
[106] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[121] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[136] 1 0 0 0 1 1 0 0 0 0 0 0 0 0 0
[151] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[166] 0 0 2 0 0 0 0 0 0 1 0 0 0 0 0
[181] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
[196] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[211] 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0
[226] 1
>SNP[,2] == 0
[1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[13] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[25] FALSE FALSE TRUE FALSE FALSE FALSE FALSE TRUE FALSE FALSE TRUE FALSE
[37] TRUE FALSE TRUE FALSE TRUE FALSE TRUE TRUE TRUE TRUE FALSE FALSE
[49] TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[61] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[73] TRUE TRUE FALSE TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE
[85] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[97] TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE
[109] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[121] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[133] TRUE TRUE TRUE FALSE TRUE TRUE TRUE FALSE FALSE TRUE TRUE TRUE
[145] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[157] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE
[169] TRUE TRUE TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE
[181] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[193] TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[205] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[217] TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE FALSE
> SNP[,2][SNP[,2] == 0]
[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[38] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[75] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[112] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[149] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
> SNP[,2][! SNP[,2] == 0]
[1] 1175 483 240 170 99 79 76 45 38 35 21 16 14 19 16
[16] 3 3 3 10 2 1 6 8 6 8 2 5 1 1 1
[31] 6 2 13 5 5 1 5 5 1 7 3 3 1 1 1
[46] 2 1 1 1
> SNP[,1]
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
[19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
[37] 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
[55] 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
[73] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
[91] 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
[109] 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126
[127] 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
[145] 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162
[163] 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180
[181] 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198
[199] 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216
[217] 217 218 219 220 221 222 223 224 225 226
> zz <- SNP[,2][! SNP[,2] == 0]
> zz
[1] 1175 483 240 170 99 79 76 45 38 35 21 16 14 19 16
[16] 3 3 3 10 2 1 6 8 6 8 2 5 1 1 1
[31] 6 2 13 5 5 1 5 5 1 7 3 3 1 1 1
[46] 2 1 1 1 1
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
[20] 20 21 22 23 24 25 26 28 29 30 31 33 34 36 38 40 42 47 48
[39] 50 75 80 104 136 140 141 168 175 195 219 226
> xx <- SNP[,1][! SNP[,2] == 0]
>xx
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
[20] 20 21 22 23 24 25 26 28 29 30 31 33 34 36 38 40 42 47 48
[39] 50 75 80 104 136 140 141 168 175 195 219 226
xlab="No. of Structure"
ylab="No. of Mutation Mapped to"
title="SNPs Structures"
barplot(zz,width = 4,names.arg= xx,axes = TRUE, axisnames = TRUE,
main=title,xlab=xlab,ylab=ylab,ylim=c(0,1175),xlim=c(0,226),col = c(for (i in
zz){if i>70 col= "lightblue" else col= "mistyrose"))
i tried to get red of zero then plot it
i want to color different result with different colors ?
any body can help pleeeeaaaaaaasssseee
thanks
Nouf
_________________________________________________________________
[[alternative HTML version deleted]]
On Thu, 28 May 2009 15:47:18 +0300 NOUF AL NUMAIR
<noufalnumair at hotmail.com> wrote:
NAN> barplot(zz,width = 4,names.arg= xx,axes = TRUE, axisnames = TRUE,
NAN> main=title,xlab=xlab,ylab=ylab,ylim=c(0,1175),xlim=c(0,226),col NAN>
c(for (i in zz){if i>70 col= "lightblue" else col=
"mistyrose"))
NAN>
NAN>
NAN> i tried to get red of zero then plot it
? I have no idea what this means...
NAN> i want to color different result with different colors ?
NAN>
NAN> any body can help pleeeeaaaaaaasssseee
If you want help try to be more clear about what you want. It is hard
with your mail because there is so much what does not belong to your
current problem, one has to scroll over plenty of numbers and still I am
not sure whether I understood what you want.
Hence you want your bar color conditional on the value or the height,
it is helpful to create a vector beforehand.(for control purposes)
y<-c(1175, 483, 240, 170, 99, 79, 76, 45, 38, 35,
21, 16, 14, 19, 16)
x<-c( 1, 2, 3, 4, 5, 6, 7,
8, 9, 10, 11, 12, 13, 14, 15)
d1<-data.frame(x=x,y=y)
Color<-ifelse(d1$y<70,"blue","red")
barplot(as.matrix(d1$y),names.arg=d1$x,col=Color,beside=T)
as an example.
hth
Stefan