Blog & Share Statistical Computations at FreeStatistics.org

Home » date » 2009 » Oct » 28 »

*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: Wed, 28 Oct 2009 05:59:03 -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/28/t1256731199qdz36d14uiwd7kz.htm/, Retrieved Wed, 28 Oct 2009 12:59:59 +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/28/t1256731199qdz36d14uiwd7kz.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 «

456,2400684 452,7802999 447,5377079 445,0707809 441,7985966 437,4871427 462,3505164 475,1631299 480,1864221 478,8120717 465,0032258 465,5244354 469,2611213 469,1641504 467,3799311 463,5860222 459,026143 457,9006442 486,8685654 492,0172761 491,8109393 492,5657317 488,4106878 490,2417771 494,7241656 494,5462567 491,4997457 490,4497936 485,8899052 486,8233766 514,3724332 519,9509592 519,273531 516,7562288 508,0482261 512,652904 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,87 etc...

Dataseries Y:

» Textbox « » Textfile « » CSV «

508,8044811 505,0346523 497,6926763 494,7514528 490,8166664 489,8295622 519,8528638 532,8161034 531,3642818 523,0391955 504,7335931 504,335206 505,8092526 504,4700189 501,9910358 496,325498 493,9827932 495,9364879 528,9669177 532,0817606 530,1858165 525,3227579 520,2826155 521,4412719 523,8434881 521,9511471 518,0357131 518,252834 514,5561194 519,5902232 552,0362307 556,2058971 555,2900143 546,2847243 537,8020082 540,3471107 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,9 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	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

Model: Y[t] = c + b X[t] + e[t]
c	129.451547698360
b	0.800848986570733

Descriptive Statistics about e[t]
# observations	103
minimum	-37.4702359710975
Q1	-4.74706714872426
median	1.35252634136970
mean	6.70815883962682e-17
Q3	9.64625534774884
maximum	22.8306446654471

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/1glz61256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/1glz61256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/2hl6s1256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/2hl6s1256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/3up1y1256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/3up1y1256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/4kmib1256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/4kmib1256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/51un11256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/51un11256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/6f0fc1256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/6f0fc1256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/7aw9m1256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/7aw9m1256731137.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/8olnl1256731137.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Oct/28/t1256731199qdz36d14uiwd7kz/8olnl1256731137.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