dnaclassificationdlg.cpp

DNA分类

// DNAClassificationDlg.cpp : implementation file
//

#include "stdafx.h"
#include "DNAClassification.h"
#include "DNAClassificationDlg.h"
#include "Cspreadsheet.h"
#include "PNNParamDlg.h"
#include "NearestParamDlg.h"

#include "PNN.h"
//#include "FishNearEx.h"
#include "bp.h"
#include "Nearest.h"


#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CDNAClassificationDlg dialog

CDNAClassificationDlg::CDNAClassificationDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CDNAClassificationDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CDNAClassificationDlg)
	m_nTotalSam = 0;
	m_dRate = 0.0;
	m_nErro = 0;
	m_dCost = 0.0;
	m_ClassifyMethod = -1;
	m_TotalTest = 0;
	m_dTestCost = 0.0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);

	m_bReaded = FALSE;
	memset(&m_TrainSamples,0,sizeof(m_TrainSamples));
	memset(&m_TestSamples,0,sizeof(m_TestSamples));

	m_pDNAClassify = NULL;

	m_hThread = NULL;

	m_dwTime = 0;

	m_nMode = 0;
}

void CDNAClassificationDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CDNAClassificationDlg)
	DDX_Control(pDX, IDC_PROGRESS_ANALYSE, m_ProgAnalyse);
	DDX_Control(pDX, IDC_CLASSIFY_GROUP, m_ClassifyGroup);
	DDX_Control(pDX, IDC_ETHOD_GROUP, m_MethodGroup);
	DDX_Control(pDX, IDC_ERRO_GROUP, m_ErroGroup);
	DDX_Control(pDX, IDC_LIST1, m_listCtrl);
	DDX_Text(pDX, IDC_EDIT_TOTAL, m_nTotalSam);
	DDX_Text(pDX, IDC_EDIT_RATE, m_dRate);
	DDX_Text(pDX, IDC_EDIT_ERRO, m_nErro);
	DDX_Text(pDX, IDC_EDIT_COST, m_dCost);
	DDX_Radio(pDX, IDC_RADIO_NEAR, m_ClassifyMethod);
	DDX_Text(pDX, IDC_EDIT_TOTALTEST, m_TotalTest);
	DDX_Text(pDX, IDC_EDIT_TESTCOST, m_dTestCost);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CDNAClassificationDlg, CDialog)
	//{{AFX_MSG_MAP(CDNAClassificationDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BUTTON_ANALYSE, OnButtonAnalyse)
	ON_WM_DESTROY()
	ON_BN_CLICKED(IDC_BUTTON_CLASSIFY, OnButtonClassify)
	ON_BN_CLICKED(IDC_BUTTON_STOPANALYSE, OnButtonStopanalyse)
	ON_BN_CLICKED(IDC_BUTTON_STOPCLASSIFY, OnButtonStopclassify)
	ON_BN_CLICKED(IDC_RADIO_PNN, OnRadioPnn)
	ON_BN_CLICKED(IDC_RADIO_NEAR, OnRadioNear)
	ON_BN_CLICKED(IDC_RADIO_BP, OnRadioBp)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CDNAClassificationDlg message handlers

BOOL CDNAClassificationDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// Add "About..." menu item to system menu.

	// IDM_ABOUTBOX must be in the system command range.
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		CString strAboutMenu;
		strAboutMenu.LoadString(IDS_ABOUTBOX);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// Set the icon for this dialog.  The framework does this automatically
	//  when the application's main window is not a dialog
	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	

	m_ErroGroup.SetXPGroupStyle(CXPGroupBox::XPGB_WINDOW)
		    .SetBackgroundColor(RGB(201, 201, 151), GetSysColor(COLOR_BTNFACE))//RGB(255, 255, 204))  
			.SetFontBold(TRUE);

	m_MethodGroup.SetXPGroupStyle(CXPGroupBox::XPGB_WINDOW)
		    .SetBackgroundColor(RGB(201, 201, 151),GetSysColor(COLOR_BTNFACE))  
			.SetFontBold(TRUE);

	m_ClassifyGroup.SetXPGroupStyle(CXPGroupBox::XPGB_WINDOW)
		    .SetBackgroundColor(RGB(201, 201, 151),GetSysColor(COLOR_BTNFACE))// RGB(255, 255, 204))  
			.SetFontBold(TRUE);

	m_ProgAnalyse.SetRange(0,100);

	GetDlgItem(IDC_BUTTON_STOPANALYSE)->EnableWindow(FALSE);
	GetDlgItem(IDC_BUTTON_STOPCLASSIFY)->EnableWindow(FALSE);

	CListCtrl *pList = (CListCtrl *)GetDlgItem(IDC_LIST1);
	pList->InsertColumn(0,"测试样本",LVCFMT_LEFT,105);
	pList->InsertColumn(1,"样本类别",LVCFMT_LEFT,105);

	pList = (CListCtrl *)GetDlgItem(IDC_LIST2);
	pList->InsertColumn(0,"训练样本",LVCFMT_LEFT,60);
	pList->InsertColumn(1,"样本类别",LVCFMT_LEFT,60);
	pList->InsertColumn(2,"判断类别",LVCFMT_LEFT,60);
	pList->InsertColumn(3,"错误/正确",LVCFMT_LEFT,70);

	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CDNAClassificationDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CDNAClassificationDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

// The system calls this to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CDNAClassificationDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CDNAClassificationDlg::OnButtonAnalyse() 
{
	UpdateData(TRUE);
	if(m_ClassifyMethod ==-1)
	{
		AfxMessageBox("请选择一种分类算法");
		return ;
	}

	m_nTotalSam = 0;
	m_nErro = 0;
	m_dRate = 0.0;
	m_dCost = 0.0;
	UpdateData(FALSE);
	GetDlgItem(IDC_BUTTON_ANALYSE)->EnableWindow(FALSE);
	GetDlgItem(IDC_BUTTON_STOPANALYSE)->EnableWindow(TRUE);
	GetDlgItem(IDC_BUTTON_CLASSIFY)->EnableWindow(FALSE);

	GetDlgItem(IDC_RADIO_NEAR)->EnableWindow(FALSE);
	GetDlgItem(IDC_RADIO_BP)->EnableWindow(FALSE);
	GetDlgItem(IDC_RADIO_PNN)->EnableWindow(FALSE);
	//GetDlgItem(IDC_RADIO_FISHER)->EnableWindow(FALSE);
	
	CListCtrl* pList = (CListCtrl*)GetDlgItem(IDC_LIST2);
	pList->DeleteAllItems();

	m_nMode = 0;
	CWinThread* pThread = AfxBeginThread(AnalyseProc,this);	
	m_hThread = pThread->m_hThread;
}

void CDNAClassificationDlg::OnButtonStopanalyse() 
{	
	if(m_hThread)
	{
		if(TerminateThread(m_hThread,0)==FALSE)
		{
			AfxMessageBox("停止失败！");
		}
		else 
		{
			GetDlgItem(IDC_BUTTON_ANALYSE)->EnableWindow(TRUE);
			GetDlgItem(IDC_BUTTON_STOPANALYSE)->EnableWindow(FALSE);
			GetDlgItem(IDC_BUTTON_CLASSIFY)->EnableWindow(TRUE);

			GetDlgItem(IDC_RADIO_NEAR)->EnableWindow(TRUE);
			GetDlgItem(IDC_RADIO_BP)->EnableWindow(TRUE);
			GetDlgItem(IDC_RADIO_PNN)->EnableWindow(TRUE);
		//	GetDlgItem(IDC_RADIO_FISHER)->EnableWindow(TRUE);		

			CListCtrl* pList = (CListCtrl*)GetDlgItem(IDC_LIST2);
			pList->DeleteAllItems();
			m_ProgAnalyse.SetPos(0);
		}
		m_hThread = NULL;
	}
}

/*****************************************************/
/*功  能：读取样本数据
/*参  数：无
/*返回值：
/*****************************************************/
BOOL CDNAClassificationDlg::ReadTrainSamples(CString filename,CString strSheet)
{
	CSpreadSheet sheet(filename,strSheet,FALSE);
	int c = sheet.GetTotalColumns()-1;//有一列是class
	int n = sheet.GetTotalRows()-1;   //有一行是标题
	
	if(n<=0||c<=0)
		return FALSE;

	ClearTrainSamples(m_TrainSamples);
	AllocTrain(m_TrainSamples,n,c);

	for(int i=2;i<=n+1;i++)
	{
		CStringArray ar;
		sheet.ReadRow(ar,i);
		for(int j =0 ;j<ar.GetSize()-1;j++)
		{
			CString str = ar.GetAt(j);
			m_TrainSamples.ppSamples[i-2][j] = atof(str);
		}
		CString str = ar.GetAt(ar.GetSize()-1);
		if(str == "EWS")
			m_TrainSamples.pClass[i-2] = 0;
		else if(str == "NB")
			m_TrainSamples.pClass[i-2] = 1;
		else if(str == "BL")
			m_TrainSamples.pClass[i-2] = 2;
		else m_TrainSamples.pClass[i-2] = 3;

	}

	m_bReaded = TRUE;
	return TRUE;
}

BOOL CDNAClassificationDlg::ReadTrainSamples()
{	
	CString filename = GetPath();
	switch(m_ClassifyMethod)
	{
	case 0:
		{
			filename+="DNA_DATA\\trainings_62c.xls";
		}
		break;
	case 1:
		{
			filename+="DNA_DATA\\trainingsv_22_normal.xls";				
		}
		break;
	case 2:
		{
			filename+="DNA_DATA\\training_60.xls";		
			
		}		
		break;
	case 3:
		{			
			filename+="DNA_DATA\\trainings_62c.xls";	
		}
		break;	
	}
	//读取样本
	return ReadTrainSamples(filename,"training");	
}

BOOL CDNAClassificationDlg::AllocTrain(TrainingSample& samples,int n,int d)
{
	m_TrainSamples.n = n;
	m_TrainSamples.c = d;
	m_TrainSamples.nExclude = 0;
	m_TrainSamples.pExludeSample = new int[n];
	m_TrainSamples.ppSamples = new double*[n];
	m_TrainSamples.pClass = new int[n];//类别
	memset(m_TrainSamples.pClass,0,sizeof(int)*n);

	for(int jj = 0;jj<n;jj++)
	{
		m_TrainSamples.pExludeSample[jj] = -1;
		m_TrainSamples.ppSamples[jj] = new double [d];
		memset(m_TrainSamples.ppSamples[jj],0,sizeof(double)*d);
	}
	return TRUE;
}

BOOL CDNAClassificationDlg::AllocTest(TestSample& samples,int n,int d)
{
	m_TestSamples.n = n;
	m_TestSamples.c = d;		
	m_TestSamples.ppTestSamples = new double*[n];
	m_TestSamples.pClass = new int[n];//类别
	memset(m_TestSamples.pClass,0,sizeof(int)*n);

	for(int jj = 0;jj<n;jj++)
	{		
		m_TestSamples.ppTestSamples[jj] = new double [d];
		memset(m_TestSamples.ppTestSamples[jj],0,sizeof(double)*d);		
	}
	return TRUE;
}

//读取测试样本
BOOL CDNAClassificationDlg::ReadTestSamples(CString filename,CString strSheet)
{
	CSpreadSheet sheet(filename,strSheet,FALSE);
	int c = sheet.GetTotalColumns()-1;//有一列是class
	int n = sheet.GetTotalRows()-1;   //有一行是标题
	
	if(n<=0||c<=0)
		return FALSE;

	ClearTestSamples(m_TestSamples);
	m_TestSamples.n = n;
	m_TestSamples.c = c;		
	m_TestSamples.ppTestSamples = new double*[n];
	m_TestSamples.pClass = new int[n];//类别
	memset(m_TestSamples.pClass,0,sizeof(int)*n);

	for(int jj = 0;jj<n;jj++)
	{		
		m_TestSamples.ppTestSamples[jj] = new double [c];
		memset(m_TestSamples.ppTestSamples[jj],0,sizeof(double)*c);
	}

	for(int i=2;i<=n+1;i++)
	{
		CStringArray ar;
		sheet.ReadRow(ar,i);
		for(int j =0 ;j<ar.GetSize()-1;j++)
		{
			CString str = ar.GetAt(j);
			m_TestSamples.ppTestSamples[i-2][j] = atof(str);
		}
	}	
	return TRUE;
}
//读取测试样本
BOOL CDNAClassificationDlg::ReadTestSamples()
{
	CString strTest = GetPath();
	switch(m_ClassifyMethod)
	{
	case 0:
		{			
			strTest += "DNA_DATA\\testsv_62.xls";	
		}
		break;
	case 1:
		{		
			strTest += "DNA_DATA\\testsv_22_normal.xls";		
		}
		break;
	case 2:
		{		
			strTest += "DNA_DATA\\testsv_22_normal.xls";
		}		
		break;