Blog & Share Statistical Computations at FreeStatistics.org

Home » date » 2009 » Oct » 27 »

*The author of this computation has been verified*

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

Title produced by software: Bivariate Explorative Data Analysis

Date of computation: Tue, 27 Oct 2009 15:54:31 -0600

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099.htm/, Retrieved Tue, 27 Oct 2009 22:55:08 +0100

BibTeX entries for LaTeX users:

@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099.htm/},
    year = {2009},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2009},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

Dataseries X:

» Textbox « » Textfile « » CSV «

12,24603827 12,23081404 12,2075216 12,19646666 12,18170823 12,16209463 12,27264659 12,32731635 12,34834882 12,34261637 12,28408869 12,28632918 12,30231875 12,30190542 12,29428497 12,27798392 12,25821433 12,25330445 12,3759884 12,39702766 12,39618875 12,39925583 12,38231325 12,38979738 12,40800073 12,40728138 12,39492284 12,39064583 12,37196413 12,37580276 12,48589516 12,507469 12,50486156 12,4951425 12,46115275 12,47919803 12,49654956 12,49637379 12,48664324 12,46405424 12,44840687 12,44850104 12,53273622 12,56147528 12,56584977 12,55335058 12,53064218 12,53507059 12,53223477 12,53214456 12,52466091 12,50673264 12,49552814 12,48745869 12,56813953 12,58171752 12,58553595 12,57136138 12,54780602 12,55198589 12,54322322 12,54400096 12,53118055 12,52531152 12,52216723 12,51102104 12,57695321 12,58013652 12,58553595 12,53719488 12,50553259 12,49072889 12,50259304 12,48437083 12,44979499 12,44253905 12,4133185 12,36892252 12,46286758 12, etc...

Dataseries Y:

» Textbox « » Textfile « » CSV «

12,46412764 12,44925409 12,41996554 12,40811104 12,39214134 12,388115 12,50709163 12,55635269 12,55089563 12,51931281 12,4480615 12,44648228 12,45231925 12,44701682 12,43716453 12,41446392 12,40500137 12,41289574 12,54185178 12,55359433 12,54645508 12,5280257 12,50874431 12,51319331 12,52238591 12,51514799 12,50008837 12,50092644 12,48660925 12,50608094 12,62722736 12,64227709 12,63898105 12,60628063 12,57498096 12,58442346 12,59443229 12,58959203 12,57954817 12,55931771 12,54901675 12,57020869 12,66432942 12,69284735 12,6946864 12,70162827 12,6681701 12,66835592 12,65631033 12,65164269 12,64236113 12,62132305 12,61259053 12,62489811 12,71732618 12,72990728 12,72466404 12,68828581 12,65691575 12,65869166 12,65465471 12,64324312 12,62136596 12,60816206 12,60626722 12,61695969 12,70613556 12,72213354 12,69959747 12,6414816 12,59783091 12,58249678 12,58293936 12,559114 12,52987568 12,51563607 12,4874134 12,48101137 12,60083401 12,6234 etc...

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	4 seconds
R Server	'George Udny Yule' @ 72.249.76.132

Model: Y[t] = c + b X[t] + e[t]
c	3.3662805914397
b	0.73851187660556

Descriptive Statistics about e[t]
# observations	103
minimum	-0.141592955289681
Q1	-0.0184264736210185
median	0.00474741703457031
mean	4.32128249378264e-18
Q3	0.0361188689008246
maximum	0.086202567411398

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/197s51256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/197s51256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/2u54u1256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/2u54u1256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/3kyon1256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/3kyon1256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/4ip641256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/4ip641256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/50zc21256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/50zc21256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/6t4ra1256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/6t4ra1256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/703so1256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/703so1256680467.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/836uz1256680467.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566805003sdnjm1jtpfm099/836uz1256680467.ps (open in new window)

Parameters (Session):

par1 = 0 ; par2 = 36 ;

Parameters (R input):

par1 = 0 ; par2 = 36 ;

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

par1 <- as.numeric(par1)

par2 <- as.numeric(par2)

x <- as.ts(x)

y <- as.ts(y)

mylm <- lm(y~x)

cbind(mylm$resid)

library(lattice)

bitmap(file='pic1.png')

plot(y,type='l',main='Run Sequence Plot of Y[t]',xlab='time or index',ylab='value')

grid()

dev.off()

bitmap(file='pic1a.png')

plot(x,type='l',main='Run Sequence Plot of X[t]',xlab='time or index',ylab='value')

grid()

dev.off()

bitmap(file='pic1b.png')

plot(x,y,main='Scatter Plot',xlab='X[t]',ylab='Y[t]')

grid()

dev.off()

bitmap(file='pic1c.png')

plot(mylm$resid,type='l',main='Run Sequence Plot of e[t]',xlab='time or index',ylab='value')

grid()

dev.off()

bitmap(file='pic2.png')

hist(mylm$resid,main='Histogram of e[t]')

dev.off()

bitmap(file='pic3.png')

if (par1 > 0)

{

densityplot(~mylm$resid,col='black',main=paste('Density Plot of e[t]   bw = ',par1),bw=par1)

} else {

densityplot(~mylm$resid,col='black',main='Density Plot of e[t]')

}

dev.off()

bitmap(file='pic4.png')

qqnorm(mylm$resid,main='QQ plot of e[t]')

qqline(mylm$resid)

grid()

dev.off()

if (par2 > 0)

{

bitmap(file='pic5.png')

acf(mylm$resid,lag.max=par2,main='Residual Autocorrelation Function')

grid()

dev.off()

}

summary(x)

load(file='createtable')

a<-table.start()

a<-table.row.start(a)

a<-table.element(a,'Model: Y[t] = c + b X[t] + e[t]',2,TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'c',1,TRUE)

a<-table.element(a,mylm$coeff[[1]])

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'b',1,TRUE)

a<-table.element(a,mylm$coeff[[2]])

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,'Descriptive Statistics about e[t]',2,TRUE)

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'# observations',header=TRUE)

a<-table.element(a,length(mylm$resid))

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'minimum',header=TRUE)

a<-table.element(a,min(mylm$resid))

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'Q1',header=TRUE)

a<-table.element(a,quantile(mylm$resid,0.25))

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'median',header=TRUE)

a<-table.element(a,median(mylm$resid))

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'mean',header=TRUE)

a<-table.element(a,mean(mylm$resid))

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'Q3',header=TRUE)

a<-table.element(a,quantile(mylm$resid,0.75))

a<-table.row.end(a)

a<-table.row.start(a)

a<-table.element(a,'maximum',header=TRUE)

a<-table.element(a,max(mylm$resid))

a<-table.row.end(a)

a<-table.end(a)

table.save(a,file='mytable.tab')

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa

Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.

We may request personal information to be submitted to our servers in order to be able to:

personalize online software applications according to your needs
enforce strict security rules with respect to the data that you upload (e.g. statistical data)
manage user sessions of online applications
alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.

FreeStatistics.org is powered by