User: Richard Pearson
Richard Pearson • 60
- Reputation:
- 60
- Status:
- New User
- Location:
- Last seen:
- 6 years, 7 months ago
- Joined:
- 7 years, 7 months ago
- Email:
- r*******@well.ox.ac.uk
Posts by Richard Pearson
<prev
• 6 results •
page 1 of 1 •
next >
0
votes
1
answer
708
views
1
answers
Comment:
C: writeVcf bug
... Hi Valerie
That's great, thanks. Will check out the new code when available.
Sorry, I worded this poorly. I simply meant that the problem I had
previously seen with the example code ("Error in row(genoMat)") had
gone
away, but I had noticed a new problem with this same code. There are
no
other pro ...
written 6.7 years ago by
Richard Pearson • 60
0
votes
1
answer
708
views
1
answer
... Hi
I've recently upgraded to R 3.0.0, but am still having a (now slightly
different) problem with writeVcf. When there is only one element in
geno(vcf), each variant gets a separate row for each sample. For
example, in the following code, out1.vcf has 5 variants, each with
genotypes for 3 samples, ...
written 6.7 years ago by
Richard Pearson • 60
• updated
6.7 years ago by
Valerie Obenchain • 6.7k
0
votes
1
answer
694
views
1
answers
Comment:
C: writeVcf bug
... Hi Michael
My devel sessionInfo is below. It seems biocLite has given me version
1.5.17of the package, but I see that .makeVcfGenolooks different in
the
latest from SVN (1.5.19) - so I'm sure you're right and I just need to
wait for things to filter through. Apologies for the noise.
Every time I u ...
written 7.0 years ago by
Richard Pearson • 60
0
votes
1
answer
694
views
1
answer
... Hi
I've noticed a bug with writeVcf when trying to write out relatively
simple vcf files, e.g. a file containing only GT in the FORMAT fields.
The following is a (hopefully) reproducible example based on examples
in
?writeVcf
library(VariantAnnotation)
fl <- system.file("extdata", "ex2.vcf", pa ...
written 7.0 years ago by
Richard Pearson • 60
• updated
7.0 years ago by
Michael Lawrence ♦ 11k
0
votes
1
answer
801
views
1
answer
... Hi
It is wonderful that we can create subsets of VariantAnnotation VCF
objects using [ but I have found that this doesn't work for VCFs that
are "sites only", i.e. have no information in geno(vcf):
> geno(vcf)
SimpleList of length 0
> passVcf <- vcf[values(fixed(vcf))[["FILTER"]] ...
written 7.3 years ago by
Richard Pearson • 60
• updated
7.3 years ago by
Valerie Obenchain • 6.7k
0
votes
1
answer
533
views
1
answer
... Hi
Is it possible to read in data on a subset of samples with readVcf?
The
only way I've found of doing this is to read in data on all samples,
then subset, but it would be much more efficient for me to read in
data
only for my samples of interest.
Apologies if I've missed something obvious.
Rich ...
written 7.6 years ago by
Richard Pearson • 60
• updated
7.6 years ago by
Valerie Obenchain • 6.7k
Latest awards to Richard Pearson
No awards yet. Soon to come :-)
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 16.09
Traffic: 189 users visited in the last hour