The algorithm prints out a incorrect p-value for 'aov'.

Sample code:

library(Devore7)

motor_stack<-stack(ex10.37)

motor_aov<-aov(values~ind, motor_stack)

summary(motor_aov)

correct p-value is '0.000187' whereas one calculated by R in solus is '0.0855'