Home » date » 2009 » Dec » 16 »

niet werkende werkzoekende tss 20-25 jaar frederic.ledent@student.kdg.be

*Unverified author*

R Software Module: /rwasp_bootstrapplot1.wasp (opens new window with default values)

Title produced by software: Bootstrap Plot - Central Tendency

Date of computation: Wed, 16 Dec 2009 06:08:21 -0700

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr.htm/, Retrieved Fri, 24 May 2013 12:37:51 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

KDGP2W32

System-generated keywords (parent):

(pk = 0)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Estimation Results of Bootstrap
statistic	Q1	Estimate	Q3	S.D.	IQR
mean	97160.681547619	97591.5476190476	98840.7619047619	1608.59590741294	1680.08035714286
median	93909	96307.5	96815.25	2120.94586680046	2906.25
midrange	96925.5	96925.5	98131	1356.77405003823	1205.5

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/1k7iu1260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/1k7iu1260968898.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/2irn91260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/2irn91260968898.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/3ximz1260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/3ximz1260968898.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/4j6c81260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/4j6c81260968898.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/5f1k91260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/5f1k91260968898.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/6htmu1260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/6htmu1260968898.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/7b2lh1260968898.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Dec/16/t1260969040gm59zjxk7sriqlr/7b2lh1260968898.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

par1 = 50 ;

Parameters (R input):

par1 = 50 ;

R code (references can be found in the software module):

par1 <- as.numeric(par1)
if (par1 < 10) par1 = 10
if (par1 > 5000) par1 = 5000
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
c(s.mean, s.median, s.midrange)
}
(r <- boot(x,boot.stat, R=par1, stype='i'))
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='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()
z <- data.frame(cbind(r$t[,1],r$t[,2],r$t[,3]))
colnames(z) <- list('mean','median','midrange')
bitmap(file='plot7.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,q1)
a<-table.element(a,r$t0[1])
a<-table.element(a,q3)
a<-table.element(a,sqrt(var(r$t[,1])))
a<-table.element(a,q3-q1)
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,q1)
a<-table.element(a,r$t0[2])
a<-table.element(a,q3)
a<-table.element(a,sqrt(var(r$t[,2])))
a<-table.element(a,q3-q1)
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,q1)
a<-table.element(a,r$t0[3])
a<-table.element(a,q3)
a<-table.element(a,sqrt(var(r$t[,3])))
a<-table.element(a,q3-q1)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')