Displaying 2 results from an estimated 2 matches for "checkcoverage".
2008 Aug 02
3
convert for loop into apply()
Hi all,I know this topic has came up multiple times, but I've never fully
understand the apply() function.
Anyway, I'm here asking for your help again to convert this loop to apply().
I have 2 data frames with the following information: a1 is the fragment that
is need to be covered, a2 is the probes that cover the specific fragment.
I need to count the number of probes cover every given
2008 Aug 05
0
P values in non linear regression and singular gradients using nls
...step if a1 has 30000 rows and a2 has 200000,
> unless you memory size is approximately infinite, so you will need a
> loop. Suppose you can handle 1000 probes at a time. You might be
able
> to get away with something like this:
> __________________
>
[[elided Yahoo spam]]
>
> checkCoverage <- function(a1, a2)
> colSums(outer(a2$st , a1$en , "<=") *
> outer(a2$en , a1$st , ">=") *
> outer(a2$cat, a1$cat, "=="))
>
> coverage <- numeric(N <- nrow(a2))
> m2 <- 0
> while(m2 < N) {
> m1 <- m2 + 1
>...