https://abaj.ai/author/annette-j.-dobson/Recent content in Annette J. Dobson on Aayush Bajaj's Augmenting InfrastructureHugoen© 2026 Aayush BajajMon, 15 Jun 2026 06:00:15 +1000https://abaj.ai/words/library/books/an_introduction_to_generalized_linear_models/Sun, 14 Jun 2026 13:42:58 +1100https://abaj.ai/words/library/books/an_introduction_to_generalized_linear_models/<div class="highlight"><pre tabindex="0" class="chroma"><code class="language-R" data-lang="R"><span class="line"><span class="cl"><span class="nf">install.packages</span><span class="p">(</span><span class="s">&#34;glmnet&#34;</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="nf">install.packages</span><span class="p">(</span><span class="s">&#34;dobson&#34;</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="nf">library</span><span class="p">(</span><span class="n">glmnet</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="nf">library</span><span class="p">(</span><span class="n">dobson</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="nf">data</span><span class="p">(</span><span class="n">carbohydrate</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="n">y</span> <span class="o">=</span> <span class="n">carbohydrate</span><span class="o">$</span><span class="n">carbohydrate</span> </span></span><span class="line"><span class="cl"><span class="n">x</span> <span class="o">=</span> <span class="nf">as.matrix</span><span class="p">(</span><span class="n">carbohydrate[</span><span class="p">,</span><span class="nf">c</span><span class="p">(</span><span class="s">&#39;age&#39;</span><span class="p">,</span> <span class="s">&#39;weight&#39;</span><span class="p">,</span> <span class="s">&#39;protein&#39;</span><span class="p">)</span><span class="n">]</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="n">fit</span> <span class="o">=</span> <span class="nf">glmnet</span><span class="p">(</span><span class="n">x</span><span class="p">,</span> <span class="n">y</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="nf">plot</span><span class="p">(</span><span class="n">fit</span><span class="p">,</span> <span class="n">xvar</span><span class="o">=</span><span class="s">&#39;lambda&#39;</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="n">cvfit</span> <span class="o">=</span> <span class="nf">cv.glmnet</span><span class="p">(</span><span class="n">x</span><span class="p">,</span><span class="n">y</span><span class="p">)</span> </span></span><span class="line"><span class="cl"><span class="nf">coef</span><span class="p">(</span><span class="n">cvfit</span><span class="p">,</span> <span class="n">s</span> <span class="o">=</span> <span class="s">&#34;lambda.1se&#34;</span><span class="p">)</span> </span></span></code></pre></div><p>(setenv &ldquo;R_PROFILE_USER&rdquo; (expand-file-name &ldquo;~/.config/R/Rprofile&rdquo;))</p>