# Setting up the data
dat1 <- list(
gauss_1 = sim_data$continuous1_data[c(1:20, 61:80, 121:140, 181:200), 1:75],
gauss_2 = sim_data$continuous2_data[c(1:20, 61:80, 121:140, 181:200), 1:75],
poisson_1 = sim_data$count_data[c(1:20, 61:80, 121:140, 181:200), 1:75])
dist <- c(
"gaussian", "gaussian",
"poisson")
fit_iClusterVB <- iClusterVB(
mydata = dat1,
dist = dist,
K = 4,
initial_method = "VarSelLCM",
VS_method = 1,
max_iter = 25
)
# We can set the colors, turn off scaling and set titles
chmap(fit_iClusterVB,
cols = c("red", "blue", "green", "purple"),
title = c("Gene Expression", "DNA Methylation", "Copy Number"),
scale = "none"
)
Run the code above in your browser using DataLab