I've started using lmPerm in order to perform regressions in R. The
equation I want to fit has the form:

out3 <- lmp(outcome ~ bin1 + bin2 + cont1 + cont2, perm="Exact")

Where "outcome" is a non-normally distributed continuous variable, and bin*
and cont are binary and continuous regressors (similarly, they are
non-normally distributed). Each variable has a length of approx. 110 cases.

Here are my questions:

This code works fine, but every time it runs in R, different p-values
appear for each regressor. Which p-value should be reported in my results?
I've tried repeating the test several times (in a loop) and getting an
estimate from that, but I'm not sure it works...

If some of the predictors are changed and (then) several models /
hypothesis are tested, should Bonferroni corrections be used in the same
way they are applied for ordinary regressions? Is lmp somehow robust to
multiple testing procedures?

Thank you in advance,

	[[alternative HTML version deleted]]

