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

113.5
121.2
130.4
115.2
117.9
110.7
107.6
124.3
115.1
112.5
127.9
117.4
119.3
130.4
126
125.4
130.5
115.9
108.7
124
119.4
118.6
131.3
111.1
124.8
132.3
126.7
131.7
130.9
122.1
113.2
133.6
119.2
129.4
131.4
117.1
130.5
132.3
140.8
137.5
128.6
126.7
120.8
139.3
128.6
131.3
136.3
128.8
133.2
136.3
151.1
145
134.4
135.7
128.7
129.2
138.6
132.7
132.5
137.3
127.1

86.7
123.6
125.3
111.1
98.4
102.3
105
128.2
124.7
116.1
131.2
97.7
88.8
132.8
113.9
112.6
104.3
107.5
106
117.3
123.1
114.3
132
92.3
93.7
121.3
113.6
116.3
98.3
111.9
109.3
133.2
118
131.6
134.1
96.7
99.8
128.3
134.9
130.7
107.3
121.6
120.6
140.5
124.8
129.9
159.4
111
110.1
132.7
135
118.6
94
117.9
114.7
113.6
130.6
117.1
123.2
106.1
87.9

105.1
106.1
106.2
103.9
109.2
110.5
98
94.1
90.2
89.5
91.2
98.2
103.7
103.9
106.5
107.2
111
111.8
101.5
95.3
92.7
93.5
96.2
102.1
102.3
127.9
130.8
134.9
141.9
124.6
118
115.1
111.2
113.5
115.2
119.4
116
115.7
121.1
120.8
125.2
124
119.1
119.2
113.9
113.3
116.8
114.8
119.2
117.8
122.5
125.1
125
125.1
121.2
118.9
109.8
109.2
109
110.9
112.5

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