<?xml version="1.0" encoding="UTF-8"?>
<rss xmlns:content="http://purl.org/rss/1.0/modules/content/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:taxo="http://purl.org/rss/1.0/modules/taxonomy/" version="2.0">
  <channel>
    <title>topic Re: interpreting variability explained by exposure in multilevel model with binary outcome in Statistical Procedures</title>
    <link>https://communities.sas.com/t5/Statistical-Procedures/interpreting-variability-explained-by-exposure-in-multilevel/m-p/533931#M26890</link>
    <description>&lt;P&gt;Under a linear regression, the parameter estimates are chosen in such a way to minimize the residual variance.&amp;nbsp; That minimization gives you the concept of an R**2 statistic, and you can talk about the percentage of variance that a model explains.&lt;/P&gt;
&lt;P&gt;&amp;nbsp;&lt;/P&gt;
&lt;P&gt;With a logistic regression, you are calculating the parameter estimates using maximum likelihood.&amp;nbsp; You are not minimizing a residual variance, hence you don't have a traditional R**2 statistic to fall back on.&amp;nbsp; Several psuedo-R**2 measures have been proposed in the statistical literature, but I don't know that any universal acceptance has been reached on one measure over another.&amp;nbsp;&amp;nbsp;&lt;/P&gt;</description>
    <pubDate>Fri, 08 Feb 2019 14:36:50 GMT</pubDate>
    <dc:creator>StatsMan</dc:creator>
    <dc:date>2019-02-08T14:36:50Z</dc:date>
    <item>
      <title>interpreting variability explained by exposure in multilevel model with binary outcome</title>
      <link>https://communities.sas.com/t5/Statistical-Procedures/interpreting-variability-explained-by-exposure-in-multilevel/m-p/532144#M26833</link>
      <description>&lt;P&gt;Hi Sas community! I am analyzing whether a census-tract level variable (exposure), which is ordinal and has 5 levels, is associated with a binary outcome (outcome), after adjusting for race (racecat), age (agegrp) and clustering by census tract (tract10). My reference category for exposure is 5 and I found significant associations between each lower level of the variable and the reference category of 5. This is the code I used:&lt;/P&gt;&lt;P&gt;&amp;nbsp;&lt;/P&gt;&lt;P&gt;proc glimmix data=dat;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&lt;/P&gt;&lt;P&gt;class tract10 exposure racecat agegrp;&lt;/P&gt;&lt;P&gt;model outcome (event='1') = exposure racecat agegrp/ dist=binary link=logit ddfm=bw solution oddsratio cl;&lt;/P&gt;&lt;P&gt;random intercept / subject=tract10 ;&amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; &amp;nbsp; run;&amp;nbsp;&lt;/P&gt;&lt;P&gt;&amp;nbsp;&lt;/P&gt;&lt;P&gt;I have been asked to look into what % variability can be explained by the exposure variable or by the model, but because I can’t use R^2 for this kind of model, I am at a loss. Can anyone provide some guidance on how to answer this question??&lt;/P&gt;</description>
      <pubDate>Fri, 01 Feb 2019 18:26:32 GMT</pubDate>
      <guid>https://communities.sas.com/t5/Statistical-Procedures/interpreting-variability-explained-by-exposure-in-multilevel/m-p/532144#M26833</guid>
      <dc:creator>kb4</dc:creator>
      <dc:date>2019-02-01T18:26:32Z</dc:date>
    </item>
    <item>
      <title>Re: interpreting variability explained by exposure in multilevel model with binary outcome</title>
      <link>https://communities.sas.com/t5/Statistical-Procedures/interpreting-variability-explained-by-exposure-in-multilevel/m-p/533931#M26890</link>
      <description>&lt;P&gt;Under a linear regression, the parameter estimates are chosen in such a way to minimize the residual variance.&amp;nbsp; That minimization gives you the concept of an R**2 statistic, and you can talk about the percentage of variance that a model explains.&lt;/P&gt;
&lt;P&gt;&amp;nbsp;&lt;/P&gt;
&lt;P&gt;With a logistic regression, you are calculating the parameter estimates using maximum likelihood.&amp;nbsp; You are not minimizing a residual variance, hence you don't have a traditional R**2 statistic to fall back on.&amp;nbsp; Several psuedo-R**2 measures have been proposed in the statistical literature, but I don't know that any universal acceptance has been reached on one measure over another.&amp;nbsp;&amp;nbsp;&lt;/P&gt;</description>
      <pubDate>Fri, 08 Feb 2019 14:36:50 GMT</pubDate>
      <guid>https://communities.sas.com/t5/Statistical-Procedures/interpreting-variability-explained-by-exposure-in-multilevel/m-p/533931#M26890</guid>
      <dc:creator>StatsMan</dc:creator>
      <dc:date>2019-02-08T14:36:50Z</dc:date>
    </item>
  </channel>
</rss>

