Hi Louise-<div><br></div><div>Great question. </div><div><br></div><div>Yes, doing sensitivity analysis with GAM's as Kosuke mentioned is beyond what we can do right now. Its possible in theory, we think, but would be a major undertaking to figure out how to do. Our Stata package does not do this either (currently it supports a smaller set of models compared to our R package).</div>
<div><br></div><div>Let me make another suggestion. When you estimate things with the GAM's, are you uncovering highly non-linear relationships (graph things out like what we do in the Psychological Methods paper). If you aren't, perhaps Mr. Linear is your friend, and you can report the sensitivity analysis based on the linear model as the robustness check?</div>
<div><br></div><div>best,</div><div>Dustin</div><div><br></div><div>PS: as an aside, we're so psyched someone from the WWF is using our software! Really cool!</div><div><br clear="all">Dustin Tingley<br>Government Department<br>
Harvard University<br><a href="http://scholar.harvard.edu/dtingley" target="_blank">http://scholar.harvard.edu/dtingley</a><br>
<br><br><div class="gmail_quote">On Sun, Dec 11, 2011 at 10:35 PM, Kosuke Imai <span dir="ltr"><<a href="mailto:kimai@princeton.edu">kimai@princeton.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Louise,<br>
<br>
Thanks for your email. Unfortunately, our software, mediation, does not have sensitivity analysis for GAM. We hope to expand the capability of the software in the future, though. However, you should still be able to use GAM to calculate mediation effects.<br>
<br>
Best,<br>
Kosuke<br>
<br>
Department of Politics<br>
Princeton University<br>
<a href="http://imai.princeton.edu" target="_blank">http://imai.princeton.edu</a><br>
<br>
<br>
On Dec 10, 2011, at 11:54 AM, Glew, Louise wrote:<br>
<br>
> Dr. Imai,<br>
><br>
> I’m working on a series of quasi-experimental assessments of the impact of biodiversity conservation interventions. I was wondering if I could ask your advice on how best to conduct a sensitivity analysis on a generalized additive mixed model.<br>
><br>
> I’m analyzing a time-series of remotely-sensed imagery, to explore the impact of a series of conservation interventions against matched controls in northern Kenya. The properties of the data have precluded using a GLMM, and are particularly challenging as they are both spatially auto-correlated and serially-dependent.<br>
><br>
> When working with less complex data, I have used the ‘rbounds’ and the ‘mediation’ package to compute Rosenbaum’s sensitivity bounds. Reading the documentation for the mediation package, it became apparent that the current function can only handle linear models. The GAMM I’ve fit to the data is technically a penalized linear model, but the errors are not normally distributed (due to the spatial and temporal correlation). I was wondering if you could advise me of any alternative sensitivity analysis methods which may be able to cope with non-linear data please?<br>
><br>
> Thank you,<br>
><br>
> Best wishes,<br>
><br>
> Louise<br>
><br>
> Louise Glew (MEnvSci)<br>
> Monitoring and Evaluation Specialist<br>
> Conservation Science Program<br>
> World Wildlife Fund<br>
> 1250 24th Street, NW<br>
> Washington, DC 20037-1193<br>
> Office: <a href="tel:202-495-4184" value="+12024954184">202-495-4184</a><br>
> Skype: louise.glew<br>
> <a href="mailto:Louise.Glew@wwfus.org">Louise.Glew@wwfus.org</a><br>
> <a href="http://worldwildlife.org" target="_blank">worldwildlife.org</a><br>
><br>
> <image001.gif><br>
><br>
<br>
<br>
_______________________________________________<br>
Mediation-information mailing list<br>
<a href="mailto:Mediation-information@lists.r-forge.r-project.org">Mediation-information@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/mediation-information" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/mediation-information</a><br>
</blockquote></div><br></div>