Qatar, unbalanced but stable

Doha Airport, Qatar.

Doha Airport, Qatar.

I visited Qatar twice this month, just passing through and only visited in the airport, but there were several things that so impressed me that I had to write. What impressed me most was not so much the size and richness of the airport, but the clothes of the locals. All of the local men wore the same, very sharp robes: blindingly white, long sleeved, and floor-length. They’re called Thobes. While other nations wear something similar. Here, every one was unwrinkled, and unstained. They all looked new, with no signs they’d ever been washed. Some were worn with cuff-links (gold), and most had a pen sticking out of the breast pocket (gold). White pants peak from underneath and a headress usually sits on the head. It’s a really dramatic look, like seeing dozens of Ricardo Montaubans of Fantasy Island in one place. Local women and children were these too, but I found the thobes so dramatic that the women and children disappeared from my mind-space almost immediately. There is a local woman in the picture above, but you hardly notice.

Not everyone wears the thobes. There are lots of stores filled with gold and technology, beer and coffee, and these are maned by non-locals, Moslems mostly, almost all men. The non-locals wear western garb, not particularly sharp; none wear thobes of any sort. Some months ago, I wrote that China had severe imbalance and speculated that it was ripe for revolution. As it happens the large number of foreign worker means that Qatar is far more unbalanced. To some extent this is shown by the male-female population pyramid below.

Qatar demographic pyramid. Vastly more males than females, mostly foreign workers.

Qatar demographic pyramid. The imbalance is caused by the presence of vastly more male than female foreign workers.

Qatar is a country of 2,500,000 residents, of whom 310,000 are locals — citizens and permanent residents. The rest are foreign workers; long term inhabitants without permanent residency or citizenship. They make up 85% of the population. They are  recruited from poor, English-speaking Muslim countries mostly: Egypt, Malaysia, Tunisia. They do all the work, as best I could tell. I saw no one who looked like a local working, male or female.

Foreign workers have very few rights, but don’t seem unhappy. There is no right to unionize, and not even the right to roam around the country. For the most part, they live in employer-owned housing, and are transported back and forth to work in employer vans. They sign up for year-long contracts, and at the end of the year, they have the choice to re-up or leave. Up a year ago, foreign workers could not become permanent residents. As of last year, the Emir’s order 10 authorized permanent residency status for as many as 100 foreign workers who had sufficient means, had been in Qatar for 10 to 20 years (depending on whether they were born there), had stayed out of trouble, and who otherwise were considered desirable. It’s a step.

I suspect that the foreign workers feel lucky to have good pay, decent hours, and a clean bed. Then again, the workers are recruited for positive outlook. And the ones I saw might have had more rights than most. The airport is part of the Umm Al Houl, free enterprise zone. These are areas of Qatar where westerners and their vices like alcohol are tolerated and welcome.

Qatar natural gas production. Natural gas provides 90% of the country's income as best I can tell.

Qatar natural gas production. Natural gas provides 90% of the country’s income as best I can tell. That’s half the GDP almost, the rest of the GDP is Qataris spending the money

There are three “free enterprise zones” in Qatar; the name for the one near the airport, “Al Houl” interestingly enough means “bird trap”. What’s going on with them, as best I can tell, is diversification. Qatar is the worlds second largest exporter of natural gas, with most going to Europe, and a significant portion to India and China. But the gas will run out eventually. They are trying to supplant this income with tourism, industry and transport: running a major airline, a bustling, air hub, and tourist hotels. The airline is only marginally profitable, and though I didn’t see the hotels, I imagine they are luxurious and marginally profitable too. Saudi Arabia, next door, is trying to diversify the same ways, aiming to control west-east, air-traffic via Emirates air.

The GDP of Qatar is $191 B as of last year at the going exchange, and over $450 B at price parity. That suggests a few things. For one that the Qatari currency is undervalued. It also suggests a per-capita GDP of at least $76,400, or perhaps of $616,000 or higher depending on how you count buying power and foreign workers. This money buys a nice lifestyle, if not republican freedoms.  In terms of government, Qatar is a real monarchy, Emir Hamad bin Khalifa al-Thani’s is an absolute ruler who came to power the traditional way: he overthrew his father. Similar to this, his father, Khalifa al-Thani, came to power by overthrowing his cousin. Supporting the Emir’s rule, there is an Advisory Council. The 35 ministers are mostly relatives, and as in North Korea, it has only advisory power. The Prime Minister and Minister of Foreign Affairs is Sheikh Hamad bin Jasim bin Jabir al-Thani; the Deputy Prime Minister is Abdallah Al-Thani. The Economy and Commerce minister is Fahd Al-Thani, and the Communications and Transport minister is Ahmad Al-Thani. Nasir al-Thani heads Cabinet Affairs; Hamad al-Thani is the Secretary of State, and the Governor of the Central Bank is Abdallah bin Saud al-Thani.

Qatar main mosque. Residents stand out from the foreign workers.

Qatar main mosque. Residents stand out from the foreign workers.

My sense was that Qatar was the Disneyland version of Islam. Life in the Qatari free zones resembled normal Islamic life the way that Main Street of Disneyland resembles an actual main street in the US. Every citizen is well dress and rich without having to work. Western visitors are welcome, and not forced to follow the local customs with vices in their own zones. And the state supports all ecological and left-wing causes except for unionization. It’s anti Israel, pro revolution (elsewhere of course) and virulently against petroleum production in all counties outside of Qatar. Al Jazeera, the Emir’s left-leaning news agency, spreads money and influence world-wide. Left-flavored news is presented with high-quality graphics, and different versions of the news story published in different languages. The Emir acknowledges that Al Jazeera is a money-losing propaganda agency, but as with Disneyland, most people seem happy to live the fiction.

Qatari woman and shop. They blend into the scenery compared to the resplendent men

Qatari woman and shop. They blend into the scenery compared to the resplendent men

The local Qataris seem happy with their lot, as best I can tell. The next world soccer tournament will be held in Qatar, 2022, and Qatari’s are excited, as best I can tell. There is a lot of building going on, some for the world cup, the rest for general tourism and the free enterprise zones. The free enterprise zones may catch on, but there is a cold war going on with Saudi Arabia, and the Saudi’s are doing what they can to pour cold water on the programs. So far Qatar seems to be winning the propaganda war at home and abroad. Its people are happy, it shows a beautiful, progressive face to the west, and it seems to have the majority of the middle east travel. Stable but for how long?

Robert Buxbaum April 15, 2019. As I side note, I just bought a Qatari Thobe.

Speed traps penalize the poor

On a street corner about 1/4 mile from my house, at the intersection of the two busiest of the local streets, in the center-median of the street, is parked a police car. He’s there, about 18 hours a day, looking to give out tickets. The cross-street that this officer watches is where drivers get off the highway. In theory, they should instantly go from 65 mph on the highway to 35 mph now. Very few people do. The officer does not ticket every car, by the way, but seems to target those of poor people from outside the city limits. The only time ai was ticketed, I was driving a broken-down car while mine was in the shop. As best I can tell, he choose cars for revenue, not for safety. It’s a speed trap. It’s appalling. And our city isn’t alone in having one.

Speed traps are an annoyance to rich, local folk who sometimes get ticketed, but they’re a disaster for the poor. Poor people are targeted, and these people don’t have any savings. They don’t have the means to pay a suddenly imposed bill of $150 or more. Meanwhile, the speed-trap officer is incentivized to increase revenue and look for other violations: expired registrations or insurance, seat-belt violations, open alcohol, unpaid tickets. Double and triple fines are handed out, and sometimes the car is impounded. A poor driver is often left without any legal way to get to work, to earn money to pay the fines. Police officers behave this way because they are evaluated based on the revenue they generate, based on the number of tickets they write. It’s a horrible situation, especially for the poor

Speed traps to little and cost much.

An article on the effect of speed traps. It appears they do little good and cause much pain, especially to the poor. Here is a link to the whole article.

The article above looks at the impact of speed traps on poor people. The damage is extreme. The folks targeted are often black, barely holding it together financially. They are generally not in a position to pay $150 for “impeding traffic,” and even less in a position to deal with having their car impounded. How are they supposed to pay the bill? And yet they are told they are lucky to have been given this ticket — impeding traffic, a ticket with no “points.” But they are not lucky. They are victims. Tickets with no points is are money generators, and many poor people realize it. If they were to get a speeding ticket, they would have the opportunity to void the penalty by going to traffic school. With a ticket for impeding traffic, there is no school option. Revenue stays local, mostly in that police precinct. Poor people know it, and they don’t like it. I don’t either. After a while, poor people cease to trust the police, or to even speak to them.

In what world should you pay $150 for impeding traffic, by the way? In what world should the police be taken from their main job protecting the people and turned into a revenue arm for the city? I’d like to see this crazy cycle ended. The first steps, I think, are to end speed traps, and to limit the incentive for giving minor tickets, like impeding traffic. As it is we have too many people in jail and too many harsh penalties. 

Robert Buxbaum, April 10, 2019. I ran for water commissioner in 2016, and may run again in 2020.

Let’s visit an earth-like planet: Trappist-1d

According to Star Trek, Vulcans and Humans meet for the first time on April 5, 2063, near the town of Bozeman, Montana. It seems that Vulcan is a relatively nearby, earth-like planet with strongly humanoid inhabitants. It’s worthwhile to speculate why they are humanoid (alternatively, how likely is it that they are), and also worthwhile to figure out which planets we’d like to visit assuming we’re the ones who do the visiting.

First things first: It’s always assumed that life evolved on earth from scratch, as it were, but it is reasonably plausible that life was seeded here by some space-traveling species. Perhaps they came, looked around and left behind (intentionally or not) some blue-green algae, or perhaps some more advanced cells, or an insect or two. A billion or so years later, we’ve evolved into something that is reasonably similar to the visiting life-form. Alternately, perhaps we’d like to do the exploring, and even perhaps the settling. The Israelis are in the process of showing that low-cost space travel is a thing. Where do we want to go this century?

As it happens we know there are thousands of stars with planets nearby, but only one that we know that has reasonably earth-like planets reasonably near. This one planet circling star is Trappist-1, or more properly Trappist 1A. We don’t know which of the seven planets that orbit Trappist-1A is most earth-like, but we do know that there are at least seven planets, that they are all roughly earth size, that several have earth-like temperatures, and that all of these have water. We know all of this because the planetary paths of this star are aligned so that seven planets cross the star as seen from earth. We know their distances from their orbital times, and we know the latter from the shadows made as the planets transit. The radiation spectrum tells us there is water.

Trappist 1A is smaller than the sun, and colder than the sun, and 1 billion years older. It’s what is known as an ultra-cool dwarf. I’d be an ultra cool dwarf too, but I’m too tall. We can estimate the mass of the star and can measure its brightness. We then can calculate the temperatures on the planets based their distance from the star, something we determine as follows:

The gravitational force of a star, mass M, on a planet of mass, m,  is MmG/r2, where G is the gravitational constant, and r is the distance from the star to the planet. Since force = mass times acceleration, and the acceleration of a circular orbit is v2/r, we can say that, for these orbits (they look circular),

MmG/r2 = mv2/r = mω2r.

Here, v is the velocity of the planet and ω is its rotational velocity, ω = v/r. Eliminating m, we find that

r3 = MG/ω2.

Since we know G and ω, and we can estimate M (it’s 0.006 solar masses, we think), we have a can make good estimates of the distances of all seven planets from their various rotation speeds around the star, ω. We find that all of these planets are much closer to their star than we are to ours, so the their years are only a few days or weeks long.

We know that three planets have a temperatures reasonably close to earths, and we know that these three also have water based on observation of the absorption of light from their atmosphere as they pass in front of their star. To tell the temperature, we use our knowledge of how bright the star is (0.0052 times Sol), and our knowledge of the distance. As best we can tell, the following three of the Trappist-1 planets should have liquid surface water: Trappist 1c, d and e, the 2nd, 3rd and 4th planets from the star. With three planets to choose from, we can be fairly sure that at least one will be inhabitable by man somewhere in the planet.

The seven orbital times are in small-number ratios, suggesting that the orbits are linked into a so-called Laplace resonance-chain. For every two orbits of the outermost planet, the next one in completes three orbits, the next one completes four, followed by 6, 9 ,15, and 24. The simple whole number relationships between the periods are similar to the ratios between musical notes that produce pleasant and harmonic sounds as I discussed here. In the case of planets, resonant ratios keep the system stable. The most earth-like of the Trappist-1 planets is likely Trappist-1d, the third planet from the star. It’s iron-core, like earth, with water and a radius 1.043 times earth’s. It has an estimated average temperature of 19°C or 66°F. If there is oxygen, and if there is life there could well be, this planet will be very, very earth-like.

The temperature of the planet one in from this, Trappist-1c, is much warmer, we think on average, 62°C (143°F). Still, this is cool enough to have liquid water, and some plants live in volcanic pools on earth that are warmer than this. Besides this is an average, and we might the planet quite comfortable at the poles. The average temperature of the planet one out from this, Trappist-1e, is ice cold, -27°C (-17°F), an ice planet, it seems. Still, life can find a way. There is life on the poles of earth, and perhaps the plant was once warmer. Thus, any of these three might be the home to life, even humanoid life, or three-eyed, green men.

Visiting Trappist-1A won’t be easy, but it won’t be out-of hand impossible. The system is located about 39 light years away, which is far, but we already have a space ship heading out of the solar system, and we are developing better, and cheaper options all the time. The Israeli’s have a low cost, rocket heading to the moon. That is part of the minimal technology we’d want to visit a nearby star. You’d want to add enough rocket power to reach relativistic speeds. For a typical rocket this requires a fuel whose latent energy is on the order mc2. That turns out to be about 1 GeV/atomic mass. The only fuel that has such high power density is matter-antimatter annihilation, a propulsion system that might have time-reversal issues. A better option, I’d suggest is ion-propulsion with hydrogen atoms taken in during the journey, and ejected behind the rocket at 100 MeV energies by a cyclotron or bevatron. This system should work if the energy for the cyclotron comes from solar power. Perhaps this is the ion-drive of Star-Trek fame. To meet the Star-Trek’s made-up history, we’d have to meet up by April, 2063: forty-four years from now. If we leave today and reach near light speed by constant acceleration for a few of years, we could get there by then, but only as time is measured on the space-ship. At high speeds, time moves slower and space shrinks.

This planetary system is named Trappist-1 after the telescope used to discover it. It was the first system discovered by the 24 inch, 60 cm aperture, TRAnsiting Planets and PlanetesImals Small Telescope. This telescope is operated by The University of Liége, Belgium, and is located in Morocco. The reason most people have not heard of this work, I think, has to do with it being European science. Our news media does an awful job covering science, in my opinion, and a worse job covering Europe, or most anything outside the US. Finally, like the Israeli moon shot, this is a low-budget project, the work to date cost less than €2 million, or about US $2.3 million. Our media seems committed to the idea that only billions of dollars (or trillions) will do anything, and that the only people worth discussing are politicians. NASA’s budget today is about $6 billion, and its existence is barely mentioned.

The Trappist system appears to be about 1 billion years older than ours, by the way, so life there might be more advanced than ours, or it might have died out. And, for all we know, we’ll discover that the Trappist folks discover space travel, went on to colonize earth, and then died out. The star is located, just about exactly on the ecliptic, in the constellation Aquarius. This is an astrological sign associated with an expansion of human consciousness, and a revelation of truths. Let us hope that, in visiting Trappist, “peace will guide the planets and love will steer the stars”.

Robert Buxbaum, April 3, 2019. Science sources are: http://www.trappist.one. I was alerted to this star’s existence by an article in the Irish Times.

A hydrogen permeation tester

Over the years I’ve done a fair amount of research on hydrogen permeation in metals — this is the process of the gas dissolving in the metal and diffusing to the other side. I’ve described some of that, but never the devices that measure the permeation rate. Besides, my company, REB Research, sells permeation testing devices, though they are not listed on our site. We recently shipped one designed to test hydrogen permeation through plastics for use in light weight hydrogen tanks, for operation at temperatures from -40°C to 85°C. Shortly thereafter we got another order for a permeation tester. With all the orders, I thought I’d describe the device a bit — this is the device for low permeation materials. We have a similar, but less complex design for high permeation rate material.

Shown below is the central part of the device. It is a small volume that can be connected to a high vacuum, or disconnected by a valve. There is an accurate pressure sensor, accurate to 0.01 Torr, and so configured that you do not get H2 + O2 reactions (something that would severely throw off results). There is also a chamber for holding a membrane so one side is help in vacuum, in connection to the gauge, and the other is exposed to hydrogen, or other gas at pressures up to 100 psig (∆P =115 psia). I’d tested to 200 psig, but currently feel like sticking to 100 psig or less. This device gives amazingly fast readings for plastics with permeabilities as low as 0.01 Barrer.

REB Research hydrogen permeation tester cell with valve and pressure sensor.

REB Research hydrogen permeation tester cell with valve and pressure sensor.

To control the temperature in this range of interest, the core device shown in the picture is put inside an environmental chamber, set up as shown below, with he control box outside the chamber. I include a nitrogen flush device as a safety measure so that any hydrogen that leaks from the high pressure chamber will not build up to reach explosive limits within the environmental chamber. If this device is used to measure permeation of a non-flammable gas, you won’t need to flush the environmental chamber.

I suggest one set up the vacuum pump right next to the entrance of the chamber; in the case of the chamber provided, that’s on the left as shown with the hydrogen tank and a nitrogen tank to the left of the pump. I’ve decided to provide a pressure sensor for the N2 (nitrogen) and a solenoidal shutoff valve for the H2 (hydrogen) line. These work together as a safety feature for long experiments. Their purpose is to automatically turn off the hydrogen if the nitrogen runs out. The nitrogen flush part of this process is a small gauge copper line that goes from the sensor into the environmental chamber with a small, N2 flow bleed valve at the end. I suggest setting the N2 pressure to 25-35 psig. This should give a good inert flow into the environmental chamber. You’ll want a nitrogen flush, even for short experiments, and most experiments will be short. You may not need an automatic N2 sensor, but you’ll be able to do this visually.

Basic setup for REB permeation tester and environmental chamber

Basic setup for REB permeation tester and environmental chamber

I shipped the permeation cell comes with some test, rubbery plastic. I’d recommend the customer leave it in for now, so he/she can use it for some basic testing. For actual experiments, you replace mutest plastic with the sample you want to check. Connect the permeation cell as shown above, using VCR gaskets (included), and connect the far end to the multi-temperature vacuum hose, provided. Do this outside of the chamber first, as a preliminary test to see if everything is working.

For a first test live the connections to the high pressure top section unconnected. The pressure then will be 1 atm, and the chamber will be full of air. eave the top, Connect the power to the vacuum pressure gauge reader and connect the gauge reader to the gauge head. Open the valve and turn on the pump. If there are no leaks the pressure should fall precipitously, and you should see little to no vapor coming out the out port on the vacuum pump. If there is vapor, you’ve got a leak, and you should find it; perhaps you didn’t tighten a VCR connection, or you didn’t do a good job with the vacuum hose. When things are going well, you should see the pressure drop to the single-digit, milliTorr range. If you close the valve, you’ll see the pressure rise in the gauge. This is mostly water and air degassing from the plastic sample. After 30 minutes, the rate of degassing should slow and you should be able to measure the rate of gas permeation in the polymer. With my test plastic, it took a minute or so for the pressure to rise by 10 milliTorr after I closed the valve.

If you like, you can now repeat this preliminary experiment with hydrogen connect the hydrogen line to one of the two ports on the top of the permeation cell and connect the other port to the rest of the copper tubing. Attach the H2 bleed restrictor (provided) at the end of this tubing. Now turn on the H2 pressure to some reasonable value — 45 psig, say. With 45 psi (3 barg upstream) you will have a ∆P of 60 psia or 4 atm across the membrane; vacuum equals -15 psig. Repeat the experiment above; pump everything down, close the valve and note that the pressure rises faster. The restrictor allows you to maintain a H2 pressure with a small, cleansing flow of gas through the cell.

If you like to do these experiments with a computer record, this might be a good time to connect your computer to the vacuum reader/ controller, and to the thermocouple, and to the N2 pressure sensor. 

Here’s how I calculate the permeability of the test polymer from the time it takes for a pressure rise assuming air as the permeating gas. The volume of the vacuumed out area after the valve is 32 cc; there is an open area in the cell of 13.0 cm2 and, as it happens, the  thickness of the test plastic is 2 mm. To calculate the permeation rate, measure the time to rise 10 millitorr. Next calculate the millitorr per hour: that’s 360 divided by the time to rise ten milliTorr. To calculate ncc/day, multiply the millitorr/hour by 24 and by the volume of the chamber, 32 cc, and divide by 760,000, the number of milliTorr in an atmosphere. I found that, for air permeation at ∆P = one atm, I was getting 1 minute per milliTorr, which translates to about 0.5 ncc/day of permeation through my test polymer sheet. To find the specific permeability in cc.mm/m2.day.atm, I multiply this last number by the thickness of the plastic (2 mm in this case), divide by the area, 0.0013 m2, and divide by ∆P, 1 atm, for this first test. Calculated this way, I got an air permeance of 771 cc.mm/m2.day.atm.

The complete setup for permeation testing.

The complete setup for permeation testing.

Now repeat the experiment with hydrogen and your own plastic. Disconnect the cell from both the vacuum line and from the hydrogen in line. Open the cell; take out my test plastic and replace it with your own sample, 1.87” diameter, or so. Replace the gasket, or reuse it. Center the top on the bottom and retighten the bolts. I used 25 Nt-m of torque, but part of that was using a very soft rubbery plastic. You might want to use a little more — perhaps 40-50 Nt-m. Seal everything up. Check that it is leak tight, and you are good to go.

The experimental method is the same as before and the only signficant change when working with hydrogen, besides the need for a nitrogen flush, is that you should multiply the time to reach 10 milliTorr by the square-root of seven, 2.646. Alternatively, you can multiply the calculated permeability by 0.378. The pressure sensor provided measures heat transfer and hydrogen is a better heat transfer material than nitrogen by a factor of √7. The vacuum gauge is thus more sensitive to H2 than to N2. When the gauge says that a pressure change of 10 milliTorr has occurred, in actuality, it’s only 3.78 milliTorr.  The pressure gauge reads 3.78 milliTorr oh hydrogen as 10 milliTorr.

You can speed experiments by a factor of ten, by testing the time to rise 1 millitorr instead of ten. At these low pressures, the gauge I provided reads in hundredths of a milliTorr. Alternately, for higher permeation plastics (or metals) you want to test the time to rise 100 milliTorr or more, otherwise the experiment is over too fast. Even at a ten millTorr change, this device gives good accuracy in under 1 hour with even the most permeation-resistant polymers.

Dr. Robert E. Buxbaum, March 27, 2019; If you’d like one of these, just ask. Here’s a link to our web site, REB Research,

Statistics for psychologists, sociologists, and political scientists

In terms of mathematical structure, psychologists, sociologists, and poly-sci folks all do the same experiment, over and over, and all use the same simple statistical calculation, the ANOVA, to determine its significance. I thought I’d explain that experiment and the calculation below, walking you through an actual paper (one I find interesting) in psychology / poly-sci. The results are true at the 95% level (that’s the same as saying p >0.05) — a significant achievement in poly-sci, but that doesn’t mean the experiment means what the researchers think. I’ll then suggest another statistic measure, r-squared, that deserves to be used along with ANOVA.

The standard psychological or poly-sci research experiments involves taking a group of people (often students) and giving them a questionnaire or test to measure their feelings about something — the war in Iraq, their fear of flying, their degree of racism, etc. This is scored on some scale to get an average. Another, near-identical group of subjects is now brought in and given a prompt: shown a movie, or a picture, or asked to visualize something, and then given the same questionnaire or test as the first group. The prompt is shown to have changed to average score, up or down, an ANOVA (analysis of variation) is used to show if this change is one the researcher can have confidence in. If the confidence exceeds 95% the researcher goes on to discuss the significance, and submits the study for publication. I’ll now walk you through the analysis the old fashioned way: the way it would have been done in the days of hand calculators and slide-rules so you understand it. Even when done this way, it only takes 20 minutes or so: far less time than the experiment.

I’ll call the “off the street score” for the ith subject, Xi°. It would be nice if papers would publish these, but usually they do not. Instead, researchers publish the survey and the average score, something I’ll call X°-bar, or X°. they also publish a standard deviation, calculated from the above, something I’ll call, SD°. In older papers, it’s called sigma, σ. Sigma and SD are the same thing. Now, moving to the group that’s been given the prompt, I’ll call the score for the ith subject, Xi*. Similar to the above, the average for this prompted group is X*, or X°-bar, and the standard deviation SD*.

I have assumed that there is only one prompt, identified by an asterix, *, one particular movie, picture, or challenge. For some studies there will be different concentrations of the prompt (show half the movie, for example), and some researchers throw in completely different prompts. The more prompts, the more likely you get false positives with an ANOVA, and the more likely you are to need to go to r-squared. Warning: very few researchers do this, intentionally (and crookedly) or by complete obliviousness to the math. Either way, if you have a study with ten prompt variations, and you are testing to 95% confidence your result is meaningless. Random variation will give you this result 50% of the time. A crooked researcher used ANOVA and 20 prompt variations “to prove to 95% confidence” that genetic modified food caused cancer; I’ll assume (trust) you won’t fall into that mistake, and that you won’t use the ANOVA knowledge I provide to get notoriety and easy publication of total, un-reproducible nonsense. If you have more than one or two prompts, you’ve got to add r-squared (and it’s probably a good idea with one or two). I’d discuss r-squared at the end.

I’ll now show how you calculate X°-bar the old-fashioned way, as would be done with a hand calculator. I do this, not because I think social-scientists can’t calculate an average, nor because I don’t trust the ANOVA function on your laptop or calculator, but because this is a good way to familiarize yourself with the notation:

X°-bar = X° = 1/n° ∑ Xi°.

Here, n° is the total number of subjects who take the test but who have not seen the prompt. Typically, for professional studies, there are 30 to 50 of these. ∑ means sum, and Xi° is the score of the ith subject, as I’d mentioned. Thus, ∑ Xi° indicates the sum of all the scores in this group, and 1/n° is the average, X°-bar. Convince yourself that this is, indeed the formula. The same formula is used for X*-bar. For a hand calculation, you’d write numbers 1 to n° on the left column of some paper, and each Xi° value next to its number, leaving room for more work to follow. This used to be done in a note-book, nowadays a spreadsheet will make that easier. Write the value of X°-bar on a separate line on the bottom.

T-table

T-table

In virtually all cases you’ll find that X°-bar is different from X*-bar, but there will be a lot of variation among the scores in both groups. The ANOVA (analysis of variation) is a simple way to determine whether the difference is significant enough to mean anything. Statistics books make this calculation seem far too complicated — they go into too much math-theory, or consider too many types of ANOVA tests, most of which make no sense in psychology or poly-sci but were developed for ball-bearings and cement. The only ANOVA approach used involves the T-table shown and the 95% confidence (column this is the same as two-tailed p<0.05 column). Though 99% is nice, it isn’t necessary. Other significances are on the chart, but they’re not really useful for publication. If you do this on a calculator, the table is buried in there someplace. The confidence level is written across the bottom line of the cart. 95% here is seen to be the same as a two-tailed P value of 0.05 = 5% seen on the third from the top line of the chart. For about 60 subjects (two groups of 30, say) and 95% certainty, T= 2.000. This is a very useful number to carry about in your head. It allows you to eyeball your results.

In order to use this T value, you will have to calculate the standard deviation, SD for both groups and the standard variation between them, SV. Typically, the SDs will be similar, but large, and the SV will be much smaller. First lets calculate SD° by hand. To do this, you first calculate its square, SD°2; once you have that, you’ll take the square-root. Take each of the X°i scores, each of the scores of the first group, and calculate the difference between each score and the average, X°-bar. Square each number and divide by (n°-1). These numbers go into their own column, each in line with its own Xi. The sum of this column will be SD°2. Put in mathematical terms, for the original group (the ones that didn’t see the movie),

SD°2 = 1/(n°-1) ∑ (Xi°- X°)2

SD° = √SD°2.

Similarly for the group that saw the movie, SD*2 = 1/(n*-1) ∑ (Xi*- X*)2

SD* = √SD*2.

As before, n° and n* are the number of subjects in each of the two groups. Usually you’ll aim for these to be the same, but often they’ll be different. Some students will end up only seeing half the move, some will see it twice, even if you don’t plan it that way; these students’ scares can not be used with the above, but be sure to write them down; save them. They might have tremendous value later on.

Write down the standard deviations, SD for each group calculated above, and check that the SDs are similar, differing by less than a factor of 2. If so, you can take a weighted average and call it SD-bar, and move on with your work. There are formulas for this average, and in some cases you’ll need an F-table to help choose the formula, but for my purposes, I’ll assume that the SDs are similar enough that any weighted average will do. If they are not, it’s likely a sign that something very significant is going on, and you may want to re-think your study.

Once you calculate SD-bar, the weighted average of the SD’s above, you can move on to calculate the standard variation, the SV between the two groups. This is the average difference that you’d expect to see if there were no real differences. That is, if there were no movie, no prompt, no nothing, just random chance of who showed up for the test. SV is calculated as:

SV = SD-bar √(1/n° + 1/n*).

Now, go to your T-table and look up the T value for two tailed tests at 95% certainty and N = n° + n*. You probably learned that you should be using degrees of freedom where, in this case, df = N-2, but for normal group sizes used, the T value will be nearly the same. As an example, I’ll assume that N is 80, two groups of 40 subjects the degrees of freedom is N-2, or 78. I you look at the T-table for 95% confidence, you’ll notice that the T value for 80 df is 1.99. You can use this. The value for  62 subjects would be 2.000, and the true value for 80 is 1.991; the least of your problems is the difference between 1.991 and 1.990; it’s unlikely your test is ideal, or your data is normally distributed. Such things cause far more problems for your results. If you want to see how to deal with these, go here.

Assuming random variation, and 80 subjects tested, we can say that, so long as X°-bar differs from X*-bar by at least 1.99 times the SV calculated above, you’ve demonstrated a difference with enough confidence that you can go for a publication. In math terms, you can publish if and only if: |X°-X*| ≥ 1.99 SV where the vertical lines represent absolute value. This is all the statistics you need. Do the above, and you’re good to publish. The reviewers will look at your average score values, and your value for SV. If the difference between the two averages is more than 2 times the SV, most people will accept that you’ve found something.

If you want any of this to sink in, you should now do a worked problem with actual numbers, in this case two groups, 11 and 10 students. It’s not difficult, but you should try at least with these real numbers. When you are done, go here. I will grind through to the answer. I’ll also introduce r-squared.

The worked problem: Assume you have two groups of people tested for racism, or political views, or some allergic reaction. One group was given nothing more than the test, the other group is given some prompt: an advertisement, a drug, a lecture… We want to know if we had a significant effect at 95% confidence. Here are the test scores for both groups assuming a scale from 0 to 3.

Control group: 0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 3.  These are the Xi° s; there are 11 of them

Prompted group: 0, 1, 1, 1, 2, 2, 2, 2, 3, 3.  These are the Xi* s; there are 10 of them.

On a semi-humorous side: Here is the relationship between marriage and a PhD.

Robert Buxbaum, March 18, 2019. I also have an explanation of loaded dice and flipped coins; statistics for high school students.

C-Pap and Apnea

A month of so ago, I went to see a sleep doctor for my snoring. I got a take-home breathing test that gave me the worst night’s sleep in recent memory. A few days later, I got a somber diagnosis: “You are a walking zombie.” Apparently, I hold my breath for ten seconds or more every minute and a half while sleeping. Normal is supposed to be every 4 to 10 minutes. But by this standard, more than half of all middle-aged men are sub-normal (how is this possible?). As a result of my breath-holding, the wrinkled, unsmiling DO claimed I’m brain-dead now and will soon be physically dead unless I change my ways. Without spending 3 minutes with me, the sleep expert told me that I need to lose weight, and that I need a C-Pap (continuous positive airway pressure) device as soon as possible. It’s supposed to help me lose that weight and get back the energy. With that he was gone. The office staff gave me the rest of the dope: I was prescribed  a “ResMed” brand C-Pap, supplied by a distributor right across the hall from the doctor (how convenient).

I picked up the C-Pap three months later. Though I was diagnosed as needing one “as soon as possible,” no one would release the device until they were sure it was covered by my insurance company. The device when I got it, was something of a horror. The first version I tried fit over the whole face and forces air into my mouth and nose simultaneously, supposedly making it easier to inhale, but harder to exhale. I found it more than a bit uncomfortable. The next version was nose only and marginally more comfortable. I found there was a major air-flow restriction when I breath in and a similar pressure penalty when I breathed out. And it’s loud. And, if you open your mouth, there is a wind blowing through. As for what happens if the pump fails or the poor goes out, I notice that there are the tiniest of air-holes to prevent me from suffocating, barely. A far better design would have given me a 0-psi flapper valve for breathing in, and a 1/10 psi flapper for breathing out. That would also reduce the pressure restriction I was feeling every time I took a deep breath. One of my first blog essays was about engineering design aesthetics; you want your designs to improve things under normal conditions and fail safe, not like here. Using this device while awake was anything but pleasant, and I found I still hold my breath, even while awake, about every 5 minutes.

Since I have a lab, and the ability to test these things, I checked the pressure of the delivered air, and found it was 3 cm of water, about 1/20 psi. The prescription was for 5 cm or water (1/14 psi). The machine registers this, but it is wrong. I used a very simple water manometer, a column of water, similar to the one I used to check the pressure drop in furnace air filters. Is 1/20 psi enough?How did he decide on 1/14 psi by the way? I’ve no idea. !/14 psi is about 1/200 atm. Is this enough to do anything? While the C-Pap should get me to breathe more, I guess, about half of all users stop after a few tries, and my guess is that they find it as uncomfortable as I have. There is no research evidence that treatment with it reduces stroke or heart attack, or extends life, or helps with weight loss. The assumption is that, if you force middle-aged men to hold their breath less, they will be healthier, but I’ve no clear logic or evidence to back the assumption. At best, anything you gain on the ease of breathing in, you lose on the difficulty of breathing out. The majority of middle-aged men are prescribed a C-Pap, if they go for a sleep study, and it’s virtually 100% for overweight men with an apple-shaped body.

I’d have asked my doctor about alternatives or for a second opinion but he was out the door too fast. Besides, I was afraid I’d get the same answer that Rodney Dangerfield got: “You want a second opinion? OK. You’re ugly, too.” Mr. Dangerfield was not a skinny comic, by the way, but he was funny, and I assume he’d have been prescribed a C-Pap (maybe he was). He died at 82, considerably older than Jim Fixx, “the running doctor,” Adelle Davis, the “eat right for health” doctor, Euell Gibbons “in search of the wild asparagus,” or Ethan Pritkin, the diet doctor. God seems to prefer fat comedians to diet experts; I expect that most-everyone does.

Benjamin Franklin and his apple-shaped body

Benjamin Franklin and his apple-shaped body; I don’t think of him as a zombie.

What really got my goat, besides my dislike of the C-Pap, is that I object to being called a walking zombie. True, I’m not as energetic as I used to be, but I manage to run a company, and to write research papers, and I get patents (this one was approved just today). And I write these blogs — I trust that any of you who’ve read this far find them amusing. Pretty good for a zombie — and I ran for water commissioner. People who use the C-Pap self-report that they have more energy, but self-reporting is poor evidence. A significant fraction of those people who start with the C-Pap, stop, and those people, presumably were not happy. Besides, a review of the internet suggests that a similarly large fraction of those who buy a “MyPillow.com” claim they have more energy. And I’ve seen the same claims from people who take a daily run, or who pray, or smoke medical marijuana (available for sleep apnea, but not from this fellow), or Mirtazapine (study results here), or  for electro-shock therapy, a device called “Inspire.” With so many different products providing the same self-reported results, I wonder if there isn’t something more fundamental going on. I’d wish the doc had spent a minute or two to speak to this, or to the alternatives.

As for weight loss, statistical analysis of lifespan suggests that there is a health advantage to being medium weight: not obese, but not skinny. I present some of this evidence here, along with evidence that extra weight helps ward off Alzheimer’s. For all I know this protection is caused by holding your breath every few minutes. It helps to do light exercise, but not necessary for mental health. In terms of mental health, the evidence suggests that weight loss is worse than nothing.

Jared Gray, author of the Alien movies, was diagnosed with apnea, so he designed his own sleep-mask.

Jared Gray, author of the Alien movies, was diagnosed with apnea, so he designed his own sleep-mask.

Benjamin Franklin was over-weight and apple-shaped, and no zombie, The same is true of John Adams, Otto Von Bismarck, and Alfred Hitchcock. All lived long, productive lives. Hitchcock was sort of morbid, it will be admitted, but I would not want him otherwise. Ed McMahon, Johnny Carson’s side-kick, apologized to America for being overweight and smoking, bu the outlived Johnny Carson by nine years, dying at 89. Henry Kissinger is still alive and writing at 95. He was always fatter than any of the people he served. He almost certainly had sleep apnea, back in the day, and still has more on the ball, in my opinion, than most of the talking-head on TV. The claim that overweight, middle-aged men are all zombies without a breath assisting machine doesn’t make no sense to me. But then, I’m not a sleep doctor. (Do sleep doctors get commissions? Why did he choose, this supplier or this brand device? With so little care about patients, I wonder who runs the doctor’s office.)

I looked up my doctor on this list provided by the American Board of Sleep Medicine. I found my doctor was not certified in sleep medicine. I suppose certified doctors would prescribe something similar  but was disappointed that you don’t need sleep certification to operate as a sleep specialist. In terms of masks, I figure, if you’ve got to wear something, you might as well wear something cool. Author Jared Gray, shown above (not the author of the Alien) was diagnosed with Apnea 6 months ago and made his own C-Pap mask to make it look like the alien was attacking him. Very cool for an ex-zombie, but I’m waiting to see a burst of creative energy.

What do we zombies want? Brains.

When do we want them? Brains.

What do vegetarian zombies want? Grains.

Robert Buxbaum, March 15, 2019. In case real zombies should attack, here’s what to do.  An odd legal/insurance issue: in order to get the device, I had to sign that, if I didn’t use it for 20 days in the first month of 4 hours per night, and thus if the insurance did not pay, I would be stuck with the full fee. I signed. This might cost me $1000 though normally in US law, companies can only charge a reasonable restock fee, but it can’t be unreasonable, like the full  price. I also had to sign that I would go back to the same, quick-take doctor, but again there has to be limits. We’ll see how the machine pans out, but one difference I see already: unlike my pillow.com, there is no money back guarantee with the C-Pap treatment.

Seize the day

It is forbidden knowledge what our term of years, mine and yours.
Don’t scan the tables of your Babylonian seers.
Better far to bear the future, my Leuconoe, like the past.
Whether Jupiter has many years yet to give,
Or this one is our last:

This, that makes the Tyrrhene waves spent against the shore.
Strain your wine and strain your wisdom.
Life is short; should we hope for long?
In the moment of our talking, precious time has slipped away.
Seize the moment. Trust tomorrow little as you may.

by Horace (23 BC Roman poet) Odes, 1.11

This poem by Horace, in 23 BC is the first appearance of the phrase “carpe diem,” translated as seize he day. I’d decided to look over the translation from Wikipedia, and to correct and update the translation as I saw fit, to some extent to extract the meanings better, to some extent to make the grammar less-clunky, and to some extent to make it rhyme. Seen in context, the whole poem looks  romantic, and the intent of the famous phrase seems more like ‘seize the moment’ when read in context. Either translation is acceptable from the Latin, as I understand it.

The phrase, “seize the day” appears in several important movies. Robin Williams speaks it to a class of literature students in the sense that I read it here — seize the moment — in this scene of “The Dead Poets Society,” He’s trying to get the boys to appreciate the purpose of poetry, and the preciousness of their years in prep-school. A well-done movie, IMHO. The newspaper sellers sing the phrase for different intent in this song in “Newsies.” For them, the intent is more like seize the opportunity, or maybe even seize power. It’s not Horace’s intent, but it’s sung in front of the statue of Horace Greeley, and it works.

In either context, there is a certain young masculinity about this phrase. In both movie, the cast experiencing the idea is male and young. I don’t think either movie would work as well with women dancing or singing to this idea.

Robert E. Buxbaum, March 9, 2019. In case you should wonder what happens to Frank Kelly (Sullivan) after the movie ends, I’ve written about that.  Also, a friend of mine notes that the grammar used in these movies is wrong:  “Carpe diem” is singular, for this 3rd declension noun. The equivalent Latin plural is “Carpite diem:” That’s the equivalent of you-all, should seize the moment. Unlike in the movies, much of classic education is spent on pedantic, uninspiring, minutia, like Latin grammar, but that’s what’s necessary to permit distinction of meaning. Thank you, David Hoenig for grammar help.

Why concrete cracks and why sealing is worthwhile

The oil tanker Palo Alto is one of several major ships made with concrete hulls.

The oil tanker Palo Alto is one of several major ships made with concrete hulls.

Modern concrete is a wonderful construction material. Major buildings are constructed of it, and major dams, and even some ships. But under the wrong circumstances, concrete has a surprising tendency to crack and fail. I thought I’d explain why that happens and what you can do about it. Concrete does not have to crack easily; ancient concrete didn’t and military or ship concrete doesn’t today. A lot of the fault lies in the use of cheap concrete — concrete with lots of filler — and with the cheap way that concrete is laid. First off, the major components of modern concrete are pretty uniform: sand and rock, Portland cement powder (made from cooked limestone, mostly), water, air, and sometimes ash. The cement component is what holds it all together — cements it together as it were — but it is not the majority of even the strongest concretes. The formula of cement has changed too, but the cement is not generally the problem. It doesn’t necessarily stick well to the rock or sand component of concrete (It sticks far better to itself) but it sticks well enough that spoliation, isn’t usually a problem by itself.

What causes problem is that the strength of concrete is strongly affected (decreased) by having lots of sand, aggregate and water. The concrete used in sidewalks is as cheap as possible, with lots of sand and aggregate. Highway and wall concrete has less sand and aggregate, and is stronger. Military and ship concrete has little sand, and is quite a lot stronger. The lowest grade, used in sidewalks, is M5, a term that refers to its compressive strength: 5 Mega Pascals. Pascals are European (Standard International) units of pressure and of strength. One Pascal is one Newton per square meter (Here’ a joke about Pascal units). In US (English) units, 5 MPa is 50 atm or 750 psi.

Ratios for concrete mixes of different strength.

Ratios for concrete mixes of different strength; the numbers I use are double these because these numbers don’t include water; that’s my “1”.

The ratio of dry ingredients in various concretes is shown at right. For M5, and including water, the ratio is 1 2 10 20. That is to say there is one part water, two parts cement, 10 parts sand, and 20 parts stone-aggregate (all these by weight). Added to this is 2-3% air, by volume, or nearly as much air as water. At least these are the target ratios; it sometimes happens that extra air and water are added to a concrete mix by greedy or rushed contractors. It’s sometimes done to save money, but more often because the job ran late. The more the mixer turns the more air gets added. If it turns too long there is extra air. It the job runs late, workers will have to add extra water too because the concrete starts hardening. I you see workers hosing down wet concrete as it comes from the truck, this is why. As you might expect, extra air and water decrease the strength of the product. M-10 and M-20 concrete have less sand, stone, and water as a proportion to cement. The result is 10 MPa or 20 MPa strength respectively.

A good on-site inspector is needed to keep the crew from adding too much water. Some water is needed for the polymerization (setting) of the concrete. The rest is excess, and when it evaporates, it leaves voids that are similar to the voids created by having air mix in. It is not uncommon to find 6% voids, in commercial concrete. This is to say that, after the water evaporates, the concrete contains about as much void as cement by volume. To get a sense of how much void space is in the normal concrete outside your house, go outside to a piece of old concrete (10 years old at least) on a hot, dry day, and pour out a cup of water. You will hear a hiss as the water absorbs, and you will see bubbles come out as the water goes in. It used to be common for cities to send inspectors to measuring the void content of the wet (and dry) concrete by a technique called “pycnometry” (that’s Greek for density measurement). I’ve not seen a local city do this in years, but don’t know why. An industrial pycnometer is shown below.

Pyncnometer used for concrete. I don't see these in use much any more.

Pycnometer used for concrete. I don’t see these in use much any more.

One of the main reason that concrete fails has to do with differential expansion, thermal stress, a concept I dealt with some years ago when figuring out how cold it had to be to freeze the balls off of a brass monkey. As an example of the temperature change to destroy M5, consider that the thermal expansion of cement is roughly 1 x 10-5/ °F or 1.8 x10-5/°C. This is to say that a 1 meter slab of cement that is heated or cooled by 100°F will expand or shrink by 10-3 m respectively; 100 x 1×10-5 = 10-3. This is a fairly large thermal expansion coefficient, as these things go. It would not cause stress-failure except that sand and rock have a smaller thermal expansion coefficients, about 0.6×10-5 — barely more than half the value for cement. Consider now what happens to concrete that s poured in the summer when it is 80°F out, and where the concrete heats up 100°F on setting (cement setting releases heat). Now lets come back in winter when it’s 0°F. This is a total of 100°F of temperature change. The differential expansion is 0.4 x 10-5/°F x 100°F =  4 x10-4 meter/meter = 4 x10-4 inch/inch.

The force created by this differential expansion is the elastic modulus of the cement times the relative change in expansion. The elastic modulus for typical cement is 20 GPa or, in English units, 3 million psi. This is to say that, if you had a column of cement (not concrete), one psi of force would compress it by 1/3,000,000. The differential expansion we calculated, cement vs sand and stone is 4×10-4 ; this much expansion times the elastic modulus, 3,000,000 = 1200 psi. Now look at the strength of the M-5 cement; it’s only 750 psi. When M-5 concrete is exposed to these conditions it will not survive. M-10 will fail on its own, from the temperature change, without any help needed from heavy traffic. You’d really like to see cities check the concrete, but I’ve seen little evidence that they do.

Water makes things worse, and not only because it creates voids when it evaporates. Water also messes up the polymerization reaction of the cement. Basic, fast setting cement is mostly Ca3SiO5

2Ca3SiO5 + 6 H2O –> 3Ca0SiO2•H2O +3Ca(OH)2•H2O.

The former of these, 3Ca0SiO2•H2O, forms something of a polymer. Monomer units of SiO4 are linked directly or by partially hydrated CaO linkages. Add too much water and the polymeric linkages are weakened or do not form at all. Over time the Ca(OH)2 can drain away or react with  CO2 in the air to form chalk.

concrete  strength versus-curing time. Slow curing of damp concrete helps; fast dry hurts. Carbonate formation adds little or no strength. Jehan Elsamni 2011.

Portland limestone cement strength versus curing time. Slow curing and damp helps; fast dry hurts. Carbonate formation adds little or no strength. Jehan Elsamni 2011.

Ca(OH)2 + CO2 → CaCO3 + H2O

Sorry to say, the chalk adds little or no strength, as the graph at right shows. Concrete made with too much water isn’t very strong at all, and it gets no stronger when dried in air. Hardening goes on for some weeks after pouring, and this is the reason you don’t drive on 1 too 2 day old concrete. Driving on weak concrete can cause cracks that would not form if you waited.

You might think to make better concrete by pouring concrete in the cold, but pouring in the cold makes things worse. Cold poured cement will expand the summer and the cement will detach from the sand and stone. Ideally, pouring should be in spring or fall, when the temperature is moderate, 40-60°F. Any crack that develops grows by a mechanism called Rayleigh crack growth, described here. Basically, once a crack starts, it concentrates the fracture forces, and any wiggling of the concrete makes the crack grow faster.

Based on the above, I’ve come to suspect that putting on a surface coat can (could) help strengthen old concrete, even long after it’s hardened. Mostly this would happen by filling in voids and cracks, but also by extending the polymer chains. I imagine it would be especially helpful to apply the surface coat somewhat watery on a dry day in the summer. In that case, I imagine that Ca3SiO5 and Ca(OH)2 from the surface coat will penetrate and fill the pores of the concrete below — the sales pores that hiss when you pour water on them. I imagine this would fill cracks and voids, and extend existing CaOSiO2•H2O chains. The coat should add strength, and should be attractive as well. At least that was my thought.

I should note that, while Portland cement is mostly Ca3SiO5, there is also a fair amount (25%) of Ca2SiO4. This component reacts with water to form the same calcium-silicate polymer as above, but does so at a slower rate using less water per gram. My hope was that this component would be the main one to diffuse into deep pores of the concrete, reacting there to strengthen the concrete long after surface drying had occurred.

Trump tower: 664', concrete and glass. What grade of concrete would you use?

Trump tower: 664′, concrete and glass. What grade of concrete would you use?

As it happened, I had a chance to test my ideas this summer and also about 3 years ago. The city inspector came by to say the concrete flags outside my house were rough, and thus needed replacing, and that I was to pay or do it myself. Not that I understand the need for smooth concrete, quite, but that’s our fair city. I applied for a building permit to apply a surface coat, and applied it watery. I used “Quickrete” brand concrete patch, and so far it’s sticking OK. Pock-holes in the old concrete have been filled in, and so far surface is smooth. We’ll have to see if my patch lasts 10-20 years like fresh cement. Otherwise, no matter how strong the concrete becomes underneath, the city will be upset, and I’ll have to fix it. I’ve noticed that there is already some crumbling at the sides of flags, something I attribute to the extra water. It’s not a problem yet, but hope this is not the beginning of something worse. If I’m wrong here, and the whole seal-coat flakes off, I’ll be stuck replacing the flags, or continuing to re-coat just to preserve my reputation. But that’s the cost of experimentation. I tried something new, and am blogging about it in the hope that you and I benefit. “Education is what you get when you don’t get what you want.” (It’s one of my wise sayings). At the worst, I’ll have spent 90 lb of patching cement to get an education. And, I’m happy to say that some of the relatively new concrete flags that the city put in are already cracked. I attribute this to: too much sand, air, water or air (they don’t look like they have much rock): Poor oversight.

Dr. Robert E. Buxbaum. March 5, 2019. As an aside, the 664 foot Trump Tower, NY is virtually the only skyscraper in the city to be built of concrete and glass. The others are mostly steel and glass. Concrete and glass is supposed to be stiffer and quieter. The engineer overseeing the project was Barbara Res, the first woman to oversee a major, NY building project. Thought question: if you built the Trump Tower, which quality of concrete would you use, and why.

A probability paradox

Here is a classic math paradox for your amusement, and perhaps your edification: (edification is a fancy word for: beware, I’m trying to learn you something).

You are on a TV game show where you will be asked to choose between two, identical-looking envelopes. All you know about the envelopes is that one of them has twice as much money as the other. The envelopes are shuffled, and you pick one. You peak in and see that your envelope contains $400, and you feel pretty good. But then you are given a choice: you can switch your envelope with the other one; the one you didn’t take. You reason that the other envelope either has $800 or $200 with equal probability. That is, a switch will either net you a $400 gain, or loose you $200. Since $400 is bigger than $200, you switch. Did that decision make sense. It seems that, at this game, every contestant should switch envelopes. Hmm.

The solution follows: The problem with this analysis is an error common in children and politicians — the confusion between your lack of knowledge of a thing, and actual variability in the system. In this case, the contestant is confusing his (or her) lack of knowledge of whether he/she has the big envelope or the smaller, with the fixed fact that the total between the two envelopes has already been set. It is some known total, in this case it is either $600 or $1200. Lets call this unknown sum y. There is a 50% chance that you now are holding 2/3 y and a 50% chance you are holding only 1/3y. therefore, the value of your current envelope is 1/3 y + 1/6y = 1/2 y. Similarly, the other envelope has a value 1/2y; there is no advantage is switching once it is accepted that the total, y had already been set before you got to choose an envelope.

And here, unfortunately is the lesson:The same issue applies in reverse when it comes to government taxation. If you assume that the total amount of goods produced by the economy is always fixed to some amount, then there is no fundamental problem with high taxes. You can print money, or redistribute it to anyone you think is worthy — more worthy than the person who has it now – and you won’t affect the usable wealth of the society. Some will gain others will lose, and likely you’ll find you have more friends than before. On the other hand, if you assume that government redistribution will affect the total: that there is some relationship between reward and the amount produced, then to the extent that you diminish the relation between work and income, or savings and wealth, you diminish the total output and wealth of your society. While some balance is needed, a redistribution that aims at identical outcomes will result in total poverty.

This is a variant of the “two-envelopes problem,” originally posed in 1912 by German, Jewish mathematician, Edmund Landau. It is described, with related problems, by Prakash Gorroochurn, Classic Problems of Probability. Wiley, 314pp. ISBN: 978-1-118-06325-5. Wikipedia article: Two Envelopes Problem.

Robert Buxbaum, February 27, 2019

Great waves, small circles, and the spread of ideas.

Simplified wave motion, GIf by Dan Russel (maybe? I think?).

The scientific method involves looking closely at things. Sometimes we look closely for a purpose — to make a better mouse-trap, say. But sometimes it’s just to understand what’s happening: to satisfy curiosity, to understand the way the world works, or to answer a child. Both motivations bring positive results, but there is a difference in how people honor the product of these motivations. Scientific knowledge developed for curiosity is considered better; it tends to become the model for social understanding, and for art and literature. Meanwhile, science developed for a purpose is considered suspect, and often that suspicion is valid. A surprising amount of our knowledge was developed for war: for the purpose of killing people, destroying things, and occupying lands.

Waves provide a wonderful example of science exploration that was developed mostly for curiosity, and so they have become models of social understanding and culture — far more so than the atom bomb and plague work discussed previously.

Waves appear magical: You poke a pond surface with a stick, and the influence of that poke travels, as if by magic, to all corners of the pond. Apparently the initial poke set off something, and that sets off something else, and we’ve come to use this as a model for cultural ideas. Any major change in music, art, or cultural thought is described as a wave (and not as a disease). The sense of wave is  that a small push occurs, and the impact travels across a continent and across an ocean. The Gifs above and below shows how this happens for the ordinary wave — the one with a peaked top. As shown, the bits of water do not move with the wave. Instead they just circulate in a small circle. The powerful waves that crosses an ocean are composed of many small circles of water rolling in the general direction of the wave. With ideas too, I think, one person can push a second, and that second a third, each acting in his or her own circle, and a powerful transmission of ideas results. Of course, for a big wave, you need a big circle, but maybe not in cases of reflection (reflected waves can add, sometimes very destructively).

simplified wave movement

In the figures I’ve shown, you will notice that the top of the circle always moves in the same direction as the top of the wave. If the wave moves to the right, the circle is clockwise. There are also Rayleigh waves. In these, the top of the wave is not peaked, but broad, with little indents between ripples. For Rayleigh wave the motion is not circular, but elliptical, and the top of the ellipse moves in the opposite direction to that of the wave. These waves go slower than the normal waves, but they are more destructive. Most of the damage of earthquakes is by the late-arriving Rayleigh waves.

If regular waves are related to fast-moving ideas, like rock n roll, Rayleigh waves might be related to slower-traveling, counter-intuitive ideas, paradigm shifts: Religions, chaos, entropyfeminism, or communism. Rayleigh waves are mostly seen in solids, and the destructive power of counter-intuitive ideas is mostly seen in rigid societies.

Then there are also pressure waves, like sound, and wiggle waves (transverse waves). Pressure waves travel the fastest, and work in both solids and liquids. Wiggle waves travel slower (and don’t travel in liquids). Both of these involve no circles at all, but just one bit of material pushing on its neighbor. I think the economy works this way: bouncing springs, for the most part. Life is made up of all of these, and life is good. The alternative to vibration, I should mention, is status. Status is a form of death. There is a certain sort of person who longs for nothing more than an unchanging, no-conflict world: one government and one leadership. Avoid such people.

Robert Buxbaum, February 10, 2019