Thanks. I tried to, but unfortunately I don’t know how to do this (I got bogged down with ‘branch comparisons’). I take your point about the name … maybe call it ‘fstat’ because it is a standard diagnostic, or else make it an allowable method for ftest)