Displaying 1 result from an estimated 1 matches for "fooqtil".
Did you mean:
fooqtile
2007 Aug 21
1
quantile() returns a value outside the data range
...qual' for "<=")?
Small example:
> foo <- h2[,"17"][h2$plate==222] # some data I have, details not
interesting
> foo <- sort(foo)[1:3] # can demo with the 3 smallest
values
> yo <- data.frame(matrix(foo,nrow=length(foo),ncol=10))
> fooQtile <- rep(0,9)
> for(i in 1:9) { # compute 0.01 quantile, for all
'types'
+ fooQtile[i] <- quantile(foo,probs=0.01,type=i)
+ yo[,i+1] <- foo <= fooQtile[i]
+ }
> names(yo) <- c("myData",paste("qType",1:9,sep=""))
> y...