Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
/* MELODIC - Multivariate exploratory linear optimized decomposition into
independent components
melpca.h - PCA and whitening
Christian F. Beckmann, FMRIB Image Analysis Group
Copyright (C) 1999-2002 University of Oxford */
/* CCOPYRIGHT */
#ifndef __MELODICPCA_h
#define __MELODICPCA_h
#include "newimageall.h"
#include "log.h"
#include "meloptions.h"
#include "meldata.h"
#include "melodic.h"
//#include "melreport.h"
#include "newmatap.h"
#include "newmatio.h"
using namespace Utilities;
using namespace NEWIMAGE;
namespace Melodic{
class MelodicReport;
class MelodicPCA
{
public:
MelodicPCA(MelodicData &pmelodat, MelodicOptions &popts, Log &plogger, MelodicReport &preport):
melodat(pmelodat),
opts(popts),
logger(plogger),
report(preport)
{
}
void perf_pca(const Matrix &Data);
void perf_white(const Matrix &Data);
private:
MelodicData &melodat;
MelodicOptions &opts;
Log &logger;
MelodicReport &report;
int pcadim();
RowVector cumsum(const RowVector& Inp);
Matrix ppca_est(const RowVector& eigenvalues, const int N);
RowVector Feta(int n1,int n2);
};
}
#endif