A simple panel for statistical comparison of two 1D waves

 This does the statistical comparisons on two 1D waves with the outputs to a notebook. It first computes the mean, sem, and number of points of the waves selected, then do the normality test, variance test then the statistical test using t-test (pair and unpaired), K-S test, or Wilcoxon Signed Rank Test based on the normality and variance tests. Finally it will calculate the effect size using Cohen's d. It is wrote in a rush, and the waves selection panel is just a copy-and-paste of the waveselector example procedures. But my colleagues like it and suggest to share this, I guess their main reason is that now all the statistical analyses go to one notebook, and they no longer need to frantically click the close button to close.



#include <WaveSelectorWidget>
Function Make2WaveSelectorPanel()


	String panelName = "WaveSelectorExample"
	// figure out what to show in the Wave Selector, and make an appropriate name for the panel

			panelName+="Waves"


	if (WinType(panelName) == 7)
		// if the panel already exists, show it
		DoWindow/F $panelName
	else
		// doesn't exist, make it
		NewPanel/N=$panelName/W=(181,179,471,510) as "Wave Selector Example"
		// list box control doesn't have any attributes set on it
		ListBox ExampleWaveSelectorList,pos={9,13},size={273,241}
		// This function does all the work of making the listbox control into a
		// Wave Selector widget. Note the optional parameter that says what type of objects to
		// display in the list.
		MakeListIntoWaveSelector(panelName, "ExampleWaveSelectorList", content = WMWS_Waves)


			PopupMenu sortKind, pos={9,270},title="Sort Waves By"
			MakePopupIntoWaveSelectorSort(panelName, "ExampleWaveSelectorList", "sortKind")

		// This is an extra bonus- you can create your own function to be notified of certain events,
		// such as a change in the selection in the list.
		WS_SetNotificationProc(panelName, "ExampleWaveSelectorList", "ExtExampleNotification", isExtendedProc=1)

		// To support this demo, provide a button that displays this code
		Button Gr2Comp,pos={9,300},size={110,20},proc=gr2stat_0,title="2 Group"
		Button GrPairComp,pos={130,300},size={110,20},proc=gr2stat_1,title="Pair"
	endif
End

Function gr2stat_0(ba) : ButtonControl
	STRUCT WMButtonAction &ba

	switch( ba.eventCode )
		case 2: // mouse up
			comp2(1)// click code here
			break
		case -1: // control being killed
			break
	endswitch

	return 0

End

Function gr2stat_1(ba) : ButtonControl
	STRUCT WMButtonAction &ba

	switch( ba.eventCode )
		case 2: // mouse up
			comp2(3)// click code here
			break
		case -1: // control being killed
			break
	endswitch

	return 0

End

Function comp2(p2)
	variable p2
	string w_sed = WS_SelectedObjectsList("Waveselectorexamplewaves","ExampleWaveSelectorList")
	string x1_n = StringFromList(0, w_sed)
	string x2_n = StringFromlist(1, w_sed)
	wave x1 = $x1_n
	wave x2 = $x2_n
	sn_nb();Append2SN_TS()
	Append2SN("*" + x1_n + ":\r	mean = " + num2str(mean(x1)) + " sem = " + num2str(sqrt(variance(x1)/numpnts(x1))) + " n = " + num2str(numpnts(x1)))
	Append2SN("*" + x2_n + ":\r	mean = " + num2str(mean(x2)) + " sem = " + num2str(sqrt(variance(x2)/numpnts(x2))) + " n = " + num2str(numpnts(x2)))

	if (p2 == 1)
		Append2SN("*** Compare 2 waves ***")
		cmp2(1, x1, x2)
	elseif (p2 == 3)
		Append2SN("*** Compare paired waves ***")
		cmp2(3, x1, x2)
	ENDIF
END


Function Append2SN_TS() //Time Stamp
	Variable stampDateTime = 1 // nonzero if we want to include stamp
  	Variable tnow
   String stamp
	Notebook StatNote selection={endOfFile, endOfFile}
	tnow = datetime
	stamp = Secs2Date(tnow,0) + ", " + Secs2Time(tnow,0)
	Notebook StatNote text="\r" + "****" + stamp + "****\r"
End

Function Append2SN(str) //Add Log
	String str
	Notebook StatNote selection={endOfFile, endOfFile}
	Notebook StatNote text= str+"\r"
End

Function sn_nb()
	if (WinType("StatNote") == 5)
		// if StatNote already exists, show it
		DoWindow/F StatNote
	else
		newnotebook /n=StatNote /F=1 as "Stat Results Note"
	ENDIF
END

Function effect_size(x1, x2)
	wave x1, x2
	variable cohen_s_d = abs(mean(x1) - mean (x2))/sqrt((variance(x1) + variance(x2))/2)
	Append2SN("*Effect Size:\r	Cohen's d = " + num2str(cohen_s_d) + " (d = 0.2, small; 0.5, medium; 0.8, large)")
END

Function cmp2(p2, x1, x2)
	variable p2
	wave x1, x2
	statsjbtest/Q x1
	wave W_JBResults
	variable nd_x1
		if (W_JBResults[3] < W_JBResults[5])
			Append2SN("x1 is Normal distribution")
			nd_x1 = 1
		else
			Append2SN("x1 is not normal")
			nd_x1 = 0
		ENDIF
	statsjbtest/Q x2
	wave W_JBResults
	variable nd_x2
		if (W_JBResults[3] < W_JBResults[5])
			Append2SN("x1 is Normal distribution")
			nd_x2 = 1
		else
			Append2SN("x1 is not normal")
			nd_x1 = 0
		ENDIF
	statsvariancestest/Q x1, x2
	wave W_statsVariancesTest
	variable nd_xs
		if (W_statsVariancesTest[2] < W_statsVariancesTest[3])
			Append2SN("x1, x2 variance same")
			nd_xs = 1
		else
			Append2SN("x1, x2 variance NOT same")
			nd_xs = 0
		ENDIF
	variable p2w_con = (nd_x1==1 && nd_x2==1 && nd_xs==1)
	variable p2_final = p2 + p2w_con
	variable hs_rn = (CaptureHistoryStart())
	switch(p2_final)
		case 2:
			Append2SN("USE T-TEST\r" + capturehistory(hs_rn, 0) )
			statsttest x1, x2
			break
		case 1:
			Append2SN("USE K-S TEST\r" + capturehistory(hs_rn, 0))
			statskstest x1, x2
			break
		case 4:
			Append2SN("USE paired T-TEST\r" + capturehistory(hs_rn, 0))
			statsttest /PAIR x1, x2
			break
		case 3:
			Append2SN("USE Wilcoxon Signed Rank Test\r" +capturehistory(hs_rn, 0))
			StatsWilcoxonRankTest /WSRT  x1, x2
			break
	endswitch
	Append2SN(CaptureHistory(hs_rn, 1))
	effect_size(x1, x2)
END

Menu "Macros"
	"2 waves analysis", Make2WaveSelectorPanel()
END

Forum

Support

Gallery

Igor Pro 10

Learn More

Igor XOP Toolkit

Learn More

Igor NIDAQ Tools MX

Learn More