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
|
#
# Copyright (C) 2004-2016 Friedrich Leisch and Bettina Gruen
# $Id: utils.R 5079 2016-01-31 12:21:12Z gruen $
#
list2object = function(from, to){
n = names(from)
s = slotNames(to)
p = pmatch(n, s)
if(any(is.na(p)))
stop(paste("\nInvalid slot name(s) for class",
to, ":", paste(n[is.na(p)], collapse=" ")))
names(from) = s[p]
do.call("new", c(from, Class=to))
}
printIter = function(iter, logLik, label="Log-likelihood")
cat(formatC(iter, width=4),
label, ":", formatC(logLik, width=12, format="f"),"\n")
## library(colorspace)
## dput(x[c(1,3,5,7,2,4,6,8)])
## x = hcl(seq(0, 360*7/8, length.out = 8), c=30)
LightColors <- c("#F9C3CD", "#D0D4A8", "#9DDDD5", "#D1CCF5",
"#EDCAB2", "#AFDCB8", "#ACD7ED", "#EFC4E8")
## x = hcl(seq(0, 360*7/8, length.out = 8), c=100, l=65)
FullColors <- c("#FF648A", "#96A100", "#00BCA3", "#9885FF",
"#DC8400", "#00B430", "#00AEEF", "#F45BE1")
###**********************************************************
## similar defaults to silhouette plots in flexclust
unipolarCols <- function(n, hue=0, chr=50, lum = c(55, 90))
{
lum <- seq(lum[1], lum[2], length=n)
hcl(hue, chr, lum)
}
bipolarCols <- function(n, hue=c(10, 130), ...)
{
if(n%%2){ # n odd
n2 <- (n-1)/2
c1 <- unipolarCols(n2, hue[1])
c2 <- rev(unipolarCols(n2, hue[2]))
return(c(c1, "white", c2))
}
else{ # n even
n2 <- n/2
c1 <- unipolarCols(n2, hue[1])
c2 <- rev(unipolarCols(n2, hue[2]))
return(c(c1, c2))
}
}
###**********************************************************
|
