FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_fitdistrnorm.wasp
Title produced by softwareMaximum-likelihood Fitting - Normal Distribution
Date of computationThu, 13 Nov 2008 02:59:11 -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/2008/Nov/13/t12265704105s6tv7vm3qt4c2x.htm/, Retrieved Sun, 19 May 2024 01:38:17 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24533, Retrieved Sun, 19 May 2024 01:38:17 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Maximum-likelihood Fitting - Normal Distribution] [Stefan Temmerman] [2008-11-13 09:18:12] [672890602bb42231e2a9172e5a234a7f]
F    D    [Maximum-likelihood Fitting - Normal Distribution] [] [2008-11-13 09:59:11] [8767719db498704e1fee27044c098ad0] [Current]
Feedback Forum
2008-11-20 14:13:15 [6066575aa30c0611e452e930b1dff53d] [reply] 
Bij deze vraag werd ook geen uitleg gegeven. Op de x-as vinden we de waarde van de data series, op de y-as vinden we de dichtheid. We zien duidelijk dat de meeste waarden zich bevinden tussen 40 en 50. Er bevinden zich zeer weinig waarden tussen 20 en 30.
2008-11-20 16:39:21 [Gert-Jan Geudens] [reply] 
Geen conclusie gegeven. We zien duidelijk dat de gegevens niet normaal verdeeld zijn. Deze bewerking kan je ook toepassen op de andere variabelen.
2008-11-20 17:19:01 [Gert-Jan Geudens] [reply] 
Graag zou ik nog even willen vermelden dat de gegevens meer naar rechts verdeeld zijn.
2008-11-24 20:30:14 [4679c4d03f1d346a85e79d87ba60ec2b] [reply] 
Er wordt geen conclusie gegeven. De histogram bevestigt duidelijk dat data niet normaal verdeeld is.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
45
24
18
20
22
39
55
35
38
47
1
57
50
33
19
2
7
15
56
53
24
48
2
49
46
32
37
10
8
16
55
46
46
45
6
45
52
44
35
15
44
51
58
23
44
43
6
51
53
47
19
18
38
43
23
43
18
43
6
31
49

Tables (Output of Computation)





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

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






ParameterEstimated ValueStandard Deviation
mean33.57377049180332.15876917269542
standard deviation16.86052623140171.52648032102941

\begin{tabular}{lllllllll}
\hline
Parameter & Estimated Value & Standard Deviation \tabularnewline
mean & 33.5737704918033 & 2.15876917269542 \tabularnewline
standard deviation & 16.8605262314017 & 1.52648032102941 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24533&T=1

[TABLE]
[ROW][C]Parameter[/C][C]Estimated Value[/C][C]Standard Deviation[/C][/ROW]
[ROW][C]mean[/C][C]33.5737704918033[/C][C]2.15876917269542[/C][/ROW]
[ROW][C]standard deviation[/C][C]16.8605262314017[/C][C]1.52648032102941[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24533&T=1

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

As an alternative you can also use a QR Code:  
QR Code


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

ParameterEstimated ValueStandard Deviation
mean33.57377049180332.15876917269542
standard deviation16.86052623140171.52648032102941


Figures (Output of Computation)

Input Parameters & R Code

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')