par3 <- '0' par2 <- '5' par1 <- '50' par1 <- as#numeric#par1# par2 <- as#numeric#par2# if #par3 == '0'# bw <- NULL if #par3 != '0'# bw <- as#numeric#par3# if #par1 < 10# par1 = 10 if #par1 > 5000# par1 = 5000 library#modeest# library#lattice# library#boot# boot#stat <- function#s,i# { s#mean <- mean#s#i## s#median <- median#s#i## s#midrange <- #max#s#i## + min#s#i### / 2 s#mode <- mlv#s#i#, method='mfv'#$M s#kernelmode <- mlv#s#i#, method='kernel', bw=bw#$M c#s#mean, s#median, s#midrange, s#mode, s#kernelmode# } bitmap#file='plot1#png'# plot#r$t#,1#,type='p',ylab='simulated values',main='Simulation of Mean'# grid## dev#off## bitmap#file='plot2#png'# plot#r$t#,2#,type='p',ylab='simulated values',main='Simulation of Median'# grid## dev#off## bitmap#file='plot3#png'# plot#r$t#,3#,type='p',ylab='simulated values',main='Simulation of Midrange'# grid## dev#off## bitmap#file='plot7#png'# plot#r$t#,4#,type='p',ylab='simulated values',main='Simulation of Mode'# grid## dev#off## bitmap#file='plot8#png'# plot#r$t#,5#,type='p',ylab='simulated values',main='Simulation of Mode of Kernel Density'# grid## dev#off## bitmap#file='plot4#png'# densityplot#~r$t#,1#,col='black',main='Density Plot',xlab='mean'# dev#off## bitmap#file='plot5#png'# densityplot#~r$t#,2#,col='black',main='Density Plot',xlab='median'# dev#off## bitmap#file='plot6#png'# densityplot#~r$t#,3#,col='black',main='Density Plot',xlab='midrange'# dev#off## bitmap#file='plot9#png'# densityplot#~r$t#,4#,col='black',main='Density Plot',xlab='mode'# dev#off## bitmap#file='plot10#png'# densityplot#~r$t#,5#,col='black',main='Density Plot',xlab='mode of kernel dens#'# dev#off## z <- data#frame#cbind#r$t#,1#,r$t#,2#,r$t#,3#,r$t#,4#,r$t#,5### colna#es#z# <- list#'mean','median','midrange','mode','mode k#dens'# bitmap#file='plot11#png'# boxplot#z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency'# grid## dev#off## load#file='createtable'# a<-table#start## a<-table#row#start#a# a<-table#element#a,'Estimation Results of Bootstrap',6,TRUE# a<-table#row#end#a# a<-table#row#start#a# a<-table#element#a,'statistic',header=TRUE# a<-table#element#a,'Q1',header=TRUE# a<-table#element#a,'Estimate',header=TRUE# a<-table#element#a,'Q3',header=TRUE# a<-table#element#a,'S#D#',header=TRUE# a<-table#element#a,'IQR',header=TRUE# a<-table#row#end#a# a<-table#row#start#a# a<-table#element#a,'mean',header=TRUE# q1 <- quantile#r$t#,1#,0#25###1## q3 <- quantile#r$t#,1#,0#75###1## a<-table#element#a,signif#q1,par2## a<-table#element#a,signif#r$t0#1#,par2## a<-table#element#a,signif#q3,par2## a<-table#element# a,signif# sqrt#var#r$t#,1###,par2 # # a<-table#element#a,signif#q3-q1,par2## a<-table#row#end#a# a<-table#row#start#a# a<-table#element#a,'median',header=TRUE# q1 <- quantile#r$t#,2#,0#25###1## q3 <- quantile#r$t#,2#,0#75###1## a<-table#element#a,signif#q1,par2## a<-table#element#a,signif#r$t0#2#,par2## a<-table#element#a,signif#q3,par2## a<-table#element#a,signif#sqrt#var#r$t#,2###,par2## a<-table#element#a,signif#q3-q1,par2## a<-table#row#end#a# a<-table#row#start#a# a<-table#element#a,'midrange',header=TRUE# q1 <- quantile#r$t#,3#,0#25###1## q3 <- quantile#r$t#,3#,0#75###1## a<-table#element#a,signif#q1,par2## a<-table#element#a,signif#r$t0#3#,par2## a<-table#element#a,signif#q3,par2## a<-table#element#a,signif#sqrt#var#r$t#,3###,par2## a<-table#element#a,signif#q3-q1,par2## a<-table#row#end#a# a<-table#row#start#a# a<-table#element#a,'mode',header=TRUE# q1 <- quantile#r$t#,4#,0#25###1## q3 <- quantile#r$t#,4#,0#75###1## a<-table#element#a,signif#q1,par2## a<-table#element#a,signif#r$t0#4#,par2## a<-table#element#a,signif#q3,par2## a<-table#element#a,signif#sqrt#var#r$t#,4###,par2## a<-table#element#a,signif#q3-q1,par2## a<-table#row#end#a# a<-table#row#start#a# a<-table#element#a,'mode k#dens',header=TRUE# q1 <- quantile#r$t#,5#,0#25###1## q3 <- quantile#r$t#,5#,0#75###1## a<-table#element#a,signif#q1,par2## a<-table#element#a,signif#r$t0#5#,par2## a<-table#element#a,signif#q3,par2## a<-table#element#a,signif#sqrt#var#r$t#,5###,par2## a<-table#element#a,signif#q3-q1,par2## a<-table#row#end#a# a<-table#end#a# table#save#a,file='mytable.tab'#
|