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

356,2
359,5
368,4
371
397,5
416,7
413,2
424,3
415
421,7
422,1
429,2
452,1
471,5
488,3
506,2
517,3
538,6
545,3
546,7
540,3
549,2
563,9
581,7
590,7
594,1
604
628,1
662,4
688,6
705,9
701,5
686,2
645,7
668,7
696,7
715,5
741,4
754,3
771,3
797,7
809,9
790,1
830,3
847,7
834,8
824,5
764,6
780
803,2
751,1
755,2
708,2
685,4
680
710,6
702,8
656,3
575,6
567,2
545,2

152823,6
123780,5
159987,1
139603,7
177831,2
173656,9
252392
228029
197300
214088
160275
186851
227777
246899
295338
243847
324602
347066
407916
312914
326127
394369
310078
422770
417974
402347
360809
298289
375873
407210
413968
457532
695731
544623
292833
534403
517030
455714
471401
451493
480615
568272
650780
553643
780711
650724
586345
725173
701257
859063
789842
512707
780845
637804
640694
553416
554622
616736
536994
407237
618796

3
2,9
2,9
2,8
2,6
2,5
2,8
2,8
2,8
2,7
2,7
2,7
2,5
2,4
2,3
2,3
2,2
2,3
2,6
2,7
2,7
2,6
2,5
2,5
2,4
2,5
2,4
2,3
2,2
2,2
2,3
2,4
2,5
2,5
2,4
2,4
2,3
2,2
2,2
2
2
2,2
2,4
2,4
2,3
2,4
2,5
2,5
2,6
2,5
2,7
2,7
3,1
3
3,4
3,3
3,5
4,1
4,7
4,4
5,4

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