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

97.3
101
113.2
101
105.7
113.9
86.4
96.5
103.3
114.9
105.8
94.2
98.4
99.4
108.8
112.6
104.4
112.2
81.1
97.1
112.6
113.8
107.8
103.2
103.3
101.2
107.7
110.4
101.9
115.9
89.9
88.6
117.2
123.9
100
103.6
94.1
98.7
119.5
112.7
104.4
124.7
89.1
97
121.6
118.8
114
111.5
97.2
102.5
113.4
109.8
104.9
126.1
80
96.8
117.2
112.3
117.3
111.1
102.2
104.3
122.9
107.6
121.3
131.5
89
104.4
128.9
135.9
133.3
121.3

104.8
105.6
118.3
89.9
90.2
107
64.5
92.6
95.8
94.3
91.2
86.3
77.6
82.5
97.7
83.3
84.2
92.8
77.4
72.5
88.8
93.4
92.6
90.7
81.6
84.1
88.1
85.3
82.9
84.8
71.2
68.9
94.3
97.6
85.6
91.9
75.8
79.8
99
88.5
86.7
97.9
94.3
72.9
91.8
93.2
86.5
98.9
77.2
79.4
90.4
81.4
85.8
103.6
73.6
75.7
99.2
88.7
94.6
98.7
84.2
87.7
103.3
88.2
93.4
106.3
73.1
78.6
101.6
101.4
98.5
99

124.9
132
151.4
108.9
121.3
123.4
90.3
79.3
117.2
116.9
120.8
96.1
100.8
105.3
116.1
112.8
114.5
117.2
77.1
80.1
120.3
133.4
109.4
93.2
91.2
99.2
108.2
101.5
106.9
104.4
77.9
60
99.5
95
105.6
102.5
93.3
97.3
127
111.7
96.4
133
72.2
95.8
124.1
127.6
110.7
104.6
112.7
115.3
139.4
119
97.4
154
81.5
88.8
127.7
105.1
114.9
106.4
104.5
121.6
141.4
99
126.7
134.1
81.3
88.6
132.7
132.9
134.4
103.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()