KR can sometimes give unusual results, especially if you have an overparameterized model (in this case, you could get ddf=1, which is really a warning of a problem), or if you have a lot of 0s for the estimated variances and covariances (where you could get somethihg like the residual df). You have several terms in your model, and I don't know which tests are giving you the "surprising" ddf calculation. I can't really comment too much on the appropriateness of your model (I can't really study the experimental layout). It is possible that the output is appropriate given your design and parameter estimates.
You should try ddfm=kr(firstorder) to see if this has an effect on the results. I am guessing not (because of the UN structure), but try this. I think you will probably need a less complex var-cov structure. Start with type=cs and work up to more complicated versions, keeping track of AIC.
... View more