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]]
______________________________________________
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.
