airplane flying

I plan to sell N79PT in October.

favorite recent movies

All are on netflix and amazon. Don't look them up first. Instead, watch them first and let yourself be surprised. If need be, watch the linked trailer:
  1. The King (2019).
  2. The Two Popes (2019).
  3. Icarus (2017)
  4. AlphaGo (2017)
  5. The Handmaiden (2016).
  6. The People Vs. Fritz Bauer (2015) (Der Staat gegen Fritz Bauer)
  7. Whiplash (2014).
  8. ... The more complete list since 2008.


My thoughts on What Makes a Great Movie? Obviously, I am an amateur — I wish I weren't.


Buy Adrienne Howley's The Naked Buddha. It's a guide to living, not a guide to theology. There is (almost) no religion in there. For more information, also see the Dalai Lama — and some of his advice on dying.


national pride

Student: Auf was kann man als Deutscher noch stolz sein? Fritz Bauer: Sehen Sie, ich glaube wir Deutschen können nicht stolz sein auf unsere Berge oder die Wälder, denn da können wir nichts für. Wir können auch nicht auf Goethe und Schiller stolz sein und auch nicht auf Einstein, denn für all das können wir ja nichts. Stolz können wir doch nur auf das Gute sein, das wir selbst gemacht haben.

(Well worth translating and contemplating, especially for Americans proud of or embarrased by this or that about their country.)

Here is Fritz Bauer's 3-minute televised Keller discussion from the 1950s.

free speech

Als die Nazis die Kommunisten holten, habe ich geschwiegen; ich war ja kein Kommunist.
Als sie die Sozialdemokraten einsperrten, habe ich geschwiegen; ich war ja kein Sozialdemokrat.
Als sie die Gewerkschafter holten, habe ich nicht protestiert; ich war ja kein Gewerkschafter.
Als sie die Juden holten, habe ich geschwiegen; ich war ja kein Jude.
Als sie mich holten, gab es keinen mehr, der protestierte. --- Martin Niemöller

bds and israel

As an expert, I often opine on whether boycotts and divestments have much chance of success. (A: Typically no—and there are often many more effective mechanisms.) I usually do not opine about the merits of causes. But BDS (Boycott, Divestment, Sanctions) deserves an exception. Please read this if you are interested in BDS and the Israeli-Palestinian conflict.

For what it's worth, I detest the right-wing political coalition in Israel. I detest Abbas. Israel is not an Apartheit state, but the West-Bank is! And, sadly, demographics imply it will only get worse. There is nothing that prevents Jews and Arabs from living happily together other than the religious Verrueckte on both sides. Both sides need better leaders.

apple reviews

Frankly, I don't like Apple, although I use their products.

  • iMac 5k (and iMac Pro) Hardware, 2021: Excellent, 4 stars. Still class-leading and affordable though now aged hardware. Great displays, good performance, nice designs. The best reason nowadays to stay in the macos eco system. Slow single-core performance.
  • MacOS Software, 2022: Good, 3.5 stars. Positives: Nice GUI, some nice built-in apps. Some good software: iTerm2 (my son claims that the new Microsoft Terminal is even better), Mathematica, Skim, syncthing. Negatives: Inconsistent. Mediocre Unix underpinnings (mediocre APFS), no built-in package manager (but great brew), emacs, ls, grep, etc.)
  • Macbook Hardware, 2022: Finally pretty good. Lots of stupid ergonomics problems fixed.
  • Macbook Pro 16, M1 Max: Superfast, supercool. Very heavy schlepper, feels heavier than its objective weight. I would instead advise to buy the lighter Macbook Pro 14 or even Macbook Pro/Air 13 for travel. The best portable to get is probably the 2022-updated Macbook Air.

It says more about their lousy competitors than about Apple that I don't have much of an alternative. Unfortunately, Dell's support for ubuntu remains half-hearted. Windows has become better. I should try it again with the WSL 2. (It's the little stuff that kills you. Does it now have a good standard terminal built in? Can I now use '/' instead of '\' for directories everywhere?)

software tools

I am a heavy user of ubuntu linux, emacs, texlive, ssh, apache2/nginx, and the Google infrastructure (gmail, docs, etc.). I program primarily in perl, C, and R; occasionally in python, php, C++, mysql, javascript, and a number of other languages. If I did not know perl already, I would learn python. R is both amazing and enraging. Terrible error diagnostics. It never figured out whether it wanted to be an interactive analysis system or a programming language. I use Mathematica. This website was created in jekyll.

I liked vanilla Android (Google Nexus). I like iPads. The lower-priced iPhones are good and worth it. The high-ends are too expensive for their performance.


  • Here is a perl parser for latex called iawltxhtml (download here). It is a good start for a clean parser of latex code into useful html/epub code. It comes in a very large zip file, primarily because it contains the source and results for four example fame jagazines on which I used it (many large images!). iawltxhtml is not documented or supported, but it contains a lot of working code based on days of experimentation of how to successfully parse code into useful html, with css, math, images, flexibility as to whether the resulting file is pdf or html, etc.
  • Here is a hack for pandoc to work with multicolumn tables and to format numeric input.
  • sylspace is the software driving
  • cpan (perl) weighted least squares regression.
  • Here is a <100-line one-pass OLS implementation in C++ for two variables. biols.hh, with standard errors, univariate means and sds, etc. For large data sets, biols could become numerically unstable, but it is very fast. The usage could not be simpler. The code in test() illustrates this. Start with
             #include <biols.h>
             int main() { BIOLS bcalc; bcalc.test(); }


glycohemoglobin (A1C) by age (diabetes and prediabetes)

Glycohemoglobin (A1C) blood levels are often classified as normal below 5.7, as pre-diabetes between 5.7 and 6.5, and as full diabetes above 6.5. From what I have read, coronary diseases are the principal problem caused by high glucose.

Because my own A1C levels are classified as "pre-diabetes," I was interested in how unusual my level was for my age. Fortunately, the CDC offers population statistics on A1C in its National Health and Nutrition Examination Survey (NHANES). There is data on about 60,000 individuals with reported A1C levels. The following plot shows the unconditional distribution of A1C levels by age. Dots are individuals (possibly many), and the lines show the more informative smoothed important percentiles, by age.

plot of a1c levels by age
  • Warning: I am not a medical expert, much less a physician. I am an economist and data scientiest. Do not consider any of the material here as medical advice.
  • By the "fixed doctor cutoffs," one quarter of the population is diagnosed with pre-diabetes by age 40, half just before age 60, and two-thirds by age 70. The fact that this is so common does not negate the potential harmful effects of A1C.
  • Just before age 70, A1C reaches its highest levels among surviving participants in the study. At this age, one in five individuals has full-blown diabetes, and two more have pre-diabetes.
  • With many additional wild assumptions on the nature of these not-so-random samples and how mortality influences dropping out from this sample, we can take a wild guess at the correlation (not causation!) of A1C with life expectancy. For individuals above the 90% level (who often have other health issues, too!), it could be on the order of years. For individuals below the 90% level, it seems more modest, less than a year. But don't get carried away. There is tremendous variation in mortality within all ages regardless of A1C! This is easier to see if you click on the figure to download the a1c.pdf file, which can give you a better idea of where the 60,000 individuals are located in the plot.
  • A simple regression can further show that caucasians have A1C levels of about 0.3 lower than Mexicans and Blacks; and men have A1C levels about 0.1 higher than women. Because of the way the CDC samples, if you want to know where you are in "your" distribution, you can add 0.15 (to your own A1C values) if you are Mexican or Black and subtract 0.15 if you are Caucasian; and add 0.05 if you are male and subtract 0.05 if you are female. If you are all of the above, good for you.
  • The program that generates this plot is a1c.R. The data is in file a1c.csv.gz, which itself was based on NHANES files DEMO.XPT, DEMO_C.XPT, DEMO_E.XPT, DEMO_G.XPT, DEMO_I.XPT, GHB_E.XPT, GHB_G.XPT, GHB_I.XPT, L10_C.XPT, DEMO_B.XPT, DEMO_D.XPT, DEMO_F.XPT, DEMO_H.XPT, GHB_D.XPT, GHB_F.XPT, GHB_H.XPT, L10_B.XPT, LAB10.XPT. Thanks to Bryan Tysinger from USC to help me locate the right data.

avery labels ruler

If you ever need to photograph your moles in order to compare their changes over time, the avery ruler stickers will help you. You need to purchase Avery 58160 stickers for your laser printer, each page of which has 3 times 10 labels. Then just print the avery ruler (or avery ruler w/o copyright) stickers over them, write the location where you attach them, affix them to your skin, and snap your photograph. If you are curious, you can also download the avery-ruler.tex file. You may want to download Avery's copyrighted 58160 template file, too, but it's not necessary. Just remove the overlay pdf.

/fun/ Last modified: