Blog & Share Statistical Computations at FreeStatistics.org

Home » date » 2009 » Dec » 08 »

*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, 08 Dec 2009 13:28:08 -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/08/t12603041224ymh1mtqgis6340.htm/, Retrieved Tue, 08 Dec 2009 21:28:51 +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/Dec/08/t12603041224ymh1mtqgis6340.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 «

543,058008 541,7453276 539,031539 533,6065967 530,8653313 536,5202699 562,3726523 570,448946 570,9737297 572,9589863 563,4536361 563,5059893 560,1223081 558,8166068 556,2292693 550,4089389 548,0109488 551,3936888 577,4738089 581,1178882 579,5964113 569,1493653 560,2918882 560,7896219 559,6588246 556,4746176 550,4207482 546,7988661 546,2810632 549,2094318 574,251687 578,8635418 572,3774978 555,984712 543,9816173 539,8268241 539,9462936 533,5522467 525,8089006 522,0785382 514,7630523 513,1179202 544,7990455 550,9972777 535,7182095 525,7632547 516,162765 516,8684552 518,0550164 512,3690076 505,511622 503,1480895 493,7641137 500,7404517 529,5705808 534,0945609 520,5535515 510,9559668 505,5719533 510,2705165 513,7149015 514,0924041 513,1042779 512,4177983 506,1294301 514,0243185 542,0885537

Dataseries Y:

» Textbox « » Textfile « » CSV «

517,1199087 517,0744627 514,5648647 508,7858096 504,8207603 504,8445305 526,5614874 534,1825531 535,3522205 532,016917 526,0104562 527,1764411 526,4294825 526,4057371 524,439701 519,7595598 516,8558793 514,774708 535,965484 539,6165305 540,6477596 536,8295819 530,5440604 531,6540228 529,3297649 529,5356456 526,1520693 524,6103316 523,7862159 520,875225 538,3326109 539,1901334 540,6477596 527,7366768 519,4477837 515,617106 518,6848754 513,9805444 505,1712581 503,3418322 496,0413289 485,1515227 508,4840214 510,87474 504,6414569 500,6425871 493,3781511 497,0965701 498,5388651 495,0141412 486,9086157 486,9137501 474,9252573 476,2761804 497,9297139 498,3302118 496,9597569 495,0737319 496,2489295 505,7321425 514,1196359 518,0222003 522,5380369 523,4147495 517,1363457 521,5294814 541,0268016

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

Model: Y[t] = c + b X[t] + e[t]
c	194.106221428094
b	0.595386919836408

Descriptive Statistics about e[t]
# observations	67
minimum	-15.9643562032483
Q1	-4.64722475640127
median	-0.479680445329828
mean	3.30633116242147e-16
Q3	2.72590383263867
maximum	24.2216734727156

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/12j9w1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/12j9w1260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/2kcus1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/2kcus1260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/3nttk1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/3nttk1260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/40w6q1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/40w6q1260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/5n6yr1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/5n6yr1260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/65no41260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/65no41260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/7etfq1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/7etfq1260304085.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/8chyx1260304085.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Dec/08/t12603041224ymh1mtqgis6340/8chyx1260304085.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