FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_qqplot.wasp
Title produced by softwareQQPlot and Histogram
Date of computationMon, 28 Oct 2024 18:33:42 +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/2024/Oct/28/t1730136859dr6nckg6la9h0ha.htm/, Retrieved Thu, 01 Oct 2026 07:49:44 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=320148, Retrieved Thu, 01 Oct 2026 07:49:44 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact346

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [QQPlot and Histogram] [] [2024-10-28 17:33:42] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
Length(cm):
21.0
21.1
22.2
22.7
22.9
23.2
23.7
23.9
24.0
24.8
28.5
29.1
Dataseries Y:
Download CSV
Histogram 
-1.41
-1.03
-0.812
-0.549
-0.319
-0.105
0.105
0.319
0.549
0.812
1.03
1.41

Tables (Output of Computation)





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=320148&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=320148&T=0

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

As an alternative you can also use a QR Code:  
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


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 0 ;
Parameters (R input):
par1 = 0 ;
R code (references can be found in the software module):
par1 <- '0'
library(car)
plot.qqline <- function(x,y,a=0.25, ...) {
y <- quantile(y[!is.na(y)],c(a, 1-a))
x <- quantile(x[!is.na(x)],c(a, 1-a))
points(x,y,...)
slope <- diff(y)/diff(x)
int <- y[1]-slope*x[1]
abline(int, slope, ...)
}
if(par1 == '') par1 <- 0
par1 <- as.numeric(par1)
if (par1 < 0) par1 <- 0
bitmap(file='test1.png')
if(par1 > 0) {
myhist<-hist(x, breaks=par1, col=2,main='Histogram (series X)')
} else {
myhist <- hist(x, col=2)
}
dev.off()
bitmap(file='test1a.png')
if(par1 > 0) {
myhist<-hist(y, breaks=par1, col=2,main='Histogram (series Y)')
} else {
myhist <- hist(y, col=2)
}
dev.off()
bitmap(file='test2.png')
qqPlot(x,main='Normal Q-Q Plot (series X)',ylab='Sample Quantiles',xlab='Normal Quantiles')
dev.off()
bitmap(file='test2a.png')
qqPlot(y,main='Normal Q-Q Plot (series Y)',ylab='Sample Quantiles',xlab='Normal Quantiles')
dev.off()
bitmap(file='test3.png')
qqplot(x,y,main='Q-Q Plot (series X vs Y)',ylab='Sample Quantiles of series Y',xlab='Sample Quantiles of series X')
plot.qqline(x,y)
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Descriptive Statistics',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean of X',header=TRUE)
a<-table.element(a,mean(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'standard deviation of X',header=TRUE)
a<-table.element(a,sd(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean of Y',header=TRUE)
a<-table.element(a,mean(y))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'standard deviation of Y',header=TRUE)
a<-table.element(a,sd(y))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')