Blog & Share Statistical Computations at FreeStatistics.org

Home » date » 2009 » Nov » 07 »

*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: Sat, 07 Nov 2009 06:36:29 -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/Nov/07/t1257601026zcn7cn0hgfto6wp.htm/, Retrieved Sat, 07 Nov 2009 14:37:07 +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/Nov/07/t1257601026zcn7cn0hgfto6wp.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 «

-4757,94575 -3501,961939 4675,799026 -7958,854315 -15749,64956 -6124,716213 -184,4590898 14345,88187 27460,89236 21294,92474 10550,92664 10201,63045 4614,163773 5480,679012 9449,551401 3259,567589 -2973,602891 5714,73617 5349,964714 19426,52282 32541,6895 23421,48854 18701,25331 15422,67711 10693,59997 12077,42759 18654,9857 6885,427588 12330,61236 12512,78284 10978,979 22258,2714 33638,93903 21010,76665 9861,955217 -1578,057173 4338,06664 -999,9800276 -719,9371594 -11274,35526 -16375,85621 -21643,18668 -15397,97623 -4937,339076 -5575,357165 -2806,125728 -16267,7324 -21705,61908 -12931,76288 -14993,62288 -25009,96574 -19861,60859 -37016,01241 -31067,70192 -23619,74289 -20614,77527 -10320,49907 -12945,85811 -22959,08765 -16954,57051 -10663,64766 -5126,929561 6361,301875 1106,273297 -7004,022896 5033,989494 11966,79233

Dataseries Y:

» Textbox « » Textfile « » CSV «

-1814,938761 -2204,024389 3414,636476 -7742,342918 -13656,14768 5,537168254 15341,93457 29611,03726 39633,20496 39834,37622 25546,34886 24352,55777 16737,71535 15984,81466 17893,10228 10417,1879 5223,136563 17548,26661 29014,43804 41297,8763 49063,98236 33086,73241 24662,33855 21679,84202 18615,35916 15987,33518 17347,77345 5312,335178 9757,499506 15525,55085 27754,55102 41297,39682 41619,4892 24048,46522 8884,57483 -1555,565515 752,2015842 -6493,082658 -7342,743698 -18169,32924 -24011,31556 -21399,25054 -1337,137376 11827,75639 -255,3018794 -5429,853309 -20389,54846 -26904,94906 -19848,74707 -24563,89435 -33177,96967 -31476,78136 -45298,25391 -34656,34985 -15791,71616 -8909,88742 -13931,12725 -24430,2882 -38781,76413 -36786,69218 -35057,17504 -33470,85313 -29093,40456 -34699,63054 -42357,42162 -28366,28127 -9675,585947

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	'RServer@AstonUniversity' @ vre.aston.ac.uk

Model: Y[t] = c + b X[t] + e[t]
c	-2.1950784821643e-07
b	1.34058739846871

Descriptive Statistics about e[t]
# observations	67
minimum	-37621.2856912609
Q1	-4399.93525870691
median	1418.76859732061
mean	2.22442296754748e-13
Q3	8426.86999224934
maximum	21842.3427623789

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/1e7231257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/1e7231257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/2gpak1257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/2gpak1257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/39ns31257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/39ns31257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/4ybzg1257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/4ybzg1257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/5si641257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/5si641257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/6p2491257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/6p2491257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/7xdb01257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/7xdb01257600982.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/88uoh1257600982.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/07/t1257601026zcn7cn0hgfto6wp/88uoh1257600982.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