inference
https://robweiss.faculty.biostat.ucla.edu/
enKilling the P-Value Messenger
https://robweiss.faculty.biostat.ucla.edu/blog/killing-p-value-messenger
<span>Killing the P-Value Messenger</span>
<div><p>Does the argument about p-values miss the big problem? One might argue that blaming p-values for our mis-uses of p-values is akin to blaming the messenger for bad news.<br />
How do we use p-values? Prior to publishing, we use p-values as gate keepers. A small p-value indicates a 'significant' and therefore 'real' result, while large p-values indicate that a particular effect is not there or not worth mentioning. Researchers emphasize their results that have small p-values, while ignoring or not mentioning results with large p-values. Often results with large p-values are omitted from a paper. Having at least one small p-value is considered de rigueur to bother submitting a paper, and a paper with no significant results will either be reworked until<br />
some p-value is small (torturing the data until the data speaks) or the paper will be laid to rest without being submitted (the file drawer problem).</p>
<p>It seems to me the problem isn't the p-value. The problem isn't even that a lot of people don't quite understand the underlying logic of what a p-value is. The problem is that p-values are used to determine which results are presented to the broader scientific community and which results are ignored. But when we read a result<br />
in an article in the scientific literature, we typically assume the canonical (but false) idea that the authors wanted to test that specific result or hypothesis, and then reported that specific result or hypothesis.</p>
<p>The truth is that there are quite a few results that authors could report in most studies. Authors choose the results that are actually reported in the paper. And this choosing is closely correlated with the p-value. This is commonly called selection bias.</p>
<p>To over-simplify slightly, there are two ways to get a significant p-value. One way is to have an underlying effect that is strong with a study that has the power to identify that strong effect. Another way is to have an unusually strong result that isn't warranted by the underlying truth. The first way is that there is an actual effect, and the second way is 'something unusual' occurred. Testing lots of hypotheses gets you lots of chances to get a significant result. With each test you may have identified something real, or you may have gotten lucky.</p>
<p>Unfortunately, when researchers get lucky, society gets unlikely. It's a situation where society's utility function and researchers' utility functions may be at odds with each other. Researchers with few exceptions do want to actually identify real effects. However, researchers operate in a milieu that rewards significant results. And especially, journal publishing rewards splashy and unusual results more than steady straightforward routine results or worse, non-results which means not-significant results. There are stories (of real researchers) getting wealthy by marketing modest (and probably lucky) but splashy findings into big businesses.</p>
<p>Scientists test a lot of hypotheses. Some of these tests are due to direct interest as when we assess a treatment's effect on a sample of subjects. Much of it is exploratory, as when we try to find the demographic characteristics that lead to higher values of some outcome. In addition, some tests of hypothesis are not straightforward, as when it requires exploratory analysis to build a model. When we build a complicated model to identify the effects of a given hypothesis, then we are subject to selection bias that occurs when we add model terms that distinguish a treatment's effects and we omit model terms that do not seem to distinguish treatment from control.</p>
<p>The problem is that, due to selection bias, our scientific literature is likely filled with results whose strength has been<br />
over-inflated. Now how over-inflated the typical result is depends on the underlying ground truth that is being explored in the literature. But when we have a large literature, and many researchers looking for many different effects, and few people checking results (it's less prestigious after all!) we're likely to have lots of candidate results that are not all correct but are all treated as real.</p>
<p>There have been a few attempts at validation studies, where researchers have gone back and tried to duplicate a number of past results in the literature. These validation studies don't validate all the past findings. Is this due to the original studies identifying false effects? Or is it because the follow-up study got unlucky and didn't manage to identify an actual effect? Perhaps a combination: the original study over estimated the effect size, and the validation study was under powered due to the original over estimation.</p>
<p>So should we toss out p-values? As a Bayesian, I made my peace with p-values a while ago. Now, two sided p-values are a bit weird and possibly a touch creepy from a Bayesian perspective. However, a one-sided p-value has a simple Bayesian interpretation that is actually useful. Suppose our estimate of a treatment effect is that treatment has a positive effect on the outcome. Then the Bayesian interpretation of a one-sided p-value (the one-side smaller than 0.5) is the probability, given the data, that the treatment effect in this study is actually harmful. It's a useful summary of the<br />
results of an analysis. Its not a sufficient statistic in either the English or statistical senses of the word with regard to how we should interpret the particular result.</p>
<p>It seems to me that the problem is more of selection. We select for strong results. We pay attention to strong results. We really pay attention to unusual and weird results. The popular press picks up unusual results and magnifies their impact immensely, without evaluating whether the reported result is likely to be true or not. Thus we propagate ideas based on something distinct from their likelihood of being true.</p>
<p>So don't shoot the messenger, neither messenger RNA nor messenger pidgeon. Oops too late, we already shot that last one. Consider how we can solve the crisis of scientific selection.</p>
</div>
<span><span lang="" about="https://robweiss.faculty.biostat.ucla.edu/user/10" typeof="schema:Person" property="schema:name" datatype="" xml:lang="">robweiss</span></span>
<span>Tue, 06/25/2019 - 16:48</span>
<div>
<div>Filed Under</div>
<div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/p-value" hreflang="en">p-value</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/inference" hreflang="en">inference</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/replication-crisis" hreflang="en">replication crisis</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/hypothesis-testing" hreflang="en">hypothesis testing</a></div>
</div>
</div>
<section></section>Tue, 25 Jun 2019 23:48:34 +0000robweiss194 at https://robweiss.faculty.biostat.ucla.eduWhy Be Bayesian? Let Me Count the Ways
https://robweiss.faculty.biostat.ucla.edu/blog/why-be-bayesian-let-me-count-ways
<span>Why Be Bayesian? Let Me Count the Ways</span>
<div><p>In answer to an old friend's question. </p>
<ol><li> Bayesians have more fun. </li>
<ol><li> Our conferences are in better places too. </li>
</ol><li> It's the model not the estimator. </li>
<li> Life's too short to be a frequentist: In an infinite number of replications ... </li>
<li> Software works better. </li>
<ol><li> Rather surprisingly, Bayesian software is a lot more general than frequentist software. </li>
</ol><li> Small sample inference comes standard with most Bayesian model fitting these days. </li>
<ol><li> But if you like your inference asymptotic, that's available, just not high on anyone's priority list. </li>
<li> We can handle the no-data problem, all the way up to very large problems. </li>
<li> Don't need a large enough sample to allow for a bootstrap. </li>
</ol><li> Hierarchical random effects models are better fit with Bayesian models and software. </li>
<ol><li> If a variance component is small, the natural Bayes model doesn't allow zero as an estimate, while the natural maximum likelihood algorithms do allow zero. If you get a zero estimate, then you're going to get poor estimates of standard errors of fixed effects. [More discussion omitted.] </li>
<li> Can handle problems where there are more parameters than data. </li>
</ol><li> Logistic regression models fit better with Bayes </li>
<ol><li> If there's perfect separation on a particular variable, the maximum likelihood estimate of the coefficient is plus or minus infinity which isn't a good estimate. </li>
<li> Bayesian modeling offers (doesn't guarantee it, there's no insurance against stupidity) the opportunity to do the estimation correctly. </li>
<li> Same thing if you're trying to estimate a very tiny (or very large) probability. Suppose you observe 20 out of 20 successes on something that you know doesn't have 100% successes. </li>
<li> To rephrase a bit: In small samples or with rare events, Bayesian estimates shrink towards sensible point estimates, (if your prior is sensible) thus avoiding the large variance of point estimates. </li>
</ol><li> Variance bias trade-off is working in your favor. </li>
<li> Frequentists keep reinventing Bayesian methods </li>
<ol><li> Shrinkage estimates </li>
<li> Empirical Bayes </li>
<li> Lasso </li>
<li> Penalized likelihood </li>
<li> Ridge regression </li>
<li> James-Stein estimators </li>
<li> Regularization </li>
<li> Pittman estimation </li>
<li> Integrated likelihood </li>
<li> In other words, it's just not possible to analyze complex data structures without Bayesian ideas. </li>
</ol><li> Your answers are admissible if you're Bayesian but usually not if you're a frequentist. </li>
<ol><li> Admissibility means never having to say you're sorry. </li>
<li> Alternatively, admissibility means that someone else can't prove that they can do a better job than you. </li>
<li> And if you're a frequentist, someone is clogging our journals with proofs that the latest idiocy is admissible or not. </li>
<li> Unless they are clogging it with yet more ways to estimate the smoothing parameter for a nonparametric estimator. </li>
</ol><li> Bayesian models are generalizations of classical models. That's what the prior buys you: more models </li>
<li> Can handle discrete, categorical, ordered categorical, trees, densities, matrices, missing data and other odd parameter types. </li>
<li> Data and parameters are treated on an equal playing field. </li>
<li> I would argue that cross-validation works because it approximates Bayesian model selection tools. </li>
<li> Bayesian Hypothesis Testing </li>
<ol><li> Treats the null and alternative hypotheses on equal terms </li>
<li> Can handle two or more than two hypotheses </li>
<li> Can handle hypotheses that are </li>
<ol><li> Disjoint </li>
<li> Nested </li>
<li> Overlapping but neither disjoint nor nested </li>
</ol><li> Gives you the probability the alternative hypothesis is true. </li>
<li> Classical inference can only handle the nested null hypothesis problem. </li>
<li> We're all probably misusing p-values anyway. </li>
</ol><li> Provides a language for talking about modeling and uncertainty that is missing in classical statistics. </li>
<ol><li> And thus provides a language for developing new models for new data sets or scientific problems. </li>
<li> Provides a language for thinking about shrinkage estimators and why we want to use them and how to specify the shrinkage. </li>
<li> Bayesian statistics permits discussion of the sampling density of the data given the unknown parameters. </li>
<li> Unfortunately this is all that frequentist statistics allows you to talk about. </li>
<li> Additionally: Bayesians can discuss the distribution of the data unconditional on the parameters. </li>
<li> Bayesian statistics also allows you to discuss the distribution of the parameters. </li>
<li> You may discuss the distribution of the parameters given the data. This is called the posterior, and is the conclusion of a Bayesian analysis. </li>
<li> You can talk about problems that classical statistics can't handle: The probability of nuclear war for example. </li>
</ol><li> Novel computing tools -- but you can often use your old tools as well. </li>
<li> Bayesian methods allow pooling of information from diverse data sources. </li>
<ol><li> Data can come from books, journal articles, older lab data, previous studies, people, experts, the horse's mouth, rats a** or it may have been collected in the traditional form of data. </li>
<li> It isn't automatic, but there is language to think about how to do this pooling. </li>
</ol><li> Less work. </li>
<ol><li> Bayesian inference is via laws of probability, not by some ad hoc procedure that you need to invent for every problem or validate every time you use it. </li>
<li> Don't need to figure out an estimator. </li>
<li> Once you have a model and data set, the conclusion is a computing problem, not a research problem. </li>
<li> Don't need to prove a theorem to show that your posterior is sensible. It is sensible if your assumptions are sensible. </li>
<li> Don't need to publish a bunch of papers to figure out sensible answers given a novel problem </li>
<li> For example, estimating a series of means $mu_1, mu_2, \ldots$ that you know are ordered $mu_j \le mu_{j+1}$ is a computing problem in Bayesian inference, but was the source of numerous papers in the frequentist literature. Finding a (good) frequentist estimator and finding standard errors and confidence intervals took lots of papers to figure out. </li>
</ol><li> Yes, you can still use SAS. </li>
<ol><li> Or R or Stata. </li>
</ol><li> Can incorporate utility functions, if you have one. </li>
<li> Odd bits of other information can be incorporated into the analysis, for example </li>
<ol><li> That a particular parameter, usually allowed to be positive or negative, must be positive. </li>
<li> That a particular parameter is probably positive, but not guaranteed to be positive. </li>
<li> That a given regression coefficient should be close to zero. </li>
<li> That group one's mean is larger than group two's mean. </li>
<li> That the data comes from a distribution that is not a Poisson, Binomial, Exponential or Normal. For example, the data may be better modeled by a t, gamma. </li>
<li> That a collection of parameters come from a distribution that is skewed, or has long tails. </li>
<li> Bayesian nonparametrics can allow you to model an unknown density as a non-parametric mixture of normals (or other density). The uncertainty in estimating this distribution is incorporated in making inferences about group means and regression coefficients. </li>
</ol><li> Bayesian modeling is about the science. </li>
<ol><li> You can calculate the probability that your hypothesis is true. </li>
<li> Bayesian modeling asks if this model describes the data, mother nature, the data generating process correctly, or sufficiently correctly. </li>
<li> Classical inference is all about the statistician and the algorithm, not the science. </li>
<li> In repeated samples, how often (or how accurately) does this algorithm/method/model/inference scheme give the right answer? </li>
<li> Classical inference is more about the robustness (in repeated sampling) of the procedure. In that way, it provides robustness results for Bayesian methods. </li>
</ol><li> Bayesian methods have had notable successes, to wit: </li>
<ol><li> Covariate selection in regression problems </li>
<li> Model selection </li>
<li> Model mixing </li>
<li> And mixture models </li>
<li> Missing data </li>
<li> Multi-level and hierarchical models </li>
<li> Phylogeny </li>
</ol></ol><p>The bottom line: <b>More tools</b>. <b>Faster progress</b>.</p>
</div>
<span><span lang="" about="https://robweiss.faculty.biostat.ucla.edu/user/10" typeof="schema:Person" property="schema:name" datatype="" xml:lang="">robweiss</span></span>
<span>Thu, 07/07/2016 - 12:18</span>
<div>
<div>Filed Under</div>
<div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/bayes" hreflang="en">Bayes</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/inference" hreflang="en">inference</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/comparative" hreflang="en">comparative</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/frequentist" hreflang="en">frequentist</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/classical" hreflang="en">classical</a></div>
<div><a href="https://robweiss.faculty.biostat.ucla.edu/why-be-bayesian" hreflang="en">Why be Bayesian?</a></div>
</div>
</div>
<ul class="links inline"><li class="comment-forbidden"><a href="https://robweiss.faculty.biostat.ucla.edu/user/login?destination=/blog/why-be-bayesian-let-me-count-ways%23comment-form">Log in</a> to post comments</li></ul><section><h2>Comments</h2>
<article data-comment-user-id="55" id="comment-1" class="comment js-comment"><mark class="hidden" data-comment-timestamp="1491322369"></mark><footer class="comment-wrap"><div class="author-details">
<article typeof="schema:Person" about="https://robweiss.faculty.biostat.ucla.edu/user/55" class="profile"></article><a href="https://robweiss.faculty.biostat.ucla.edu/comment/1#comment-1" hreflang="und">Permalink</a>
</div>
<div class="author-comments">
<p class="comment-submitted">Submitted by <span lang="" about="https://robweiss.faculty.biostat.ucla.edu/user/55" typeof="schema:Person" property="schema:name" datatype="" xml:lang="">biostat202Cuser</span> on Tue, 04/04/2017 - 09:12</p>
<div class="content">
<h3><a href="https://robweiss.faculty.biostat.ucla.edu/comment/1#comment-1" class="permalink" rel="bookmark" hreflang="und">Biostat 285</a></h3>
<div><p>Student name: Changhee Lee<br />
Department: Electrical engineering</p>
</div>
<drupal-render-placeholder callback="comment.lazy_builders:renderLinks" arguments="0=1&1=default&2=und&3=" token="R_Gu3EzvWLFCHetA-cUhY3EDBnkBpLlLEkaVICt8Hrw"></drupal-render-placeholder></div>
</div>
</footer></article></section>Thu, 07 Jul 2016 19:18:45 +0000robweiss183 at https://robweiss.faculty.biostat.ucla.edu