Thanks, Erich!
On Mon, Nov 1, 2010 at 4:55 AM, Erich Neuwirth
<erich.neuwi...@univie.ac.at>wrote:
> My guess is that what you want probably is best done by using ecdf.
> You might want to look it up in the docs.
>
> 1-ecdf(test)(x)
> Will give you the percentage of values in test larger than x.
>
>
> On 11/1/2010 2:24 AM, Changbin Du wrote:
> > Thanks Joshua! Yes, i is not going up sequentially by 1, as i here is
> the
> > raw number of reads for each DNA base. Thanks so much for the great help!
> >
> >
> > On Sun, Oct 31, 2010 at 6:03 PM, Joshua Wiley <jwiley.ps...@gmail.com
> >wrote:
> >
> >> On Sun, Oct 31, 2010 at 5:44 PM, Joshua Wiley <jwiley.ps...@gmail.com>
> >> wrote:
> >> <snip>
> >>> #cover is a vector
> >>> cover_per <- function(cover) {
> >>> ## create a vector to store the results of your for loop
> >>> output <- vector("numeric", length(min(cover):max(cover)))
> >>> for (i in min(cover):max(cover)) {
> >>> ## rather than print()ing the output, assign it to an object
> >>> output[i] <- 100*sum(ifelse(cover >= i, 1, 0))/length(cover)
> >>> }
> >>> ## have the return value from the function be
> >>> ## the object 'output'
> >>> return(output)
> >>> }
> >>
> >> I did not catch that i was not necessarily starting at 1 going
> >> sequentially up, so that would have to be done manually (or use cumsum
> >> per David rather than the function you wrote).
> >>
> >> cover_per2 <- function(cover) {
> >> output <- vector("numeric", length(min(cover):max(cover)))
> >> j <- 1
> >> for (i in min(cover):max(cover)) {
> >> output[j] <- 100*sum(ifelse(cover >= i, 1, 0))/length(cover)
> >> j <- j + 1
> >> }
> >> return(output)
> >> }
> >>
> >> Josh
> >>
> >
> >
> >
>
>
--
Sincerely,
Changbin
--
Changbin Du
DOE Joint Genome Institute
Bldg 400 Rm 457
2800 Mitchell Dr
Walnut Creet, CA 94598
Phone: 925-927-2856
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
