#pragma rtGlobals=1		// Use modern global access method.
//Jobic, H. & Theodorou, D.
//"Quasi-elastic neutron scattering and molecular dynamics simulation
// as complementary techniques for studying diffusion in zeolites",
// Microporous and Mesoporous Materials, 2007, 102, 21-50
 
//note that we are fitting FWHM, NOT HWHM, so there are factors of two difference.
Function chudleyelliot(w, y, x):fitfunc
	wave w, y, x
	//eqn 51
	//w[0] = tau
	//w[1] = d
 
	y = sin(x * w[1]) / (x * w[1])
	y = 1 - y
	y /= w[0] / 2
End
 
Function hallross(w, y, x):fitfunc
	wave w, y, x

Function/s parseDirectoryListLine(directoryline, isdirectory)
string directoryline
variable &isdirectory
 
//directoryline = " dr-xr-xr-x 2 1 1  45056 Jan 6 20:11 q02ee8dddds"
string a, b, c,d,e,f,g,h,i, j,k,l,m,n,o,p,q
//perm
string regex="([dwrx-]+)"
//u
regex+="(\\s)+"
regex+="([[:alnum:][:punct:]]+)"
//g
regex+="(\\s)+"
regex+="([[:alnum:][:punct:]]+)"
//o
regex+="(\\s)+"
regex+="([[:alnum:][:punct:]]+)"
 
//size
regex+="(\\s)+"
regex+="([[:digit:]]+)"
 
//month
regex+="(\\s)+"
regex+="([[:alpha:]]{3})"
 
//day
regex+="(\\s)+"

//*********************************************************************
// statsBHFDR
//
// Implements Benjamini and Hochberg False Discovery Rate. Tolerates
// more false positives but less false negatives. Least stringent
// correction.
//
// Parameters:
// w - input wave of p-values
// 
// Optional Parameters:
// ncol - column number of p-values if input wave w is a 2D matrix
//*********************************************************************
function statsBHFDR(w,[ncol])
 
	//-----------------------------------------------------------------
	// parameters

//*********************************************************************
// statsBSDC
//
// Implements Bonferroni Step-down correction. 
//
// Parameters:
// w - input wave of p-values
// 
// Optional Parameters:
// ncol - column number of p-values if input wave w is a 2D matrix
//*********************************************************************
function statsBSDC(w,[ncol])
 
	//-----------------------------------------------------------------
	// parameters
	//-----------------------------------------------------------------
	wave w;
	variable ncol;
 

//*********************************************************************
// statsBC
//
// Shawn Driscoll
// Salk Institute for Biological Studies
//
// Implements Bonferroni correction. With a Family-wise error rate of
// 0.05, the expected number of false positives will be 0.05. 
//
// Parameters:
// w - input wave of p-values
// 
// Optional Parameters:
// ncol - column number of p-values if input wave w is a 2D matrix
//*********************************************************************
function statsBC(w,[ncol])

Function Test(in)
       Wave in //in is an unsigned 32 bit integer wave
 
       variable numPoints = numpnts(in)
 
       Make/O/W/N=(numPoints) out1
       MultiThread out1 = (in & 0xFFFF0000) / 65536
       Make/O/W/N=(numPoints) out2
       MultiThread out2 = in & 0xFFFF
End

function tri(x, [v])
	variable x, v
	if (ParamIsDefault(v))
		v = 0.5