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
|
/* quick.cc
Mark Woolrich, FMRIB Image Analysis Group
Copyright (C) 1999-2000 University of Oxford */
/* CCOPYRIGHT */
#include <iostream>
#include <fstream>
#include <iomanip>
#include <sstream>
#define WANT_STREAM
#define WANT_MATH
#include "miscmaths.h"
#include "t2z.h"
//#include "libprob.h"
using namespace MISCMATHS;
int main(int argc, char *argv[])
{
try{
Matrix X = read_vest("/usr/people/woolrich/matlab/vbbabe/data/design2.mat").t();
ColumnVector Y = read_vest("/usr/people/woolrich/matlab/vbbabe/data/sdf2.mat").t();
ColumnVector m_B;
SymmetricMatrix ilambda_B;
glm_vb(X, Y, m_B, ilambda_B, 30);
write_ascii_matrix(m_B,"/usr/people/woolrich/matlab/vbbabe/data/m_B");
}
catch(Exception p_excp)
{
cerr << p_excp.what() << endl;
}
return 0;
}
|
