code: def outputpairedstats(fname,writemode,name1,n1,m1,se1,min1,max1,name2,n2,m2,se2,min2,max2,statname,stat,prob):
      suffix = ''
      try:
          x = prob.shape
          prob = prob[0]
      except:
          pass
      if  prob < 0.001:  suffix = '  ***'
      elif prob < 0.01:  suffix = '  **'
      elif prob < 0.05:  suffix = '  *'
      title = [['Name','N','Mean','SD','Min','Max']]
      lofl = title+[[name1,n1,round(m1,3),round(math.sqrt(se1),3),min1,max1],
                    [name2,n2,round(m2,3),round(math.sqrt(se2),3),min2,max2]]
      if type(fname)!=StringType or len(fname)==0:
          print()
          print(statname)
          print()
          pstat.printcc(lofl)
          print()
          try:
              if stat.shape == ():
                  stat = stat[0]
              if prob.shape == ():
                  prob = prob[0]
          except:
              pass
          print('Test statistic = ',round(stat,3),'   p = ',round(prob,3),suffix)
          print()
      else:
          file = open(fname,writemode)
          file.write('\n'+statname+'\n\n')
          file.close()
          writecc(lofl,fname,'a')
          file = open(fname,'a')
          try:
              if stat.shape == ():
                  stat = stat[0]
              if prob.shape == ():
                  prob = prob[0]
          except:
              pass
          file.write(pstat.list2string(['\nTest statistic = ',round(stat,4),'   p = ',round(prob,4),suffix,'\n\n']))
          file.close()
      return None