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

28.9
23.7
21.2
26.1
22.2
16.8
19.7
21.1
22.9
24.7
20.8
23
25.5
23.5
23.2
25.8
27
23.4
23.9
18.4
25.4
22.6
29.5
24.9
19.7
23.2
21.1
21.6
21.9
21.5
22.1
24.2
24.6
25.3
24.9
21.2
24.3
19
25.1
23.8
27.2
25.8
20.4
18.7
23.1
22.2
27.6
23.2
20.8

30.6
32.4
36.4
27.6
26
33.4
23.2
35.4
33.8
29.5
30.4
31.3
28.3
34.4
30
36.8
37.1
31.2
26.1
27.8
29.3
31.3
35.1
33.4
41.5
36.6
44
43.4
24.1
37.8
23
30.2
41.7
33.5
27
30.6
29.6
28
37
39
30.2
37.4
19.5
31.4
29.3
27.9
31.3
38.2
36.4

1.72
3.83
3.29
1.34
4.53
3.66
4.62
5.14
4.09
2.38
3.06
3.1
2
1.33
1.34
1.39
2.35
3.57
2.95
4.89
2.47
2.66
1.06
2.36
3.06
1.71
5.75
6.26
4.46
2.25
3.61
2.43
1.56
2.01
1.36
4.27
2.09
4.36
2.06
1.5
1.6
2.13
1.37
4.75
3.09
4.1
0.88
2.63
2

par7 <- 'Variable Z'
par6 <- 'Variable Y'
par5 <- 'Variable X'
par4 <- 'Y'
par3 <- 'Y'
par2 <- '50'
par1 <- '50'
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()