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

100,70
97,90
96,50
96,60
96,60
95,50
91,80
89,30
87,00
85,90
88,00
87,90
89,20
90,90
91,60
90,20
89,10
87,50
86,30
86,00
84,40
86,10
91,00
92,70
88,00
84,30
82,20
80,80
79,40
80,20
82,20
82,20
81,20
82,10
88,10
88,50
92,10
98,60
100,90
100,60
101,10
102,10
103,60
102,80
108,30
104,00
106,10
106,30
109,00
111,00
113,70
112,70
110,30
114,50
119,30
121,80
125,40
129,70
129,40
134,50
141,20
141,40
152,20
167,70
173,30
168,70
172,60
169,80
172,00
179,40
174,60
172,50
172,60
176,30
178,90
179,60
179,90
180,30
180,90
177,70

101,30
97,60
96,40
97,00
96,40
94,70
89,30
85,90
83,30
81,50
85,00
84,80
87,50
89,00
90,00
89,60
87,40
84,80
81,90
81,10
79,10
80,50
88,50
90,90
84,90
80,00
76,50
75,40
73,50
74,30
77,70
77,90
76,70
77,20
86,00
86,90
92,00
101,70
104,50
101,70
100,60
100,30
102,50
101,00
108,60
103,40
106,40
106,60
108,90
110,50
114,00
112,80
109,60
116,00
124,60
129,00
131,50
138,60
138,10
146,30
157,60
158,40
176,30
199,90
210,40
202,60
207,10
202,00
203,40
216,30
207,30
203,50
204,40
203,70
205,70
208,00
209,30
208,70
206,50
204,50

99,80
103,50
103,10
105,60
105,70
106,60
107,00
105,20
105,70
105,00
105,10
105,90
105,30
104,90
103,20
103,40
104,40
104,50
105,90
110,60
112,40
111,80
111,00
111,00
109,10
107,80
107,20
108,40
107,50
106,40
106,20
104,90
106,20
107,60
107,00
104,50
105,10
104,70
103,70
104,90
105,90
106,10
106,10
106,80
106,40
107,80
107,60
107,60
108,40
109,50
109,20
109,10
110,00
109,00
109,00
111,90
109,30
112,10
112,10
112,50
113,60
112,90
114,00
116,10
116,50
117,10
117,10
117,10
116,50
116,50
116,30
116,50
119,20
118,60
117,50
117,10
117,60
118,30
118,60
116,70

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