{"id":615,"date":"2018-05-21T08:30:00","date_gmt":"2018-05-21T12:30:00","guid":{"rendered":"http:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/?p=615"},"modified":"2018-05-20T12:53:13","modified_gmt":"2018-05-20T16:53:13","slug":"causal-inference-in-ecology-setting-the-stage-for-the-rubin-causal-model","status":"publish","type":"post","link":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/blog\/2018\/05\/21\/causal-inference-in-ecology-setting-the-stage-for-the-rubin-causal-model\/","title":{"rendered":"Causal inference in ecology &#8211; Setting the stage for the Rubin causal model"},"content":{"rendered":"<p><a href=\"http:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/causal-inference-in-ecology\/\">Causal inference in ecology &#8211; links to the series<\/a><\/p>\n<p>If you\u2019ve been following along, you realize by now that it\u2019s not easy to infer the cause of a phenomenon, even in a well-controlled experiment. What about observational experiments, which are what many ecologists and evolutionary biologists have? Take one paper of mine that I\u2019m reasonably happy with (<em>PLoS One<\/em> e52035; 2012. doi: <a href=\"https:\/\/dx.doi.org\/10.1371\/journal.pone.0052035\">10.1371\/journal.pone.0052035<\/a>). One part of that paper included an experiment of sorts. We<sup><a id=\"ffn1\" class=\"footnote\" href=\"#fn1\">1<\/a><\/sup> established experimental gardens at the <a href=\"https:\/\/www.sanbi.org\/gardens\/kirstenbosch\/\">Kirstenbosch National Botanical Garden<\/a> and at a mid-elevation site on <a href=\"https:\/\/mapcarta.com\/19078592\">Jonaskop mountain<\/a> about 100km due east of Kirstenbosch. Among other things we were interested both in how certain traits in newly formed leaves (specific leaf area, stomatal pore index, and leaf area) differed depending on the age of the plant and on the garden in which they were grown.<\/p>\n<p>This part of the paper is analogous to the thought experiment on corn that we\u2019ve discussed so far. Since the plants were grown from wild-collected seedlings, we obviously couldn\u2019t duplicate genotypes across gardens. We also didn\u2019t have enough seed from individual maternal plants to replicate families across the gardens. So we did the best we could. We randomized seedlings within populations and split populations across gardens. You\u2019ll see the results from this part of the analysis in the figure below.<\/p>\n<div id=\"attachment_617\" style=\"width: 700px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-617\" class=\"size-large wp-image-617\" src=\"http:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-9-1024x1013.png\" alt=\"\" width=\"690\" height=\"683\" srcset=\"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-9-1024x1013.png 1024w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-9-300x297.png 300w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-9-768x760.png 768w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-9-980x970.png 980w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-9.png 1357w\" sizes=\"auto, (max-width: 690px) 100vw, 690px\" \/><p id=\"caption-attachment-617\" class=\"wp-caption-text\">Figure 4 from <em>PLoS One<\/em> e52035; 2012. doi: <a href=\"https:\/\/dx.doi.org\/10.1371\/journal.pone.0052035\"> 10.1371\/journal.pone.0052035<\/a><\/p><\/div>\n<p>The trends with plant age are clear.<sup><a id=\"ffn2\" class=\"footnote\" href=\"#fn2\">2<\/a><\/sup> Specific leaf area (SLA) declines with plant age, stomatal pore index increases with plant age, and leaf area increases with plant age. Given that the trends are consistent across species and gardens, I\u2019m reasonably confident that plant age influences these traits <em><strong>in this group of Protea<\/strong><\/em>.<sup><a id=\"ffn3\" class=\"footnote\" href=\"#fn3\">3<\/a><\/sup> Notice that I wrote \u201cinfluences\u201d, which is a short way (for me) of writing that plant age is <em><strong>a<\/strong><\/em> causal factor that influences the traits but that I am not claiming that it is <em><strong>the<\/strong><\/em> causal factor.<sup><a id=\"ffn4\" class=\"footnote\" href=\"#fn4\">4<\/a><\/sup><\/p>\n<div id=\"attachment_619\" style=\"width: 700px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-619\" class=\"size-large wp-image-619\" src=\"http:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-6-905x1024.png\" alt=\"\" width=\"690\" height=\"781\" srcset=\"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-6-905x1024.png 905w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-6-265x300.png 265w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-6-768x869.png 768w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-6-980x1109.png 980w, https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-content\/uploads\/2018\/05\/pone-2012-7-e52035-6.png 1546w\" sizes=\"auto, (max-width: 690px) 100vw, 690px\" \/><p id=\"caption-attachment-619\" class=\"wp-caption-text\">Figure 3 from <em>PLoS One<\/em> e52035; 2012. doi: <a href=\"https:\/\/dx.doi.org\/10.1371\/journal.pone.0052035\"> 10.1371\/journal.pone.0052035<\/a><\/p><\/div>\n<p>Similarly, the figure above makes it clear that which garden the plants are grown in influences these (and other) traits. These results won\u2019t surprise anyone whose worked with plants. The traits plants have depend both on how old the plant is and on where its grown. So far a reasonably straightforward experiment, but how about this? We also wanted to know whether the amount of change in leaf traits depended on measures of resource availability and rainfall seasonality in the places that seeds were collected from. Here we\u2019re asking a more complicated question.<\/p>\n<p>Take SLA, for example, and imagine that we\u2019re asking the question just about changes in SLA for plants grown at Jonaskop. Now the fully fleshed out the question is something like this:<\/p>\n<ul>\n<li>I know that plants are often adapted to the local circumstances in which they are growing.<\/li>\n<li>If SLA reflects plant characteristics that are important in local adaptation to nutrients or water availability, then plants that grow in places that differ in nutrients or water availability should also differ in SLA in ways that make them well-suited to the place where they occur.<\/li>\n<li>Do we have evidence that there is an association between changes in SLA and nutrient availability or precipitation patterns in the site from which they are derived?<\/li>\n<\/ul>\n<p>It\u2019s that middle step that\u2019s tricky. We don\u2019t need to do anything special to run a regression on changes in SLA and home site characteristics, but to <em><strong>interpret<\/strong><\/em> that regression as evidence for the causal story in that middle step we need to do something more. Unlike the nicely randomized experiment with which we began this post, we aren\u2019t randomizing plants across sites and allowing them to adapt. What we have is purely observational data to address this question. To what extent can we make a causal inference from these data? That\u2019s the question I\u2019ll turn to in the next installment.<\/p>\n<ol id=\"footnotes\">\n<li id=\"fn1\">By \u201cwe\u201d I should be clear that Jane Carlson, Rachel Prunier, and Ann Marie Gawel collected all of the seeds that \u201cwe\u201d used to establish the gardens, and they did all of the work of germinating seedlings and establishing the gardens. They also collected nearly all of the data. I helped collect a little, but my help mostly consisted of standing there with a clipboard and data sheet and writing down numbers in the appropriate columns. <a href=\"#ffn1\">&#x21a9;<\/a><\/li>\n<li id=\"fn2\">The same plants were measured in 2009 and 2010. In both cases, measurements were made on newly formed, but fully expanded leaves. I\u2019m not reporting P-values, but you can find them in Table 1. <a href=\"#ffn2\">&#x21a9;<\/a><\/li>\n<li id=\"fn3\">The species are all members of a small, recently evolved monophyletic clade, <em>Protea<\/em> sect. <em>Exsertae<\/em>. <a href=\"#ffn3\">&#x21a9;<\/a><\/li>\n<li id=\"fn4\">Notice also that I am discounting the possibility that it is weather in the year the plants were growing that influences their traits rather than their age. <a href=\"#ffn4\">&#x21a9;<\/a><\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Causal inference in ecology &#8211; links to the series If you\u2019ve been following along, you realize by now that it\u2019s not easy to infer the cause of a phenomenon, even&#8230; <a class=\"read-more-button\" href=\"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/blog\/2018\/05\/21\/causal-inference-in-ecology-setting-the-stage-for-the-rubin-causal-model\/\">Read more &gt;<\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[10],"tags":[],"class_list":["post-615","post","type-post","status-publish","format-standard","hentry","category-statistics"],"_links":{"self":[{"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/posts\/615","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/comments?post=615"}],"version-history":[{"count":4,"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/posts\/615\/revisions"}],"predecessor-version":[{"id":621,"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/posts\/615\/revisions\/621"}],"wp:attachment":[{"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/media?parent=615"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/categories?post=615"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/darwin.eeb.uconn.edu\/uncommon-ground\/wp-json\/wp\/v2\/tags?post=615"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}