-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathSIMULATIONSRandIndex.R
More file actions
61 lines (39 loc) · 1.11 KB
/
SIMULATIONSRandIndex.R
File metadata and controls
61 lines (39 loc) · 1.11 KB
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
59
60
### This simulation Checks if the model can predict well with 10 Dimensions
### Also CHECK THE TIME REQUIRED FOR THE MODEL
### Remove The Past
rm(list = ls())
RIground <- c(0)
CIground <- c(0)
for ( v in 1:20){
## Number of points
N.test = 100
N.train = 100
## Number of Clusters
F = 2
## Distribution of the points within three clusters
p.dist = c(0.5,0.5)
## Total Number of features D
D = 20
## Total Percentage of irrelevant feature
prob.noise.feature = 0.5
## Overlap between Cluster of molecular Data of the relevant features
prob.overlap = 0.10
#################################### SIMULATED DATA PROPERTIES ####################################################
############################# PARAMETERS for GIBB's SAMPLING ####
iter = 100
iter.burnin = 100
iter.thin = 5
k = 2
F = k
## Initialize the Training Data
source('simulateDPMM.R')
simulateDPMM()
######### OR Don't preprocess use this ########3
Y <- Y.dat
Y.new <- Y.new.dat
########### Check Prediction Ground Truth
source('predictionGroundTruth.R')
predictionGroundTruth()
RIground[v] <- predRandIndex
CIground[v] <- predCIndex
}