Sandbox: Difference between revisions

From ChanceWiki
Jump to navigation Jump to search
No edit summary
 
Line 1: Line 1:
==Autism Statistics Lesson==


Autism is a devastating disease.  Recent attempts to pinpoint a cause have been recently in the news in the United States: [http://www.nytimes.com/2009/02/13/health/13vaccine.html?scp=1&sq=donald%20g.%20mcneil%20autism&st=cse New York Times article] by Donald G. McNeil, Jr.] and [http://www.nytimes.com/2009/02/13/opinion/13fri2.html?scp=1&sq=editorial%20vaccines&st=cse a follow-on New York Times editorial].  The focus of these articles is on law suits regarding the measles, mumps and rubella (MMR) vaccine “or its combination with thimerosal, a mercury-based preservative that was used in most childhood vaccines until 2001,” as a cause of autism.  After “5000 pages of testimony from experts and 939 medical articles,” judges concluded the plaintiffs failed to prove their assertions.  One judge “ruled that the evidence was ‘overwhelmingly contrary’ to their argument.” 


Coincidentally, in England autism was also in the news in February, [http://briandeer.com/mmr/lancet-summary.htm here] and [http://briandeer.com/mmr-lancet.htm here].  In this instance, the story begins back in 1998 [http://download.videohelp.com/vitualis/downloads/Wakefield_%20LancetVolume%20351(9103)February28-1998.pdf The Lancet, February, 1998] and ignores thimerosal but introduces a problem additional to autism due to the MMR vaccine, Crohn’s disease (inflammatory bowel disease).  The Lancet article had an extraordinary impact on the general public in England as the following graph indicates:
==Forsooth==
 
==Quotations==
“We know that people tend to overestimate the frequency of well-publicized, spectacular
events compared with more commonplace ones; this is a well-understood phenomenon in
the literature of risk assessment and leads to the truism that when statistics plays folklore,
folklore always wins in a rout.”
<div align=right>-- Donald Kennedy (former president of Stanford University), ''Academic Duty'', Harvard University Press, 1997, p.17</div>
 
----
 
"Using scientific language and measurement doesn’t prevent a researcher from conducting flawed experiments and drawing wrong conclusions — especially when they confirm preconceptions."
 
<div align=right>-- Blaise Agüera y Arcas, Margaret Mitchell and Alexander Todoorov, quoted in: The racist history behind facial recognition, ''New York Times'', 10 July 2019</div>
 
==In progress==
[https://www.nytimes.com/2018/11/07/magazine/placebo-effect-medicine.html What if the Placebo Effect Isn’t a Trick?]<br>
by Gary Greenberg, ''New York Times Magazine'', 7 November 2018
 
[https://www.nytimes.com/2019/07/17/opinion/pretrial-ai.html The Problems With Risk Assessment Tools]<br>
by Chelsea Barabas, Karthik Dinakar and Colin Doyle, ''New York Times'', 17 July 2019
 
==Hurricane Maria deaths==
Laura Kapitula sent the following to the Isolated Statisticians e-mail list:
 
:[Why counting casualties after a hurricane is so hard]<br>
:by Jo Craven McGinty, Wall Street Journal, 7 September 2018
 
The article is subtitled: Indirect deaths—such as those caused by gaps in medication—can occur months after a storm, complicating tallies
   
   
Laura noted that
:[https://www.washingtonpost.com/news/fact-checker/wp/2018/06/02/did-4645-people-die-in-hurricane-maria-nope/?utm_term=.0a5e6e48bf11 Did 4,645 people die in Hurricane Maria? Nope.]<br>
:by Glenn Kessler, ''Washington Post'', 1 June 2018


MMR inoculation rates fall off sharply after the Lancet article and start to rise in 2004 because of a (London) Sunday Times investigation which revealed serious deficiencies in the Lancet studyThese deficiencies often fall under the rubric of “follow the money,” a concept not given enough attention when discussing what constitutes statistical literacy.
The source of the 4645 figure is a [https://www.nejm.org/doi/full/10.1056/NEJMsa1803972 NEJM article]Point estimate, the 95% confidence interval ran from 793 to 8498.


The phrase, “follow the money,” is often thought to have originated in the book, All the President’s Men.  According to [http://newsmine.org/content.php?ol=deceptions/deep-throat/follow-the-money-phrase-written-by-princess-bride-author.txt Frank Rich] the book never uses that phrase.  It is however, from the film of the same name.  Obviously, the pharmaceutical industry has a vested financial interest in vaccines and [http://www.huffingtonpost.com/deirdre-imus/on-vaccinations-consider_b_165347.html?show_comment_id=20767201 Deirdre Imus] is suspicious of any “big pharma” vaccine and any doctor who sides with it. The main author of the Lancet article, Andrew Wakefield, unbeknownst to the twelve other authors of the Lancet study, had been paid “about $780,000 plus expenses, for his role in backing the generic case against MMR.”  Further, he had a patent on “a single vaccine against measles—a potential competitor to MMR” which he claimed would cure “both inflammatory bowel disease and autism.”  As cited by [http://www.washingtonpost.com/wp-dyn/articles/A41450-2004Jul10.htm Glenn Frankel] ten of the twelve other authors in 2004 issued a “Retraction of an interpretation” because “no causal link was established between MMR vaccine and autism. Wakefield has since moved to the U.S. and according to a supporter of Wakefield, “The United States, with its privatized health care system and entrepreneurial spirit is much more fertile ground than Britain for a medical pioneer like Wakefield.”  According to [http://www.timesonline.co.uk/tol/life_and_style/health/article5683671.ece Brian Deer] official figures showed that 1,348 confirmed cases of measles in England and Wales were reported last year [2008], compared with 56 in 1998. Two children have died of the disease.”
President Trump has asserted that the actual number is
[https://twitter.com/realDonaldTrump/status/1040217897703026689 6 to 18].
The ''Post'' article notes that Puerto Rican official had asked researchers at George Washington University to do an estimate of the death tollThat work is not complete.
[https://prstudy.publichealth.gwu.edu/ George Washington University study]


Discussion
:[https://fivethirtyeight.com/features/we-still-dont-know-how-many-people-died-because-of-katrina/?ex_cid=538twitter We sttill don’t know how many people died because of Katrina]<br>
:by Carl Bialik, FiveThirtyEight, 26 August 2015


1.  If the medical profession overwhelmingly believes the MMR vaccine to be safe, why are parents of autistic children actively seeking litigation? That is, what element of emotional guilt might there be?
----
[https://www.nytimes.com/2018/09/11/climate/hurricane-evacuation-path-forecasts.html These 3 Hurricane Misconceptions Can Be Dangerous. Scientists Want to Clear Them Up.]<br>
[https://journals.ametsoc.org/doi/abs/10.1175/BAMS-88-5-651 Misinterpretations of the “Cone of Uncertainty” in Florida during the 2004 Hurricane Season]<br>
[https://www.nhc.noaa.gov/aboutcone.shtml Definition of the NHC Track Forecast Cone]
----
[https://www.popsci.com/moderate-drinking-benefits-risks Remember when a glass of wine a day was good for you? Here's why that changed.]
''Popular Science'', 10 September 2018
----
[https://www.economist.com/united-states/2018/08/30/googling-the-news Googling the news]<br>
''Economist'', 1 September 2018
 
[https://www.cnbc.com/2018/09/17/google-tests-changes-to-its-search-algorithm-how-search-works.html We sat in on an internal Google meeting where they talked about changing the search algorithm — here's what we learned]
----
[http://www.wyso.org/post/stats-stories-reading-writing-and-risk-literacy Reading , Writing and Risk Literacy]
 
[http://www.riskliteracy.org/]
-----
[https://twitter.com/i/moments/1025000711539572737?cn=ZmxleGlibGVfcmVjc18y&refsrc=email Today is the deadliest day of the year for car wrecks in the U.S.]
 
==Some math doodles==
<math>P \left({A_1 \cup A_2}\right) = P\left({A_1}\right) + P\left({A_2}\right) -P \left({A_1 \cap A_2}\right)</math>
 
<math>P(E)  = {n \choose k} p^k (1-p)^{ n-k}</math>
 
<math>\hat{p}(H|H)</math>
 
<math>\hat{p}(H|HH)</math>
 
==Accidental insights==
 
My collective understanding of Power Laws would fit beneath the shallow end of the long tail. Curiosity, however, easily fills the fat endI long have been intrigued by the concept and the surprisingly common appearance of power laws in varied natural, social and organizational dynamics.  But, am I just seeing a statistical novelty or is there meaning and utility in Power Law relationships? Here’s a case in point.
 
While carrying a pair of 10 lb. hand weights one, by chance, slipped from my grasp and fell onto a piece of ceramic tile I had left on the carpeted floor. The fractured tile was inconsequential, meant for the trash.
<center>[[File:BrokenTile.jpg | 400px]]</center>
As I stared, slightly annoyed, at the mess, a favorite maxim of the Greek philosopher, Epictetus, came to mind: “On the occasion of every accident that befalls you, turn to yourself and ask what power you have to put it to use.”  Could this array of large and small polygons form a Power Law? With curiosity piqued, I collected all the fragments and measured the area of each piece.
 
<center>
{| class="wikitable"
|-
! Piece !! Sq. Inches !! % of Total
|-
| 1 || 43.25 || 31.9%
|-
| 2 || 35.25 ||26.0%
|-
|  3 || 23.25 || 17.2%
|-
| 4 || 14.10 || 10.4%
|-
| 5 || 7.10 || 5.2%
|-
| 6 || 4.70 || 3.5%
|-
| 7 || 3.60 || 2.7%
|-
| 8 || 3.03 || 2.2%
|-
| 9 || 0.66 || 0.5%
|-
| 10 || 0.61 || 0.5%
|}
</center>
<center>[[File:Montante_plot1.png | 500px]]</center>
The data and plot look like a Power Law distribution. The first plot is an exponential fit of percent total area. The second plot is same data on a log normal format. Clue: Ok, data fits a straight line.  I found myself again in the shallow end of the knowledge curve. Does the data reflect a Power Law or something else, and if it does what does it reflect?  What insights can I gain from this accident? Favorite maxims of Epictetus and Pasteur echoed in my head:
“On the occasion of every accident that befalls you, remember to turn to yourself and inquire what power you have to turn it to use” and “Chance favors only the prepared mind.”
 
<center>[[File:Montante_plot2.png | 500px]]</center>
My “prepared” mind searched for answers, leading me down varied learning paths. Tapping the power of networks, I dropped a note to Chance News editor Bill Peterson. His quick web search surfaced a story from ''Nature News'' on research by Hans Herrmann, et. al. [http://www.nature.com/news/2004/040227/full/news040223-11.html Shattered eggs reveal secrets of explosions].  As described there, researchers have found power-law relationships for the fragments produced by shattering a pane of glass or breaking a solid object, such as a stone. Seems there is a science underpinning how things break and explode; potentially useful in Forensic reconstructions.
Bill also provided a link to [http://cran.r-project.org/web/packages/poweRlaw/vignettes/poweRlaw.pdf a vignette from CRAN] describing a maximum likelihood procedure for fitting a Power Law relationship. I am now learning my way through that.


2.  Although there were 13 authors of the Lancet article, there were only twelve children in the study.  A multiplicity of authors is a common phenomenon in medical journals.  Why is this so?
Submitted by William Montante


3.  The Lancet article claims that in eight of the twelve children, “the average exposure to first behavioral symptoms was 6.3 days (range 1-14)” after receiving the MMR vaccine.  An earlier version of the paper, not unearthed until 2005, puts the average at 14 days with the maximum time as 56 days.  Further, it was later revealed that there was “no trace of measles virus [or mumps and rubella viruses] in any of the children.”  Subsequent investigation indicated that instead of Crohn’s disease, the children were suffering from a benign condition, severe constipation.  Moreover, the children were not randomly referred by general practitioners but were recruited from a lawyer “who had been attempting to raise a speculative lawsuit.”  In 2007, Wakefield abandoned a libel claim, and agreed to pay costs, “estimated at about £500,000.”  Assuming all of this is factually correct, explain why some parents still view Wakefield as a hero.
----

Latest revision as of 20:58, 17 July 2019


Forsooth

Quotations

“We know that people tend to overestimate the frequency of well-publicized, spectacular events compared with more commonplace ones; this is a well-understood phenomenon in the literature of risk assessment and leads to the truism that when statistics plays folklore, folklore always wins in a rout.”

-- Donald Kennedy (former president of Stanford University), Academic Duty, Harvard University Press, 1997, p.17

"Using scientific language and measurement doesn’t prevent a researcher from conducting flawed experiments and drawing wrong conclusions — especially when they confirm preconceptions."

-- Blaise Agüera y Arcas, Margaret Mitchell and Alexander Todoorov, quoted in: The racist history behind facial recognition, New York Times, 10 July 2019

In progress

What if the Placebo Effect Isn’t a Trick?
by Gary Greenberg, New York Times Magazine, 7 November 2018

The Problems With Risk Assessment Tools
by Chelsea Barabas, Karthik Dinakar and Colin Doyle, New York Times, 17 July 2019

Hurricane Maria deaths

Laura Kapitula sent the following to the Isolated Statisticians e-mail list:

[Why counting casualties after a hurricane is so hard]
by Jo Craven McGinty, Wall Street Journal, 7 September 2018

The article is subtitled: Indirect deaths—such as those caused by gaps in medication—can occur months after a storm, complicating tallies

Laura noted that

Did 4,645 people die in Hurricane Maria? Nope.
by Glenn Kessler, Washington Post, 1 June 2018

The source of the 4645 figure is a NEJM article. Point estimate, the 95% confidence interval ran from 793 to 8498.

President Trump has asserted that the actual number is 6 to 18. The Post article notes that Puerto Rican official had asked researchers at George Washington University to do an estimate of the death toll. That work is not complete. George Washington University study

We sttill don’t know how many people died because of Katrina
by Carl Bialik, FiveThirtyEight, 26 August 2015

These 3 Hurricane Misconceptions Can Be Dangerous. Scientists Want to Clear Them Up.
Misinterpretations of the “Cone of Uncertainty” in Florida during the 2004 Hurricane Season
Definition of the NHC Track Forecast Cone


Remember when a glass of wine a day was good for you? Here's why that changed. Popular Science, 10 September 2018


Googling the news
Economist, 1 September 2018

We sat in on an internal Google meeting where they talked about changing the search algorithm — here's what we learned


Reading , Writing and Risk Literacy

[1]


Today is the deadliest day of the year for car wrecks in the U.S.

Some math doodles

<math>P \left({A_1 \cup A_2}\right) = P\left({A_1}\right) + P\left({A_2}\right) -P \left({A_1 \cap A_2}\right)</math>

<math>P(E) = {n \choose k} p^k (1-p)^{ n-k}</math>

<math>\hat{p}(H|H)</math>

<math>\hat{p}(H|HH)</math>

Accidental insights

My collective understanding of Power Laws would fit beneath the shallow end of the long tail. Curiosity, however, easily fills the fat end. I long have been intrigued by the concept and the surprisingly common appearance of power laws in varied natural, social and organizational dynamics. But, am I just seeing a statistical novelty or is there meaning and utility in Power Law relationships? Here’s a case in point.

While carrying a pair of 10 lb. hand weights one, by chance, slipped from my grasp and fell onto a piece of ceramic tile I had left on the carpeted floor. The fractured tile was inconsequential, meant for the trash.

BrokenTile.jpg

As I stared, slightly annoyed, at the mess, a favorite maxim of the Greek philosopher, Epictetus, came to mind: “On the occasion of every accident that befalls you, turn to yourself and ask what power you have to put it to use.” Could this array of large and small polygons form a Power Law? With curiosity piqued, I collected all the fragments and measured the area of each piece.

Piece Sq. Inches % of Total
1 43.25 31.9%
2 35.25 26.0%
3 23.25 17.2%
4 14.10 10.4%
5 7.10 5.2%
6 4.70 3.5%
7 3.60 2.7%
8 3.03 2.2%
9 0.66 0.5%
10 0.61 0.5%
Montante plot1.png

The data and plot look like a Power Law distribution. The first plot is an exponential fit of percent total area. The second plot is same data on a log normal format. Clue: Ok, data fits a straight line. I found myself again in the shallow end of the knowledge curve. Does the data reflect a Power Law or something else, and if it does what does it reflect? What insights can I gain from this accident? Favorite maxims of Epictetus and Pasteur echoed in my head: “On the occasion of every accident that befalls you, remember to turn to yourself and inquire what power you have to turn it to use” and “Chance favors only the prepared mind.”

Montante plot2.png

My “prepared” mind searched for answers, leading me down varied learning paths. Tapping the power of networks, I dropped a note to Chance News editor Bill Peterson. His quick web search surfaced a story from Nature News on research by Hans Herrmann, et. al. Shattered eggs reveal secrets of explosions. As described there, researchers have found power-law relationships for the fragments produced by shattering a pane of glass or breaking a solid object, such as a stone. Seems there is a science underpinning how things break and explode; potentially useful in Forensic reconstructions. Bill also provided a link to a vignette from CRAN describing a maximum likelihood procedure for fitting a Power Law relationship. I am now learning my way through that.

Submitted by William Montante