[Phylobase-devl] characters not being read
François Michonneau
francois.michonneau at gmail.com
Thu Aug 5 15:17:56 CEST 2010
Hi again Conrad,
In the file you sent me there were no state labels (just
charstatelabels), so the warning is normal.
For the underscore issue, there isn't much we can do about. NCL
converts underscores to spaces because there is no way to
differentiate them in newick strings for instance. Then when the name
is assigned to the column name in R, the space is converted to a dot.
Cheers,
-- François
On Wed, Aug 4, 2010 at 22:31, Conrad Stack <stack at psu.edu> wrote:
> Looks great. Two things: it looks like the warning message about 'state
> labels' fires even when labels are found (only when DATATYPE=STANDARD) and,
> perhaps relatedly, underscores in CHARSTATELABELS are automatically
> converted to periods.
> Thanks!
> -Conrad
>
>
>
> On Wed, Aug 4, 2010 at 12:53 PM, François Michonneau
> <francois.michonneau at gmail.com> wrote:
>>
>> Hi,
>>
>> My last commit should fix the issue Conrad mentioned.
>>
>> The problem came from the fact that the nexus file didn't have state
>> labels and the default option in readNexus was to return them. It
>> ended up attaching empty labels to the the 0 and 1 which resulted in
>> the empty column observed. I modified readNexus to (1) return a
>> warning when state labels are missing and return.labels is TRUE which
>> leads to ignoring the return.labels option.
>>
>> Let me know if you find any other bugs.
>>
>> Cheers,
>> -- François
>>
>> On Tue, Aug 3, 2010 at 17:42, Conrad Stack <stack at psu.edu> wrote:
>> > I just upgraded to phylobase 0.6.1 and found that the data in my
>> > CHARACTERS
>> > blocks are no longer being read in (they show up as blanks in the
>> > returned
>> > phylo4d object).
>> > this is from the geospiza.nex dataset:
>> > BEGIN CHARACTERS;
>> > DIMENSIONS NCHAR=1;
>> > FORMAT DATATYPE = STANDARD GAP = - MISSING = ? SYMBOLS = " 0 1";
>> > MATRIX
>> > fuliginosa 1
>> > fortis 1
>> > magnirostris 0
>> > conirostris 0
>> > scandens 0
>> > difficilis 1
>> > pallida 0
>> > parvulus 0
>> > psittacula 0
>> > pauper 0
>> > Platyspiza 0
>> > fusca 0
>> > Pinaroloxias 0
>> > ;
>> > END;
>> > When it reads in, it shows up like this:
>> > readNexus("geospiza.nex")
>> > label node ancestor edge.length node.type standard_char
>> > 1 fuliginosa 1 22 0.05500 tip
>> > 2 fortis 2 22 0.05500 tip
>> > 3 magnirostris 3 21 0.11000 tip
>> > 4 conirostris 4 20 0.18333 tip
>> > 5 scandens 5 19 0.19250 tip
>> > 6 difficilis 6 18 0.22800 tip
>> > 7 pallida 7 23 0.08667 tip
>> > 8 parvulus 8 25 0.02000 tip
>> > 9 psittacula 9 25 0.02000 tip
>> > 10 pauper 10 24 0.03500 tip
>> > 11 Platyspiza 11 16 0.46550 tip
>> > 12 fusca 12 15 0.53409 tip
>> > 13 Pinaroloxias 13 14 0.58333 tip
>> > 14 <NA> 14 0 0.29744 root <NA>
>> > 15 <NA> 15 14 0.04924 internal <NA>
>> > 16 <NA> 16 15 0.06859 internal <NA>
>> > .......cont....
>> > Any thoughts? (I'm using windows 7 32-bit with R 2.11, but the error
>> > occurs
>> > on my Mac leopard install as well.
>> >
>> >
>> > --
>> > Conrad Stack
>> > -----------------------
>> > PSU Department of Biology
>> > 208 Mueller Lab
>> > University Park, PA 16802
>> > cell: 814.409.8310
>> > email: conrad.stack at gmail.com
>> >
>> > _______________________________________________
>> > Phylobase-devl mailing list
>> > Phylobase-devl at lists.r-forge.r-project.org
>> >
>> > https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/phylobase-devl
>> >
>> >
>>
>
>
>
> --
> Conrad Stack
> -----------------------
> PSU Department of Biology
> 208 Mueller Lab
> University Park, PA 16802
> cell: 814.409.8310
> email: conrad.stack at gmail.com
>
More information about the Phylobase-devl
mailing list