Hi, Thanks, I should have thought of that. If I do tst<-(10*seq(100,125,by=.2))%in%(10*seq(0,800,by=.1)) > sum(tst) [1] 122 > tst<-(seq(100,125,by=.2))%in%(seq(0,800,by=.1)) > sum(tst) [1] 76 The problem is corrected. However in my code, where I am comparing much longer sequences it doesn't always work, even with rounding and multiplying to an integer. Is there a more robust way to do this? Here is the function where things are going wrong
spectots<-function(t,x,begin=min(t,na.rm=TRUE),end=max(t,na.rm=TRUE), by.inter=.5,precis=1){ #browser() j<-!is.na(t) k<-!((t>end) | (t<begin)) j<-j&k t<-t[j] x<-x[j] get.mult<-function(pre)return(10^pre) new.t<-seq(begin,end,by=by.inter) new.x<-rep(NA,length(new.t)) #new.n<-length(new.t) t.err<-rep(NA,length(new.t)) ######### Here is the matching ############ tst<-(get.mult(precis)*new.t)%in%(get.mult(precis)*round(t,precis)) check<-sum(tst) if(check!=length(t))warning("not all values matched!") new.x[tst]<-x t.err[tst]<-t-new.t[tst] x.impute<-ifelse(is.na(new.x),.5*min(new.x,na.rm=TRUE),new.x) ts(cbind(x.impute,new.x,t.err),start=begin,end=end,freq=1/by.inter) } On Sat, 2003-04-05 at 03:09, Peter Dalgaard BSA wrote: > Nicholas Lewin-Koh <[EMAIL PROTECTED]> writes: > > > Hi, > > Am I hitting some limit in match? Consider the following example: > > > > > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > > > sum(tst) > > [1] 76 > > > Gives the correct answer. Did I miss something? > > The fact that 1/5 and 1/10 are not represented exactly in a binary > computer? > > > sum(tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1)) > [1] 76 > > sum(tst<-seq(100,125,by=.2)%in%round(seq(0,800,by=.1),1)) > [1] 126 > > And, just to be sure: > > > sum(round(seq(100,125,by=.2),1)%in%round(seq(0,800,by=.1),1)) > [1] 126 > > > -- > O__ ---- Peter Dalgaard Blegdamsvej 3 > c/ /'_ --- Dept. of Biostatistics 2200 Cph. N > (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 > ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 > ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help