Sunday, August 14, 2011

martial arts, global warming, and chance.

Scientific knowledge is produced piece meal through the slow process of research. One particular challenge in this process occurs when a researcher observes a difference between the experimental group and the control group. Did the difference arise because of chance or because of the experimental condition? Consider a child psychologist who is attempting to discern the most effective method for modifying explosive behaviors in children (say, ages 5 to 8). She designs a study in which the control group receives standard family counseling intervention while the treatment group receives standard family counseling plus Aikido martial arts training.
Suppose that the counseling + Aikido group exhibits a greater average reduction in explosive behaviors than the counseling alone group. What should the psychologist conclude? Was the difference really due to martial arts training? Or was there simply random variation between the two groups?

The observable world (a.k.a. the phenomenal or that which is accessible by our five senses) is characterized by variation. The world does have unifying undercurrents, however on the whole a vast multitude of divergent atomic structures abounds. In short, all things are not the same. There are alkali metals, alkaline earth metals, inert transition metals, transition metals and post-transition metals. There are tall people, short people, fat people, skinny people, brown skinned people, olive skinned people, red haired, brown eyed, toothy and toothless. Botanically, there are at least 40 different types of daffodils; each type of daffodil itself has a wide range of variable characteristics. This variation is what statisticians refer to when they use the terms "random error" or "random variation." On the other hand, variation due to some factor (e.g.: medication, exercise, psychotherapy, air pollution, etc) is referred to as explained variation. For example, consider an experiment which showed that lower carbon pollution levels curiously cause increased average temperatures in a simulated atmosphere.
Researchers want to know if this temperature variation is explained by reduced carbon pollution or by error (random variation). A high degree of variation in temperature recordings collected during the experiment may indicate that the average increase in temperature was a chance finding. In other words, if the variation in temperatures is great, then it may be that the researcher simply collected data during a time when temperatures were varying on the high end of the thermometer. Below I will show how the statisticians would answer the question, "Is the temperature rise due to the decreased carbon pollution or due to chance variation?"

When analyzing data sets from an experiment/study, researchers want to know if the observed change is a chance finding or a true reflection of how things really are. To do this they use a statistic called standard error (a quotient of standard deviation divided by the square root of the study sample size) to predict how variable the sample average would be if the experiment were repeated a large number of times. Essentially, the standard error statistic is used to calculate a range of values (termed "confidence interval") in which the researcher can be 95% confident that the true population average lies. Consider the atmospheric temperature study above. If the average temperature increased by 3.0 degrees Celsius, but the 95% confidence interval was -0.3 to 8.2 then it is entirely possible that no temperature change happened whatsoever (note that the value "0" is included in the range...indicating that the true average change in temperature may be 0 degrees Celsius). Such a scenario would indicate that the increase in temperature may truly be due to random variations or chance, not due to reduced carbon pollution.
Further testing can be done viz a viz Fisher's all powerful p-value to determine how insignificant the temperature rise truly is. In fact, selecting levels of significance and calculation of p-values will be the topic of next post. From there I will switch out of statistics mode and into philosophy mode.

Saturday, August 13, 2011

descriptive stats revisited

Last post I attempted a brief and somewhat unclear introduction to descriptive statistics. So by way of review I will summarize last post with a graph to help depict things. The three most basic ways of describing a data set are its average, spread and shape. These three concepts are very important as performing statistical analysis depends on these three descriptions of the data being "normal." The average of the data set represents where the majority of the data lie. A normally distributed set of data (referred to as a bell curve when plotted) will have an obvious mean at the peak of that curve. See how on the normal curve below, the line representing the mean divides the curve into two symmetric halves.

The spread of the data indicates to what degree all the data points vary away from the average. Standard deviation is a number that indicates the average distance of all data points from the average. Higher standard deviations indicate that the range of data spread further on either side of the average. See the lines on either side of the central average line on the graph below? The first line on either side of the average line indicates a point that is 1 standard deviation from the mean. The next line indicates a point that is 2 standard deviations away from the mean. Particular to the normal bell curve is the reality that 67% of all data will be only one standard deviation from the mean and 95% of the data only 2 standard deviations away from the mean.

The shape indicates whether the bell curve has a long tail to the left [negative skew], long tail to the right [positive skew], tall curve [leptokurtic], flat cruve platykurtic). Any distortion of the bell curve impacts statistical calculations as basic statistics assume a normal bell curve. See below for depictions of skewed data.

Next post I will answer the following question: how do we know if the difference between two groups or two measurements is due to chance or due to skill or treatment effect?

Thursday, August 4, 2011

redlegs and descriptive stats

Honestly, my appreciation of baseball is pathetic. That is a veritable reality (as those acquainted with me can attest). However, as a youth I followed Redlegs baseball compulsively. Sunday mornings were great because the Cincinnati Enquirer would have an expanded sports section with all the breakdown of player and team stats (this was before the dawn of ESPN.com/mlb/statistics....in a former age when all fans waited with baited breath for our information to arrive through the paper delivery man). Batting averages, ERA's, standings, etc. Typical of my analytic approach to things, I would perform my own calculations to see how many hits Hal Morris would need to break .400 over the next 50 at bats. Much has changed since my youth, but I still gain a good deal of satisfaction in performing calculations. For those who might find such a thought foreign or even perverse, I will frame this post on summary statistics around baseball. Perhaps no one will even notice that math is being conducted.

Data collection is the process of compiling bits of information (e.g.: number of hits, innings pitched) in an organized system (I personally prefer Excel for its ubiquitous use). Without doing at least some basic calculations the data are not very meaningful. Summary statistics are employed to describe several qualities of the data: location, spread, and shape. Location captures where the majority of the data lie (e.g.: the team batting average represents a central number around which most of the players on a baseball team are averaging). Spread captures how widely the data varies (e.g.: once the location is known we want to know how big the difference is between the worst batting average and the best batting average). Shape captures how closely a graph of the data resembles a bell shaped curve (e.g.: if the batting average is relatively high compared to the spread of individual batting averages then the shape will be abnormally tall and departs from a normal bell curve). In statistical jargon these three qualities are referred to as: mean, variance and skewness/kurtosis. (WARNING: those who experience anxiety when exposed to mathematical formulas should skip over the next paragraph...you won't miss much).



Mean is the sum of the data divided by the number of data points and is denote by the upper case "Y" with a line over it in the equation:
    • YBAR = SUM[Y(i)/N] where the summation is for 1 to N.
  • Variance is the square of each data points distance from the mean, all summed and divided by total number of data points minus 1. The square root of variance can then be taken to reflect the average distance that data points vary from the mean

      s = SQRT[SUM(Y(i) - YBAR)**2/(N-1)] where the summation is from 1 to N.

    Skewness is obtained by the cubed value of (the mean subtracted from each data point), summed for all data points and then divided by the product of (total number of data points minus 1) and (the cube of standard deviation). skewness = SUM(Y(i) - YBAR)**3/(N-1)*s**3 where the  summation is from 1 to N When this number is significantly positive then the data is not normally shaped and is lopsided to the right; when this number is negative then the data is lopsided to the left. Kurtosis is obtained similarly to skewness only the numerator is taken to the 4th power and standard deviation in the denominator is taken to the 4th power
    kurtosis = SUM(Y(i) - YBAR)**4/(N-1)*s**4 where the  summation is from 1 to N. A significantly negative kurtosis number indicates a flat curve and a positive kurtosis number indicates a tall, lanky curve.


    Why are descriptive statistics important? There is a phenomenon out there in the sports world known as fanaticism. A true baseball fanatic will go to extreme lengths to stand by her team. It is even conceivable that during pleasant interchange over brews she may overstate or even understate the skill level of her team. In the world of statistics we call this bias. Descriptive statistics help us to reduce the influence of bias and see things more like they really are. No matter how devoted I am to the Redlegs offense, if I compare a variety of reliable offensive indicators to the same indicators for the Cardinals, I will find out the truth of how well they compare. If we compare the two teams and find that the Redlegs have better numbers and that the difference between the two teams is not due to chance, then I may not only brag as a fanatic but as a man who knows the statistical truth ("statistical truth," is that an oxymoron?).

    So that raises a question, how would I know if the better Redlegs numbers were due to chance or skill? Tune in to next post for a breathtaking answer.

    Saturday, July 30, 2011

    ents

    Last post I mentioned that this post would be a technical one on statistics. However, after studying further on the subject of statistics I quickly realized that I am not qualified to author a technical post on the matter. Instead I will here offer a post that works to hold together statistics and scientific epistemology. Think of scientific epistemology being somewhat like Tolkien's tree Ents that he depicted as gathering around in a counsel to decide on a course of action. They refuse to be rushed in their decision making by the immediacy of danger: rather they thoroughly work the problem over to ensure that their decision places them on the side of good. Likewise with scientific epistemology, establishment of true knowledge occurs slowly and without any deference to the urgent need of human beings. So here ensues an elaboration of Ent-like knowledge building interspersed with pictures of past heroes of statistics (everybody should choose at least one statistical hero...such a thing ought to be as ubiquitous as having a NASCAR driver to root for).

    I believe that natural human knowledge is capable of more than a mere approximation of reality. I believe that natural human knowledge can achieve a one-to-one relationship with reality; that is to say that theories of reality (i.e.: how things are) can represent specified slices of reality with perfect accuracy (e.g.: Copernicus and heliocentrism; Newton and gravitational forces, etc.). However, I also believe that the vast majority of natural knowledge produced in a given calendar year is specious at worst (i.e.: false knowledge) and an approximation of reality at best (i.e.: analogous knowledge). Here I refer not to the sort of assertions that pass off as knowledge among the popular pundits, but rather I am referring to critically tested, peer-reviewed scientific knowledge. It is this high level of knowledge that ranges from false to analogous. To characterize this range of knowledge quality I have generated a list of four (4) descriptors beginning with the most basic and contingent sort of knowledge progressing up to true knowledge or "fact" (as science would have us call it).

    Ronald Fisher: what's not to like about this statistician? The very incarnation of an Ent.

    Possible. Definition: that can be; capable of existing. For example an oval circle is quite possible; an oblong circle is quite possible. However, a square circle could never be a circle (i.e.: essentialistically impossible). An MRI machine that produces only rare cheese could never be an MRI machine (i.e.: nominalistically impossible). Certain things are impossible. Importantly, though, many possible things are not true. For instance, although it is possible for a live cow to go over the moon, if I claim to know that a live cow actually went over the moon last night this would be a false knowledge claim. Although it is possible for the sun to rise in the West and set in the East, this is not the truth of what happens on planet Earth. That which is conceivably possible is not necessarily existentially so. This being held true, there are three possible descriptions of possibilities and existence: existent possibility, nonexistent possibility, nonexistent impossibility. There is no such thing as an existent impossibility.

    Karl Pearson: note with what alacrity he wields his pen in noble service of statistics.

    Probable. Definition: likely to occur or to be so; that can reasonably be expected or believed on the basis of the available evidence, though not proved or certain. For example, a 2002 study published in Spine journal deduced a clinical prediction rule to identify which patients with low back pain are most likely to be successfully treated with spinal manipulation. The study deduced these 5 parts for a prediction rule: 1) pain onset less than 16 days prior to treatment; 2) no symptoms distal to the knee; 3) Fear Avoidance Behavior Questionnaire score less than or equal to 19; 4) one or more hypomobile segments in the lumbar spine; 5) at least one hip with more than 35 degrees of internal rotation motion. A patient must exhibit at least 4 of the above 5 traits to be considered positive on the rule. A subsequent randomized controlled trial published in 2004 compared patients who satisfy this rule with those who do not. They found that patients positive on the rule and receiving spinal manipulation have an adjusted odds ratio for successful treatment of 60.8 when compared to those negative on the rule and receiving exercise. That is to say, a person who satisfies the prediction rule and is treated with spinal manipulation is 60.8 times more likely to have a successful outcome than someone negative on the rule and receiving exercise. Although adjusted odds ratios are not exactly probabilities, they approximate mathematical probabilities closely when baseline risks are low. Consequently we can say that a person with low back pain who satisfies the prediction rule is very likely to gain significant function after only 2 treatments of spinal manipulation and this functional gain is likely to last for 6 months. But, this is a probability. This means that a small number of people who are positive on the rule and receive spinal manipulation will feel no improvement or worsening symptoms after treatment. So this is a practical example of probabilistic knowledge. The above studies did not elucidate fact, they elucidated pragmatic probabilities to help clinicians, insurance companies and patients make a decision about treatment options.

    Valid. Definition: sound; well grounded on principles or evidence; able to withstand criticism or objection, as an argument. The move from probable knowledge up to valid knowledge is like passing over the Natural Light beer for Kentucky Bourbon Barrel Ale. Like declining the Hostess Snack Cake to save room for Graeter's ice cream. Like operating on the spine with a high precision Medtronic drill instead of a high speed DeWalt house drill. Valid knowledge within medical diagnostics is produced when a new test for diagnosis is measured against a gold standard. Assessing the new test (lets say for detecting prostate cancer) against the gold standard can result in four possible outcomes: true positive, false positive, true negative, and false negative. The box below illustrates this well:

    Gold Standard Negative (no prostate cancer present)

    Gold Standard Positive (prostate cancer present)

    New Test Positive

    False Positive (invalid test result)

    True Positive (valid test result)

    New Test Negative

    True Negative (valid test result)

    False Negative (invalid test result)

    A variety of statistics are available to describe these results, however I will spare the reader details of these at this time. Let the reader note, though, that validity does not measure probability of departure from truth but rather it measures reality of departure from truth. In this way validity is not a guess at the truth but a true measure of the truth. The greatest mistake in popular conceptions of research occurs when studies analyzed by probabilistic statistics (i.e.: p-values and confidence intervals) are interpreted as if they utilize validity statistics (i.e.: likelihood ratios). In this manner, many studies are presented as valid when they ought truthfully be presented as some grade (e.g.: low, moderate, high grade) of probabilistic evidence for or against a hypothesis.

    Thomas Bayes: the right Reverend had a penchant for more than divine truth.

    Veritable. Definition: true; real; actual. I am taking the name of this descriptor from the Latin word veritas: truth. How does the descriptor "veritable" differ from the descriptor "valid"? Knowledge that is veritable is perfectly true (1-for-1 correspondence with reality or logic) whereas knowledge that is described by validity is expressed as a degree of departure from truth. Given the above diagnostic example: the new prostate cancer test would be characterized as being valid if it deviates from the gold standard to only a small degree and invalid if it deviates from the gold standard to a large degree. However, the gold standard test is veritable: a container of 100% true knowledge, finding prostate cancer whenever it exists and ruling out prostate cancer whenever it does not exist. Veritable knowledge is not simply produced by a mathematical formula, but is rather arrived at through years of practice, research (statistically analyzed), critical review (e.g.: statistical assessment, construct assessment, etc.), technical reformulation, and even large scale disciplinary enactment (i.e.: use within the guild at large).

    Jerzy Neyman: a fierce opponent to Fisher; his calculations were as killer as his Hitler-ish look.

    One final cautionary note is due here. Having moved through this post, and engaged with the proposed hierarchical categories of knowledge (low quality "possibility" knowledge up to highest quality "veritable" knowledge) the reader might be tempted to scorn knowledge that is beneath "veritable" in my hierarchical scheme. The purpose of these hierarchies is not to cast aspersion upon "inferior" knowledge, but rather to appropriately characterize the levels of knowledge and their relation to truth. This has very pragmatic implications. Possible knowledge should be utilized with extreme caution and proposed lightly with a ready willingness to believe alternative possibilities; probable knowledge should be utilized discriminantly and defended publicly with openness to evidence to the contrary; valid knowledge should be ubiquitously employed and publicly defended with great vigor and energy; veritable knowledge should be published and re-published with regular frequency as well as being proclaimed with certainty to winsomely disabuse falsehoods consistently where falsehoods appear.

    Next post will be devoted to a description of basic statistics and some of their mathematical formulations.

    Saturday, July 23, 2011

    monochrome

    Two and one half years have ticked off the mantle top clock since the last post on this blog. Even more embarrassing is a swift look at the goals mentioned in that last post: publishing a case study in the national physical therapy journal (that didn't happen), become a researcher (not much traction on that over 2.5 years), becoming a voice that shapes physical therapy knowledge (I haven't even started talking yet). Alas and alack, my life was actually interesting over the last couple years: job changes, selling/buying a house, first born daughter arrived on the scene, theological growth, starting graduate classes again. No time for the monochrome world of research. I had some life to live. Life is settling down a little now. Hopefully a little more tranquility will invade my work life over the coming months and then I can really begin in earnest tackling scientific epistemology viz a vi statistics and clinical research.

    Interestingly enough, my introduction to biostatistics course this summer is a refreshing reminder that learning is fun. My instructors Drs. Shukla and Dwivedi are doing a superb job of instructing the class at a philosophical level. They are engaging questions of statistical certainty within the constraints of probabilistic mathematics (a world they refer to as fuzzy or gray) counterbalancing the clinical drive for a definite answer for ailing patients. The most unique concept I have learned to date is this: hypothesis testing (i.e.: level of significance and p-value) does not result in a mathematical acceptance or rejection of the hypothesis. Rather, hypothesis testing adds or removes weight to an a priori belief that the hypothesis is true. In this manner, statistics are legal procedures. The test statistic is the prosecuting trial attorney arguing the state's case in front of the judge and jury. The data is considered insignificant unless otherwise proven beyond reasonable doubt by the test statistic.

    Next post will be a technical one exploring some aspect of statistics and how they relate to contemporary medical knowledge.

    Saturday, February 14, 2009

    Towards publishing...

    Adventures excite me. Being a thirty year old white male in America means that any adventure I experience is usually between my ears. My adventure is an intellectual one, no doubt. It is one that I find difficult to share with others. While I might take the nuances of epistemology to be fascinating, most others have heralded it it as the great cure for insomnia. Nonetheless, this blog page is my feeble attempt at bringing my adventure to others in my world.

    I am closer to publishing than I was one year ago. I have written over ten papers and am positioning myself to write a case study and attempt to publish it in the national PT journal. This in itself is encouraging.

    Recently I have writen several essays for graduate school. The one recurring theme in my essays is my desire to become a researcher. In this career I hope to be able to contribute to the body of knowledge in physical therapy. Research is not an end point in my desire, though. Perhaps more motivating for me is a desire to become a voice that shapes the way in which physical therapy knowledge is intuited, refined, and published as truth.

    I would like to follow in the foot steps of Michael Polanyi. Polanyi was a Hungarian physician who went on to become a physical chemist and then philospher. During his career he worked brilliantly to shape scientific epistemology. Through his work he outlined the personal scope of knowledge; he showed that learning knowledge is not some mechanical process, but is rather a personal process. Gaining knowledge involves everything from epiphany to mentorship to critique. And that's what I want to bring to the table in my writing.

    Saturday, September 13, 2008

    My how time flies...

    It has been over a month since I last blogged. Time to catch up briefly.

    I never did read von Bertalanffy's book. Instead I read and analyzed a bunch of California and Ohio supreme court rulings about duty to warn. For my term paper I chose a California ruling dubbed the Tarasoff rule. This rule holds a mental health professional liable in a situation where a patient confides his intention to harm a third party (that is legal speak for someone with whom the professional has no relationship). The professional is liable only if: 1) the patient ends up harming the third party and; 2) the professional did not warn the third party. I supplemented my case law (case law is court established law) analysis with some statutory law (statutory law is legislatively established law) analysis. I then made some applications for physical therapists. I hope to publish the paper.

    I'm curious what my vast readership thinks. Consider a scenario where Bob Jones is being treated by Jack Handy (psychologist) for depression after a rough break-up with his girlfriend. What is the psychologist's duty when Bob makes specific claims that he is "going to bring a knife over to her house and end it all?" According to traditional Anglo-American common law there is no duty to protect between strangers. However, the California Supreme Court made an exception to the common law, stating that when a psychologist enters into a relationship with a patient he then has a duty to protect the public from the patient.

    What do you think?