[R] ANCOVA with random factor

[David Semmens]

I would like to know if there is a way of directly calculating the 
F-ratio of a random effect using the "aov" function. I have 2 factors 
in my model, "population" which is random and "length" which is the 
length of female fish within each population. The dependent variable is 
"diam" which is the average diameter of eggs produced by each female. 
At present I set up the model like this:

 >model <- aov(diam~population*length, data=data)
 >anova(model)

then using the output:

 >ratio=MS-population/MS-interaction
 >1-pf(ratio, df-population,df-interaction)

Is there a way of getting r to calculate the correct F-ratio for a 
random factor and present it in the original output?

Thanks,
David Semmens.