Home » date » 2009 » Nov » 04 »

*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, 04 Nov 2009 12:12:31 -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/04/t1257362010js645ntqbpeniap.htm/, Retrieved Wed, 04 Nov 2009 20:13:32 +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/04/t1257362010js645ntqbpeniap.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:
Rob_WS5
 
Dataseries X:
» Textbox « » Textfile « » CSV «
103784.5647 106789.1029 113889.3931 120709.6235 112224.8795 115336.2211 123876.6481 128645.1607 117955.5279 121200.0661 127221.878 130871.7754 120395.7155 122826.724 127608.0832 130425.8272 118864.802 119456.0505 125726.0767 128969.7266 117441.3338 120415.4962 127525.8041 129939.5224 121157.9842 122346.3257 131061.1889 133972.0435 124923.4364 126359.2311 134604.5898 137128.9057 126495.632 129016.3242 136825.2729 140536.7256 131359.965 134527.7854 142275.3921 146687.7849 136205.9041 140667.6305 147443.0576 152116.5872 142296.2281 146169.9972 153908.1597 158496.0314 150180.0485 155418.8604 160725.5782 167376.3729 159676.2361 163171.9967 168006.6033 172087.5947 163770.3981 167646.3382 173556.5007 179462.6033 169624.6717 175262.6119 179659.9197 185496.4408 174269.7487 180679.4066 185375.475 192056.3211 180441.979 186922.0218 191210.7995 197904.9619 185547.9192 191137.9534 194532.0646 202546.723 189157.7658 194107.1591 197191.1591 20476 etc...
 
Dataseries Y:
» Textbox « » Textfile « » CSV «
11079.57694 10688.9842 10086.36942 9433.814847 8802.448758 8245.710955 7760.601438 7348.820207 6813.17089 6424.178145 6265.495915 6180.641972 6131.902172 6144.765 6451.903802 7085.769891 7914.432353 8977.60533 9849.744763 10403.01974 10919.38057 11060.27431 11321.63614 11932.31375 12323.60649 11766.56869 11650.76343 11527.79534 11012.33451 10840.52663 10295.91448 9715.788191 9484.677675 9225.915845 9182.738874 8718.117843 8922.277044 8771.657643 8437.867527 7923.106695 7175.935349 6945.524833 6460.415316 6062.259743 5757.720941 5564.624568 5703.118312 5594.875883 5608.70154 5450.719311 5211.145967 5038.03808 4922.532822 4722.273621 4388.483505 4382.320676 4215.075618 4166.335818 4304.429562 4391.383505 4449.586134 4400.846334 4662.908164 4257.789761 4520.251591 4231.838447 4290.441076 4160.810161 3871.197016 3910.011159 3722.277616 3862.071359 3826.257216 3857.508531 3952.625303 4512.963106 4549.777249 4887.267365 4888.967365 4862. etc...
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time5 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk


Model: Y[t] = c + b X[t] + e[t]
c16328.3126027675
b-0.0611443483368808


Descriptive Statistics about e[t]
# observations90
minimum-2834.8928639675
Q1-1213.72214071740
median-572.548849863763
mean1.53852239590277e-14
Q3952.756486111555
maximum5473.35624961697
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/10zlj1257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/10zlj1257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/25vb61257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/25vb61257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/3zmp81257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/3zmp81257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/4foxg1257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/4foxg1257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/5b9y61257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/5b9y61257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/66zrn1257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/66zrn1257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/7hk691257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/7hk691257361944.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/88tqj1257361944.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Nov/04/t1257362010js645ntqbpeniap/88tqj1257361944.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')
 




Copyright

Creative Commons License

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.


Privacy Policy

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