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

 106 474.50 
 106 559.75 
 106 636.75 
 106 722.00 
 106 804.50 
 106 889.75 
 106 972.25 
 107 057.50 
 107 142.75 
 107 225.25 
 107 310.50 
 107 393.00 
 107 478.25 
 107 563.50 
 107 640.50 
 107 725.75 
 107 808.25 
 107 893.50 
 107 976.00 
 108 061.25 
 108 146.50 
 108 229.00 
 108 314.25 
 108 396.75 
 108 482.00 
 108 567.25 
 108 647.00 
 108 732.25 
 108 814.75 
 108 900.00 
 108 982.50 
 109 067.75 
 109 153.00 
 109 235.50 
 109 320.75 
 109 403.25 
 109 488.50 
 109 573.75 
 109 650.75 
 109 736.00 
 109 818.50 
 109 903.75 
 109 986.25 
 110 071.50 
 110 156.75 
 110 239.25 
 110 324.50 
 110 407.00 
 110 492.25 
 110 577.50 
 110 654.50 
 110 739.75 
 110 822.25 
 110 907.50 
 110 990.00 
 111 075.25 
 111 160.50 
 111 243.00 
 111 328.25 
 111 410.75 
 111 496.00 
 111 581.25 
 111 658.25 
 111 743.50 
 111 826.00 
 111 911.25 
 111 993.75 
 112 079.00 
 112 164.25 
 112 246.75 
 112 332.00 
 112 414.50 
 112 499.75 
 112 585.00 
 112 664.75 
 112 750.00 
 112 832.50 
 112 917.75 
 113 000.25 
 113 085.50 
 113 170.75 
 113 253.25 
 113 338.50 
 113 421.00 
 113 506.25 
 113 591.50 
 113 668.50 
 113 753.75 
 113 836.25 
 113 921.50 
 114 004.00 
 114 089.25 
 114 174.50 
 114 257.00 
 114 342.25 
 114 424.75 
 114 510.00 
 114 595.25 
 114 672.25 
 114 757.50 
 114 840.00 
 114 925.25 
 115 007.75 
 115 093.00 
 115 178.25 
 115 260.75 
 115 346.00 
 115 428.50 
 115 513.75 
 115 599.00 
 115 676.00 
 115 761.25 
 115 843.75 
 115 929.00 
 116 011.50 
 116 096.75 
 116 182.00 
 116 264.50 
 116 349.75 
 116 432.25 

par1 <- as.numeric(par1)
(n <- length(x))
(np <- floor(n / par1))
arr <- array(NA,dim=c(par1,np))
j <- 0
k <- 1
for (i in 1:(np*par1))
{
j = j + 1
arr[j,k] <- x[i]
if (j == par1) {
j = 0
k=k+1
}
}
arr
arr.mean <- array(NA,dim=np)
arr.sd <- array(NA,dim=np)
arr.range <- array(NA,dim=np)
for (j in 1:np)
{
arr.mean[j] <- mean(arr[,j],na.rm=TRUE)
arr.sd[j] <- sd(arr[,j],na.rm=TRUE)
arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE)
}
arr.mean
arr.sd
arr.range
(lm1 <- lm(arr.sd~arr.mean))
(lnlm1 <- lm(log(arr.sd)~log(arr.mean)))
(lm2 <- lm(arr.range~arr.mean))
bitmap(file='test1.png')
plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation')
dev.off()
bitmap(file='test2.png')
plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Section',header=TRUE)
a<-table.element(a,'Mean',header=TRUE)
a<-table.element(a,'Standard Deviation',header=TRUE)
a<-table.element(a,'Range',header=TRUE)
a<-table.row.end(a)
for (j in 1:np) {
a<-table.row.start(a)
a<-table.element(a,j,header=TRUE)
a<-table.element(a,arr.mean[j])
a<-table.element(a,arr.sd[j] )
a<-table.element(a,arr.range[j] )
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,lm1$coefficients[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,lm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,4])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,lnlm1$coefficients[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,lnlm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,4])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Lambda',header=TRUE)
a<-table.element(a,1-lnlm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable2.tab')