R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > y <- c(1081 + ,309 + ,458 + ,588 + ,302 + ,156 + ,481 + ,353 + ,452 + ,109 + ,115 + ,110 + ,239 + ,247 + ,505 + ,159 + ,109 + ,519 + ,248 + ,373 + ,119 + ,84 + ,102 + ,295 + ,105 + ,64 + ,282 + ,182 + ,37 + ,361 + ,28 + ,85 + ,45 + ,49 + ,22 + ,155 + ,91 + ,81 + ,79 + ,145 + ,855 + ,61 + ,226 + ,105 + ,62 + ,25 + ,217 + ,322 + ,84 + ,33 + ,108 + ,150 + ,115 + ,162 + ,158 + ,97 + ,9 + ,66 + ,107 + ,101 + ,47 + ,38 + ,34 + ,87 + ,79 + ,947 + ,74 + ,53 + ,94 + ,63 + ,58 + ,49 + ,34 + ,11 + ,35 + ,20 + ,47 + ,43 + ,117 + ,171 + ,26 + ,75 + ,59 + ,18 + ,15 + ,72 + ,86 + ,14 + ,64 + ,11 + ,52 + ,41 + ,99 + ,75 + ,45 + ,43 + ,8 + ,198 + ,22 + ,11 + ,33 + ,23 + ,80 + ,18 + ,40 + ,23 + ,60 + ,20 + ,61 + ,36 + ,30 + ,47 + ,71 + ,14 + ,9 + ,39 + ,26 + ,21 + ,16 + ,69 + ,92 + ,14 + ,107 + ,29 + ,37 + ,23 + ,0 + ,7 + ,28 + ,0 + ,8 + ,63 + ,0 + ,3 + ,5 + ,9 + ,13 + ,2 + ,5 + ,0 + ,14 + ,0 + ,15 + ,3 + ,15 + ,11 + ,0 + ,6 + ,2 + ,1 + ,10 + ,73 + ,0 + ,11 + ,3 + ,0 + ,0 + ,2 + ,7 + ,0 + ,0 + ,0 + ,0 + ,27 + ,51 + ,3 + ,0 + ,19 + ,393 + ,0 + ,0 + ,4 + ,0 + ,9 + ,0 + ,10 + ,152 + ,1 + ,0 + ,0 + ,34 + ,10 + ,57 + ,52 + ,5 + ,14 + ,29 + ,5 + ,5 + ,0 + ,4 + ,0 + ,6 + ,0 + ,2 + ,0 + ,91 + ,0 + ,0 + ,20 + ,0 + ,0 + ,0 + ,27 + ,17 + ,2 + ,4 + ,0 + ,32 + ,31 + ,0 + ,0 + ,32 + ,20 + ,7 + ,0 + ,8 + ,28 + ,0 + ,29 + ,4 + ,0 + ,2 + ,21 + ,2 + ,26 + ,14 + ,0 + ,4 + ,0 + ,9 + ,10 + ,0 + ,17 + ,0 + ,1 + ,0 + ,6 + ,0 + ,0 + ,0 + ,3 + ,0 + ,8 + ,4 + ,0 + ,0 + ,11 + ,0 + ,0 + ,9 + ,0 + ,2 + ,73 + ,94 + ,0 + ,8 + ,35 + ,0 + ,0 + ,0 + ,0 + ,12 + ,15 + ,0 + ,0 + ,0 + ,11 + ,0 + ,6 + ,0 + ,0 + ,0 + ,12 + ,0 + ,30 + ,33 + ,0 + ,117 + ,28 + ,0 + ,0 + ,72 + ,0 + ,13 + ,0 + ,6 + ,4 + ,0 + ,62 + ,0 + ,24 + ,21 + ,0 + ,14 + ,21 + ,0 + ,0 + ,0 + ,0 + ,4 + ,2 + ,0 + ,0 + ,53 + ,9 + ,0 + ,13 + ,22 + ,0 + ,0 + ,83 + ,0 + ,8 + ,4 + ,14 + ,1 + ,17 + ,6 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,2 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,5 + ,2 + ,5 + ,78 + ,0 + ,1 + ,13 + ,15 + ,0 + ,0 + ,0 + ,48 + ,0 + ,6 + ,0 + ,0 + ,17 + ,14 + ,10 + ,12 + ,2 + ,0 + ,0 + ,52 + ,0 + ,0 + ,0 + ,4 + ,0 + ,0 + ,24 + ,11 + ,0 + ,0 + ,0 + ,0 + ,21 + ,0 + ,40 + ,9 + ,1 + ,0 + ,0 + ,24 + ,11 + ,14 + ,0 + ,0 + ,60 + ,80 + ,0 + ,16 + ,40 + ,6 + ,8 + ,3 + ,16 + ,10 + ,8 + ,7 + ,8 + ,12 + ,13 + ,42 + ,118 + ,9 + ,138 + ,5 + ,9 + ,8 + ,25 + ,7 + ,13 + ,16 + ,11 + ,11 + ,3 + ,61 + ,29 + ,17 + ,33 + ,15 + ,3 + ,66 + ,17 + ,26 + ,3 + ,2 + ,67 + ,70 + ,26 + ,24 + ,97 + ,30 + ,223 + ,48 + ,90 + ,180) > x <- c(213118 + ,81767 + ,153198 + ,-26007 + ,126942 + ,157214 + ,129352 + ,234817 + ,60448 + ,47818 + ,245546 + ,48020 + ,-1710 + ,32648 + ,95350 + ,151352 + ,288170 + ,114337 + ,37884 + ,122844 + ,82340 + ,79801 + ,165548 + ,116384 + ,134028 + ,63838 + ,74996 + ,31080 + ,32168 + ,49857 + ,87161 + ,106113 + ,80570 + ,102129 + ,301670 + ,102313 + ,88577 + ,112477 + ,191778 + ,79804 + ,128294 + ,96448 + ,93811 + ,117520 + ,69159 + ,101792 + ,210568 + ,136996 + ,121920 + ,76403 + ,108094 + ,134759 + ,188873 + ,146216 + ,156608 + ,61348 + ,50350 + ,87720 + ,99489 + ,87419 + ,94355 + ,60326 + ,94670 + ,82425 + ,59017 + ,90829 + ,80791 + ,100423 + ,131116 + ,100269 + ,27330 + ,39039 + ,106885 + ,79285 + ,118881 + ,77623 + ,114768 + ,74015 + ,69465 + ,117869 + ,60982 + ,90131 + ,138971 + ,39625 + ,102725 + ,64239 + ,90262 + ,103960 + ,106611 + ,103345 + ,95551 + ,82903 + ,63593 + ,126910 + ,37527 + ,60247 + ,112995 + ,70184 + ,130140 + ,73221 + ,76114 + ,90534 + ,108479 + ,113761 + ,68696 + ,71561 + ,59831 + ,97890 + ,101481 + ,72954 + ,67939 + ,48022 + ,86111 + ,74020 + ,57530 + ,56364 + ,84990 + ,88590 + ,77200 + ,61262 + ,110309 + ,67000 + ,93099 + ,107577 + ,62920 + ,75832 + ,60720 + ,60793 + ,57935 + ,60720 + ,60630 + ,55637 + ,60720 + ,60887 + ,60720 + ,60505 + ,60945 + ,60720 + ,60720 + ,60720 + ,58990 + ,60720 + ,56750 + ,60894 + ,63346 + ,56535 + ,60720 + ,60835 + ,60720 + ,61016 + ,58650 + ,60438 + ,60720 + ,58625 + ,60938 + ,60720 + ,60720 + ,61490 + ,60845 + ,60720 + ,60720 + ,60720 + ,60720 + ,60830 + ,63261 + ,60720 + ,60720 + ,45689 + ,60720 + ,60720 + ,60720 + ,61564 + ,60720 + ,61938 + ,60720 + ,60951 + ,60720 + ,60745 + ,60720 + ,60720 + ,71642 + ,71641 + ,55792 + ,71873 + ,62555 + ,60370 + ,64873 + ,62041 + ,65745 + ,60720 + ,59500 + ,60720 + ,61630 + ,60720 + ,60890 + ,60720 + ,113521 + ,60720 + ,60720 + ,80045 + ,60720 + ,60720 + ,60720 + ,50804 + ,87390 + ,61656 + ,65688 + ,60720 + ,48522 + ,60720 + ,60720 + ,60720 + ,57640 + ,61977 + ,62620 + ,60720 + ,60831 + ,60646 + ,60720 + ,56225 + ,60510 + ,60720 + ,60698 + ,60720 + ,60805 + ,61404 + ,60720 + ,60720 + ,65276 + ,60720 + ,63915 + ,60720 + ,60720 + ,61686 + ,60720 + ,60743 + ,60720 + ,60349 + ,60720 + ,60720 + ,60720 + ,61360 + ,60720 + ,59818 + ,72680 + ,60720 + ,60720 + ,61808 + ,60720 + ,60720 + ,53110 + ,60720 + ,64245 + ,73007 + ,82732 + ,60720 + ,54820 + ,47705 + ,60720 + ,60720 + ,60720 + ,60720 + ,72835 + ,58856 + ,60720 + ,60720 + ,60720 + ,77655 + ,60720 + ,69817 + ,60720 + ,60720 + ,60720 + ,60798 + ,60720 + ,62452 + ,64175 + ,60720 + ,67440 + ,68136 + ,60720 + ,60720 + ,56726 + ,60720 + ,70811 + ,60720 + ,60720 + ,62045 + ,60720 + ,54323 + ,60720 + ,62841 + ,81125 + ,60720 + ,59506 + ,59365 + ,60720 + ,60720 + ,60720 + ,60720 + ,60798 + ,58790 + ,60720 + ,60720 + ,61808 + ,60735 + ,60720 + ,64016 + ,54683 + ,60720 + ,60720 + ,87192 + ,60720 + ,64107 + ,60761 + ,65990 + ,59988 + ,61167 + ,60719 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,60722 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,60379 + ,60727 + ,60720 + ,60925 + ,60720 + ,60896 + ,59734 + ,62969 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,59118 + ,60720 + ,60720 + ,60720 + ,58598 + ,61124 + ,59595 + ,62065 + ,60720 + ,60720 + ,78780 + ,60720 + ,60720 + ,60720 + ,60722 + ,60720 + ,60720 + ,61600 + ,59635 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,60720 + ,59781 + ,76644 + ,64820 + ,60720 + ,60720 + ,56178 + ,60436 + ,60720 + ,60720 + ,60720 + ,73433 + ,41477 + ,60720 + ,62700 + ,67804 + ,59661 + ,58620 + ,60398 + ,58580 + ,62710 + ,59325 + ,60950 + ,68060 + ,83620 + ,58456 + ,52811 + ,121173 + ,63870 + ,21001 + ,70415 + ,64230 + ,59190 + ,69351 + ,64270 + ,70694 + ,68005 + ,58930 + ,58320 + ,69980 + ,69863 + ,63255 + ,57320 + ,75230 + ,79420 + ,73490 + ,35250 + ,62285 + ,69206 + ,65920 + ,69770 + ,72683 + ,-14545 + ,55830 + ,55174 + ,67038 + ,51252 + ,157278 + ,79510 + ,77440 + ,27284) > ylab = 'trades' > xlab = 'dividens' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > postscript(file="/var/www/rcomp/tmp/198ph1292774902.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > histx <- hist(x, plot=FALSE) > histy <- hist(y, plot=FALSE) > maxcounts <- max(c(histx$counts, histx$counts)) > xrange <- c(min(x),max(x)) > yrange <- c(min(y),max(y)) > nf <- layout(matrix(c(2,0,1,3),2,2,byrow=TRUE), c(3,1), c(1,3), TRUE) > par(mar=c(4,4,1,1)) > plot(x, y, xlim=xrange, ylim=yrange, xlab=xlab, ylab=ylab) > par(mar=c(0,4,1,1)) > barplot(histx$counts, axes=FALSE, ylim=c(0, maxcounts), space=0) > par(mar=c(4,0,1,1)) > barplot(histy$counts, axes=FALSE, xlim=c(0, maxcounts), space=0, horiz=TRUE) > dev.off() null device 1 > lx = length(x) > makebiased = (lx-1)/lx > varx = var(x)*makebiased > vary = var(y)*makebiased > corxy <- cor.test(x,y,method='pearson') > cxy <- as.matrix(corxy$estimate)[1,1] > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Pearson Product Moment Correlation - Ungrouped Data',3,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Statistic',1,TRUE) > a<-table.element(a,'Variable X',1,TRUE) > a<-table.element(a,'Variable Y',1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm','Mean',''),header=TRUE) > a<-table.element(a,mean(x)) > a<-table.element(a,mean(y)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/biased.htm','Biased Variance',''),header=TRUE) > a<-table.element(a,varx) > a<-table.element(a,vary) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/biased1.htm','Biased Standard Deviation',''),header=TRUE) > a<-table.element(a,sqrt(varx)) > a<-table.element(a,sqrt(vary)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/covariance.htm','Covariance',''),header=TRUE) > a<-table.element(a,cov(x,y),2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/pearson_correlation.htm','Correlation',''),header=TRUE) > a<-table.element(a,cxy,2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/coeff_of_determination.htm','Determination',''),header=TRUE) > a<-table.element(a,cxy*cxy,2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ttest_statistic.htm','T-Test',''),header=TRUE) > a<-table.element(a,as.matrix(corxy$statistic)[1,1],2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'p-value (2 sided)',header=TRUE) > a<-table.element(a,(p2 <- as.matrix(corxy$p.value)[1,1]),2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'p-value (1 sided)',header=TRUE) > a<-table.element(a,p2/2,2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Degrees of Freedom',header=TRUE) > a<-table.element(a,lx-2,2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Number of Observations',header=TRUE) > a<-table.element(a,lx,2) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/2292z1292774902.tab") > > try(system("convert tmp/198ph1292774902.ps tmp/198ph1292774902.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.580 0.190 0.771