FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_fitdistrchisq1.wasp
Title produced by softwareMaximum-likelihood Fitting - Chi-squared Distribution
Date of computationWed, 16 Dec 2015 14:59:35 +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/t14502779886zcj9ra4lihbrbu.htm/, Retrieved Thu, 16 May 2024 04:15:40 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=286712, Retrieved Thu, 16 May 2024 04:15:40 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Maximum-likelihood Fitting - Chi-squared Distribution] [] [2015-12-16 14:59:35] [c7d580dabf700704a3d6c2773d62a8e4] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
11.28213147
19.80976738
14.68725698
18.24650513
13.41384973
11.06389324
9.448028228
6.626878796
12.67032425
13.85249866
5.827720375
13.42563838
13.87352764
7.702996539
21.45831683
11.41456151
14.90524342
15.69843882
9.940968653
13.47435426
10.09514188
26.77186892
3.616551467
8.787921821
14.43049419
14.34610876
17.27412249
7.838995367
19.23249398
13.99180483
6.07853677
13.77494623
8.552342548
9.222865946
20.96687335
18.92515854
6.751124117
13.99710458
13.668061
12.29080775
7.49616521
14.20439884
12.42500634
11.34422843
9.805490973
17.65550643
24.14277747
8.223136808
10.8113172
8.889826106
6.938388402
17.59084112
10.14508619
15.25291132
7.845452531
15.61462373
7.299101692
9.780458111
10.98151097
15.38552714
15.90571297
6.389202182
12.27389423
20.70376358
12.64326228
11.06932013
6.018691925
20.49527356
11.09215162
14.99159978
20.08656489
10.77515701
6.856662284
9.460471204
20.15421393
16.26630398
20.33885739
9.04785309
10.60225122
11.08780214
11.90266828
3.222265826
9.87068255
4.986454153
9.736563134
4.98329832
12.00445042
8.229331664
13.75474856
15.8550942
15.41071084
15.48823785
10.31550008
12.10455324
15.46498584
7.62487599
22.08332131
5.359110996
5.117908595
11.30406574
17.40429208
16.71522774
9.89339975
14.89752515
15.51859733
14.63531003
9.783288851
4.947282293
15.70655178
8.657379929
12.12106752
11.98961255
7.759935801
11.14105246
11.50460439
6.238063139
13.83594229
15.53427595
6.571811424
13.50164688
4.63393299
9.28348361
7.712119379
21.96700508
3.732706175
4.980194458
15.68292657
8.884450075
11.77849441
6.461078911
3.786166947
14.21040511
18.41954551
7.372216977
6.938859713
6.894730972
4.016769402
15.97575731
9.353299427
18.09360447
6.938790204
5.905288571
19.57651319
7.53499932
14.79370398
10.67050305
6.788939275
14.38916753
13.16471132

Tables (Output of Computation)





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=286712&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=286712&T=0

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






ParameterEstimated ValueStandard Deviation
df11.91936577542750.383373000251455

\begin{tabular}{lllllllll}
\hline
Parameter & Estimated Value & Standard Deviation \tabularnewline
df & 11.9193657754275 & 0.383373000251455 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=286712&T=1

[TABLE]
[ROW][C]Parameter[/C][C]Estimated Value[/C][C]Standard Deviation[/C][/ROW]
[ROW][C]df[/C][C]11.9193657754275[/C][C]0.383373000251455[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=286712&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=286712&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
df11.91936577542750.383373000251455


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
library(MASS)
library(car)
chii <- x
chi_df <- fitdistr(chii,'chi-squared',start=list(df=3), method='Brent',lower=0.1,upper=100)
chi_k <- chi_df[[1]][1]
bitmap(file='test1.png')
chi_hist <- hist(chii,breaks=50,col='gray',freq=FALSE,main='Histogram scaled to density',xlab='x')
curve(dchisq(x,df=chi_k),add=TRUE,col='red')
dev.off()
f<-chi_df
f$estimate
f$sd
xlab <- paste('Chisq(df = ',round(f$estimate[[1]],2),sep='')
xlab <- paste(xlab,')',sep='')
xlab
bitmap(file='test2.png')
qqPlot(x, dist='chisq', df=f$estimate[[1]], ncp=0, main='QQ plot (Chi-squared 1 param.)', xlab=xlab )
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,'df',header=TRUE)
a<-table.element(a,f$estimate[1])
a<-table.element(a,f$sd[1])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')