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

96,8
91,2
97,1
104,9
110,9
104,8
94,1
95,8
99,3
101,1
104,0
99,0
105,4
107,1
110,7
117,1
118,7
126,5
127,5
134,6
131,8
135,9
142,7
141,7
153,4
145,0
137,7
148,3
152,2
169,4
168,6
161,1
174,1
179,0
190,6
190,0
181,6
174,8
180,5
196,8
193,8
197,0
216,3
221,4
217,9
229,7
227,4
204,2
196,6
198,8
207,5
190,7
201,6
210,5
223,5
223,8
231,2
244,0
234,7
250,2

96,9
98,0
97,9
100,9
103,9
103,1
102,5
104,3
102,6
101,7
102,8
105,4
110,9
113,5
116,3
124,0
128,8
133,5
132,6
128,4
127,3
126,7
123,3
123,2
124,4
128,2
128,7
135,7
139,0
145,4
142,4
137,7
137,0
137,1
139,3
139,6
140,4
142,3
148,3
157,7
161,6
161,7
171,8
185,1
176,7
184,4
183,0
180,9
187,0
189,9
193,8
194,5
198,7
204,7
213,2
214,7
211,0
213,2
206,3
210,8

96,7
88,0
96,7
106,8
114,3
105,7
90,1
91,6
97,7
100,8
104,6
95,9
102,7
104,0
107,9
113,8
113,8
123,1
125,1
137,6
134,0
140,3
152,1
150,6
167,3
153,2
142,0
154,4
158,5
180,9
181,3
172,4
192,0
199,3
215,4
214,3
201,5
190,5
196,0
215,7
209,4
214,1
237,8
239,0
237,8
251,5
248,8
215,4
201,2
203,1
214,2
188,9
203,0
213,3
228,5
228,2
240,9
258,8
248,5
269,2

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