## User: Harris A. Jaffee

Reputation:
590
Status:
Trusted
Location:
United States
Last seen:
3 years, 10 months ago
Joined:
8 years, 6 months ago
Email:
h*@jhu.edu

#### Posts by Harris A. Jaffee

<prev • 58 results • page 1 of 6 • next >
1
319
views
1
Answer: A: mismatch nb in Biostrings
... First, see ?lowlevel-matching, and if you are really adventurous, the function nedit_at() in Biostrings/src/lowlevel_matching.c. The degree to which your pattern and subjects disagree, in a sense dictated by your with.indels=FALSE and fixed=TRUE arguments, is determined first, although actually o ...
written 4.2 years ago by Harris A. Jaffee590
2
394
views
2
... If X and Y are your data frames, the rows you want are given by library(IRanges) nearest(IRanges(X[,"Position"], width=1), IRanges(start=Y[,"Start"], end=Y[,"Stop"])) On May 6, 2014, at 3:34 PM, Elena Sorokin wrote: > I have two data frames as input: > >> small Pos ...
written 4.4 years ago by Harris A. Jaffee590
2
1.8k
views
2
... I would like to know the answer! This is all I have. > width(tr_by_gene) IntegerList of length 25529 [["1"]] 6694 [["10"]] 9969 [["100"]] 32214 [["1000"]] 226516 [["10000"]] 355050 343564 343866 355040 343554 343856 [["100008586"]] 187588 6118 6109 [["100008587"]] 156 156 156 156 156 156 156 15 ...
written 4.5 years ago by Harris A. Jaffee590
1
2.5k
views
1
... Thanks, Val, for the details. Now I have what is either a shameless plug or a request for explanation. matchGenes() in _bumphunter_ is similar, although definitely not the same. I am not sure how to translate between the results, except to note that matchGenes uses nearest() with either select="a ...
written 4.7 years ago by Harris A. Jaffee590
3
874
views
3
Comment: C: rhdf5 segfault
... On Dec 23, 2013, at 10:56 PM, Ben Bolstad wrote: > You might try changing your stack size limits in your shell, before launching R. e.g. > > ulimit -s unlimited I don?t think that is possible, at least as you put it. The soft limit tends to be smaller than the hard limit and is the onl ...
written 4.8 years ago by Harris A. Jaffee590
3
544
views
3
... I needed essentially your "last_hits" for bumphunter:::fuzzy.match2 In your notation, it comes out: last_hits_rows <- cumsum(rle(queryHits(all_hits))\$lengths) last_hits <- all_hits[last_hits_rows,] Having expended that effort, you could do first_hits_rows <- c(1, ...
written 4.9 years ago by Harris A. Jaffee590
1
921
views
1
... I believe gzfile only applies when reading. With gzcon, you may have to flush or close the connection (not sure which one - I think your fout2), or your whole session, to get the file written. On Sep 25, 2013, at 4:07 PM, Marcus Davy wrote: > Thanks for the responses, Martin I will give your ...
written 5.0 years ago by Harris A. Jaffee590
1
776
views
1
... Thank you again for the very precise bug report and your data (off list). The annotation handlers for the 10000-nearest-transcript chunks were not being sent the user's value of 'all', so it defaulted to FALSE. Having fixed that, I found out that these handlers also needed a 'queries' and 'Tx' argu ...
written 5.2 years ago by Harris A. Jaffee590
1
776
views
1
... Anne, Can you run it with 'verbose=TRUE'? Also, tell us what your mc.cores value is (This will also appear in the verbose output.)? If it's >1, please try 'mc.cores=1' to see if that makes any difference. On Jul 26, 2013, at 9:45 AM, Kasper Daniel Hansen wrote: > This sounds weird and it ...
written 5.2 years ago by Harris A. Jaffee590
6
1.3k
views
6
... On Jul 24, 2013, at 12:26 PM, Taylor, Sean D wrote: > Hi Harris, > > No you are right, I will mostly do 3' trimming. But for my test data it just seemed more convenient to test using Lpattern and I assumed (correctly I hope) that all the same things apply to Rpattern so long as you keep y ...
written 5.2 years ago by Harris A. Jaffee590

#### Latest awards to Harris A. Jaffee

No awards yet. Soon to come :-)

Content
Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.