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

132.7
132.7
132.8
132.4
131.5
131
130.7
130.1
129.1
128
127
126
125
123.7
122.1
120.2
118.5
117.2
116.2
115.1
113.8
112.7
111.8
111.2
110.6
109.7
108.6
107.8
107.5
107.5
107.4
107.2
107.1
107.2
107.5
107.9
108.4
108.8
109
109
109
108.8
108.4
108.1
107.8
107.3
106.6
105.9
105.4
105.4
105.6
105.6
105.6
105.7
105.6
105.4
105.4
105.5
105.7
105.8

151.3
151.2
150.8
149.9
148.9
148
146.8
145.3
144.1
143.1
142
140.8
138.9
136.2
133.4
130.7
128.5
126.9
125.8
124.6
122.7
120.9
119.3
118
117
115.9
115.1
114.7
114.1
113.3
112.7
112.3
112.1
112.1
111.9
111.6
111.3
110.9
110.7
110.8
110.6
110
109.1
108.9
108.9
108.2
107
106
105.2
104.4
103.9
103.4
103
102.6
101.8
100.7
99.9
99.9
100.6
101

98.5
97.9
96.2
93.7
93
93.6
93
92.1
91.4
90.3
90.3
91.6
92.4
93
93
91.9
91.6
93.2
94.4
94.2
93.6
93
92.3
91.4
90.5
89.1
87.4
86.7
87.4
89.4
90.5
88.5
85.3
83.6
83.9
85
85
84.2
83.8
83.6
84.2
86
88
89.4
89.6
88.4
86.1
84.5
84.7
86.1
87.3
87.4
86.6
85.3
84.1
84.8
86
86.8
88.2
89.6

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