Tuesday, February 4, 2014

Fulbright meeting, Paris

 Today and tomorrow are the midyear meeting in Paris for France Fulbright and Chateaubriand Scholars. The train left out of Gare St. Jean (Bordeaux) at 8:23, exactly as the sun was rising. As we went out of town to the north, the sun was lighting up the towers over the bridge beside the CAP Sciences Museum. This side of town I don’t know at all (we live on the west side, in Caudéran, and my commute to work takes me over to Pessac to catch the train, so I have missed north Bordeaux altogether so far). Like the rest of Bordeaux, this side grades into vineyards pretty quickly; there’s a chateau now off to the east, active, another we’ve just passed that appears to be an individual’s home, no longer active in grape-growing or wine-making. Fourteen minutes out and we’ve just passed another village, all I could see of which is the cathedral surrounded by its cemetery. There are horses and little round haystacks in the sodden fields, and the ditches are brilliant with duckweed or something else bright green. The train is full of people going about their business, one knitting, another reading the paper. I asked my son David two nights ago what his favorite part of the trip has been so far, from the minute we left Downers Grove. He said it was the flight over. He didn’t know why. It certainly wasn’t for lack of enjoying the rest of the trip: he’s been happy, even exuberant, over so much on this trip. But I understand. The time between destinations transcends what we do at either end. Neither here nor there, a time when plans are evolving, interruptions are minimal, you are with family or strangers or by yourself, and you have time to think and watch out the window, and nothing else you can do. “Old men ought to be explorers / Here and there does not matter / We must be still and still moving” (T.S. Eliot, East Coker). When you travel, constraints and desires come into a different kind of balance.

Last time we took this trip our trip from Paris to Bordeaux to begin this project. The trees were all full of round masses of bright green stuff, which we took to be some kind of squirrel’s nest. Anyone reading this who knows better will recognize immediately that this was in fact mistletoe. It’s all over here, today perhaps a bit greener than it was at the end of December. The fields are much wetter. We’ve had loads of rain since we’ve been here. I wonder if the horses get colds from standing in it all day?

Back to work. I’m trying to finish up documenting the codebase for RADami today. That has been a ridiculously slow process with little to write about, but I think I can finish it up on the train and the manuscript revised and resubmitted this week, so I can get on with new RAD work next week. One realization: there may be a bias in the partitioned RAD visualization the way I’ve been doing it, because the expectation is linear with a slope of 1.0 when (1) the tree likelihood is calculated as the sum of locus log-likelihoods for all loci based on just the loci used for the partitioned analysis, rather than the global likelihood for all loci; and (2) locus log-likelihoods are assigned based on topological identity with pruned trees that are voted on by each locus. I’m not sure this introduces a bias. It should reduce noise, as there is noise associated with requiring a locus to vote on trees that are topologically identical to one another when pruned down to just those taxa that are in the locus, but not topologically identical when all taxa are included. Because the optimization runs only until it stops improving by epsilon, trees that all lie at one point in the likelihood surface may appear to be at different positions. So in this sense the visualization as written biases us toward finding a tighter fit plot; but does it bias us toward a linear relationship? I’ll set this up to run the global likelihood on each tree as well, and see whether that plot is also so nice.

No comments:

Post a Comment