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

2,86
2,55
2,27
2,26
2,57
3,07
2,76
2,51
2,87
3,14
3,11
3,16
2,47
2,57
2,89
2,63
2,38
1,69
1,96
2,19
1,87
1,6
1,63
1,22
1,21
1,49
1,64
1,66
1,77
1,82
1,78
1,28
1,29
1,37
1,12
1,51
2,24
2,94
3,09
3,46
3,64
4,39
4,15
5,21
5,8
5,91
5,39
5,46
4,72
3,14
2,63
2,32
1,93
0,62
0,6
-0,37
-1,1
-1,68
-0,78
-1,19

8,9
8,6
8,3
8,3
8,3
8,4
8,5
8,4
8,6
8,5
8,5
8,5
8,5
8,5
8,5
8,5
8,5
8,5
8,5
8,5
8,6
8,4
8,1
8
8
8
8
7,9
7,8
7,8
7,9
8,1
8
7,6
7,3
7
6,8
7
7,1
7,2
7,1
6,9
6,7
6,7
6,6
6,9
7,3
7,5
7,3
7,1
6,9
7,1
7,5
7,7
7,8
7,8
7,7
7,7
7,8
7,9

100,64
100,63
100,43
100,8
101,33
101,88
101,85
102,04
102,22
102,63
102,65
102,54
102,37
102,68
102,76
102,82
103,31
103,23
103,6
103,95
103,93
104,25
104,38
104,36
104,32
104,58
104,68
104,92
105,46
105,23
105,58
105,34
105,28
105,7
105,67
105,71
106,19
106,93
107,44
107,85
108,71
109,32
109,49
110,2
110,62
111,22
110,88
111,15
111,29
111,09
111,24
111,45
111,75
111,07
111,17
110,96
110,5
110,48
110,66
110,46

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