diff --git a/melhlprfns.cc b/melhlprfns.cc index 7d5a013f134347379716e127e5bbe71b02b6b461..9035a8f762041784106671af05a1b90976880410 100644 --- a/melhlprfns.cc +++ b/melhlprfns.cc @@ -908,8 +908,7 @@ namespace Melodic{ pf_fmf = f_fmf.Row(1); for(int ctr1=1;ctr1<=f_fmf.Ncols();ctr1++) - pf_fmf(1,ctr1) = 1.0-MISCMATHS::fdtr(design.Ncols(), - int(design.Nrows() -1 -design.Ncols()),f_fmf.Column(ctr1).AsScalar()); + pf_fmf(1,ctr1) = 1.0-MISCMATHS::fdtr(design.Ncols(),dof,f_fmf.Column(ctr1).AsScalar()); if(contrasts.Storage()>0 && contrasts.Ncols()==beta.Nrows()){ cbeta = contrasts*beta;