First linear model


\begin{Schunk}
\begin{Sinput}
> vs <- mincLm(jacobians ~ Genotype, gf)
\end{Sinp...
...lm.mnc
Sizes: 141 274 210
Range: 19.596058 -18.206707
\end{Soutput}\end{Schunk}

mincLm is the command to run linear models in RMINC. Its basic use is to provide a formula (same syntax as the R lm command) with the left side containing the filenames, the right side the variables to be regressed. The output of mincLm depends on the formula. There will always be a column of F-statistics, representing the significance of the entire model. Then there is one column for each of the terms in the mode. The above linear model, relating the Jacobian determinant to Genotype, will thus have three columns:

F-statistic
representing the significance of the entire model.
(Intercept)
the intercept term - this term is rarely interesting, as it tests for whether the intercept is 0. There's no reason to believe it should be in most cases, so this value will be highly significant but meaningless.
Genotype+
the term testing whether the ``+'' level of the Genotype factor is significant. In this case this term is the most interesting and therefore the one written to file.

The output is placed into a variable that can be written to file in the same way as described in the descriptive statistics section.

Jason Lerch 2008-02-17