I am trying to return the p value for a stat from the ECDF. That is the
index of the first occurrence,
on an ordered vector, of a value either greater than or equal to a given
value.
Ideally I would not have to order the vector beforehand.
Currently I use:
PValue<-function(stat, ECDF){
###Get the length of the ECFD
L<-length(ECDF)
###Loop through the ECDF until the p value is found
for(i in 1:L){
if(ECDF[i]>=stat){
break
}
}
###Return the 3 values that bracket the p value
c((i-1)/L, i/L, (i+1)/L)
}
Is there a way of doing this that avoids the explicit loop?
platform i386-pc-mingw32
arch i386
os mingw32
system i386, mingw32
status
major 1
minor 9.1
year 2004
month 06
day 21
language R
Phineas Campbell
______________________________________________
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
