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

153,4
159,5
157,4
169,1
172,6
161,7
159,2
157,4
153,9
144,8
142,2
140,1
143,4
153,3
166,9
170,6
182,8
170,3
156,6
155,2
154,7
151,6
152,1
153,2
149,5
149,7
144,3
140
137,8
132,2
128,9
123,1
120,4
122,8
126
124,5
120,6
114,7
111,7
109,1
108
107,7
99,9
103,7
103,4
103,4
104,7
105,8
105,3
103
103,8
103,4
105,8
101,4
97
94,3
96,6
97,1
95,7
96,9
97,4
95,3
93,6
91,5
93,1
91,7
94,3
93,9
90,9
88,3
91,3
91,7
92,4
92
95,6
95,8
96,4
99
107
109,7
116,2
115,9
113,8
112,6
113,7
115,9
110,3
111,3
113,4
108,2
104,8
106
110,9
115
118,4
121,4
128,8
131,7
141,7
142,9
139,4
134,7
125
113,6
111,5
108,5
112,3
116,6
115,5
120,1
132,9
128,1
129,3
132,5
131
124,9
120,8
122
122,1
127,4
135,2
137,3
135
136
138,4
134,7
138,4
133,9
133,6
141,2
151,8
155,4
156,6
161,6
160,7
156
159,5
168,7
169,9
169,9
185,9

174,1
180,4
182,6
207,1
213,7
186,5
179,1
168,3
156,5
144,3
138,9
137,8
136,3
140,3
149,1
149,2
140,4
129
124,7
130,8
130,1
133,2
130,1
126,6
124,8
125,3
126,9
120,1
118,7
117,7
113,4
107,5
107,6
114,3
114,9
111,2
109,9
108,6
109,2
106,4
103,7
103
96,9
104,7
102,2
99
95,8
94,5
102,7
103,2
105,6
103,9
107,2
100,7
92,1
90,3
93,4
98,5
100,8
102,3
104,7
101,1
101,4
99,5
98,4
96,3
100,7
101,2
100,3
97,8
97,4
98,6
99,7
99
98,1
97
98,5
103,8
114,4
124,5
134,2
131,8
125,6
119,9
114,9
115,5
112,5
111,4
115,3
110,8
103,7
111,1
113
111,2
117,6
121,7
127,3
129,8
137,1
141,4
137,4
130,7
117,2
110,8
111,4
108,2
108,8
110,2
109,5
109,5
116
111,2
112,1
114
119,1
114,1
115,1
115,4
110,8
116
119,2
126,5
127,8
131,3
140,3
137,3
143
134,5
139,9
159,3
170,4
175
175,8
180,9
180,3
169,6
172,3
184,8
177,7
184,6
211,4

145,3
143,6
142,8
155,9
156,2
149,8
152,7
155,5
159,3
143
141,4
142,8
146,4
152,3
164,3
168
171,3
162,7
150,2
142,5
138,2
138
145,1
138,4
131,8
130,8
126,3
123
124
120,8
122,1
106,5
104,3
108,7
113,8
112,5
106,1
98,4
96
99,3
97,5
95,3
88
94,7
99,4
98,9
96,4
95,3
99,5
101,6
103,9
106,6
108,3
102
93,8
91,6
97,7
94,8
98
103,8
97,8
91,2
89,3
87,5
90,4
94,2
102,2
101,3
96
90,8
93,2
90,9
91,1
90,2
94,3
96
99
103,3
113,1
112,8
112,1
107,4
111
110,5
110,8
112,4
111,5
116,2
122,5
121,3
113,9
110,7
120,8
141,1
147,4
148
158,1
165
187
190,3
182,4
168,8
151,2
120,1
112,5
106,2
107,1
108,5
106,5
108,3
125,6
124
127,2
136,9
135,8
124,3
115,4
113,6
114,4
118,4
117
116,5
115,4
113,6
117,4
116,9
116,4
111,1
110,2
118,9
131,8
130,6
138,3
148,4
148,7
144,3
152,5
162,9
167,2
166,5
185,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()