a quick look at it shows you would be trying to access y[n+1] in the
last part of that loop and that is greater than the number of entries
in 'y' so you will get an NA and this is not legal for comparisons.
On Feb 9, 2008 6:07 PM, mohamed nur anisah <nuranisah_mohamed at
yahoo.com> wrote:> Dear lists,
>
> i want to find the non-overlapping interval values with this code:
>
> mysetdiff=function(x,y){
> m=length(x)
> n=length(y)
> bx = logical(m)
> by = logical(n)
> for(i in 1:m){
> for(j in 1:n){
> if(x[i]<=y[j+1]){
> bx[i] = T
> by[j] = T
> NA = NA
> }
> }
> }
> sx = x[!bx]
> sy = y[!by]
> s=c(sx,sy)
> return(s)
> }
>
> Below is my dataset. When i called back my function with the
code;>mysetdiff(f,e). An error had occur: Error in if (x[i] <= y[j + 1]) {
: missing value where TRUE/FALSE needed. How am i going to fix my function so
that i can get the values of my non-overlapping interval. Any suggestion??
Thanks a bunch!!
>
> e
> [1] 17130612 17712302 21225764 25012714 33852816 36012944 36252300
> [8] 36737468 43693832 44148616 45318876 45852632 53258208 58530988
> [15] 60437872 72516480 79673224 93128744 94269896 95868704 99651504
> [22] 113688560 131101008 132955984 135891280 141318144 148257888 156158176
> [29] 157797616 162055856 168221296 173125232 176267104 182826240 183742528
> [36] 184401728 190671888 196639616 17587118 18221688 21387314 30748348
> [43] 34480192 36209144 36280276 36971144 43878548 44496056 45740012
> [50] 46752088 53700056 58603536 60691012 72757696 80077728 93181480
> [57] 94474624 97418088 106596368 120128352 132462320 132980744 135998880
> [64] 142259520 151591840 156920960 157838176 162743136 168466848 173167936
> [71] 176338384 182930096 184149776 185735712 190910576
> > f
> [1] 17712302 21203780 25012714 33852816 34794536 36012944 37891284
> [8] 43693832 44148616 45852632 53289188 61573112 63664928 72516480
> [15] 79673224 94474624 95868704 99651504 113688560 125159688 127388568
> [22] 131101008 154599216 176267104 181504912 182562720 182826240 183742528
> [29] 196841904 18100404 21387314 30748348 34384588 35996440 36252300
> [36] 37942556 43878548 44496056 46752088 53700056 62637560 63969952
> [43] 72757696 80077728 94617360 97144032 106596368 120128352 127220456
> [50] 127504536 132462320 154717312 176338384 181836032 182687824 182930096
> [57] 184149776
>
>
>
>
> ---------------------------------
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem you are trying to solve?