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

97,57
97,74
97,92
98,19
98,23
98,41
98,59
98,71
99,14
99,62
100,18
100,66
101,19
101,75
102,2
102,87
98,81
97,6
96,68
95,96
98,89
99,05
99,2
99,11
99,19
99,77
100,70
100,78
100,53
101,01
100,92
101,10
103,11
102,99
102,31
102,61
103,68
104,72
107,66
108,87
108,12
107,61
106,42
105,61
105,71
105,49
105,57
105,18
106,09
106,34
108,47
116,87
121,08
123,27
124,18
125,60
126,57
127,18
128,04
128,55
129,67

127,96
127,47
126,47
125,75
125,42
125,14
125,15
125,51
125,63
126,22
126,88
127,96
128,74
129,6
131,2
132,72
134,67
135,94
136,39
136,74
137,2
137,36
138,63
141,07
143,32
147,91
152,56
151,61
156,56
157,45
158,13
159,18
159,47
159,79
161,65
162,77
163,48
166,16
163,86
162,12
149,08
145,32
141,21
134,68
133,65
139,17
138,61
144,96
157,99
167,18
174,48
182,77
190,00
189,70
188,90
198,28
201,18
204,14
221,02
221,12
220,68

20,72
21,45
22,09
21,53
23,35
23,57
26,42
25,21
26,44
29,34
29,40
33,05
28,38
26,01
29,31
30,36
35,75
36,15
34,21
37,91
38,70
42,12
42,16
39,80
37,36
38,35
42,60
41,25
42,16
46,94
47,43
47,06
50,18
50,13
43,23
40,04
40,37
42,21
37,00
39,74
42,68
46,29
46,97
48,73
52,37
50,05
54,04
57,78
64,72
63,41
64,36
66,03
72,14
76,60
86,97
93,48
95,59
81,89
70,55
50,38
36,25

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