

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


