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 softwareML Fitting and QQ Plot- Normal Distribution
Date of computationThu, 26 Nov 2015 09:28:36 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2015/Nov/26/t14485301757ka9gcejatpmbbv.htm/, Retrieved Tue, 14 May 2024 10:23:34 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=284173, Retrieved Tue, 14 May 2024 10:23:34 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact113
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [ML Fitting and QQ Plot- Normal Distribution] [Normaal verdeling...] [2015-11-26 09:28:36] [20fcaaf1d4bc4a12bf87c6c50d624c14] [Current]
Feedback Forum

Post a new message
Dataseries X:
63966
60410
59440
59445
57614
55396
53030
50090
48764
61658
63943
64878
60634
57905
57224
60953
60621
57258
54903
53278
53042
63753
69210
71446
68408
65427
64630
66086
65058
62689
60841
57346
56222
68202
70745
73690
68992
65925
65546
67221
65315
62038
58774
55320
53900
65544
67906
70911
66544
63657
61720
62140
60837
58632
54642
51125
51000
60486
62685
66142
61895




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net

\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 & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=284173&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=284173&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=284173&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 time3 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net







ParameterEstimated ValueStandard Deviation
mean61427.9016393443732.456795755516
standard deviation5720.67045166514517.925167204895

\begin{tabular}{lllllllll}
\hline
Parameter & Estimated Value & Standard Deviation \tabularnewline
mean & 61427.9016393443 & 732.456795755516 \tabularnewline
standard deviation & 5720.67045166514 & 517.925167204895 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=284173&T=1

[TABLE]
[ROW][C]Parameter[/C][C]Estimated Value[/C][C]Standard Deviation[/C][/ROW]
[ROW][C]mean[/C][C]61427.9016393443[/C][C]732.456795755516[/C][/ROW]
[ROW][C]standard deviation[/C][C]5720.67045166514[/C][C]517.925167204895[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=284173&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=284173&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
mean61427.9016393443732.456795755516
standard deviation5720.67045166514517.925167204895



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)
library(car)
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()
bitmap(file='test3.png')
qqPlot(x,dist='norm',main='QQ plot (Normal) with confidence intervals')
grid()
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')