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

97.33
97.89
98.69
99.01
99.18
98.45
98.13
98.29
99.1
99.26
98.85
98.05
98.53
99.34
100.14
100.3
100.22
99.9
99.58
99.9
100.78
100.78
100.46
100.06
100.28
100.78
101.58
102.06
102.02
101.68
101.32
101.81
102.3
102.12
102.1
101.75
101.5
102.16
103.47
104.05
104.09
103.55
102.77
102.89
103.6
103.76
103.92
103.35
103.32
104.2
105.44
105.81
106.25
105.94
105.82
105.96
106.49
106.32
105.88
105.07

5.7
6.1
6
5.9
5.8
5.7
5.6
5.4
5.4
5.5
5.6
5.7
5.9
6.1
6
5.8
5.8
5.7
5.5
5.3
5.2
5.2
5
5.1
5.1
5.2
4.9
4.8
4.5
4.5
4.4
4.4
4.2
4.1
3.9
3.8
3.9
4.2
4.1
3.8
3.6
3.7
3.5
3.4
3.1
3.1
3.1
3.2
3.3
3.5
3.6
3.5
3.3
3.2
3.1
3.2
3
3
3.1
3.4

91.4
91.1
104.4
97.6
93.7
104.5
95.4
86.5
102.9
101.9
103.7
100.7
94.2
93.6
104.7
101
97.6
105.8
93.7
91.2
106.3
103.4
107.4
101.2
96.9
96.3
109.8
97.9
105.1
107.9
95
95.2
105.8
110.1
112.2
102.5
103.7
102
112.3
103.3
106.9
104.6
100.7
99
106.5
114.9
114.1
102.2
107
107.4
107.4
110.1
105.6
110.9
101.9
93.2
110.5
113.1
101.7
96.7

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