• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

9
8
10
8
9
10
9
9
9
9
10
9
13
11
10
8
11
6
8
9
9
10
9
10
8
10
10
11
9
9
7
11
6
11
10
9
10
10
9
8
9
9
10
11
7
9
10
8
3
10
10
10
4
10
8
9
13
10
8
9
11
10
9
10
7
10
10
11
12
8
10
6
9
11
10
10
8
10
9
9
10
10
11
9
12
7
9
9
11
8
9
9
9
9
11
9
7
15
9
9
12
10
9
10
10
9
10
10
9
9
9
9
11
9
7
11
9
7
12
8
9
9
9
9
11
9
10
10
9
8
10
9
8
10
9
13
13
11
9
8
9
10
10
10
9
11
7
9
10
11
8
8
7
9
9
10
9
11
8
11
9
9
9
7
10
9
9

10
13
14
12
12
13
13
13
13
14
14
12
12
11
12
14
12
11
13
13
12
13
12
13
11
12
12
13
13
10
12
13
13
10
14
12
10
10
14
12
14
10
13
12
12
13
12
10
9
14
15
14
8
11
10
12
14
12
12
14
13
13
13
12
10
14
11
10
13
12
12
11
10
14
12
13
11
10
14
13
7
13
13
13
15
13
14
12
13
11
12
14
13
14
12
12
13
14
13
12
13
12
10
12
13
12
13
12
12
12
11
12
9
14
12
13
13
13
11
12
11
12
12
13
12
13
13
12
12
8
12
13
10
8
12
13
12
15
14
10
11
12
10
14
10
15
11
12
9
12
13
12
9
12
14
11
12
14
12
15
11
12
12
10
12
11
11

18
21
21
18
21
21
17
19
23
23
21
20
20
19
19
20
20
21
21
19
22
21
20
21
19
20
20
20
21
19
23
22
21
18
21
19
19
19
20
21
23
18
23
20
20
20
18
21
15
21
23
22
21
19
20
17
18
20
21
24
20
22
17
22
18
19
20
20
20
19
23
20
20
22
19
21
20
19
19
21
18
21
22
22
21
23
22
20
19
17
19
21
21
21
20
21
22
23
23
20
22
21
22
18
24
21
18
20
20
19
20
20
19
19
19
19
20
23
19
19
19
20
21
21
19
19
21
18
18
15
20
21
17
18
19
23
21
24
23
20
20
20
20
22
22
21
19
20
17
20
19
21
20
21
22
19
17
21
20
19
22
21
20
20
21
21
19

x <- array(x,dim=c(length(x),1))
colnames(x) <- par5
y <- array(y,dim=c(length(y),1))
colnames(y) <- par6
z <- array(z,dim=c(length(z),1))
colnames(z) <- par7
d <- data.frame(cbind(z,y,x))
colnames(d) <- list(par7,par6,par5)
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
if (par1>500) par1 <- 500
if (par2>500) par2 <- 500
if (par1<10) par1 <- 10
if (par2<10) par2 <- 10
library(GenKern)
library(lattice)
panel.hist <- function(x, ...)
{
usr <- par('usr'); on.exit(par(usr))
par(usr = c(usr[1:2], 0, 1.5) )
h <- hist(x, plot = FALSE)
breaks <- h$breaks; nB <- length(breaks)
y <- h$counts; y <- y/max(y)
rect(breaks[-nB], 0, breaks[-1], y, col='black', ...)
}
bitmap(file='cloud1.png')
cloud(z~x*y, screen = list(x=-45, y=45, z=35),xlab=par5,ylab=par6,zlab=par7)
dev.off()
bitmap(file='cloud2.png')
cloud(z~x*y, screen = list(x=35, y=45, z=25),xlab=par5,ylab=par6,zlab=par7)
dev.off()
bitmap(file='cloud3.png')
cloud(z~x*y, screen = list(x=35, y=-25, z=90),xlab=par5,ylab=par6,zlab=par7)
dev.off()
bitmap(file='pairs.png')
pairs(d,diag.panel=panel.hist)
dev.off()
x <- as.vector(x)
y <- as.vector(y)
z <- as.vector(z)
bitmap(file='bidensity1.png')
op <- KernSur(x,y, xgridsize=par1, ygridsize=par2, correlation=cor(x,y), xbandwidth=dpik(x), ybandwidth=dpik(y))
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (x,y)',xlab=par5,ylab=par6)
if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par4=='Y') points(x,y)
(r<-lm(y ~ x))
abline(r)
box()
dev.off()
bitmap(file='bidensity2.png')
op <- KernSur(y,z, xgridsize=par1, ygridsize=par2, correlation=cor(y,z), xbandwidth=dpik(y), ybandwidth=dpik(z))
print(op)
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (y,z)',xlab=par6,ylab=par7)
if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par4=='Y') points(y,z)
(r<-lm(z ~ y))
abline(r)
box()
dev.off()
bitmap(file='bidensity3.png')
op <- KernSur(x,z, xgridsize=par1, ygridsize=par2, correlation=cor(x,z), xbandwidth=dpik(x), ybandwidth=dpik(z))
print(op)
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (x,z)',xlab=par5,ylab=par7)
if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par4=='Y') points(x,z)
(r<-lm(z ~ x))
abline(r)
box()
dev.off()