Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_fitdistrnorm.wasp
Title produced by softwareML Fitting and QQ Plot- Normal Distribution
Date of computationTue, 15 Jun 2021 04:38:42 +0200
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2021/Jun/15/t1623724999g4uu5en97fbmjfb.htm/, Retrieved Fri, 26 Apr 2024 23:50:08 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=319452, Retrieved Fri, 26 Apr 2024 23:50:08 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact84
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] [] [2021-06-15 02:38:42] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
hgt
161.2
167.5
159.5
157
155.8
170
159.1
166
176.2
160.2
172.5
170.9
172.9
153.4
160
147.2
168.2
175
157
167.6
159.5
175
166.8
176.5
170.2
174
173
179.9
170.5
160
154.4
162
176.5
160
152
162.1
170
160.2
161.3
166.4
168.9
163.8
167.6
160
161.3
167.6
165.1
160
170
157.5
167.6
160.7
163.2
152.4
157.5
168.3
180.3
165.5
165
164.5
156
160
163
165.7
161
162
166
174
172.7
167.6
151.1
164.5
163.5
152
169
164
161.2
155
170
176.2
170
162.5
170.3
164.1
169.5
163.2
154.5
159.8
173.2
170
161.4
169
166.2
159.4
162.5
159
162.8
159
179.8
162.9
161
151.1
168.2
168.9
173.2
171.8
178
164.3
163
168.5
166.8
172.7
163.5
169.4
167.8
159.5
167.6
161.2
160
163.2
162.2
161.3
149.5
157.5
163.2
172.7
155
156.5
164
160.9
162.8
167
160
160
168.9
158.2
156
160
167.1
158
167.6
156
162.1
173.4
159.8
170.5
159.2
157.5
161.3
162.6
160
168.9
165.1
162.6
165.1
166.4
160
152.4
170.2
162.6
170.2
158.8
172.7
167.6
162.6
167.6
156.2
175.2
172.1
162.6
160
165.1
182.9
166.4
165.1
177.8
165.1
175.3
154.9
158.8
172.7
168.9
161.3
167.6
165.1
175.3
157.5
163.8
167.6
165.1
165.1
168.9
162.6
164.5
176.5
168.9
175.3
159.4
160
170.2
162.6
167.6
162.6
160.7
160
157.5
162.6
152.4
170.2
165.1
172.7
165.1
170.2
170.2
170.2
161.3
167.6
167.6
165.1
162.6
152.4
168.9
170.2
175.2
175.2
160
165.1
174
170.2
160
167.6
167.6
167.6
154.9
162.6
175.3
171.4
157.5
165.1
160
174
162.6
174
162.6
161.3
156.2
149.9
169.5
160
175.3
169.5
160
172.7
162.6
157.5
176.5
164.4
160.7
174
163.8




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center

\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 time1 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=319452&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]1 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=319452&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=319452&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 time1 seconds
R ServerBig Analytics Cloud Computing Center



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 <- na.omit(x)
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')