#include <origin.h>
#include <GetNBox.h>
#include <xfutils.h>
enum{
	tail_two,
	tail_upper,
	tail_lower,
};
#include <event_utils.h>

//put additional include files here
//put your own support static functions here


void ttestpair(const Range& irng, const double& mdiff, const int& tail, const double& alpha, double& stat, double& df, double& prob, double& lcl, double& ucl)
{
	
	/// *****		Terminate if error setting		***** /// 
	if (!irng.IsValid())
	{
		XF_THROW(CER_NO_DATA);
	}
	
	vector vData1, vData2;
	int nRet = get_2_data_in_var(irng, vData1, vData2, NULL, 0, true);
	if (nRet != STATS_NO_ERROR)			
		XF_THROW(nRet);
	
	if (tail < 0 || tail > 2)
	{
		XF_THROW(CER_TAIL_SETTING);
	}
	
	if (!check_sig_level(alpha))
		XF_THROW(CER_INVALID_SIG_LEV);
	
	///Echo 10/23/06 Different size of data will be trimmed automatically
	//Size of two datastes should be the same for paired t-test
	//if( vData1.GetSize() != vData2.GetSize())
	//{
		//XF_THROW(ERR_DATA_LENGTH);
	//}
	//
	if (0 == vData1.GetSize() || 0 == vData2.GetSize())
	{
		XF_THROW_EX(CER_DATA_TOO_FEW, 1);
	}
	
	/// *****		Paired T-test Computation		***** ///
	HypotTestOptions HTO;
	tTestMean2SampleResults tRes;
	
	HTO.HypotValue = mdiff;
	HTO.TailType=tail;
	nRet = ocmath_two_sample_t_test_pair(vData1, vData1.GetSize(), vData2, vData2.GetSize(), &HTO, &tRes);
	if (nRet != STATS_NO_ERROR)
	{
		XF_THROW_EX(CER_UNKNOWN_ERROR_EX, nRet);
	}
	
	if (stat != NULL)	stat = tRes.tValue;
	if (df != NULL)	df = tRes.DOF;
	if (prob != NULL)	prob = tRes.Prob;
	if (lcl != NULL)	lcl = NANUM;
	if (ucl != NULL)	ucl = NANUM;
	
	//confidence interval computation
	vector vLevels(1);
	vLevels[0] = 1 - alpha;
	vector vLCLs(1); 
	vector vUCLs(1);
	
	ocmath_two_sample_confidence_levels_for_mean_pair(vLevels, vLevels.GetSize(), &tRes, &HTO, vLCLs, vUCLs);	
	if (nRet != STATS_NO_ERROR)
	{
		XF_THROW_EX(CER_UNKNOWN_ERROR_EX, nRet);
	}
	
	if (lcl != NULL)	lcl = vLCLs[0];
	if (ucl != NULL)	ucl = vUCLs[0];
	
}

//nEventID < 0 on event, like GETNE_ON_INIT, GETNE_ON_OK etc, >0 on change from control
//return 0 to do nothing, return PEVENT_GETN_RET_TO_CLOSE to close dialog
//dwEnables = bOKEnable, bBTN1Enable, bBTN2Enable, bBTN3Enable, , can set to false to disable OK/Apply buttons
//lpcszNodeName = the trGetN TreeNode tag name that generates the event, or is  if GETNE_ON_OK etc
//strAux = optional string for XF Wizard Bar, like DLG_NEED_CHANGE_GRID etc, it will be invalid if from internal node value change event
//strErrMsg = allowing messages to be shown
int ttestpair_event1(TreeNode& trGetN, int nRow, int nEventID, DWORD& dwEnables, LPCSTR lpcszNodeName, WndContainer& DynaCntrlContainer, string& strAux, string& strErrMsg)
{
	if( xf_is_extended_event(nEventID) )
		return ttestpair_event1_ex(trGetN, nRow, nEventID, dwEnables, lpcszNodeName, DynaCntrlContainer, strAux, strErrMsg);
	DECLARE_BUTTON_ENABLES
	return false;

}

int ttestpair_make_tree(TreeNode& tr, LPCSTR lpcszVarName)
{
	return false;
}
//nGetNDialog = 0,1,2, corresponding to the XF GetN dropdown, 0 = dialog is not used
//you can change nRet to be one of XFEVT_PROCEED (defualt), XFEVT_ABORT, XFEVT_PROCEED_NO_DLG or XFEVT_PROCEED_WITH_DLG(when nGetNDialog = 0)
//nGetNDialog = -1 when Ok on GetN dialog to allow handler to preocess <auto> etc variables
//dwCntrl is combination of bits
//LTXF_FROM_GUI_MENU, LTXF_FROM_GUI_PROMPT, LTXF_FROM_AUTO_UPDATE, LTXF_FROM_FUNCTION, LTXF_THEME_USED, LTXF_CHANGE_PARAM
void ttestpair_before_execute(TreeNode& trGetN, int nGetNDialog, int& nRet, int dwCntrl)
{

}
//nEventID < 0 on event, like GETNE_ON_INIT, GETNE_ON_OK etc, >0 on change from control
//return 0 to do nothing, return PEVENT_GETN_RET_TO_CLOSE to close dialog
//dwEnables = bOKEnable, bBTN1Enable, bBTN2Enable, bBTN3Enable, , can set to false to disable OK/Apply buttons
//lpcszNodeName = the trGetN TreeNode tag name that generates the event, or is  if GETNE_ON_OK etc
//strAux = optional string for XF Wizard Bar, like DLG_NEED_CHANGE_GRID etc, it will be invalid if from internal node value change event
//strErrMsg = allowing messages to be shown
int ttestpair_event1_ex(TreeNode& trGetN, int nRow, int nEventID, DWORD& dwEnables, LPCSTR lpcszNodeName, WndContainer& DynaCntrlContainer, string& strAux, string& strErrMsg)
{
	return false;
}
