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