R help - Oct 2006 - ANOVA with 3 factors

Hi,

I am trying to perform ANOVA with 3 factors:
individuals (3), time(3) and treatment(2). I tried
something like,

fit <- lm(dat ~ ind*time*treat)
res <- anova(fit)

Is this the correct approach? Is there any R pacakge
to calculate P values using permutation for ANOVA with
3 factors? 

Thanks in advance.

Best regards,
Ezhil