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

15.9
15.5
15.3
14.5
14.4
14.7
19.1
21.6
20.2
17.9
15.7
14.5
14.1
13.9
14.2
15.3
15.4
15.2
16.5
18.2
18.6
21
19.2
18.7
18.4
17.8
17.2
16.2
15.5
15.3
18.3
19.2
19
18.7
18.1
18.5
21.1
21
20.4
19.5
18.6
18.8
23.7
24.8
25
23.6
22.3
21.8
20.8
19.7
18.3
17.4
17
18.1
23.9
25.6
25.3
23.6
21.9
21.4
20.6
20.5
20.2
20.6
19.7
19.3
22.8
23.5
23.8
22.6
22
21.7
20.7
20.2
19.1
19.5
18.7
18.6
22.2
23.2
23.5
21.3
20
18.7
18.9
18.3
18.4
19.9
19.2
18.5
20.9
20.5
19.4
18.1
17
17
17.3
16.7
15.5
15.3
13.7
14.1
17.3
18.1
18.1

15.5
15.5
15.5
12.2
12.2
12.2
17.3
17.3
17.3
12.7
12.7
12.7
14.1
14.1
14.1
14.3
14.3
14.3
16.8
16.8
16.8
18.7
18.7
18.7
16.8
16.8
16.8
16
16
16
18
18
18
18
18
18
21.6
21.6
21.6
20.1
20.1
20.1
22.9
22.9
22.9
24
24
24
20.3
20.3
20.3
15.8
15.8
15.8
23.2
23.2
23.2
20.9
20.9
20.9
19.8
19.8
19.8
20.6
20.6
20.6
21.1
21.1
21.1
22.4
22.4
22.4
20.5
20.5
20.5
18.4
18.4
18.4
17.6
17.6
17.6
18.5
18.5
18.5
17.3
17.3
17.3
16.2
16.2
16.2
18.5
18.5
18.5
16.3
16.3
16.3
16.8
16.8
16.8
14.8
14.8
14.8
17.6
17.6
17.6

15.6
15.6
15.6
17.4
17.4
17.4
24
24
24
20.2
20.2
20.2
14.1
14.1
14.1
16.6
16.6
16.6
19.1
19.1
19.1
20.8
20.8
20.8
19
19
19
15.2
15.2
15.2
19.9
19.9
19.9
19
19
19
19.9
19.9
19.9
17.5
17.5
17.5
26.6
26.6
26.6
20.8
20.8
20.8
18.9
18.9
18.9
19.5
19.5
19.5
27.1
27.1
27.1
24.1
24.1
24.1
21.3
21.3
21.3
19.1
19.1
19.1
26.2
26.2
26.2
21.7
21.7
21.7
19.4
19.4
19.4
19.5
19.5
19.5
28.7
28.7
28.7
21.8
21.8
21.8
20
20
20
22.6
22.6
22.6
22.4
22.4
22.4
18.6
18.6
18.6
16.2
16.2
16.2
13.8
13.8
13.8
18
18
18

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))
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))
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()