Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_fitdistrnorm.wasp
Title produced by softwareMaximum-likelihood Fitting - Normal Distribution
Date of computationMon, 27 Oct 2008 15:44:06 -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/2008/Oct/27/t1225143886ysobe6nke56fma3.htm/, Retrieved Tue, 28 May 2024 20:12:04 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=19656, Retrieved Tue, 28 May 2024 20:12:04 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact147
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Maximum-likelihood Fitting - Normal Distribution] [Taak 2, model] [2008-10-25 12:48:23] [87cabf13a90315c7085b765dcebb7412]
F    D    [Maximum-likelihood Fitting - Normal Distribution] [Q7] [2008-10-27 21:44:06] [f6a332ba2d530c028d935c5a5bbb53af] [Current]
Feedback Forum
2008-10-28 20:36:18 [Ken Van den Heuvel] [reply
Task 2:

Je gebruikt hier 2 nieuwe methoden die buiten de vraagstelling staan. De bedoeling was om de linkerkant van de vergelijking aan te passen en daaruit conclusies te trekken.

Zo kon je bijvoorveeld:

Prijs kledij / Totale Productie = constant + random component
en
Investeringen / Totale Productie = constant + random component
Wederom de assumpties nagaan per model, de resultaten vergelijken, en conclusies hieruit trekken.

2008-11-03 06:57:01 [Jeroen Michel] [reply
Je maakt hier inderdaad niet echt gebruik van twee nieuwe modellen. Enerzijds kon je best andere software gebruiken en anderzijds was het dus aan te raden geweest om de linkerhelft van de vergelijking te veranderen.

Zie opmerkingen van vorige student om deze linkerhelft te wijzigen.
2008-11-03 07:14:29 [Jeroen Michel] [reply
Voor deze vraag heb ik toch een herberekening gemaakt. Dit aangezien er enkele foute redeneringen zijn gemaakt. Je geeft vaak het juiste antwoord alleen de redenering erachter zit niet altijd juist.

Assumptie 1:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/Nov/03/t1225696290mpt0qxyer9725k9.htm, Retrieved Mon, 03 Nov 2008 07:11:34 +0000

Assumptie 2:
Zie ook hier naar de vorige link!

Assumptie 3:
zie ook hier de eerste link. Jouw redenering blijft me vreemd.

Assumptie 4:
Ook hier bekom je de juiste resultaten, maar hier is de redenering weer het tegenovergestelde.
2008-11-03 19:50:36 [Kevin Truyts] [reply
De student heeft de vraag verkeerd begrepen. Het was namelijk de bedoling om een nieuwe tijdreeks te maken van de reeks van kledij door op een creatieve manier de reeks te veranderen/beïnvloeden.

Post a new message
Dataseries X:
109.20
88.60
94.30
98.30
86.40
80.60
104.10
108.20
93.40
71.90
94.10
94.90
96.40
91.10
84.40
86.40
88.00
75.10
109.70
103.00
82.10
68.00
96.40
94.30
90.00
88.00
76.10
82.50
81.40
66.50
97.20
94.10
80.70
70.50
87.80
89.50
99.60
84.20
75.10
92.00
80.80
73.10
99.80
90.00
83.10
72.40
78.80
87.30
91.00
80.10
73.60
86.40
74.50
71.20
92.40
81.50
85.30
69.90
84.20
90.70
100.30




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 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 & 2 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=19656&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]2 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=19656&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=19656&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 time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24







ParameterEstimated ValueStandard Deviation
mean86.89344262295081.34220127658893
standard deviation10.48292708548020.949079624393276

\begin{tabular}{lllllllll}
\hline
Parameter & Estimated Value & Standard Deviation \tabularnewline
mean & 86.8934426229508 & 1.34220127658893 \tabularnewline
standard deviation & 10.4829270854802 & 0.949079624393276 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=19656&T=1

[TABLE]
[ROW][C]Parameter[/C][C]Estimated Value[/C][C]Standard Deviation[/C][/ROW]
[ROW][C]mean[/C][C]86.8934426229508[/C][C]1.34220127658893[/C][/ROW]
[ROW][C]standard deviation[/C][C]10.4829270854802[/C][C]0.949079624393276[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=19656&T=1

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

As an alternative you can also use a QR Code:  

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

ParameterEstimated ValueStandard Deviation
mean86.89344262295081.34220127658893
standard deviation10.48292708548020.949079624393276



Parameters (Session):
par1 = 8 ; par2 = 0 ;
Parameters (R input):
par1 = 8 ; par2 = 0 ;
R code (references can be found in the software module):
library(MASS)
par1 <- as.numeric(par1)
if (par2 == '0') par2 = 'Sturges' else par2 <- as.numeric(par2)
x <- as.ts(x) #otherwise the fitdistr function does not work properly
r <- fitdistr(x,'normal')
r
bitmap(file='test1.png')
myhist<-hist(x,col=par1,breaks=par2,main=main,ylab=ylab,xlab=xlab,freq=F)
curve(1/(r$estimate[2]*sqrt(2*pi))*exp(-1/2*((x-r$estimate[1])/r$estimate[2])^2),min(x),max(x),add=T)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Parameter',1,TRUE)
a<-table.element(a,'Estimated Value',1,TRUE)
a<-table.element(a,'Standard Deviation',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean',header=TRUE)
a<-table.element(a,r$estimate[1])
a<-table.element(a,r$sd[1])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'standard deviation',header=TRUE)
a<-table.element(a,r$estimate[2])
a<-table.element(a,r$sd[2])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')