Blog & Share Statistical Computations at FreeStatistics.org

Home » date » 2009 » Nov » 03 »

*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, 03 Nov 2009 15:53:47 -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/03/t1257288895gelbq1uiton6b7h.htm/, Retrieved Tue, 03 Nov 2009 23:55:04 +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/03/t1257288895gelbq1uiton6b7h.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 «

258881,9875 255059,9883 247697,9873 244778,9879 240900,9888 239932,9876 270246,9882 283892,9885 282347,986 273569,9859 254755,9873 254353,9863 255842,9886 254489,9882 251994,9864 246338,9879 244018,9884 245952,9859 279805,9879 283110,9877 281096,9855 275963,9864 270693,9865 271900,9855 274411,9877 272432,9881 268360,9865 268585,9875 264767,9875 269973,9865 304743,989 309364,9879 308346,9857 298426,9868 289230,9865 291974,9869 294911,9875 293487,9873 290554,9853 284735,9877 281817,9866 287853,986 316262,9885 325411,9873 326010,9856 328281,9854 317479,9858 317538,9865 313736,987 312275,9872 309390,9854 302949,9874 300315,9868 304034,9857 333475,9874 337697,9871 335931,9864 323930,9852 313926,9859 314484,9864 313217,9867 309663,9868 302962,9864 298988,9856 298422,9873 301630,986 329764,987 335082,9882 327615,9853 309118,9855 295915,9872 291412,9859 291458,4609 284594,461 276391,4593 272482,4605 264897,461 263206,4597 296722,4612

Dataseries Y:

» Textbox « » Textfile « » CSV «

208069,33 204924,3308 200204,3298 198002,3304 195100,3314 191309,3301 213682,3307 225694,331 230493,3285 229175,3283 216142,3298 216627,3287 220120,3312 220029,3308 218358,3289 214826,3304 210619,331 209587,3283 236955,3304 241995,3303 241792,3279 242535,3288 238459,329 240251,3279 244666,3303 244490,3306 241486,3289 240455,33 236003,33 236911,3289 264493,3316 270263,3304 269559,3281 266951,3293 258027,3289 262727,3294 267327,33 267280,3298 264691,3277 258777,3302 254758,3291 254782,3285 277181,3311 285265,3298 286516,328 282956,3278 276601,3282 277829,329 277042,3294 277017,3297 274951,3278 270064,3299 267054,3293 264907,3281 287173,3299 291100,3296 292214,3288 288100,3275 281391,3283 282570,3289 280104,3292 280322,3292 276750,3289 275130,328 274266,3298 271225,3284 289716,3295 290640,3307 292214,3276 278420,3278 269740,3297 265775,3283 268948,3304 264090,3305 255112,3287 253267,33 245971,3306 235286,3291 258470,3308

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	-32896.0583342823
b	0.991869352784941

Descriptive Statistics about e[t]
# observations	79
minimum	-22994.3654293957
Q1	-4538.09274498655
median	1308.28316866395
mean	-2.4891200212096e-13
Q3	6178.15265032349
maximum	16123.7157378378

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/1j6ru1257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/1j6ru1257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/2i0711257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/2i0711257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/35akc1257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/35akc1257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/4t1e51257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/4t1e51257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/5jj3n1257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/5jj3n1257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/668ex1257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/668ex1257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/7vprt1257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/7vprt1257288823.ps (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/851cf1257288823.png (open in new window)

http://www.freestatistics.org/blog/date/2009/Nov/03/t1257288895gelbq1uiton6b7h/851cf1257288823.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