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

106.48
106.83
107.14
107.94
108.46
108.81
108.92
108.99
109.16
109.22
109.43
109.23
109.93
110.09
110.33
110.11
110.35
110.09
110.44
110.39
110.62
110.43
110.46
110.55
110.94
111.56
111.82
111.73
111.57
111.85
112.06
112.2
112.47
112.15
112.36
112.32
112.67
113.02
113.05
113.5
113.67
113.65
114
114.03
114.08
114.49
114.48
114.25
114.68
115.28
115.9
115.87
116.09
116.29
116.76
116.78
116.65
116.46
116.82
116.91

106.22
106.31
107.38
109.31
110.82
111.22
110.66
110.76
110.69
111.08
110.97
110.24
112.51
111.52
112.13
112.23
112.92
111.89
111.99
111.51
112.33
112.04
112.09
111.41
112.61
113.14
113.65
114.26
114.4
114.93
114.86
114.95
116.17
114.6
114.62
113.82
115.02
115.18
115.59
116.6
117.07
116.96
116.66
116.07
116.04
115.81
116.22
115.85
116.43
117.39
119.17
119.24
120.03
119.34
118.49
118.59
117.5
117.56
118.25
118.01

107,97
108,13
108,54
109,86
109,75
109,99
112,01
111,96
111,41
112,11
111,67
111,95
112,31
113,26
113,5
114,43
115,02
115,1
117,11
117,52
116,1
116,39
116,01
116,74
116,68
117,45
117,8
119,37
118,9
119,05
120,46
120,99
119,86
120,18
119,81
120,15
119,8
120,27
120,71
121,87
121,87
121,92
123,72
124,38
123,21
123,17
122,95
123,46
123,24
123,86
124,28
124,78
125,19
125,46
127,6
127,8
126,63
127,06
126,77
127,05

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