wPath = sys.argv[4]
perms = int(sys.argv[5])
damp = float(sys.argv[6])
#outPath = sys.argv[7]
t = np.loadtxt(tPath)
X = np.loadtxt(xPath)
if wPath == 'NA':
wV = []
else:
wV = np.loadtxt(wPath)
x = X[index,:]
p,w,g = estSigGStat(x,t,wV = wV, nPerm=perms,permEst=False,damp=damp)
if np.isnan(p) or p==0:
p,w,g = estSigGStat(x,t,wV = wV, nPerm=perms*100,permEst=False,damp=damp)
#f = open(outPath,'a')
#f.write('>>\t'+str(index)+'\t'+str(p)+'\t'+str(w)+'\n')
#f.flush()
#os.fsync(f)
#f.close()
if ~np.isnan(p): out = '>>\t'+str(index)+'\t%05.4E\t%05.4E\t%05.4E\n' % (p,w,g)
else: out = '>>\t'+str(index)+'\tnan\tnan\tnan\n'
print(out)
wV = (2*np.pi/n)*np.arange(1,int(n/2))
for i in range(100):
x[i,:] = e*np.random.randn(n)
b = np.sqrt(2)
for i in range(100,200):
w = np.random.rand(1)*(np.max(wV)-np.min(wV))+np.min(wV)
s = np.random.rand(1)*2.0*n - n
x[i,:] = b*np.cos(w*t+s) + e*np.random.randn(n)
#g = np.zeros(200)
#for i in range(200):
# g[i],wTmp = fourierSeries.gStatMax(x[i,:],t,wV)
#
#np.savetxt('g.dat',g)
p = np.zeros(200)
for i in range(200):
p[i],wTmp = fourierSeries.estSigGStat(x[i,:],t,wV,300,False)
np.savetxt('p.dat',p)
tmp = -1*np.log10(p)
m = np.median(tmp)
plt.plot(tmp,'o');plt.plot([0,len(p)],[m,m]);plt.show()
