FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_QRW1.wasp
Title produced by softwareQRW Figure 1
Date of computationThu, 18 Oct 2007 05:46:35 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2007/Oct/18/26qptxvtt3arzpz1192711431.htm/, Retrieved Mon, 29 Apr 2024 00:10:24 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=889, Retrieved Mon, 29 Apr 2024 00:10:24 +0000
QR Codes:

Original text written by user:archived for reference purposes
IsPrivate?No (this computation is public)
User-defined keywordsquasi random walk, descriptive statistics, figure 1, article
Estimated Impact408

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [QRW Figure 1] [Quasi Random Walk...] [2007-10-18 12:46:35] [8326296a4a969e220946fc6855491489] [Current]
- RMP     [] [Alternatieve opdr...] [-0001-11-30 00:00:00] [74be16979710d4c4e7c6647856088456]
- RMP     [] [Alternatieve opdr...] [-0001-11-30 00:00:00] [74be16979710d4c4e7c6647856088456]
Feedback Forum

Post a new message

Dataset

Tables (Output of Computation)





Summary of compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time42 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001

\begin{tabular}{lllllllll}
\hline
Summary of compuational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 42 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=889&T=0

[TABLE]
[ROW][C]Summary of compuational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]42 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=889&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=889&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 compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time42 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 1995 ; par2 = 2006 ; par3 = 20 ;
Parameters (R input):
par1 = 1995 ; par2 = 2006 ; par3 = 20 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
if (par3 > 50) par3 = 50
if (par3 < 5) par3 = 5
load('Q2.RData')
maxnumsubsets <- 5
startyear <- par1
endyear <- par2
numseries <- 66
numsim <- par3
simarr <- array(NA,dim=c(numseries,length(tsarr[1,]),numsim))
lenser <- array(NA,dim=c(numseries,2,numsim))
minser <- array(NA,dim=c(numseries,2,numsim))
maxser <- array(NA,dim=c(numseries,2,numsim))
medser <- array(NA,dim=c(numseries,2,numsim))
q1ser <- array(NA,dim=c(numseries,2,numsim))
q3ser <- array(NA,dim=c(numseries,2,numsim))
for (i in 1:numseries)
{
dum <- tsarr[i,!is.na(tsarr[i,])]
lenser[i,1,1] <- length(dum)
minser[i,1,1] <- min(dum)
maxser[i,1,1] <- max(dum)
medser[i,1,1] <- median(dum)
q1ser[i,1,1] <- quantile(dum,0.25)
q3ser[i,1,1] <- quantile(dum,0.75)
omean <- mean(tsarr[i,],na.rm=TRUE)
osd <- sd(tsarr[i,],na.rm=TRUE)
for (jj in 1:numsim)
{
dum <- rnorm(lenser[i,1,1],omean,osd)
for (j in 1:lenser[i,1,1])
{
simarr[i,j,jj] <- dum[j]
}
dum <- simarr[i,!is.na(simarr[i,,jj]),jj]
lenser[i,2,jj] <- length(dum)
minser[i,2,jj] <- min(dum)
maxser[i,2,jj] <- max(dum)
medser[i,2,jj] <- median(dum)
q1ser[i,2,jj] <- quantile(dum,0.25)
q3ser[i,2,jj] <- quantile(dum,0.75)
}
}
iqr1<-(rep(q3ser[,1,1],numsim)-rep(q1ser[,1,1],numsim))
iqr2<-(c(q3ser[,2,])-c(q1ser[,2,]))
bitmap(file='figure1.png')
op <- par(mfrow=c(3,3))
hist(lenser[,1,1],main='',xlab='length of QRW')
hist(medser[,1,1],main='',xlab='median of QRW')
hist(maxser[,1,1]-minser[,1,1],main='',xlab='range of QRW')
hist(minser[,1,1],main='',xlab='min of QRW')
hist(maxser[,1,1],main='',xlab='max of QRW')
hist(q3ser[,1,1]-q1ser[,1,1],main='',xlab='IQR of QRW')
plot(rep(minser[,1,1],numsim),c(minser[,2,]),xlab='min of QRW',ylab='min of RW')
plot(rep(maxser[,1,1],numsim),c(maxser[,2,]),xlab='max of QRW',ylab='max of RW')
plot(iqr1,iqr2,xlab='IQR of QRW',ylab='IQR of RW')
par(op)
dev.off()