02.09.2014 Views

multivariate poisson hidden markov models for analysis of spatial ...

multivariate poisson hidden markov models for analysis of spatial ...

multivariate poisson hidden markov models for analysis of spatial ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

The loglinear model is one special case <strong>of</strong> Generalized Linear Models (GLM) <strong>for</strong><br />

Poisson distributed data (Agresti, 2002 and Brijs, 2002). Further, loglinear <strong>analysis</strong> can<br />

be considered as an extension <strong>of</strong> the two-way contingency table to where the<br />

conditional relationship between two or more discrete categorical variables is analyzed<br />

by taking the natural logarithm <strong>of</strong> the cell frequencies within the contingency table.<br />

Loglinear <strong>models</strong> are generally used to summarize multi-way contingency tables that<br />

involve three or more variables. There<strong>for</strong>e, loglinear <strong>models</strong> are very useful to evaluate<br />

the association between variables. PROC CATMOD procedure in SAS s<strong>of</strong>tware<br />

(SAS/STAT, 2003) can be used to fit the <strong>models</strong>.<br />

The fundamental strategy in loglinear <strong>analysis</strong> involves fitting <strong>models</strong> to the observed<br />

frequencies in the cross-tabulation <strong>of</strong> categorical variables. The <strong>models</strong> can then be<br />

represented by a set <strong>of</strong> expected frequencies that may or may not look like the observed<br />

frequencies. Different <strong>models</strong> can be described in terms <strong>of</strong> marginal <strong>models</strong> that they fit<br />

and in terms <strong>of</strong> the constraints they impose on the associations that are present in the<br />

data. Using expected frequencies, different <strong>models</strong> can be fitted and compared that are<br />

hierarchical to one another. The idea <strong>of</strong> modelling is then to choose a preferred model,<br />

which is the most suitable model that fits the data. The choice <strong>of</strong> the preferred model is<br />

based on a <strong>for</strong>mal comparison <strong>of</strong> goodness-<strong>of</strong>-fit statistics (likelihood ratio test)<br />

associated with <strong>models</strong> that are related hierarchically (i.e. <strong>models</strong> containing higher<br />

order terms also implicitly include all lower order terms).<br />

104

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!