Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_edauni.wasp
Title produced by softwareUnivariate Explorative Data Analysis
Date of computationWed, 21 Oct 2009 11:34:32 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Oct/21/t1256146522vcj8g7t2tasq2h5.htm/, Retrieved Thu, 02 May 2024 00:03:28 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=49527, Retrieved Thu, 02 May 2024 00:03:28 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact115
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Data Series] [Bivariate dataset] [2008-01-05 23:51:08] [74be16979710d4c4e7c6647856088456]
F RMPD  [Univariate Explorative Data Analysis] [Colombia Coffee] [2008-01-07 14:21:11] [74be16979710d4c4e7c6647856088456]
F RMPD    [Univariate Data Series] [] [2009-10-14 08:30:28] [74be16979710d4c4e7c6647856088456]
F RMP         [Univariate Explorative Data Analysis] [] [2009-10-21 17:34:32] [3e9f70e60513fc8919624add68d96eca] [Current]
Feedback Forum
2009-10-23 18:14:28 [Joris Van Mol] [reply
Beste medestudent,

Om één of andere reden is iets misgegaan op het forum van Moodel, daardoor zal het komen dat jij een review ziet die totaal niet over jouw workshop gaat. blijkbaar als je een vorige review opent om daar iets uit te kopieëren flipt het systeem en beetje en ben je één van de twee kwijt...
Ik zou dan via deze weg even willen zeggen wat ik van de review vond:

In grote lijnen was deze WS zeer goed, zeker bij de beteren die ik gereviewd heb, maar ik heb toch enkele kleine bemerkingen.
Part 1 is eigenlijk volledig in orde voor mij. Het begint enkel bij part 2 dat ik wat kwijt wil.
Jouw reeksen zijn de werkloosheid van de mannen en van de vrouwen, beide uitgedrukt in percentages. Het model dat je gebruikt (Y(t) + X(t))/2 is zeker en vast interessant. Wat ik ook zou kunnen voorstellen is bijvoorbeeld Y(t) / (Y(t) + X(t)) of X(t) / (Y(t) + X(t) omdat je dan de mannen, respectievelijk de vrouwen vergelijkt met het totaal. Dat is misschien ook een interessante die je kon proberen.
Verder ook een puntje van aandacht bij één van de laatste dingen van deze WS: op een bepaald moment lees ik het volgende ''Ik vermoed dat de volgende verhouding tussen het werkloosheidspercentage bij vrouwen en mannen ronde de 0,815% zal schommelen'' Hier wel opletten voor de plaatsing van de komma, en dat je niet vergelijkt met het totaal aantal maar enkel kijkt naar de verhouding tussen de 2 deelgroepen. Je deelt hier twee deelgroepen door elkaar, al wat we uit deze breuk kunnen afleiden is dat het percentage werkloze vrouwen groter is dan dat van de mannen. En dat de waarde van de mannen slechts 81,5% bedraagt in vergelijking met de waarde van de vrouwen. En hierdoor wordt het dus duidelijk dat je de komma verkeerd geplaatst hebt en dat de breuk ook niet met het totaal # werklozen vergelijkt. Jou gegevens leenden zich er misschien toe om wel met het totaal te vergelijken.....
Voor de rest was de WS heel goed, volledig en correct, blijkbaar wel wat tijd in gestoken en je snapt de dingen wel goed blijkbaar.

Post a new message
Dataseries X:
87.28
87.28
87.09
86.92
87.59
90.72
90.69
90.3
89.55
88.94
88.41
87.82
87.07
86.82
86.4
86.02
85.66
85.32
85
84.67
83.94
82.83
81.95
81.19
80.48
78.86
69.47
68.77
70.06
73.95
75.8
77.79
81.57
83.07
84.34
85.1
85.25
84.26
83.63
86.44
85.3
84.1
83.36
82.48
81.58
80.47
79.34
82.13
81.69
80.7
79.88
79.16
78.38
77.42
76.47
75.46
74.48
78.27
80.7
79.91
78.75
77.78
81.14
81.08
80.03
78.91
78.01
76.9
75.97
81.93
80.27
78.67
77.42
76.16
74.7
76.39
76.04
74.65
73.29
71.79
74.39
74.91
74.54
73.08
72.75
71.32
70.38
70.35
70.01
69.36
67.77
69.26
69.8
68.38
67.62
68.39
66.95
65.21
66.64
63.45
60.66
62.34
60.32
58.64
60.46
58.59
61.87
61.85
67.44
77.06
91.74
93.15
94.15
93.11
91.51
89.96
88.16
86.98
88.03
86.24
84.65
83.23
81.7
80.25
78.8
77.51
76.2
75.04
74
75.49
77.14
76.15
76.27
78.19
76.49
77.31
76.65
74.99
73.51
72.07
70.59
71.96
76.29
74.86
74.93
71.9
71.01
77.47
75.78
76.6
76.07
74.57
73.02
72.65
73.16
71.53
69.78
67.98
69.96
72.16
70.47
68.86
67.37
65.87
72.16
71.34
69.93
68.44
67.16
66.01
67.25
70.91
69.75
68.59
67.48
66.31
64.81
66.58
65.97
64.7
64.7
60.94
59.08
58.42
57.77
57.11
53.31
49.96
49.4
48.84
48.3
47.74
47.24
46.76
46.29
48.9
49.23
48.53
48.03
54.34
53.79
53.24
52.96
52.17
51.7
58.55
78.2
77.03
76.19
77.15
75.87
95.47
109.67
112.28
112.01
107.93
105.96
105.06
102.98
102.2
105.23
101.85
99.89
96.23
94.76
91.51
91.63
91.54
85.23
87.83
87.38
84.44
85.19
84.03
86.73
102.52
104.45
106.98
107.02
99.26
94.45
113.44
157.33
147.38
171.89
171.95
132.71
126.02
121.18
115.45
110.48
117.85
117.63
124.65
109.59
111.27
99.78
98.21
99.2
97.97
89.55
87.91
93.34
94.42
93.2
90.29
91.46
89.98
88.35
88.41
82.44
79.89
75.69
75.66
84.5
96.73
87.48
82.39
83.48
79.31
78.16
72.77
72.45
68.46
67.62
68.76
70.07
68.55
65.3
58.96
59.17
62.37
66.28
55.62
55.23
55.85
56.75
50.89
53.88
52.95
55.08
53.61
58.78
61.85
55.91
53.32
46.41
44.57
50
50
53.36
46.23
50.45
49.07
45.85
48.45
49.96
46.53
50.51
47.58
48.05
46.84
47.67
49.16
55.54
55.82
58.22
56.19
57.77
63.19
54.76
55.74
62.54
61.39
69.6
79.23
80
93.68
107.63
100.18
97.3
90.45
80.64
80.58
75.82
85.59
89.35
89.42
104.73
95.32
89.27
90.44
86.97
79.98
81.22
87.35
83.64
82.22
94.4
102.18




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 4 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=49527&T=0

[TABLE]
[ROW][C]Summary of computational 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]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=49527&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=49527&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24







Descriptive Statistics
# observations360
minimum44.57
Q166.8725
median76.775
mean77.5381666666667
Q387.0025
maximum171.95

\begin{tabular}{lllllllll}
\hline
Descriptive Statistics \tabularnewline
# observations & 360 \tabularnewline
minimum & 44.57 \tabularnewline
Q1 & 66.8725 \tabularnewline
median & 76.775 \tabularnewline
mean & 77.5381666666667 \tabularnewline
Q3 & 87.0025 \tabularnewline
maximum & 171.95 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=49527&T=1

[TABLE]
[ROW][C]Descriptive Statistics[/C][/ROW]
[ROW][C]# observations[/C][C]360[/C][/ROW]
[ROW][C]minimum[/C][C]44.57[/C][/ROW]
[ROW][C]Q1[/C][C]66.8725[/C][/ROW]
[ROW][C]median[/C][C]76.775[/C][/ROW]
[ROW][C]mean[/C][C]77.5381666666667[/C][/ROW]
[ROW][C]Q3[/C][C]87.0025[/C][/ROW]
[ROW][C]maximum[/C][C]171.95[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=49527&T=1

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

As an alternative you can also use a QR Code:  

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

Descriptive Statistics
# observations360
minimum44.57
Q166.8725
median76.775
mean77.5381666666667
Q387.0025
maximum171.95



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)
library(lattice)
bitmap(file='pic1.png')
plot(x,type='l',main='Run Sequence Plot',xlab='time or index',ylab='value')
grid()
dev.off()
bitmap(file='pic2.png')
hist(x)
grid()
dev.off()
bitmap(file='pic3.png')
if (par1 > 0)
{
densityplot(~x,col='black',main=paste('Density Plot bw = ',par1),bw=par1)
} else {
densityplot(~x,col='black',main='Density Plot')
}
dev.off()
bitmap(file='pic4.png')
qqnorm(x)
qqline(x)
grid()
dev.off()
if (par2 > 0)
{
bitmap(file='lagplot1.png')
dum <- cbind(lag(x,k=1),x)
dum
dum1 <- dum[2:length(x),]
dum1
z <- as.data.frame(dum1)
z
plot(z,main='Lag plot (k=1), lowess, and regression line')
lines(lowess(z))
abline(lm(z))
dev.off()
if (par2 > 1) {
bitmap(file='lagplotpar2.png')
dum <- cbind(lag(x,k=par2),x)
dum
dum1 <- dum[(par2+1):length(x),]
dum1
z <- as.data.frame(dum1)
z
mylagtitle <- 'Lag plot (k='
mylagtitle <- paste(mylagtitle,par2,sep='')
mylagtitle <- paste(mylagtitle,'), and lowess',sep='')
plot(z,main=mylagtitle)
lines(lowess(z))
dev.off()
}
bitmap(file='pic5.png')
acf(x,lag.max=par2,main='Autocorrelation Function')
grid()
dev.off()
}
summary(x)
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,'# observations',header=TRUE)
a<-table.element(a,length(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'minimum',header=TRUE)
a<-table.element(a,min(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q1',header=TRUE)
a<-table.element(a,quantile(x,0.25))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'median',header=TRUE)
a<-table.element(a,median(x))
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.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Q3',header=TRUE)
a<-table.element(a,quantile(x,0.75))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'maximum',header=TRUE)
a<-table.element(a,max(x))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')