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

91,25
91,5
91,68
91,81
91,84
91,93
92,08
92,11
92,26
92,28
92,39
92,46
92,82
93,16
93,33
93,51
93,56
93,67
93,76
93,88
94,01
94,21
94,31
94,4
94,9
95,31
95,52
95,68
95,91
95,97
96,15
96,34
96,42
96,54
96,72
96,81
97,19
97,5
97,71
97,86
98,04
98,2
98,25
98,41
98,56
98,62
98,75
98,71
99,05
99,52
99,71
99,8
100,01
99,99
100,12
100,15
100,27
100,42
100,43
100,5
100,95
101,26
101,42
101,68
101,75
101,89
102,07
102,22
102,45
102,62
102,67
102,86
104,78
104,87
105,06
105,14
105,32
105,54
105,68
105,77
106,07
106,03
106,13
106,28
106,61
106,74
107,01
107,1
107,28
107,4
107,59
107,69
107,78

1,79
1,95
2,26
2,04
2,16
2,75
2,79
2,88
3,36
2,97
3,1
2,49
2,2
2,25
2,09
2,79
3,14
2,93
2,65
2,67
2,26
2,35
2,13
2,18
2,9
2,63
2,67
1,81
1,33
0,88
1,28
1,26
1,26
1,29
1,1
1,37
1,21
1,74
1,76
1,48
1,04
1,62
1,49
1,79
1,8
1,58
1,86
1,74
1,59
1,26
1,13
1,92
2,61
2,26
2,41
2,26
2,03
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

90,49
90,65
90,95
91,19
91,07
91,15
91,81
91,95
91,62
91,27
91,4
91,76
91,99
92,34
92,3
92,85
92,94
93,26
94,21
94,08
93,98
94,23
94,93
95,09
95,37
96,23
96,2
95,43
95,63
95,96
96,51
96,65
96,21
95,54
95,96
96,41
96,32
96,94
96,97
97,63
97,33
97,66
98,18
98,22
97,91
97,93
98,4
98,78
98,73
99,4
99,04
99,68
99,62
99,8
100,65
100,59
100,46
100,57
100,75
100,7
101,44
101,77
101,79
101,52
101,83
102,23
103,04
102,81
102,48
102,81
103,21
103,21
102,92
103,48
103,18
103,39
103,5
103,73
104,42
104,53
104,09
104,23
104,23
104,54
104,65
105,48
105,61
105,74
105,86
105,81
106,49
106,43
105,73

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