<?xml version='1.0' encoding='UTF-8'?><?xml-stylesheet href="http://www.blogger.com/styles/atom.css" type="text/css"?><feed xmlns='http://www.w3.org/2005/Atom' xmlns:openSearch='http://a9.com/-/spec/opensearchrss/1.0/' xmlns:blogger='http://schemas.google.com/blogger/2008' xmlns:georss='http://www.georss.org/georss' xmlns:gd="http://schemas.google.com/g/2005" xmlns:thr='http://purl.org/syndication/thread/1.0'><id>tag:blogger.com,1999:blog-7931813</id><updated>2024-09-07T22:20:47.326-04:00</updated><category term="actuaries"/><category term="blog"/><category term="rstats"/><category term="twitter"/><category term="insurance"/><category term="operational risk"/><category term="reinsurance"/><category term="social media"/><category term="Bayesian"/><category term="CAS"/><category term="CIA"/><category term="FEM"/><category term="Gibbs"/><category term="HTML"/><category term="Hastings"/><category term="IE6"/><category term="JAGS"/><category term="LaTeX"/><category term="MCMC"/><category term="MathJax"/><category term="Metropolis"/><category term="SOA"/><category term="accounting"/><category term="actuary"/><category term="bit.ly"/><category term="blogger"/><category term="cartoon"/><category term="currency"/><category term="curve fitting"/><category term="economy"/><category term="facebook"/><category term="firefox"/><category term="identi.ca"/><category term="internet"/><category term="linkedin"/><category term="loss distributions"/><category term="online behavior"/><category term="pingfm"/><category term="solvency"/><category term="twitterfeed"/><category term="web design"/><title type='text'>Random Fluctuations</title><subtitle type='html'>Frequency and Severity&lt;br&gt;&#xa;Random Thoughts for a Random World</subtitle><link rel='http://schemas.google.com/g/2005#feed' type='application/atom+xml' href='http://severity.blogspot.com/feeds/posts/default'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default?alt=atom&amp;redirect=false'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/'/><link rel='hub' href='http://pubsubhubbub.appspot.com/'/><link rel='next' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default?alt=atom&amp;start-index=26&amp;max-results=25&amp;redirect=false'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><generator version='7.00' uri='http://www.blogger.com'>Blogger</generator><openSearch:totalResults>28</openSearch:totalResults><openSearch:startIndex>1</openSearch:startIndex><openSearch:itemsPerPage>25</openSearch:itemsPerPage><entry><id>tag:blogger.com,1999:blog-7931813.post-2807682338129789808</id><published>2012-01-19T20:44:00.000-05:00</published><updated>2012-01-19T21:33:43.218-05:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="HTML"/><category scheme="http://www.blogger.com/atom/ns#" term="LaTeX"/><category scheme="http://www.blogger.com/atom/ns#" term="MathJax"/><title type='text'>MathJax</title><content type='html'>In the latest issue of the &lt;i&gt;Notices&lt;/i&gt; of the American Mathematical Society, there is an article on &lt;a href=&quot;http://www.mathjax.org/&quot;&gt;MathJax&lt;/a&gt;, a project to help nicely render mathematics on the web, using TeX or MathML. I decided to see if I can get it working here on blogspot, so here goes nuthin&#39;:

$$\LARGE TVaR_\alpha =\frac{\int_{F^{-1}(\alpha)}^\infty x\;f(x)\;dx}{1-\alpha}$$

NICE!

I&#39;m using the MathJax &lt;a href=&quot;http://www.mathjax.org/docs/1.1/start.html#mathjax-cdn&quot;&gt;content delivery network&lt;/a&gt;; well, actually, the &lt;a href=&quot;http://www.mathjax.org/docs/1.1/start.html#secure-access-to-the-cdn&quot;&gt;secure connection&lt;/a&gt; (of course). Once that is done, then &lt;code&gt;$$...$$&lt;/code&gt; or  &lt;code&gt;\[...\]&lt;/code&gt; delimit displayed mathematics, and &lt;code&gt;\(...\)&lt;/code&gt; delimits in-line math. What should also work is right-clicking the image to see the source, although that may only function if MathJax is installed locally.

Regardless, this should make it easier to display mathematics on the web, and should obviate the need to have images rendered as gifs and loaded (as can be seen in previous posts. Maybe I&#39;ll post a bit more often now!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/2807682338129789808/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2012/01/mathjax.html#comment-form' title='1 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/2807682338129789808'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/2807682338129789808'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2012/01/mathjax.html' title='MathJax'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>1</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-4711985487589400360</id><published>2011-08-22T18:46:00.001-04:00</published><updated>2011-08-22T18:46:46.066-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="Bayesian"/><category scheme="http://www.blogger.com/atom/ns#" term="Gibbs"/><category scheme="http://www.blogger.com/atom/ns#" term="Hastings"/><category scheme="http://www.blogger.com/atom/ns#" term="JAGS"/><category scheme="http://www.blogger.com/atom/ns#" term="MCMC"/><category scheme="http://www.blogger.com/atom/ns#" term="Metropolis"/><category scheme="http://www.blogger.com/atom/ns#" term="rstats"/><title type='text'>Bayesian analysis: Comparing algorithms Part 1?</title><content type='html'>I recently had the opportunity to engage in some Bayesian analysis at work. I was able to state the problem in terms of the lognormal distribution, and took advantage of &lt;a href=&quot;http://calvin.iarc.fr/%7Emartyn/software/jags/&quot;&gt;JAGS&lt;/a&gt; and its integration with &quot;R&quot; using the &lt;a href=&quot;http://cran.r-project.org/web/packages/R2jags/index.html&quot;&gt;R2jags&lt;/a&gt; package. The client was very happy, and life moved on. However, as the underlying problem may be better described by a transformed beta than a lognormal undergoing transforms (say that five times quickly!), it got me to thinking about how to teach myself at least rudimentary Bayesian analysis using MCMC methods.As my main text, I have been using Dr. Scott M. Lynch&#39;s excellent text,&lt;span style=&quot;font-size: small;&quot;&gt;&lt;span id=&quot;btAsinTitle&quot;&gt; &quot;&lt;a href=&quot;http://www.amazon.com/Introduction-Statistics-Estimation-Scientists-Behavioral/dp/038771264X&quot;&gt;&lt;i&gt;Introduction to Applied Bayesian Statistics and Estimation for Social Scientist&lt;/i&gt;&lt;/a&gt;&quot;. Starting small (I haven&#39;t worked on the four-parameter distribution yet), I decided to play with data for which I knew the answer, and see if I could write my own Gibbs or Metropolis-Hastings algorithms. I picked a gamma distribution as my test case.&lt;/span&gt;&lt;/span&gt;&lt;br /&gt;
&lt;span style=&quot;font-size: small;&quot;&gt;&lt;span id=&quot;btAsinTitle&quot;&gt;&lt;br /&gt;
&lt;/span&gt;&lt;/span&gt;&lt;br /&gt;
&lt;span style=&quot;font-size: small;&quot;&gt;&lt;span id=&quot;btAsinTitle&quot;&gt;I&#39;m no expert on MCMC methods, so I&#39;m not going to even attempt to explain how this works, but I did see how the different algorithms behave, and it is rather insightful. For the gamma distribution, while I always parameterize it using shape and scale, for the purposes of Gibbs sampling, I used shape/rate. The nice thging about Gibbs sampling is that as long as your can sample from a distribution which is &lt;i&gt;proportional&lt;/i&gt; to the actual &lt;/span&gt;&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;&lt;span id=&quot;btAsinTitle&quot;&gt;conditional distribution, you&#39;re fine, and do not have to calculate the constants of proportionality that bring the distribution back to an integral of 1.&lt;/span&gt;&lt;/span&gt;&lt;br /&gt;
&lt;span style=&quot;font-size: small;&quot;&gt;&lt;span id=&quot;btAsinTitle&quot;&gt;&lt;br /&gt;
&lt;/span&gt;&lt;/span&gt;&lt;br /&gt;
&lt;span style=&quot;font-size: small;&quot;&gt;&lt;span id=&quot;btAsinTitle&quot;&gt;In the case of the gamma distribution, and using a flat prior, it is the likelihood which dominates the conditional distribution. Using the shape/rate parameterization, &lt;a href=&quot;http://www.codecogs.com/eqnedit.php?latex=%5Cfrac%7B%5Cbeta%5E%5Calpha%20x%5E%7B%5Calpha-1%7De%5E%7B-%5Cbeta%20x%7D%7D%7B%7B%5CGamma%28%5Calpha%29%7D%7D&quot; target=&quot;_blank&quot;&gt;&lt;img src=&quot;http://latex.codecogs.com/gif.latex?%5Cfrac%7B%5Cbeta%5E%5Calpha%20x%5E%7B%5Calpha-1%7De%5E%7B-%5Cbeta%20x%7D%7D%7B%7B%5CGamma%28%5Calpha%29%7D%7D&quot; title=&quot;\frac{\beta^\alpha x^{\alpha-1}e^{-\beta x}}{{\Gamma(\alpha)}}&quot; /&gt;&lt;/a&gt;, the likelihood over the observations &lt;a href=&quot;http://www.codecogs.com/eqnedit.php?latex=x_i&quot; target=&quot;_blank&quot;&gt;&lt;img src=&quot;http://latex.codecogs.com/svg.latex?x_i&quot; title=&quot;x_i&quot; /&gt;&lt;/a&gt; is:&lt;/span&gt;&lt;/span&gt;&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;text-align: center;&quot;&gt;&lt;a href=&quot;http://www.codecogs.com/eqnedit.php?latex=%5Cpropto%20%5Cfrac%7B%5Cbeta%5E%7Bn%5Calpha%7D%20%28%5Cprod_i%20x_i%5E%7B%5Calpha-1%7D%29%20e%5E%7B-%5Cbeta%5Csum_i%7B%20x%7D%7D%7D%7B%7B%5CGamma%28%5Calpha%29%5En%7D%7D&quot; target=&quot;_blank&quot;&gt;&lt;img src=&quot;http://latex.codecogs.com/gif.latex?%5Cpropto%20%5Cfrac%7B%5Cbeta%5E%7Bn%5Calpha%7D%20%28%5Cprod_i%20x_i%5E%7B%5Calpha-1%7D%29%20e%5E%7B-%5Cbeta%5Csum_i%7B%20x%7D%7D%7D%7B%7B%5CGamma%28%5Calpha%29%5En%7D%7D&quot; title=&quot;\propto \frac{\beta^{n\alpha} (\prod_i x_i^{\alpha-1}) e^{-\beta\sum_i{ x}}}{{\Gamma(\alpha)^n}}&quot; /&gt;&lt;/a&gt;&lt;/div&gt;&lt;br /&gt;
When looking for a conditional distribution for alpha, this does not resolve into anything remotely nice (although see Fink (1997) &quot;&lt;a href=&quot;https://encrypted.google.com/url?sa=t&amp;amp;source=web&amp;amp;cd=1&amp;amp;ved=0CBsQFjAA&amp;amp;url=http%3A%2F%2Fciteseerx.ist.psu.edu%2Fviewdoc%2Fdownload%3Fdoi%3D10.1.1.157.5540%26rep%3Drep1%26type%3Dpdf&amp;amp;rct=j&amp;amp;q=A%20Compendium%20of%20Conjugate%20Priors&amp;amp;ei=ksFSTrf2OKfw0gHjzOmNBw&amp;amp;usg=AFQjCNEyQcx9zQ_zWmPED0fOwg0kErsMsw&amp;amp;sig2=jozf4mOojgeqBxFMAzAvLA&amp;amp;cad=rja&quot;&gt;A Compendium of Conjugate Priors&lt;/a&gt;&quot; (pdf) for more) so Gibbs sampling will have to give way to a Metropolis algorithm. I&#39;m not sophisticated enough to use a non-symmetrical proposal distribution, so it is not a Metropolis-Hastings algorithm, maybe in a future post. However, when looking at beta, and discarding the x&#39;s and alphas, the core of the likelihood is rather gamma-like. Actually, it is proportional to a gamma distribution with the rate equal to the sum of the observed data points and the shape equal to one more than the product of the number of observations and the underlying alpha.&lt;br /&gt;
&lt;br /&gt;
With that observation, and properly writing the algorithm (which took more time than I care to admit), a comparison between using Gibbs sampling or straight Metropolis on the scale parameter can be made.&lt;br /&gt;
&lt;br /&gt;
Here is a &quot;pure&quot; Metropolis algorithm to calculate the posterior distribution, assuming the input data is gamma (this and all &quot;pretty&quot; R code &lt;a href=&quot;http://www.inside-r.org/pretty-r&quot; title=&quot;created by Pretty R at inside-R.org&quot;&gt;created by Pretty R at inside-R.org&lt;/a&gt; :&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;overflow: auto;&quot;&gt;&lt;div class=&quot;geshifilter&quot;&gt;&lt;pre class=&quot;r geshifilter-R&quot; style=&quot;font-family: monospace;&quot;&gt;MHGam &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/function&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;function&lt;/span&gt;&lt;/a&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;x&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter=&lt;span style=&quot;color: #cc66cc;&quot;&gt;10000&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.burnin=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
alpha &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/matrix&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;matrix&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;.5&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
theta &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/matrix&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;matrix&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;10&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
acca&amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;
accq &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;
lnpost &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/function&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;function&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;a&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/q&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;q&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; z&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
  NLL &amp;lt;- &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;a-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;*&lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;z&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;-&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;z/q&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;-a*&lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/q&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;q&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;-&lt;a href=&quot;http://inside-r.org/r-doc/base/lgamma&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;lgamma&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;a&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
  &lt;a href=&quot;http://inside-r.org/r-doc/base/return&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;return&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/sum&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;sum&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;NLL&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
&lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;for&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;i &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;in&lt;/span&gt; &lt;span style=&quot;color: #cc66cc;&quot;&gt;2&lt;/span&gt;:n.iter&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
&lt;span style=&quot;color: #666666; font-style: italic;&quot;&gt;#sample theta&lt;/span&gt;
  acc &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;
  theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; + &lt;a href=&quot;http://inside-r.org/r-doc/stats/rnorm&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;rnorm&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/mean&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;mean&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/stats/sd&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;sd&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;100&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
  &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&amp;lt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;acc=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;;&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
  LL &amp;lt;- lnpost&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; x&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
  &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;LL - lnpost&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; x&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &amp;lt; &lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/stats/runif&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;runif&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
    theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;
    acc &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;
  &lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
  accq &amp;lt;- accq + acc
&lt;span style=&quot;color: #666666; font-style: italic;&quot;&gt;#sample alpha&lt;/span&gt;
    acc &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;
    alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; + &lt;a href=&quot;http://inside-r.org/r-doc/stats/runif&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;runif&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/min&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;min&lt;/span&gt;&lt;/a&gt;=-&lt;span style=&quot;color: #cc66cc;&quot;&gt;.5&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/max&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;max&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;.5&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
    &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&amp;lt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;acc=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;;&lt;/span&gt; alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
    &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;lnpost&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; x&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; - LL&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &amp;lt; &lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/stats/runif&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;runif&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
      alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;
      acc &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
    acca &amp;lt;- acca+acc
    &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;i%%100==&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/print&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;print&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/c&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;c&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;i&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; acca/i&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; accq/i&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
    &lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
amean &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/mean&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;mean&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;n.burnin+&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;:n.iter&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
qmean&amp;lt;-&lt;a href=&quot;http://inside-r.org/r-doc/base/mean&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;mean&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;n.burnin+&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;:n.iter&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
Fit &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/structure&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;structure&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/list&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;list&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha=alpha&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta=theta&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; amean=amean&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; qmean=qmean&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;a href=&quot;http://inside-r.org/r-doc/base/return&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;return&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;Fit&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;&lt;/pre&gt;&lt;/div&gt;&lt;/div&gt;&lt;br /&gt;
The above code will just hand off the parameter sampling between alpha and theta (which is the &lt;i&gt;scale&lt;/i&gt; and not rate) and will, hopefully, eventually settle down. I added two extra parameters for iterations and burnin, and eventually, I could add another, outer, for-next loop for multiple chains.&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
The next code takes advantage of Gibbs sampling for the scale parameter:&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;overflow: auto;&quot;&gt;&lt;div class=&quot;geshifilter&quot;&gt;&lt;pre class=&quot;r geshifilter-R&quot; style=&quot;font-family: monospace;&quot;&gt;MHGamGib &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/function&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;function&lt;/span&gt;&lt;/a&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;x&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter=&lt;span style=&quot;color: #cc66cc;&quot;&gt;10000&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.burnin=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
alpha &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/matrix&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;matrix&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;.5&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
theta &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/matrix&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;matrix&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;10&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
acca&amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;
lnpost &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/function&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;function&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;a&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/q&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;q&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; z&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
  NLL &amp;lt;- &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;a-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;*&lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;z&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;-&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;z/q&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;-a*&lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/q&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;q&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;-&lt;a href=&quot;http://inside-r.org/r-doc/base/lgamma&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;lgamma&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;a&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
  &lt;a href=&quot;http://inside-r.org/r-doc/base/return&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;return&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/sum&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;sum&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;NLL&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
&lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;for&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;i &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;in&lt;/span&gt; &lt;span style=&quot;color: #cc66cc;&quot;&gt;2&lt;/span&gt;:n.iter&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
&lt;span style=&quot;color: #666666; font-style: italic;&quot;&gt;#sample theta&lt;/span&gt;
  theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;/rgamma&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; rate=&lt;a href=&quot;http://inside-r.org/r-doc/base/sum&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;sum&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;Test&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/packages/cran/shape&quot;&gt;shape&lt;/a&gt;=alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;*&lt;a href=&quot;http://inside-r.org/r-doc/base/length&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;length&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;x&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;+&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #666666; font-style: italic;&quot;&gt;#sample alpha&lt;/span&gt;
    acc &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;
    alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; + &lt;a href=&quot;http://inside-r.org/r-doc/stats/runif&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;runif&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/min&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;min&lt;/span&gt;&lt;/a&gt;=-&lt;span style=&quot;color: #cc66cc;&quot;&gt;.5&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/max&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;max&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;.5&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
    &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&amp;lt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;acc=&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;;&lt;/span&gt; alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
    &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;lnpost&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; x&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; - lnpost&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; x&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &amp;lt; &lt;a href=&quot;http://inside-r.org/r-doc/base/log&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;log&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/stats/runif&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;runif&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;
      alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt; &amp;lt;- alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i-&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;
      acc &amp;lt;- &lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
    acca &amp;lt;- acca+acc
    &lt;span style=&quot;color: black; font-weight: bold;&quot;&gt;if&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;i%%100==&lt;span style=&quot;color: #cc66cc;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt; &lt;span style=&quot;color: #009900;&quot;&gt;{&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/print&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;print&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/c&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;c&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;i&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;i&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; acca/i&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
    &lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;
amean &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/mean&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;mean&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;n.burnin+&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;:n.iter&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
qmean&amp;lt;-&lt;a href=&quot;http://inside-r.org/r-doc/base/mean&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;mean&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;theta&lt;span style=&quot;color: #009900;&quot;&gt;[&lt;/span&gt;n.burnin+&lt;span style=&quot;color: #cc66cc;&quot;&gt;1&lt;/span&gt;:n.iter&lt;span style=&quot;color: #009900;&quot;&gt;]&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
Fit &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/base/structure&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;structure&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/list&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;list&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;alpha=alpha&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; theta=theta&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; amean=amean&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; qmean=qmean&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;a href=&quot;http://inside-r.org/r-doc/base/return&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;return&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;Fit&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;span style=&quot;color: #009900;&quot;&gt;}&lt;/span&gt;&lt;/pre&gt;&lt;/div&gt;&lt;/div&gt;&lt;br /&gt;
Note that I&#39;m still parameterizing using the scale by sampling the rate and using its reciprocal.&lt;br /&gt;
&lt;br /&gt;
To set up, run, and compare the models, I used the following code:&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;overflow: auto;&quot;&gt;&lt;div class=&quot;geshifilter&quot;&gt;&lt;pre class=&quot;r geshifilter-R&quot; style=&quot;font-family: monospace;&quot;&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/set.seed&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;set.seed&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;22&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&amp;nbsp;
Test &amp;lt;- &lt;a href=&quot;http://inside-r.org/r-doc/stats/rgamma&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;rgamma&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;color: #cc66cc;&quot;&gt;100&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/packages/cran/shape&quot;&gt;shape&lt;/a&gt;=&lt;span style=&quot;color: #cc66cc;&quot;&gt;2.5&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/scale&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;scale&lt;/span&gt;&lt;/a&gt;=1e5&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&amp;nbsp;
MHT&amp;lt;-&lt;a href=&quot;http://inside-r.org/r-doc/base/system.time&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;system.time&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;MH&amp;lt;-MHGam&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;Test&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter=&lt;span style=&quot;color: #cc66cc;&quot;&gt;50000&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&amp;nbsp;
MHGT&amp;lt;-&lt;a href=&quot;http://inside-r.org/r-doc/base/system.time&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;system.time&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;MHG&amp;lt;-MHGamGib&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;Test&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; n.iter=&lt;span style=&quot;color: #cc66cc;&quot;&gt;50000&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;/pre&gt;&lt;/div&gt;&lt;/div&gt;&lt;br /&gt;
The Gibbs version runs about 40% faster, comparing MHT with MHGT (at least on my laptop), which stands to reason as Gibbs sampling is Metropolis sampling where every variate is accepted. What is more fascinating is how fast the parameter values converge:&lt;br /&gt;
&lt;br /&gt;
Using the following code to overlay the trace plots for theta:&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;overflow: auto;&quot;&gt;&lt;div class=&quot;geshifilter&quot;&gt;&lt;pre class=&quot;r geshifilter-R&quot; style=&quot;font-family: monospace;&quot;&gt;&lt;a href=&quot;http://inside-r.org/r-doc/graphics/plot&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;plot&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;MHG$theta&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; type=&lt;span style=&quot;color: blue;&quot;&gt;&quot;l&quot;&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/col&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;col&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: blue;&quot;&gt;&quot;red&quot;&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;a href=&quot;http://inside-r.org/r-doc/graphics/points&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;points&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;MH$theta&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; type=&lt;span style=&quot;color: blue;&quot;&gt;&quot;l&quot;&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/col&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;col&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: blue;&quot;&gt;&quot;blue&quot;&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;/pre&gt;&lt;/div&gt;&lt;/div&gt;&lt;br /&gt;
&lt;br /&gt;
We see:&lt;br /&gt;
&lt;div class=&quot;separator&quot; style=&quot;clear: both; text-align: center;&quot;&gt;&lt;a href=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi4apBDw5ZVX8JIYYu8PA-_51MMpucsRpc5WppkUg5WwesVPt5iyn_vtITpXJeW5XuF79eVYdfjo1LoaUCurwAiF-q1pN2KS-cEgUuhZI5Jbcym-bt6E7CwTuQjHtSr7LnLLBES/s1600/ThetaPlot.png&quot; imageanchor=&quot;1&quot; style=&quot;margin-left: 1em; margin-right: 1em;&quot;&gt;&lt;img border=&quot;0&quot; height=&quot;261&quot; src=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi4apBDw5ZVX8JIYYu8PA-_51MMpucsRpc5WppkUg5WwesVPt5iyn_vtITpXJeW5XuF79eVYdfjo1LoaUCurwAiF-q1pN2KS-cEgUuhZI5Jbcym-bt6E7CwTuQjHtSr7LnLLBES/s320/ThetaPlot.png&quot; width=&quot;320&quot; /&gt;&lt;/a&gt;&lt;/div&gt;The Gibbs sampling (in red) gets to the right neighborhood much more quickly for theta than the Metropolis (blue), and gives a more robust (less serially correlated) answer.&lt;br /&gt;
&lt;br /&gt;
What is also surprising is how the better estimates for theta affect the random-walk Metropolis algorithm for alpha. Using the following code:&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;overflow: auto;&quot;&gt;&lt;div class=&quot;geshifilter&quot;&gt;&lt;pre class=&quot;r geshifilter-R&quot; style=&quot;font-family: monospace;&quot;&gt;&lt;a href=&quot;http://inside-r.org/r-doc/graphics/plot&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;plot&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;MH$alpha&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; type=&lt;span style=&quot;color: blue;&quot;&gt;&quot;l&quot;&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/col&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;col&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: blue;&quot;&gt;&quot;blue&quot;&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;
&lt;a href=&quot;http://inside-r.org/r-doc/graphics/points&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;points&lt;/span&gt;&lt;/a&gt;&lt;span style=&quot;color: #009900;&quot;&gt;(&lt;/span&gt;MHG$alpha&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; type=&lt;span style=&quot;color: blue;&quot;&gt;&quot;l&quot;&lt;/span&gt;&lt;span style=&quot;color: #339933;&quot;&gt;,&lt;/span&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/col&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;col&lt;/span&gt;&lt;/a&gt;=&lt;span style=&quot;color: blue;&quot;&gt;&quot;red&quot;&lt;/span&gt;&lt;span style=&quot;color: #009900;&quot;&gt;)&lt;/span&gt;&lt;/pre&gt;&lt;/div&gt;&lt;/div&gt;&lt;br /&gt;
We get the following trace:&lt;br /&gt;
&lt;div class=&quot;separator&quot; style=&quot;clear: both; text-align: center;&quot;&gt;&lt;a href=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh0TGR6-QRKMy4XOzlpQAnV1iwLYGH60fs-SVGqsFVZD1CB0FxnR_Gj_onpZa3lBODZnCVuugYYxD-EovmvJrp8CZQYK7P0Ail9uPhXKQBNqIzc4ExmYEpYaI85Q59VnLNhJUcQ/s1600/AlphaPlot.png&quot; imageanchor=&quot;1&quot; style=&quot;margin-left: 1em; margin-right: 1em;&quot;&gt;&lt;img border=&quot;0&quot; height=&quot;261&quot; src=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh0TGR6-QRKMy4XOzlpQAnV1iwLYGH60fs-SVGqsFVZD1CB0FxnR_Gj_onpZa3lBODZnCVuugYYxD-EovmvJrp8CZQYK7P0Ail9uPhXKQBNqIzc4ExmYEpYaI85Q59VnLNhJUcQ/s320/AlphaPlot.png&quot; width=&quot;320&quot; /&gt;&lt;/a&gt;&lt;/div&gt;Personally, I found that, writing the algorithms and seeing the direct results helped me understand the concepts better than just reading, and I&#39;ll certainly try to use Gibbs sampling where I can!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/4711985487589400360/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/08/bayesian-analysis-comparing-algorithms.html#comment-form' title='5 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/4711985487589400360'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/4711985487589400360'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/08/bayesian-analysis-comparing-algorithms.html' title='Bayesian analysis: Comparing algorithms Part 1?'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><media:thumbnail xmlns:media="http://search.yahoo.com/mrss/" url="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi4apBDw5ZVX8JIYYu8PA-_51MMpucsRpc5WppkUg5WwesVPt5iyn_vtITpXJeW5XuF79eVYdfjo1LoaUCurwAiF-q1pN2KS-cEgUuhZI5Jbcym-bt6E7CwTuQjHtSr7LnLLBES/s72-c/ThetaPlot.png" height="72" width="72"/><thr:total>5</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-2190543106867924452</id><published>2011-06-07T13:24:00.000-04:00</published><updated>2011-06-07T13:24:36.383-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="blog"/><category scheme="http://www.blogger.com/atom/ns#" term="internet"/><category scheme="http://www.blogger.com/atom/ns#" term="online behavior"/><category scheme="http://www.blogger.com/atom/ns#" term="operational risk"/><category scheme="http://www.blogger.com/atom/ns#" term="social media"/><category scheme="http://www.blogger.com/atom/ns#" term="twitter"/><title type='text'>Social media redux</title><content type='html'>I&#39;m certain that by now, everyone knows about Congressman &lt;a href=&quot;http://www.anthonyweiner.com/&quot;&gt;Anthony Weiner&lt;/a&gt;, his &lt;a href=&quot;https://secure.wikimedia.org/wikipedia/en/wiki/Anthony_Weiner#Twitter_photo_scandal&quot;&gt;social media behavior&lt;/a&gt;, and &lt;a href=&quot;http://www.cbsnews.com/8301-503544_162-20069435-503544.html&quot;&gt;yesterday&#39;s press conference&lt;/a&gt;. I&#39;m not going to discuss Weiner&#39;s immoral and, dare I say it, despicable behavior in his online actions and the ensuing attempt at a coverup. What I &lt;i&gt;do&lt;/i&gt; want to mention is that this serves to reiterate one of the important points I mentioned in &lt;a href=&quot;http://severity.blogspot.com/2011/06/social-media-guidelines-specific-for.html&quot;&gt;yesterday&#39;s blog post&lt;/a&gt;: that anything you say or post is there forever. If you would not want something on the &lt;a href=&quot;http://www.nydailynews.com/news/galleries/june_2011_daily_news_front_pages/june_2011_daily_news_front_pages.html&quot;&gt;front page&lt;/a&gt; of your &lt;a href=&quot;http://www.nypost.com/rw/nypost/2011/06/07/covers/front060711.jpg&quot;&gt;local tabloid&lt;/a&gt;, don&#39;t post it &lt;b&gt;anywhere&lt;/b&gt; on the Internet. Of course, if you always attempt to behave appropriately, chances are you won&#39;t get into such a mess to begin with, but as we say in the business - that is out of scope for this discussion.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/2190543106867924452/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/06/social-media-redux.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/2190543106867924452'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/2190543106867924452'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/06/social-media-redux.html' title='Social media redux'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-4698339451386033903</id><published>2011-06-06T14:31:00.001-04:00</published><updated>2011-06-06T14:33:14.341-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="blog"/><category scheme="http://www.blogger.com/atom/ns#" term="facebook"/><category scheme="http://www.blogger.com/atom/ns#" term="insurance"/><category scheme="http://www.blogger.com/atom/ns#" term="linkedin"/><category scheme="http://www.blogger.com/atom/ns#" term="operational risk"/><category scheme="http://www.blogger.com/atom/ns#" term="reinsurance"/><category scheme="http://www.blogger.com/atom/ns#" term="social media"/><category scheme="http://www.blogger.com/atom/ns#" term="twitter"/><title type='text'>Social Media guidelines specific for (re)insurers</title><content type='html'>With social media here to stay, and many of us (actuaries, cat modelers, and other (re)insurance company personnel) actively involved in the social media on both personal and professional levels, it behooves us to remember certain social niceties. Everything that you post on the Internet is pretty much there for posterity; 50 years from now, the &lt;a href=&quot;http://wayback.archive.org/web/&quot;&gt;wayback machine &lt;/a&gt;will very likely be able to retrieve that one time you posted a picture of yourself at the holiday party—&lt;i&gt;after&lt;/i&gt; spiked eggnog #7. My personal, primary rules, whether it be on twitter, linkedin, facebook, wikipedia, or some random bulletin board are:&lt;br /&gt;
&lt;ol&gt;&lt;li&gt;Never post anything on-line that you would be uncomfortable saying to someone&#39;s face &lt;i&gt;in public&lt;/i&gt;.&lt;/li&gt;
&lt;li&gt;Never post anything on-line that you would be devastated seeing on the front page of the New York Times.&lt;/li&gt;
&lt;/ol&gt;I am all for having a certain amount of silliness and fun on-line, I mean, where are you going to let your hair down, at work? However, there is a &lt;i&gt;big&lt;/i&gt; difference between good-natured silliness and flaming. Too many people believe that there is anonymity on the Internet, and thus act selfishly, spitefully, and all-too-often just plain cruelly with the mistaken assumption that they will never be uncovered. Besides for this being awful and immature behavior, which speaks loads about the person demonstrating it, it can be a serious operational risk hazard for the person&#39;s employer, and thus potentially grounds for termination for the employee. For example, it has not even been a year yet since a Massachusetts schools&#39; supervisor &lt;a href=&quot;http://www.thebostonchannel.com/news/24670937/detail.html&quot;&gt;lost her job over a facebook post.&lt;/a&gt; To demonstrate Santayana&#39;s maxim still holds, a teacher in Brooklyn did not learn from history, and seems to be &lt;a href=&quot;http://nymag.com/daily/intel/2011/03/brooklyn_teacher.html&quot;&gt;in danger of losing her job for the same reason.&lt;/a&gt;&lt;br /&gt;
&lt;br /&gt;
Mairi Mallon, founder and managing director of &lt;a href=&quot;http://www.rein4ce.co.uk/index.php&quot;&gt;rein4ce&lt;/a&gt;, and @reinsurancegirl on twitter, has recently authored &lt;a href=&quot;http://www.rein4ce.co.uk/blog/?p=315&quot;&gt;a detailed post&lt;/a&gt; on her blog about this very topic (from which I shamelessly stole the idea). I very highly recommend anyone even tangentially involved in social media read it. Better yet, subscribe to her blog (which can be found in the blogroll at the side of this post). I think her summary is important enough that it bears repetition:&lt;br /&gt;
&lt;br /&gt;
&lt;blockquote&gt;&lt;ol&gt;&lt;li&gt;Know      and follow our &lt;i&gt;(INSERT LINK TO OWN CORPORATE GUIDELINES) &lt;/i&gt;corporate guidelines . The      same rules apply online.&lt;/li&gt;
&lt;li&gt;Users      are personally responsible for the content they publish  on blogs, Facebook      LinkedIn, Twitter or any other form of  user-generated media.&lt;/li&gt;
&lt;li&gt;Identify      yourself—name and, when relevant, role within the  organisation—when you      discuss company or company-related matters.  You must make it clear that      you are speaking for yourself and not  on behalf of the company.&lt;/li&gt;
&lt;li&gt;Respect      copyright, fair use and financial disclosure laws.&lt;/li&gt;
&lt;li&gt;Don’t      provide our or another’s confidential or other  proprietary information.      Ask permission to publish or report on  conversations that may be deemed to      be private or internal to the  company.&lt;/li&gt;
&lt;li&gt;Don’t      cite or reference clients, partners or suppliers without  their approval.      When you do make a reference, where possible link  back to the source.&lt;/li&gt;
&lt;li&gt;Don’t      use ethnic slurs, personal insults, obscenity, or engage  in any conduct      that would not be acceptable in our workplace (Mallon 2011).&lt;/li&gt;
&lt;/ol&gt;&lt;/blockquote&gt;Engaging the social media can be exciting and rewarding, both on a personal and a professional level; make certain that you don&#39;t hurt yourself, your image, your integrity, and your future when you do so. Everyone wants to be the next Richard Branson; I&#39;m less certain about who wants to be the next &lt;a href=&quot;http://dealbook.nytimes.com/2006/10/19/the-resume-mocked-round-the-world-vayner-speaks/&quot;&gt;Aleksey Vayner.&lt;/a&gt;&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;font-family: Verdana,sans-serif;&quot;&gt;&lt;span style=&quot;font-size: x-small;&quot;&gt;References:&lt;/span&gt;&lt;/div&gt;&lt;span style=&quot;font-size: x-small;&quot;&gt;&lt;span style=&quot;font-family: Verdana,sans-serif;&quot;&gt;Mallon, Mairi, “Social Media guidelines — a freebie for insurance and reinsurance bods.” Weblog entry. Reinsurance girl&#39;s blog. June 6, 2011. June 6, 2011 (&lt;a href=&quot;http://www.rein4ce.co.uk/blog/?p=315&quot;&gt;http://www.rein4ce.co.uk/blog/?p=315&lt;/a&gt;). &lt;/span&gt;&lt;/span&gt;</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/4698339451386033903/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/06/social-media-guidelines-specific-for.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/4698339451386033903'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/4698339451386033903'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/06/social-media-guidelines-specific-for.html' title='Social Media guidelines specific for (re)insurers'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-1550459023621215991</id><published>2011-06-05T15:49:00.000-04:00</published><updated>2011-06-05T15:49:29.485-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="accounting"/><category scheme="http://www.blogger.com/atom/ns#" term="actuaries"/><category scheme="http://www.blogger.com/atom/ns#" term="insurance"/><category scheme="http://www.blogger.com/atom/ns#" term="reinsurance"/><category scheme="http://www.blogger.com/atom/ns#" term="solvency"/><title type='text'>Solvency 2, what I don&#39;t understand</title><content type='html'>There has been much fanfare about &lt;a href=&quot;https://secure.wikimedia.org/wikipedia/en/wiki/Solvency_II&quot;&gt;Solvency II&lt;/a&gt;, and how it will be coming to the US now that it is pretty much &lt;i&gt;fait accompli&lt;/i&gt; to be a requirement for UK and European insurance companies starting in 2013. Perhaps I am missing something basic, but blow are some of the issues that I have with Solvency II:&lt;br /&gt;
&lt;ol&gt;&lt;li&gt;Isn&#39;t Solvency II based on the same principles as Basel II/III? Now how well has &lt;i&gt;that&lt;/i&gt; worked for Europe. Greece, Iceland, Ireland anyone?&lt;/li&gt;
&lt;li&gt;The first pillar of Solvency II (&lt;a href=&quot;http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2009:335:0001:0155:EN:PDF&quot;&gt;Article 75 (1)(b)&lt;/a&gt;) requires that &quot;...liabilities shall be valued at the amount for which they could be transferred, or settled, between knowledgeable willing par­ties in an arm’s length transaction.&quot; I don&#39;t know about the rest of you, but there really is no liquid market for (re)insurer loss reserves. I&#39;ve been involved in pricing a few loss portfolio attempts--none of which came to fruition mind you--and each one is really a bespoke transaction. Different counter-parties to the same transaction will arrive at different values for the reserves. So how is Solvency II going to handle this?&lt;/li&gt;
&lt;li&gt;At its heart, the capital requirement is still a Value at Risk (VaR) measure, albeit at the 99.5%-ile. When will people learn that VaR is a rather non-robust statistic? It is a &lt;b&gt;point&lt;/b&gt; on the cumulative frequency distribution, with no &quot;knowledge&quot; of what is above or below it. It is very susceptible to discontinuities (think step function), and its components are non-additive. At the very least, the measure should be based on TVaR, or the expected value above a given point. As an expectation, it is additive in its components (Co-TVaR measures exist and are meaningful) and it, as a first moment, reflects to some extent the entire distribution in the tail above the point, not just a point. As we all know, actually seeing any particular result from a continuous distribution &lt;a href=&quot;https://secure.wikimedia.org/wikipedia/en/wiki/Almost_surely&quot;&gt;almost never&lt;/a&gt; happens, which is why we talk intervals and not points.&lt;/li&gt;
&lt;li&gt;While &lt;a href=&quot;https://secure.wikimedia.org/wikipedia/en/wiki/Operational_risk&quot;&gt;Operational Risk&lt;/a&gt; is certainly a significant factor in a company&#39;s risk profile, I have yet to see any good measure or process to quantify the expected value of said risk.&amp;nbsp;&lt;/li&gt;
&lt;/ol&gt;Unlike the banking industry, the insurance industry, especially in the US, has been rather stable for the past century or so, and while I agree that certain elements of Statutory Accounting, such as the holding of nominal reserves, may seem somewhat extreme with respect to other businesses, it has done rather well in ensuring a viable and solvent insurance industry. Why are we rushing into something which has neither the provenance (Basel) nor the track record (not even officially implemented in Europe) to justify adding a new accounting and solvency standard.&lt;br /&gt;
&lt;ol&gt;&lt;/ol&gt;</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/1550459023621215991/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/06/solvency-2-what-i-dont-understand.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/1550459023621215991'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/1550459023621215991'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/06/solvency-2-what-i-dont-understand.html' title='Solvency 2, what I don&#39;t understand'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-6520476795162949277</id><published>2011-05-18T20:48:00.000-04:00</published><updated>2011-05-18T20:48:41.493-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="blog"/><category scheme="http://www.blogger.com/atom/ns#" term="rstats"/><title type='text'>R-bloggers</title><content type='html'>As I decided to try and blog a little more often now, and touch on &quot;R&quot; every now and then, I decided to take R-bloggers up on their standing offer to include R-related feeds at their site. So, everything I tag with &quot;rstats&quot; (you can guess where &lt;i&gt;that &lt;/i&gt;came from) should flow through to them. I&#39;ve added them to my (tiny) blogroll at the side of the blog, but if you just cannot wait to see what they are all about, and I recommend it, you can just go &lt;a href=&quot;http://www.r-bloggers.com/&quot;&gt;here&lt;/a&gt;.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/6520476795162949277/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/05/r-bloggers.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6520476795162949277'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6520476795162949277'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/05/r-bloggers.html' title='R-bloggers'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-722934309125754226</id><published>2011-05-15T20:37:00.001-04:00</published><updated>2011-05-15T20:42:44.482-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="actuaries"/><category scheme="http://www.blogger.com/atom/ns#" term="curve fitting"/><category scheme="http://www.blogger.com/atom/ns#" term="loss distributions"/><category scheme="http://www.blogger.com/atom/ns#" term="rstats"/><title type='text'>Why method of moments doesn&#39;t always work</title><content type='html'>A number of years ago, someone asked me &quot;why does my company need actuaries to fit curves, once I have the mean and standard deviation of my losses, isn&#39;t that enough?&quot; I explained to him that not every distribution is completely determined by its mean and standard deviation (as the normal and lognormal are), and as at that point, I did not have &quot;R&quot; installed on my laptop, I demonstrated it to him in Excel. Having wanted to start blogging about &quot;R&quot;, even ever so infrequently, I figured I&#39;d toss together a little code to demonstrate.&lt;br /&gt;
&lt;br /&gt;
The example I gave was to compare a gamma and a pareto distribution, each of which has mean 10,000 and a CV of 150% (making the standard deviation 15,000). I will spare all of you the algebra, but suffice to say, that using the Klugman-Panjer-Wilmot parameterization (which is used by most casualty actuaries in the past 20 years or so) the parameters of the gamma would be theta (R&#39;s scale) = 22500 and alpha (R&#39;s shape) = 4/9. The equivalent pareto would have theta (R&#39;s scale) = 26000 and alpha (R&#39;s shape) = 3.6.&lt;br /&gt;
&lt;br /&gt;
&lt;div class=&quot;separator&quot; style=&quot;clear: both; text-align: center;&quot;&gt;&lt;a href=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhm-VrzbJzNNmcJwsdACaslYjYzCSx7bDGwR7Fh20p4ldtsG6RkmlYLjnVvO8BglRUqqA14q7cVUc0Jvre3yBFaWb0bF9oFYyoGkNUlyJV5hivc3T0aF6qTsQJnXxRH2XBoUprJ/s1600/Pareto+vs.+Gamma.png&quot; imageanchor=&quot;1&quot; style=&quot;margin-left: 1em; margin-right: 1em;&quot;&gt;&lt;/a&gt;&lt;/div&gt;Graphing the two (and Hadley, please forgive me for using default R&#39; plotting, I left my ggplot book in the office; &lt;i&gt;mea culpa&lt;/i&gt;) you can easily see how the distributions are rather different.&lt;br /&gt;
&lt;br /&gt;
&lt;a href=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhm-VrzbJzNNmcJwsdACaslYjYzCSx7bDGwR7Fh20p4ldtsG6RkmlYLjnVvO8BglRUqqA14q7cVUc0Jvre3yBFaWb0bF9oFYyoGkNUlyJV5hivc3T0aF6qTsQJnXxRH2XBoUprJ/s1600/Pareto+vs.+Gamma.png&quot; style=&quot;margin-left: 1em; margin-right: 1em;&quot;&gt;&lt;img border=&quot;0&quot; height=&quot;320&quot; src=&quot;https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhm-VrzbJzNNmcJwsdACaslYjYzCSx7bDGwR7Fh20p4ldtsG6RkmlYLjnVvO8BglRUqqA14q7cVUc0Jvre3yBFaWb0bF9oFYyoGkNUlyJV5hivc3T0aF6qTsQJnXxRH2XBoUprJ/s320/Pareto+vs.+Gamma.png&quot; width=&quot;320&quot; /&gt;&amp;nbsp;&lt;/a&gt;&lt;br /&gt;
&lt;br /&gt;
To make things easier for me, I used the &lt;span style=&quot;font-family: &amp;quot;Courier New&amp;quot;,Courier,monospace;&quot;&gt;actuar&lt;/span&gt;&amp;nbsp; package to do the graphing:&lt;br /&gt;
&lt;br /&gt;
&lt;div style=&quot;overflow: auto;&quot;&gt;&lt;div class=&quot;geshifilter&quot;&gt;&lt;pre class=&quot;r geshifilter-R&quot; style=&quot;font-family: monospace;&quot;&gt;&lt;span style=&quot;font-size: small;&quot;&gt;&lt;a href=&quot;http://inside-r.org/r-doc/base/library&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;library&lt;/span&gt;&lt;/a&gt;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;&lt;a href=&quot;http://inside-r.org/packages/cran/actuar&quot;&gt;actuar&lt;/a&gt;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;
&lt;a href=&quot;http://inside-r.org/r-doc/graphics/curve&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;curve&lt;/span&gt;&lt;/a&gt;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;dpareto&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;x&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; &lt;a href=&quot;http://inside-r.org/packages/cran/shape&quot;&gt;shape&lt;/a&gt;=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;3.6&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/scale&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;scale&lt;/span&gt;&lt;/a&gt;=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;26000&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; from=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; to=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;100000&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/col&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;col&lt;/span&gt;&lt;/a&gt;=&lt;/span&gt;&lt;span style=&quot;color: blue; font-size: small;&quot;&gt;&quot;blue&quot;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;
&lt;a href=&quot;http://inside-r.org/r-doc/graphics/curve&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;curve&lt;/span&gt;&lt;/a&gt;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;&lt;a href=&quot;http://inside-r.org/r-doc/stats/dgamma&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;dgamma&lt;/span&gt;&lt;/a&gt;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;(&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;x&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; &lt;a href=&quot;http://inside-r.org/packages/cran/shape&quot;&gt;shape&lt;/a&gt;=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;4&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt;/&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;9&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/scale&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;scale&lt;/span&gt;&lt;/a&gt;=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;22500&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;)&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; from=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;0&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; to=&lt;/span&gt;&lt;span style=&quot;color: #cc66cc; font-size: small;&quot;&gt;100000&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; add=&lt;/span&gt;&lt;span style=&quot;color: black; font-size: small; font-weight: bold;&quot;&gt;TRUE&lt;/span&gt;&lt;span style=&quot;color: #339933; font-size: small;&quot;&gt;,&lt;/span&gt;&lt;span style=&quot;font-size: small;&quot;&gt; &lt;a href=&quot;http://inside-r.org/r-doc/base/col&quot;&gt;&lt;span style=&quot;color: #003399; font-weight: bold;&quot;&gt;col&lt;/span&gt;&lt;/a&gt;=&lt;/span&gt;&lt;span style=&quot;color: blue; font-size: small;&quot;&gt;&quot;green&quot;&lt;/span&gt;&lt;span style=&quot;color: #009900; font-size: small;&quot;&gt;)&lt;/span&gt;&lt;/pre&gt;&lt;/div&gt;&lt;/div&gt;&lt;a href=&quot;http://www.inside-r.org/pretty-r&quot; title=&quot;Created by Pretty R at inside-R.org&quot;&gt;Created by Pretty R at inside-R.org&lt;/a&gt;&lt;br /&gt;
&lt;br /&gt;
Obviously, the tails of the distributions, and thus the survival function at a given loss size, is different for the two, notwithstanding their sharing identical first two moments. So, this was just a brief but effective visualization as to how the first two moments do not contain all the information needed to find a &quot;best fit,&quot; and why we like to use distributional fitting methods (maximum likelihood, maximum spacing, various minimum distance metrics like Cramer-von Mises, etc.) to get a better understanding of the potential underlying loss processes.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/722934309125754226/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/05/why-method-of-moments-doesnt-always.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/722934309125754226'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/722934309125754226'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/05/why-method-of-moments-doesnt-always.html' title='Why method of moments doesn&#39;t always work'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><media:thumbnail xmlns:media="http://search.yahoo.com/mrss/" url="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhm-VrzbJzNNmcJwsdACaslYjYzCSx7bDGwR7Fh20p4ldtsG6RkmlYLjnVvO8BglRUqqA14q7cVUc0Jvre3yBFaWb0bF9oFYyoGkNUlyJV5hivc3T0aF6qTsQJnXxRH2XBoUprJ/s72-c/Pareto+vs.+Gamma.png" height="72" width="72"/><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-390697585260338937</id><published>2011-05-13T14:42:00.001-04:00</published><updated>2011-05-13T14:43:50.757-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="blogger"/><category scheme="http://www.blogger.com/atom/ns#" term="web design"/><title type='text'>New look!</title><content type='html'>I decided to change things around a bit, and applied a new template and color scheme. Maybe now my next post won&#39;t take another five months! Let me know if you like the new look.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/390697585260338937/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/05/new-look.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/390697585260338937'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/390697585260338937'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/05/new-look.html' title='New look!'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-3917573457649168280</id><published>2011-01-12T10:39:00.002-05:00</published><updated>2011-01-12T10:42:14.495-05:00</updated><title type='text'>It has been a while</title><content type='html'>It has been quite a while since I posted an update here. Possibly, because the longer I work, the less time I seem to have. Just to share something, I did recently have some fun programming a random gamma variate generator in VBA for some copula work. While I&#39;m certain there are much more elegant implementations out there, it was good to figuratively roll up my sleeves and dust off some VBA programming skills. Now I can keep that line on the resume!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/3917573457649168280/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2011/01/it-has-been-while.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3917573457649168280'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3917573457649168280'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2011/01/it-has-been-while.html' title='It has been a while'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-6452205265501710972</id><published>2009-09-22T17:07:00.003-04:00</published><updated>2009-09-22T17:10:26.421-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="firefox"/><category scheme="http://www.blogger.com/atom/ns#" term="IE6"/><title type='text'>IE6 no more</title><content type='html'>I just came across this article, entitled &lt;a class=&quot;fn url&quot; title=&quot;Permanent Link to IE6 No More! Popular Web Companies Start Project to Kill IE6&quot; href=&quot;http://mashable.com/2009/08/04/ie6-no-more/&quot; rel=&quot;bookmark&quot;&gt;&lt;span style=&quot;color: rgb(0, 0, 0);&quot;&gt;&quot;&lt;span style=&quot;color: rgb(51, 51, 255);&quot;&gt;IE6 No More!  Popular Web Companies Start Project to Kill IE6&lt;/span&gt;&quot;&lt;/span&gt;&lt;/a&gt; and decided to add the code to my blog, so if you are using IE6, you may see a note to that effect. FireFox is the way to go anyway!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/6452205265501710972/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/09/ie6-no-more.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6452205265501710972'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6452205265501710972'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/09/ie6-no-more.html' title='IE6 no more'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-6203078068833954493</id><published>2009-09-21T16:49:00.003-04:00</published><updated>2009-09-21T16:53:19.660-04:00</updated><title type='text'>Hey, that sounds like my studying routine!</title><content type='html'>In my daily perusal of the intellectual side of the news, the comics, I came across another gem, which while ostensibly applied to elementary school, may be the best actuarial study technique ever printed. Enjoy!&lt;br /&gt;&lt;br /&gt;&lt;a href=&quot;http://d.yimg.com/a/p/umedia/20090920/largeimage.7b0f57af9a54b45f31e113377789cccb.gif&quot;&gt;http://d.yimg.com/a/p/umedia/20090920/largeimage.7b0f57af9a54b45f31e113377789cccb.gif&lt;/a&gt;&lt;br /&gt;&lt;br /&gt;&lt;a onblur=&quot;try {parent.deselectBloggerImageGracefully();} catch(e) {}&quot; href=&quot;http://d.yimg.com/a/p/umedia/20090920/largeimage.7b0f57af9a54b45f31e113377789cccb.gif&quot;&gt;&lt;img style=&quot;margin: 0pt 0pt 10px 10px; float: right; cursor: pointer; width: 503px; height: 233px;&quot; src=&quot;http://d.yimg.com/a/p/umedia/20090920/largeimage.7b0f57af9a54b45f31e113377789cccb.gif&quot; alt=&quot;&quot; border=&quot;0&quot; /&gt;&lt;/a&gt;&lt;br /&gt;&lt;br /&gt;&lt;img src=&quot;file:///C:/DOCUME%7E1/aadler00/LOCALS%7E1/Temp/moz-screenshot.png&quot; alt=&quot;&quot; /&gt;</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/6203078068833954493/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/09/hey-that-sounds-like-my-studying.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6203078068833954493'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6203078068833954493'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/09/hey-that-sounds-like-my-studying.html' title='Hey, that sounds like my studying routine!'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-4007436530990338626</id><published>2009-09-10T16:49:00.004-04:00</published><updated>2009-09-10T16:51:25.963-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="actuaries"/><category scheme="http://www.blogger.com/atom/ns#" term="cartoon"/><title type='text'>Actuaries in the media…</title><content type='html'>…well at least in cartoons. Unfortunately, even we are not immune from politics.&lt;br /&gt;&lt;br /&gt;&lt;a href=&quot;http://d.yimg.com/a/p/uc/20090910/largeimagepo090909.gif&quot;&gt;http://d.yimg.com/a/p/uc/20090910/largeimagepo090909.gif&lt;/a&gt;</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/4007436530990338626/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/09/actuaries-in-media.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/4007436530990338626'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/4007436530990338626'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/09/actuaries-in-media.html' title='Actuaries in the media…'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-5461402336974051388</id><published>2009-09-09T12:09:00.004-04:00</published><updated>2009-09-09T12:18:38.221-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="actuary"/><category scheme="http://www.blogger.com/atom/ns#" term="CAS"/><category scheme="http://www.blogger.com/atom/ns#" term="CIA"/><category scheme="http://www.blogger.com/atom/ns#" term="FEM"/><category scheme="http://www.blogger.com/atom/ns#" term="SOA"/><title type='text'>Future Education Methods and the Actuarial Profession</title><content type='html'>Recently, the &quot;big three&quot; actuarial organizations in North America, the CAS, SOA, and CIA, have promoted a proposal to radically change the actuarial education methods used to credential actuaries in the United States and Canada. Their letter describing the suggestions may be found &lt;a href=&quot;http://www.soa.org/files/pdf/fem-letter.pdf&quot;&gt;here&lt;/a&gt;.&lt;br /&gt;&lt;br /&gt;Personally, I think it is a very bad idea, and I sent the following to the feedback address at &lt;a href=&quot;mailto:fem@actuarialdirectory.org&quot;&gt;fem@actuarialdirectory.org&lt;/a&gt;. The deadline for responses is tomorrow, so if you have not already commented, you have one more day.&lt;br /&gt;&lt;br /&gt;&lt;blockquote style=&quot;font-family: arial;&quot;&gt;I believe that the proposed future education method (FEM) of exempting actuaries from examinations based on passing specific university courses is a bad idea, and I am strongly opposed to it.&lt;br /&gt;&lt;br /&gt;Currently, one of the strongest elements of our profession is its being close to a true meritocracy; much more so than any given university. A candidate takes an examination simultaneously with his or her peers, has his or her exam anonymized, and then graded by anonymous referees (anonymous to the candidates) similar to a double-blind test, and has his or her results, still anonymous, compared with the results of all other test-takers. Then an pass score is set (we will leave the discussion of the score being set based on the actual distribution of results to another time) and all who have met or exceeded the score will be deemed to have net the requirements. There is very little opportunity for either &quot;gaming the system&quot; or for allowing any other external influence (pernicious, malicious, or otherwise) to effect the outcomes. Contrast the above with a university setting. Notwithstanding the various societies&#39; intentions of ensuring &quot;specific standard&quot;, it is undeniable that different teachers, and different universities for that matter, will have different standards. Furthermore, this opens the hideous potential for student-teacher relationships impacting the results, be they for good or ill. We lose the anonymity, we lose the simultaneity, we lose the ability to compare oneself to one&#39;s peers, we lose the &quot;double-blind&quot; nature of the exam process.&lt;br /&gt;&lt;br /&gt;Another distinct shortcoming of the new process is the complete loss of standardization. Currently, all who pass Exam &quot;X&quot; in a given sitting do so with the exact same exam. Symmetrically, all those who fail Exam &quot;X&quot; in a given sitting, fail the same exam that those who pass have passed. When a change is made to the syllabus, and thus the exams, it will affect all sitting candidates in the same way. There is no preferential treatment possible. Under the new system, there can be many different &quot;exams;&quot; as many as there are approved courses. The syllabus, and any changes thereto, will not affect all candidates equally, but is completely dependant on the institution, and for that matter the teacher. This also includes the potential dilution of candidates, in that their passing is no longer controlled by the respective societies, but by external institutions.&lt;br /&gt;&lt;br /&gt;A third distinct shortcoming of this new process, while less egregious, is the loss of its screening abilities. As a former employer of mine, a credentialed actuary, once told me in half-jest, &quot;you do not pass actuarial exams, you win them.&quot; One of the benefits of the exam process, and I say this in hindsight having completed it, is its ability to isolate people who can deliver good answers under pressure. I cannot speak for anyone else, but to me, this is pretty much the epitome of what we as actuaries are called upon to do. We are not supposed to always have the best answer, but we are supposed to be able to help our clients (be they actual clients, employers, etc.) to address their concerns as best we can in the ever-shortening time-frames we are given. Furthermore, in actual work, we are not always given all the information &quot;handed&quot; to us on a platter. We need to be able to quickly identify good and useful data from bad, and use it as best we can. Exam work is very similar to this. We have to learn the material on our own—exam seminars are good refreshers, but I have never met an actuary who has only used them without other studying who has had any success—which includes triaging the material, identifying the areas on which to focus (the Brosious quartering comes to mind), and quickly &quot;ingesting&quot; what is needed to achieve success. In a university setting, much of this is lost. For all we know, the final can be a one week, open book, take home affair—unthinkable under the current method.&lt;br /&gt;&lt;br /&gt;Regarding the purported benefits of such a system, I must respectfully disagree with some.&lt;br /&gt;&lt;br /&gt;I do not see the quality of actuarial education improving any more because of institutionalized exam exceptions than it would should the current syllabus and exams be appropriately updated. Some of our syllabus papers, at least when I was a candidate, were written not for candidates but for actuaries, and experienced ones at those. What we need is a centralized effort at generating better exam material, not diffusing that effort throughout 20 or so universities, each with their own agenda and intentions.&lt;br /&gt;&lt;br /&gt;The following, I feel, highlights a critical chasm between academic actuaries and practicing actuaries. I would like to know why &quot;university&quot; training in actuarial methods is considered a benefit? Many would argue that in truth it is a detriment. Universities are often known for their being &quot;out-of-touch&quot; with what is occurring in the &quot;trenches&quot; as it were. The ivory towers of academia are bastions of theory, but oft-times are highly impractical. As actuaries, most of us will practice, and not remain theoreticians, and as such, I would counter that one of the most important elements in becoming and being an actuary is the &quot;apprenticeship&quot; one undergoes in one&#39;s early actuarial career. Knowledge should never be confused with, or substituted for, experience, and it is experience that is key in our discipline, which at times is as much an art as it is a science (at least in my field of casualty reinsurance).&lt;br /&gt;&lt;br /&gt;I can understand how allowing university classes to exempt students from actuarial exams would increase the career&#39;s attractiveness. However, that can be addressed by having the university courses geared towards the specific exams and preparing their students for taking said exam.&lt;br /&gt;&lt;br /&gt;In closing, I feel that the FEM described would do far more harm than good to our profession, and I urge the societies to reject it in favor of enhancing our current process.&lt;br /&gt;&lt;br /&gt;Thank you.&lt;/blockquote&gt;</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/5461402336974051388/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/09/future-education-methods-and-actuarial.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/5461402336974051388'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/5461402336974051388'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/09/future-education-methods-and-actuarial.html' title='Future Education Methods and the Actuarial Profession'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-6595142801665889100</id><published>2009-04-21T00:10:00.002-04:00</published><updated>2009-04-21T00:13:44.063-04:00</updated><title type='text'>Cassini photos of Saturn and its environs</title><content type='html'>Cassini photos of Saturn, its rings, and its moons, as seen here &lt;a href=&quot;http://www.boston.com/bigpicture/2009/04/cassinis_continued_mission.html&quot;&gt;http://www.boston.com/bigpicture/2009/04/cassinis_continued_mission.html&lt;/a&gt;, are simultaneously amazing and awe-inspiring. It is further serendipity that I am currently reading Jack McDevitt&#39;s Academy series; the one where the first book features the monument on Iapetus, whose surface is captured in glorious detail by Cassini.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/6595142801665889100/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/04/cassini-photos-of-saturn-and-its.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6595142801665889100'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6595142801665889100'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/04/cassini-photos-of-saturn-and-its.html' title='Cassini photos of Saturn and its environs'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-3272645234902938672</id><published>2009-04-06T10:23:00.001-04:00</published><updated>2009-04-06T10:23:34.732-04:00</updated><title type='text'>Another bit.ly test</title><content type='html'>Hmm, the API key didn&#39;t work. One more try.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/3272645234902938672/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/04/another-bitly-test.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3272645234902938672'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3272645234902938672'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/04/another-bitly-test.html' title='Another bit.ly test'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-3553313149515933205</id><published>2009-04-06T02:18:00.002-04:00</published><updated>2009-04-06T02:21:37.484-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="bit.ly"/><category scheme="http://www.blogger.com/atom/ns#" term="identi.ca"/><category scheme="http://www.blogger.com/atom/ns#" term="pingfm"/><category scheme="http://www.blogger.com/atom/ns#" term="twitter"/><category scheme="http://www.blogger.com/atom/ns#" term="twitterfeed"/><title type='text'>Bit.Ly test</title><content type='html'>Now that I have twitterfeed posting to ping.fm, and from there to twitter and identi.ca, I&#39;m trying one last time to get a shortened URL that does not contain &quot;ping.fm&quot; in it, for people behind websense or the like. This post will check the activity of a &lt;a href=&quot;http://www.bit.ly/&quot;&gt;bit.ly&lt;/a&gt; account (not anonymous) as the shortener in twitterfeed. Let&#39;s hope this works!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/3553313149515933205/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/04/bitly-test.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3553313149515933205'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3553313149515933205'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/04/bitly-test.html' title='Bit.Ly test'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-226922836966940751</id><published>2009-04-05T20:37:00.001-04:00</published><updated>2009-04-05T20:40:38.055-04:00</updated><title type='text'>Twitterfeed Test</title><content type='html'>Well, twitterfeed through Ping seemed to update both Twitter and Identica, but the url shortener has a ping domain which is blocked at work. Trying something else now. This may be a many-post day trying to get this set up.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/226922836966940751/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/04/twitterfeed-test.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/226922836966940751'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/226922836966940751'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/04/twitterfeed-test.html' title='Twitterfeed Test'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-6444917069213205590</id><published>2009-04-05T12:11:00.003-04:00</published><updated>2009-04-05T12:16:23.614-04:00</updated><title type='text'>Another test</title><content type='html'>Well, it seems that &lt;a href=&quot;http://ping.fm/dashboard/&quot;&gt;Ping.fm&lt;/a&gt; is great for cross-posting my &lt;a href=&quot;http://twitter.com/timeline/home&quot;&gt;tweets&lt;/a&gt; to &lt;a href=&quot;http://identi.ca/&quot;&gt;identi.ca&lt;/a&gt; but it didn&#39;t take the previous blog post. So, let&#39;s try &lt;a href=&quot;http://twitterfeed.com/&quot;&gt;Twitterfeed&lt;/a&gt; feeding Ping.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/6444917069213205590/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/04/another-test.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6444917069213205590'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/6444917069213205590'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/04/another-test.html' title='Another test'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-3571319653444394285</id><published>2009-04-05T09:56:00.001-04:00</published><updated>2009-04-05T10:01:54.776-04:00</updated><title type='text'>Ping.Fm and the social-interaction revolution</title><content type='html'>We are witnessing a social-interaction revolution. First blogging, and then micro-blogging, have taken cyberspace by storm. With so many different options and niches, it may be difficult to manage all of one&#39;s separate service providers. There is a service, &lt;a href=&quot;http://ping.fm/dashboard/&quot;&gt;Ping.Fm&lt;/a&gt;, which promises to handle that. You can set it to accept feeds from one service and cross-post them to others. You can specify three different kinds of feeds (status, blog, and microblog) to decide which goes where; you can set up special tags to post to s defined subset of feeds, or you can even specify a particular feed. This blog entry is being made to test if Ping.fm is going to handle blogger posts the way I expect. Enjoy!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/3571319653444394285/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/04/pingfm-and-social-interaction.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3571319653444394285'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/3571319653444394285'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/04/pingfm-and-social-interaction.html' title='Ping.Fm and the social-interaction revolution'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-7676513952989233045</id><published>2009-03-31T00:39:00.005-04:00</published><updated>2009-03-31T01:20:09.304-04:00</updated><category scheme="http://www.blogger.com/atom/ns#" term="currency"/><category scheme="http://www.blogger.com/atom/ns#" term="economy"/><title type='text'>A single currency?</title><content type='html'>On &lt;a href=&quot;http://www.bloomberg.com/tvradio/radio/&quot;&gt;Bloomberg Surveillance&lt;/a&gt; this morning the concept of a single world-wide currency was discussed. Personally, I am not sure how that would even be possible. How would investors or the market be able to reflect their opinions as to the relative security and stability of various governments?&lt;br /&gt;&lt;br /&gt;Even with a single world-wide currency, governments and countries will have different levels of output, stability, production, etc. There needs to be a way for the markets to reflect this. Currently, this can be done by the relative value of one currency versus another. Now consider every country has just one currency. The governments, natural resources, and indigenous economies are still different. Unless the currency is 100% backed by some universal commodity (gold, for example) part of the value of the currency is the &quot;full faith and credit&quot; in the issuing body, and that the currency note is redeemable for equivalent goods and services across national borders. That discrepancy between entities still exists regardless of shape, color, size, and portraits on the bills or coins in use.&lt;br /&gt;&lt;br /&gt;Perhaps the discrepancies will become evident in different interest rates that need to be offered by the treasuries of the respective countries; a government with a stronger economy can offer lower rates in return for the inherent stability. Well, doesn&#39;t that make government securities in and of themselves a &quot;shadow currency&quot;? We can have the two bonds with the same par value, the same coupon, and the same duration with different prices, due to the perceived difference in risk and stability. I&#39;d guess that pretty soon thereafter a market for trading these identical yet different instruments will be created, and we are pretty much back where we started.&lt;br /&gt;&lt;br /&gt;As this is just semi-knowledgeable ramblings on my part, I&#39;m curious to know how investors bet the German economy against the French, if both use the Euro. While an &lt;a href=&quot;http://en.wikipedia.org/wiki/Optimum_currency_area&quot;&gt;optimal currency area&lt;/a&gt; may be larger than a country, I&#39;m pretty sure it is smaller than the entire world.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/7676513952989233045/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2009/03/single-currency.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/7676513952989233045'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/7676513952989233045'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2009/03/single-currency.html' title='A single currency?'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-7275746276400633589</id><published>2008-06-18T00:46:00.002-04:00</published><updated>2008-06-18T00:54:19.725-04:00</updated><title type='text'>AIG or ARRGGHH?</title><content type='html'>Interesting news today, from &lt;a href=&quot;http://www.cnbc.com/id/25177201&quot;&gt;CNBC&lt;/a&gt;:&lt;br /&gt;&lt;blockquote&gt;American International Group CEO Martin Sullivan stepped down from his post, as expected, after growing discontent among the board and shareholders with his management of the insurance giant.…Sullivan has come under intense criticism since AIG recently announced billions of dollars in writedowns from losses after Sullivan and the company assured investors the writedowns would be minimal.&lt;/blockquote&gt;&lt;br /&gt;&lt;br /&gt;This is another example of where a good ERM study would have exposed the asset weakness. As &lt;a href=&quot;http://www.economist.com/finance/displaystory.cfm?story_id=11565362&quot;&gt;The Economist&lt;/a&gt; writes:&lt;br /&gt;&lt;blockquote&gt;Mr Sullivan did himself no favours by mismanaging expectations. He and his team poo-pooed critics, insisting until this year that AIG had oodles of excess capital and that its actual (as opposed to mark-to-market) losses would be modest. The firm has since been forced to admit that its accounting models were too optimistic, after its auditors found “material weaknesses” in them&lt;/blockquote&gt;&lt;br /&gt;&lt;br /&gt;Perhaps an outside analysis of the capital model may have been a good idea?</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/7275746276400633589/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2008/06/aig-or-arrgghh.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/7275746276400633589'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/7275746276400633589'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2008/06/aig-or-arrgghh.html' title='AIG or ARRGGHH?'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-2636931629374599312</id><published>2007-09-04T13:37:00.000-04:00</published><updated>2007-09-04T13:41:41.308-04:00</updated><title type='text'>MA INSURANCE COMMISSIONER PROPOSES BANNING USE OF SOCIOECONOMIC FACTORS IN RATING AND UNDERWRITING</title><content type='html'>From &lt;a href=&quot;http://www.insurancejournal.com/news/east/2007/08/29/83069.htm&quot;&gt;Insurance Journal - August 29, 2007&lt;/a&gt;&lt;br /&gt;&lt;br /&gt;&lt;blockquote&gt;The draft regulation makes years of driving experience, driving record and a vehicle&#39;s model and safety features as the primary rating and underwriting factors in private passenger auto insurance.&lt;br /&gt;&lt;br /&gt;The draft regulation bans insurers from using gender, marital status, education, occupation, national origin, religion, homeownership and other socioeconomic factors, some of which Burnes says are prohibited by statute and others which she deemed violate public policy, for either rating or underwriting.&lt;br /&gt;&lt;br /&gt;It also forbids insurers from using information from credit reports or scores for rating for at least a year, during which time she said she would study the overall use of credit scores by insurers. It does not prohibit the use of credt scoring for underwriting.&lt;/blockquote&gt;&lt;br /&gt;&lt;br /&gt;This does mean that subsidization will still be alive and well. Let&#39;s hope that there is less of it.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/2636931629374599312/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2007/09/ma-insurance-commissioner-proposes.html#comment-form' title='1 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/2636931629374599312'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/2636931629374599312'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2007/09/ma-insurance-commissioner-proposes.html' title='MA INSURANCE COMMISSIONER PROPOSES BANNING USE OF SOCIOECONOMIC FACTORS IN RATING AND UNDERWRITING'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>1</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-9138617307841397345</id><published>2007-08-20T12:47:00.000-04:00</published><updated>2007-08-20T12:50:58.976-04:00</updated><title type='text'>Florida turbulence</title><content type='html'>At least &lt;em&gt;some&lt;/em&gt; Florida legislators are seeing the light. &lt;a href=&quot;http://www.floridaunderwriter.com/cms/fl/Monthly%20Issues/Issues/2007/06/Index/Capitol%20Line/Capitol%20Line&quot;&gt;Florida Underwriter&lt;/a&gt; writes:&lt;br /&gt;&lt;blockquote&gt;The resulting 91-page bill swept through the insurance code with all the fury of a major storm that left insurance representatives stunned and some legislator’s in disbelief. Representative Dennis Ross (R-Lakeland) said the bill was “as close to a socialist policy that we ever came to.” He was joined by Representative Don Brown (R-De Funiak Springs) who decried the changes to Citizens Property Insurance Corporation, which makes it certain the former residual market will be the state’s largest insurer. Speaking of the freeze of Citizens’ rates until January 1, 2009, Brown noted that the bill ensured that homeowners would be eventually face potentially crippling assessments. “A one-year rate freeze is playing Russian roulette and every year you continue the freeze you’re adding another bullet to the chamber,” he said.&lt;/blockquote&gt;&lt;br /&gt;&lt;br /&gt;Leave it to Governor Crist, however, to never confuse reality for political expediency:&lt;br /&gt;&lt;blockquote&gt;But no matter how dire the message, Crist and the legislature were more than happy to press forward and place their imprint on the property market. “The Florida House and the Senate have been working very hard to bring relief to our homeowners, and I am grateful for their dedication,” said Crist . “The people have pleaded for relief, and the good men and women of the Florida Legislature have answered the call. This is the right thing to do.” &lt;/blockquote&gt;&lt;br /&gt;The right thing to do? For whom?!</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/9138617307841397345/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2007/08/florida-turbulence.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/9138617307841397345'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/9138617307841397345'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2007/08/florida-turbulence.html' title='Florida turbulence'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-112291126915117931</id><published>2005-08-01T11:45:00.000-04:00</published><updated>2005-08-01T11:47:49.156-04:00</updated><title type='text'>Workers Compensation Data</title><content type='html'>This is interesting, and useful. Needing to do research on workers comp rate filings, I stumbled on this site, which links to near every rating bureaus you can think of: &lt;a href=&quot;http://www.cutcomp.com/depts.htm&quot;&gt;http://www.cutcomp.com/depts.htm&lt;/a&gt;&lt;br /&gt;&lt;br /&gt;&lt;br /&gt;Rather helpful.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/112291126915117931/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2005/08/workers-compensation-data.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/112291126915117931'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/112291126915117931'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2005/08/workers-compensation-data.html' title='Workers Compensation Data'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry><entry><id>tag:blogger.com,1999:blog-7931813.post-112290925954160065</id><published>2005-08-01T11:08:00.000-04:00</published><updated>2005-08-01T11:14:19.546-04:00</updated><title type='text'>The Further adventures of a failed behemoth</title><content type='html'>From the &lt;a href=&quot;http://www.insurancejournal.com/news/national/2005/08/01/57880.htm&quot;&gt;Insurance Journal: August 1, 2005&lt;/a&gt;&lt;br /&gt;&lt;br /&gt;&lt;blockquote&gt;&lt;strong&gt;&lt;span style=&quot;font-family:verdana;font-size:130%;&quot;&gt;Reinsurers Sue AIG Alleging Claims Fraud; AIG&lt;br /&gt;Denies All Wrongdoing&lt;/span&gt;&lt;/strong&gt;&lt;br /&gt;&lt;span style=&quot;font-family:Times New Roman;font-size:100%;&quot;&gt;&lt;p&gt;A group of 18 insurers is suing American International Group and bankrupt&lt;br /&gt;fronting company Trenwick America Reinsurance Corp. for allegedly scheming to&lt;br /&gt;collect as much as $73 million in what the insurers claim are &quot;grossly inflated&quot;&lt;br /&gt;workers compensation and other reinsurance claims.&lt;/p&gt;&lt;br /&gt;&lt;p&gt;The suit alleges that only about $15 million of the $73 million in claims for which AIG has demanded payment appear to be legitimate paid losses eligible for reinsurance coverage. The remaining amount reflects &quot;highly suspect&quot; estimates of future (incurred but not reported) losses, according to the complaint, which was filed in Suffolk Superior Court in Boston on July 6.&lt;/p&gt;&lt;/span&gt;&lt;/blockquote&gt;Click on the above link for the rest of the article.&lt;br /&gt;It is amazing to see what is happening with AIG. It is even more amazing to see just how much AIG &amp;amp; Co. was doing &lt;em&gt;sub rosa. &lt;/em&gt;I am very curious to see what else is unearthed. Stay tuned, this is going to develop for a long time.</content><link rel='replies' type='application/atom+xml' href='http://severity.blogspot.com/feeds/112290925954160065/comments/default' title='Post Comments'/><link rel='replies' type='text/html' href='http://severity.blogspot.com/2005/08/further-adventures-of-failed-behemoth.html#comment-form' title='0 Comments'/><link rel='edit' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/112290925954160065'/><link rel='self' type='application/atom+xml' href='http://www.blogger.com/feeds/7931813/posts/default/112290925954160065'/><link rel='alternate' type='text/html' href='http://severity.blogspot.com/2005/08/further-adventures-of-failed-behemoth.html' title='The Further adventures of a failed behemoth'/><author><name>Avi</name><uri>http://www.blogger.com/profile/17021838649882950034</uri><email>noreply@blogger.com</email><gd:image rel='http://schemas.google.com/g/2005#thumbnail' width='31' height='24' src='//blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh3k6rMIRfexg9TADsrG6BHfOYj_ixEjrC_Du1lLyeojoyNbILGXXwLLFkWtTUfVAi5tRnMWL7kMDnvMc5aHPPLeJAEm7R38pQPH8LVO3A9_dv0eYTkUxvfsJ5yKmY-RQ/s220/AQa2.PNG'/></author><thr:total>0</thr:total></entry></feed>