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, 16 Dec 2015 19:56:57 +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/Dec/16/t14502958393lgnqut6aomhjut.htm/, Retrieved Thu, 16 May 2024 13:49:38 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=286743, Retrieved Thu, 16 May 2024 13:49:38 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact78
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] [] [2015-12-16 19:56:57] [7b81fac622814275349f5d25cf6bd6bd] [Current]
Feedback Forum

Post a new message
Dataseries X:
191
190
220
194
144
192
162
96
265
222
182
203
198
159
196
140
163
177
123
100
227
205
133
190




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Gertrude Mary Cox' @ cox.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 & 4 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=286743&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]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=286743&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=286743&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 time4 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net







ParameterEstimated ValueStandard Deviation
mean1788.11954945383876
standard deviation39.77750620639765.74138847898892

\begin{tabular}{lllllllll}
\hline
Parameter & Estimated Value & Standard Deviation \tabularnewline
mean & 178 & 8.11954945383876 \tabularnewline
standard deviation & 39.7775062063976 & 5.74138847898892 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=286743&T=1

[TABLE]
[ROW][C]Parameter[/C][C]Estimated Value[/C][C]Standard Deviation[/C][/ROW]
[ROW][C]mean[/C][C]178[/C][C]8.11954945383876[/C][/ROW]
[ROW][C]standard deviation[/C][C]39.7775062063976[/C][C]5.74138847898892[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=286743&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=286743&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
mean1788.11954945383876
standard deviation39.77750620639765.74138847898892



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