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 computationWed, 30 Nov 2016 13:57:29 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2016/Nov/30/t148051068074ikykd3g1ct9e8.htm/, Retrieved Sat, 18 May 2024 23:42:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=297386, Retrieved Sat, 18 May 2024 23:42:19 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact65
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] [] [2016-11-30 12:57:29] [d06ec19b175650a2a09ee5879d174acf] [Current]
Feedback Forum

Post a new message
Dataseries X:
3406
3416
3426
3444
3464
3482
3496
3510
3530
3556
3580
3604
3620
3636
3648
3662
3676
3690
3700
3712
3726
3740
3760
3762
3766
3784
3800
3814
3828
3840
3850
3858
3866
3872
3878
3884
3920
3932
3946
3962
3976
3994
4012
4042
4062
4084
4106
4128
4132
4146
4154
4190
4134
4160
4174
4182
4224
4290
4330
4370
4398
4426
4458
4492
4528
4576
4624
4672
4720
4764
4808
4848
4866
4896
4926
4958
4988
5020
5052
5084
5116
5150
5182
5216
5342
5360
5380
5404
5430
5458
5490
5524
5560
5592
5626
5660
5702
5742
5784
5824
5866
5908
5956
6002
6050
6096
6140
6182
6226
6268
6312
6356
6398
6440
6480
6520
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA
NA




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center
R Engine error message
Error in fitdistr(x, "normal") : 'x' contains missing or infinite values
Execution halted

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time2 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
R Engine error message & 
Error in fitdistr(x, "normal") : 'x' contains missing or infinite values
Execution halted
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=297386&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]2 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [ROW]R Engine error message[/C][C]
Error in fitdistr(x, "normal") : 'x' contains missing or infinite values
Execution halted
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=297386&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=297386&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center
R Engine error message
Error in fitdistr(x, "normal") : 'x' contains missing or infinite values
Execution halted



Parameters (Session):
par1 = 8 ; par2 = 0 ;
Parameters (R input):
par1 = 8 ; par2 = 0 ;
R code (references can be found in the software module):
par2 <- '0'
par1 <- '8'
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')
print(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')