重复实验的方差分析与拟合

Quanc

有一个有关重复实验设计的问题向各位大神请教一下
数据包括一个因素X（四水平），观测值Y，和14个被试S
每个被试都参加了这个因素的所有水平的试验，也就是有4*14=56条记录

现在对该因素进行方差分析
两种方法
一种是ez包里的ezAnova(dv=.(Y),wid=.(S),within=.(X))
一种是aov(Y~X+Error(S/X))
到此还没有问题，自由度均是F(3,39)
问题是，我想进一步进行多项式拟合，把X因素转化成数值，进行最高次为3的多项式拟合（本质是想做Trend Analysis）
采用 fit <- lme(Y~X+I(X^2)+I(X^3)),random=~1|S/X),
然后将结果方差分析看效果
anova(fit)
得出每一项都是numDF=1,denDF=39

那么问题就来了，文献中求趋势分析最后得出的自由度为F(1,13)，而我的自由度却是F(1,39)
我的拟合过程究竟在哪出问题了呢？
求大神指点～～感激不尽

4337222liyi

新手，请吃教

Quanc

没有大神来指教一下么？..还是要附上数据？