[R] How to get Greenhouse-Geisser epsilons from anova?

[Skotara]

Dear all,

I apologize for my basic question.
I try to calculate an anova for repeated measurements with 3 factors 
(A,B,C) having 2, 2, and 7 levels.
or with an additional fourth between subjects factor D.
Everything works fine using
aov(val ~ A*B*C  + Error(subject/ (A*B*C) ) )  or
aov(val ~ (D*A*B*C)  + Error(subject/(A*B*C)) + D )
val, A, B, C, D and subject are columns in a data.frame.

How can I get the estimated Greenhouse-Geisser and Huynh-Feldt epsilons?
I know Peter Dalgaard described it in R-News Vol. 7/2, October 2007. 
However, unfortunately I am not able to apply that using my data...

Furthermore, I am still confused of how SPSS calculates the epsilons 
since it is mentioned that perhaps there are any errors in SPSS??

I would be glad if anyone could help me!
I am looking forward to hearing from you!

Thank you!
Nils