#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 vartest2(const Range& irng, const int& tail, const double& alpha, double& stat, double& numdf, double& dendf, double& prob, double& lcl, double& ucl)
{
		
	/// *****		Terminate if error setting		***** ///
	if (!irng.IsValid())
	{
		XF_THROW(CER_NO_DATA);
	}
	
	/// Trim out missing values before calculation
	vector vData1, vData2;
	int nRet = get_2_data_in_var(irng, vData1, vData2);
	if (nRet != STATS_NO_ERROR)			
		XF_THROW(nRet);
	
	int nSize1 = vData1.GetSize();
	int nSize2 = vData2.GetSize();
	
	if (tail < 0 || tail > 2)
	{
		XF_THROW(CER_TAIL_SETTING);
	}
	
	if (!check_sig_level(alpha))
		XF_THROW(CER_INVALID_SIG_LEV);
		
	/// *****		F test computation		***** ///		
	HypotTestOptions opt;
	FVarResults FVres;
	opt.HypotValue = 1;
	opt.TailType=tail;
	nRet = ocmath_f_test(nSize1, vData1, nSize2, vData2, &opt, &FVres);
	if (nRet != STATS_NO_ERROR)
	{
		XF_THROW_EX(CER_UNKNOWN_ERROR_EX, nRet);
	}
	
	//Output missing values if there are
	//get the data range including missing vaue
	//dwRules = DRR_GET_MISSING | DRR_GET_DEPENDENT | DRR_NO_FACTORS;
	//irng.GetData(dwRules, 0, &dwPlotID, NULL, &vData2Tmp, &vData1Tmp);
	//int nMissing = vData2Tmp.GetSize() - vData2.GetSize();
	//if (nMissing > 0)
	//{
		//XF_TRACE(MISSING_IN_DATA_MSG, nMissing);
	//}
	
	if (stat != NULL)	stat= FVres.FValue;
	if (dendf != NULL)	dendf = FVres.FDOFDen;
	if (numdf != NULL)	numdf= FVres.FDOFNum;
	if (prob != NULL)	prob = FVres.FProb;
	if (lcl != NULL)	lcl = NANUM;
	if (ucl != NULL)	ucl = NANUM;
	
	//Confidence computation for a variance
	vector vLevels(1);
	vLevels[0] = 1 - alpha;
	vector vLCLs(1); 
	vector vUCLs(1);
	
	vector vSD(2);	
	vector<int> vN(2);
	ocmath_basic_summary_stats(nSize1, vData1, &vN[0], NULL, &vSD[0]);
	ocmath_basic_summary_stats(nSize2, vData2, &vN[1], NULL, &vSD[1]);
															
	nRet = ocmath_two_sample_confidence_levels_for_variance(vLevels.GetSize(), vLevels, vSD, vN, 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 vartest2_event1(TreeNode& trGetN, int nRow, int nEventID, DWORD& dwEnables, LPCSTR lpcszNodeName, WndContainer& DynaCntrlContainer, string& strAux, string& strErrMsg)
{
	if( xf_is_extended_event(nEventID) )
		return vartest2_event1_ex(trGetN, nRow, nEventID, dwEnables, lpcszNodeName, DynaCntrlContainer, strAux, strErrMsg);
	DECLARE_BUTTON_ENABLES
	return false;

}

int vartest2_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 vartest2_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 vartest2_event1_ex(TreeNode& trGetN, int nRow, int nEventID, DWORD& dwEnables, LPCSTR lpcszNodeName, WndContainer& DynaCntrlContainer, string& strAux, string& strErrMsg)
{
	return false;
}
