Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_correlation.wasp
Title produced by softwarePearson Correlation
Date of computationMon, 05 Mar 2018 20:29:37 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2018/Mar/05/t1520278205dgav776uu1rg6ie.htm/, Retrieved Mon, 06 May 2024 03:40:32 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=315004, Retrieved Mon, 06 May 2024 03:40:32 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact137
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Pearson Correlation] [] [2018-03-05 19:29:37] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
length
144
123
169
121
130
165
97
131
142
136
161
102
146
112
91
98
104
143
130
130
169
115
143
165
88
136
93
106
143
142
133
117
103
130
154
106
91
102
136
131
100
146
148
95
124
103
93
124
98
108
128
162
130
158
153
94
150
100
96
106
89
108
152
92
98
92
122
150
126
124
101
87
138
144
111
169
93
139
117
108
144
108
151
154
117
104
149
126
91
187
143
157
115
116
140
120
86
140
118
115
125
115
90
127
142
124
93
163
127
201
154
129
147
143
109
100
101
138
134
179
133
161
95
106
106
88
145
142
121
81
178
116
152
92
119
92
183
90
130
131
143
127
104
108
130
104
130
123
82
155
92
128
122
107
88
95
124
136
125
136
123
95
169
119
127
151
85
88
139
120
139
194
119
91
124
98
105
129
126
86
139
103
121
145
95
115
91
136
110
113
104
81
130
127
135
140
81
121
100
131
107
123
141
101
141
142
89
116
91
117
195
125
141
125
130
154
105
127
112
117
138
90
129
120
128
126
100
118
127
94
189
142
84
99
128
84
137
143
113
118
111
181
103
127
127
124
113
118
93
119
145
108
83
93
114
126
93
108
128
127
133
107
85
104
131
116
104
104
108
92
121
102
102
119
117
118
119
100
111
110
120
94
119
80
97
137
93
113
96
110
154
161
106
102
91
116
117
114
96
112
105
107
126
105
106
105
118
106
96
107
117
132
91
99
93
131
116
114
134
95
116
96
124
114
127
115
114
113
99
119
94
109
110
97
127
106
127
95
27
115
107
110
116
138
111
109
100
88
104
124
133
124
105
94
132
122
153
117
126
117
95
119
127
182
114
143
99
109
99
101
110
116
111
137
118
98
125
123
175
96
121
93
99
96
139
134
119
114
111
132
138
103
125
100
113
88
136
109
Dataseries Y:
$952,624,490.27 
$750,140,046.17 
$671,220,455.10 
$772,158,880.00 
$706,988,165.89 
$1,088,559,994.68 
$756,677,675.77 
$746,171,720.76 
$707,732,953.69 
$712,903,691.09 
$973,402,222.02 
$1,291,514,391.87 
$876,652,703.76 
$653,534,988.15 
$726,116,532.72 
$983,938,347.10 
$753,652,179.78 
$677,113,113.86 
$799,384,854.30 
$1,231,937,550.46 
$1,045,906,537.07 
$853,264,813.24 
$1,140,066,020.14 
$1,115,258,567.81 
$902,175,843.80 
$779,470,849.49 
$768,148,576.60 
$641,761,449.87 
$1,561,752,970.70 
$710,892,844.47 
$729,242,277.12 
$747,568,480.64 
$591,768,996.92 
$1,408,187,497.58 
$1,179,649,301.10 
$587,678,383.43 
$698,778,760.90 
$615,927,877.92 
$1,097,688,249.77 
$657,258,197.48 
$640,817,168.50 
$1,039,829,831.30 
$893,915,146.90 
$588,103,656.97 
$756,351,938.78 
$1,151,241,371.27 
$814,943,703.06 
$675,617,796.27 
$535,872,792.75 
$1,110,413,235.80 
$576,744,932.84 
$3,025,614,789.35 
$781,234,707.77 
$846,695,292.79 
$1,028,416,946.70 
$975,885,765.83 
$920,434,198.70 
$514,511,873.56 
$804,914,384.79 
$642,755,355.79 
$662,286,859.38 
$668,787,232.35 
$1,093,392,302.94 
$684,753,923.91 
$571,713,515.70 
$692,823,108.65 
$862,689,853.26 
$460,239,821.32 
$641,750,240.94 
$520,842,241.26 
$666,589,609.88 
$600,223,218.85 
$1,070,332,714.40 
$808,210,413.44 
$710,271,758.12 
$1,097,133,563.54 
$909,845,580.02 
$1,014,515,715.39 
$519,365,892.24 
$672,845,512.89 
$701,617,723.36 
$450,143,548.74 
$1,248,736,559.96 
$458,044,146.70 
$541,086,305.35 
$538,013,485.95 
$888,069,619.93 
$465,971,959.24 
$774,110,496.19 
$665,579,149.09 
$900,723,586.96 
$1,083,174,496.16 
$574,240,396.06 
$715,424,296.77 
$450,608,889.05 
$578,156,074.83 
$644,014,476.99 
$1,026,150,164.19 
$445,035,858.01 
$645,787,263.70 
$453,418,828.21 
$789,282,558.61 
$459,082,108.05 
$979,683,031.61 
$994,517,314.65 
$680,323,214.86 
$1,149,769,230.81 
$621,746,515.85 
$764,854,839.34 
$1,437,468,137.47 
$586,265,893.86 
$548,509,850.95 
$350,840,740.50 
$839,770,411.82 
$556,246,757.96 
$1,113,971,831.87 
$621,991,654.79 
$952,547,456.42 
$523,311,824.52 
$1,215,888,673.56 
$567,025,727.92 
$1,153,349,012.27 
$484,030,881.28 
$535,885,534.32 
$535,885,534.32 
$579,952,267.88 
$470,417,261.03 
$852,430,975.14 
$1,078,613,527.65 
$503,081,452.61 
$1,160,945,184.82 
$600,919,721.86 
$1,299,565,844.68 
$700,446,276.21 
$482,919,671.58 
$491,677,898.59 
$598,924,902.95 
$645,839,752.37 
$577,315,986.69 
$468,894,560.29 
$589,107,501.38 
$512,978,409.67 
$473,254,604.73 
$453,101,727.42 
$399,592,809.86 
$406,337,937.14 
$450,735,586.05 
$749,202,165.88 
$479,673,643.72 
$627,512,195.02 
$687,402,755.40 
$512,818,154.00 
$504,972,371.74 
$953,555,515.37 
$635,213,710.15 
$442,215,537.18 
$515,733,894.84 
$1,310,017,631.92 
$589,494,477.14 
$652,487,687.22 
$419,100,786.51 
$526,412,918.83 
$697,986,463.76 
$535,808,551.59 
$413,490,148.56 
$802,094,487.32 
$426,544,779.85 
$440,833,092.36 
$549,033,595.61 
$506,228,509.20 
$462,201,855.32 
$2,711,618,942.05 
$489,908,039.66 
$379,335,283.22 
$463,641,222.04 
$867,079,779.93 
$440,297,359.03 
$910,108,547.07 
$388,265,227.30 
$445,331,636.78 
$411,669,626.96 
$482,605,706.88 
$465,755,910.42 
$1,230,107,590.71 
$412,284,269.25 
$364,630,884.66 
$489,602,830.84 
$504,236,144.09 
$688,787,804.91 
$744,130,795.25 
$407,161,572.42 
$560,796,160.01 
$545,667,343.27 
$507,116,462.17 
$409,363,135.75 
$547,646,277.68 
$536,194,290.82 
$521,397,145.36 
$446,097,882.15 
$567,215,110.80 
$393,970,405.08 
$356,352,655.71 
$343,962,158.23 
$560,158,721.17 
$603,652,739.56 
$1,079,245,007.52 
$1,224,609,826.48 
$558,349,969.97 
$544,302,423.76 
$391,461,772.95 
$525,027,843.81 
$337,721,152.76 
$319,076,111.52 
$721,079,965.93 
$602,758,308.48 
$441,597,053.06 
$372,234,469.27 
$1,494,643,267.07 
$416,680,933.02 
$435,658,442.89 
$409,362,990.66 
$848,295,321.88 
$691,707,059.46 
$604,173,235.83 
$363,287,833.02 
$449,050,801.17 
$389,781,727.82 
$541,460,630.25 
$593,961,676.43 
$308,145,094.86 
$356,095,013.09 
$521,567,841.06 
$609,997,249.69 
$331,992,373.90 
$316,024,944.41 
$313,951,873.33 
$901,210,651.82 
$676,967,922.64 
$310,376,502.92 
$472,832,496.33 
$364,851,165.52 
$766,363,955.48 
$861,165,331.11 
$399,800,210.00 
$913,588,289.03 
$433,633,526.28 
$472,122,382.99 
$441,756,758.76 
$364,866,028.19 
$837,176,928.80 
$306,577,096.21 
$632,094,375.18 
$350,666,921.09 
$565,542,732.63 
$433,876,252.58 
$783,284,649.89 
$424,108,135.23 
$410,151,533.40 
$449,121,409.11 
$902,911,069.85 
$344,948,173.18 
$223,421,641.89 
$157,192,136.58 
$156,127,457.44 
$165,680,195.48 
$255,783,543.19 
$312,268,614.98 
$229,470,092.55 
$218,167,988.16 
$147,115,232.32 
$257,574,173.26 
$167,612,331.26 
$348,734,793.08 
$325,590,003.39 
$136,502,617.91 
$132,539,476.89 
$158,529,078.59 
$319,394,671.09 
$139,249,933.14 
$135,533,817.08 
$298,446,153.99 
$171,957,403.89 
$236,362,728.24 
$182,134,789.87 
$376,591,505.78 
$185,621,091.22 
$154,306,521.94 
$247,976,609.65 
$196,603,537.43 
$380,853,328.77 
$206,660,859.39 
$191,070,314.62 
$166,716,691.55 
$131,860,128.89 
$280,607,936.32 
$462,165,115.18 
$187,212,356.19 
$134,714,403.83 
$330,073,804.72 
$150,771,621.56 
$533,507,308.04 
$178,168,716.07 
$206,383,818.47 
$542,305,856.47 
$347,890,557.86 
$174,015,016.98 
$408,766,010.76 
$178,103,266.34 
$184,529,371.05 
$181,885,265.49 
$160,358,685.02 
$257,036,604.41 
$153,581,249.47 
$150,636,755.46 
$153,839,280.60 
$160,731,608.43 
$214,320,190.47 
$188,628,895.22 
$598,793,082.96 
$220,285,758.34 
$433,578,101.55 
$192,977,753.79 
$317,587,527.62 
$170,547,987.84 
$139,613,624.14 
$878,893,732.43 
$187,442,468.68 
$193,086,527.57 
$303,766,229.77 
$153,186,172.93 
$258,122,068.33 
$309,852,111.31 
$110,047,509.57 
$247,968,282.64 
$142,380,021.95 
$221,564,822.26 
$281,019,809.85 
$187,491,419.29 
$359,365,648.51 
$551,261,301.66 
$130,989,277.73 
$295,981,184.10 
$202,556,128.37 
$534,814,499.80 
$290,277,362.86 
$200,149,101.49 
$189,504,484.45 
$239,137,939.02 
$168,644,860.87 
$600,034,614.87 
$163,993,499.96 
$192,513,743.73 
$345,592,354.55 
$239,668,652.38 
$267,531,465.18 
$243,487,203.97 
$271,476,717.71 
$281,106,889.67 
$228,664,986.61 
$263,216,748.21 
$212,051,774.10 
$277,040,309.05 
$213,741,861.13 
$177,376,532.95 
$308,535,411.68 
$486,064,788.22 
$282,473,650.31 
$434,904,844.66 
$179,549,193.51 
$295,656,760.75 
$579,356,607.45 
$281,522,192.14 
$397,428,019.80 
$453,517,815.51 
$367,087,834.87 
$198,952,451.13 
$182,500,010.40 
$184,478,157.71 
$197,739,921.62 
$493,813,774.95 
$1,197,208,492.59 
$149,041,059.73 
$231,139,337.43 
$191,866,398.66 
$331,977,504.39 
$218,335,151.83 
$486,492,894.10 
$211,963,752.42 
$252,813,235.76 
$179,350,846.75 
$292,970,147.63 
$152,174,306.67 
$219,441,449.81 
$161,741,515.00 
$183,599,406.33 
$153,609,014.87 
$171,164,330.86 
$216,842,095.62 




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time1 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=315004&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]1 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=315004&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=315004&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center



Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
library(psychometric)
x <- x[!is.na(y)]
y <- y[!is.na(y)]
y <- y[!is.na(x)]
x <- x[!is.na(x)]
bitmap(file='test1.png')
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, sub=main)
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()
lx = length(x)
makebiased = (lx-1)/lx
varx = var(x)*makebiased
vary = var(y)*makebiased
corxy <- cor.test(x,y,method='pearson', na.rm = T)
cxy <- as.matrix(corxy$estimate)[1,1]
load(file='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,'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,'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,'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,'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,'Correlation',header=TRUE)
a<-table.element(a,cxy,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Determination',header=TRUE)
a<-table.element(a,cxy*cxy,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'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,'95% CI of Correlation',header=TRUE)
a<-table.element(a,paste('[',CIr(r=cxy, n = lx, level = .95)[1],', ', CIr(r=cxy, n = lx, level = .95)[2],']',sep=''),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='mytable.tab')
library(moments)
library(nortest)
jarque.x <- jarque.test(x)
jarque.y <- jarque.test(y)
if(lx>7) {
ad.x <- ad.test(x)
ad.y <- ad.test(y)
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Normality Tests',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('jarque.x'),'
',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('jarque.y'),'
',sep=''))
a<-table.row.end(a)
if(lx>7) {
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('ad.x'),'
',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('ad.y'),'
',sep=''))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
library(car)
bitmap(file='test2.png')
qqPlot(x,main='QQplot of variable x')
dev.off()
bitmap(file='test3.png')
qqPlot(y,main='QQplot of variable y')
dev.off()