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

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

89.28
89.47
89.53
90.72
90.91
91.38
91.49
90.9
90.93
90.57
91.28
90.83
91.5
91.58
92.49
94.16
95.46
95.8
95.32
95.41
95.35
95.68
95.59
94.96
96.92
96.06
96.59
96.67
97.27
96.38
96.47
96.05
96.76
96.51
96.55
95.97
97
97.46
97.9
98.42
98.54
99
98.94
99.02
100.07
98.72
98.73
98.04
99.08
99.22
99.57
100.44
100.84
100.75
100.49
99.98
99.96
99.76
100.11
99.79
100.29
101.12
102.65
102.71
103.39
102.8
102.07
102.15
101.21
101.27
101.86
101.65
101.94
102.62
102.71
103.39
104.51
104.09
104.29
104.57
105.39
105.15
106.13
105.46
106.47
106.62
106.52
108.04
107.15
107.32
107.76
107.26
107.89

92.58
93.15
93.7
93.34
93.93
94.41
94.3
94.44
96.09
95.99
96.23
95.64
94.88
95.48
95.52
96.18
96.82
96.81
96.25
96.24
96.95
96.25
96.04
95.78
95.86
96.02
96.34
96.84
96.73
96.34
96.6
96.64
97.2
97.5
96.99
97.08
97.55
98.42
98.78
97.49
96.99
97.16
97.29
97.8
98.12
98.03
98.11
98.07
98.21
98.48
98.83
99.2
99.88
99.71
100.03
100.6
100.85
101.96
101.4
100.81
100.66
101.55
102.23
102.9
102.68
103.41
104.62
104.93
105.88
105.18
104.54
104.58
104.34
104.66
104.73
105.44
105.72
105.68
105.9
105.97
105.21
104.75
104.89
105.26
104.84
105.47
105.4
105.73
105.72
105.63
105.97
105.92
106.32

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