Hi, Jorge's method will be faster. #system.time(res1<-sum(apply(mat2,1,function(x) x[1]<12 & x[2]>12))) #instead of 2, it should be 1 ?#? user? system elapsed ?# 0.440?? 0.000?? 0.445 ?system.time(res1<-sum(apply(mat2,1,function(x) x[1]<=12 & x[2]>12))) # ?#? user? system elapsed ?# 0.500?? 0.000?? 0.502 ?res1 #[1] 80070 A.K.
From: Jim Silverton <jim.silverton at gmail.com>
To: arun <smartpink111 at yahoo.com>
Sent: Monday, March 18, 2013 10:08 AM
Subject: Re: [R] Counting confidence intervals
To: arun <smartpink111 at yahoo.com>
Sent: Monday, March 18, 2013 10:08 AM
Subject: Re: [R] Counting confidence intervals
thanks arun!! On Mon, Mar 18, 2013 at 10:06 AM, arun <smartpink111 at yahoo.com> wrote: Hi, >Try this: >set.seed(25) >mat1<- matrix(cbind(sample(1:15,20,replace=TRUE),sample(16:30,20,replace=TRUE)),ncol=2) >?nrow(mat1[sapply(seq_len(nrow(mat1)),function(i) any(seq(mat1[i,1],mat1[i,2])==12)),]) >#[1] 17 > > >set.seed(25) >mat2<- matrix(cbind(sample(1:15,1e5,replace=TRUE),sample(16:30,1e5,replace=TRUE)),ncol=2) > >system.time(res<-nrow(mat2[sapply(seq_len(nrow(mat2)),function(i) any(seq(mat2[i,1],mat2[i,2])==12)),])) >?#? user? system elapsed >?# 1.552?? 0.000?? 1.549 >res >#[1] 80070 >?head(mat2[sapply(seq_len(nrow(mat2)),function(i) any(seq(mat2[i,1],mat2[i,2])==12)),]) >#???? [,1] [,2] >#[1,]??? 7?? 29 >#[2,]?? 11?? 30 >#[3,]??? 3?? 30 >#[4,]??? 2?? 26 >#[5,]?? 10?? 22 >#[6,]??? 6?? 22 >A.K. > > > > > > >________________________________ >From: Jim Silverton <jim.silverton at gmail.com> >To: r-help at r-project.org >Sent: Monday, March 18, 2013 9:03 AM >Subject: Re: [R] Counting confidence intervals > > >Hi, >I have a 2 x 10000 matrix of confidence intervals. The first column is the >lower and the next column is the upper. I want to cont how many times a >number say 12 lies in the interval. Can anyone assist? > >-- >Thanks, >Jim. > > >??? [[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. >
Thanks, Jim.