515HW3
Zihao Lin
5/10/2021
R Markdown
This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
summary(cars)## speed dist
## Min. : 4.0 Min. : 2.00
## 1st Qu.:12.0 1st Qu.: 26.00
## Median :15.0 Median : 36.00
## Mean :15.4 Mean : 42.98
## 3rd Qu.:19.0 3rd Qu.: 56.00
## Max. :25.0 Max. :120.00Including Plots
You can also embed plots, for example:
Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.
Part A – Risks
## 1. Using two securities from your portfolio (from different sectors)
## a) Construct a variance-covariance matrix
library("quantmod")## Warning: package 'quantmod' was built under R version 3.6.3## Loading required package: xts## Warning: package 'xts' was built under R version 3.6.3## Loading required package: zoo## Warning: package 'zoo' was built under R version 3.6.3##
## Attaching package: 'zoo'## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric## Loading required package: TTR## Warning: package 'TTR' was built under R version 3.6.3## Registered S3 method overwritten by 'quantmod':
## method from
## as.zoo.data.frame zoodata.T = read.csv("/Users/FrankLin/Desktop/515/Communication/T.csv", header = TRUE)
date.T <- as.Date(data.T$Date,format="%Y-%m-%d")
data.T <- cbind(date.T, data.T[,-1])
data.T <- data.T[order(data.T$date),]
data.T <- xts(data.T[,2:7],order.by=data.T[,1])
names(data.T) <- paste(c("T.Open","T.High","T.Low","T.Close","T.Adjusted","T.Volume"))
T.ret <- Delt(data.T$T.Adjusted)
T.ret[c(1:3,nrow(T.ret)),]## Delt.1.arithmetic
## 2009-01-02 NA
## 2009-01-05 -0.033650103
## 2009-01-06 -0.004572869
## 2013-12-31 -0.001135968data.HD = read.csv("/Users/FrankLin/Desktop/515/ConsDiscr/HD.csv", header = TRUE)
date.HD <- as.Date(data.HD$Date,format="%Y-%m-%d")
data.HD <- cbind(date.HD, data.HD[,-1])
data.HD <- data.HD[order(data.HD$date),]
data.HD <- xts(data.HD[,2:7],order.by=data.HD[,1])
names(data.HD) <- paste(c("HD.Open","HD.High","HD.Low","HD.Close","HD.Adjusted","HD.Volume"))
HD.ret <- Delt(data.HD$HD.Adjusted)
HD.ret[c(1:3,nrow(HD.ret)),]## Delt.1.arithmetic
## 2009-01-02 NA
## 2009-01-05 0.024036682
## 2009-01-06 0.022257990
## 2013-12-31 0.004881718## Combine the two stocks' return
returns <- cbind(T.ret,HD.ret)
names(returns) <- paste(c("T.Ret","HD.Ret"))
returns <- returns[-1,]
returns[c(1:3,nrow(return)),]## T.Ret HD.Ret
## 2009-01-05 -0.033650103 0.02403668
## 2009-01-06 -0.004572869 0.02225799
## 2009-01-07 -0.024381764 -0.02731610## Create a vector of weights T: 50% , HD: 50%
WGT.2asset <- c(0.5,0.5)
WGT.2asset <- matrix(WGT.2asset,1)
WGT.2asset## [,1] [,2]
## [1,] 0.5 0.5## Create a transposed vector of weights (command: t)
tWGT.2asset <- t(WGT.2asset)
tWGT.2asset## [,1]
## [1,] 0.5
## [2,] 0.5## Construct Variance-Covariance Matrix
mat.Ret <- as.matrix(returns)
head(mat.Ret)## T.Ret HD.Ret
## 2009-01-05 -0.033650103 0.024036682
## 2009-01-06 -0.004572869 0.022257990
## 2009-01-07 -0.024381764 -0.027316100
## 2009-01-08 -0.001102208 -0.007325663
## 2009-01-09 -0.018396017 -0.021730271
## 2009-01-12 -0.027736112 -0.037301117Part B – Suitable distributions for returns
## 1. Fit the data using GHD, HYP and NIG
Part2.HD = read.csv("/Users/FrankLin/Desktop/515/ConsDiscr/HD.csv", header = TRUE)head(Part2.HD)## Date Open High Low Close Adj.Close Volume
## 1 2009-01-02 23.07 24.19 22.96 24.13 17.79792 14902500
## 2 2009-01-05 24.02 24.88 23.81 24.71 18.22573 22570900
## 3 2009-01-06 24.94 25.49 24.62 25.26 18.63139 22856600
## 4 2009-01-07 24.87 25.16 24.43 24.57 18.12246 19857300
## 5 2009-01-08 24.12 24.46 23.70 24.39 17.98970 17319800
## 6 2009-01-09 24.48 24.48 23.54 23.86 17.59878 16927400date <- Part2.HD$Date
HDPrice <- Part2.HD$Adj.Closestr(date)## Factor w/ 1258 levels "2009-01-02","2009-01-05",..: 1 2 3 4 5 6 7 8 9 10 ...str(HDPrice)## num [1:1258] 17.8 18.2 18.6 18.1 18 ...attr(HDPrice, 'time')<-date
head(HDPrice)## [1] 17.79792 18.22573 18.63139 18.12246 17.98970 17.59878install.packages("FRAPO", repos = "http://cran.us.r-project.org")## Installing package into 'C:/Users/FrankLin/Documents/R/win-library/3.6'
## (as 'lib' is unspecified)## package 'FRAPO' successfully unpacked and MD5 sums checked
##
## The downloaded binary packages are in
## C:\Users\FrankLin\AppData\Local\Temp\RtmpiQZClE\downloaded_packageslibrary(FRAPO)## Warning: package 'FRAPO' was built under R version 3.6.3## Loading required package: cccp## Warning: package 'cccp' was built under R version 3.6.3## Loading required package: Rglpk## Warning: package 'Rglpk' was built under R version 3.6.3## Loading required package: slam## Warning: package 'slam' was built under R version 3.6.3## Using the GLPK callable library version 4.47## Loading required package: timeSeries## Warning: package 'timeSeries' was built under R version 3.6.3## Loading required package: timeDate## Warning: package 'timeDate' was built under R version 3.6.3##
## Attaching package: 'timeSeries'## The following object is masked from 'package:zoo':
##
## time<-## Financial Risk Modelling and Portfolio Optimisation with R (version 0.4-1)HDRet <- returnseries(HDPrice)
summary(HDRet)## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## -6.22163 -0.64923 0.07836 0.10882 0.91000 9.96251 1attr(HDRet, 'time')<-date
str(HDRet)## num [1:1258] NA 2.375 2.201 -2.77 -0.735 ...
## - attr(*, "time")= Factor w/ 1258 levels "2009-01-02","2009-01-05",..: 1 2 3 4 5 6 7 8 9 10 ...install.packages("ghyp", repos = "http://cran.us.r-project.org")## Installing package into 'C:/Users/FrankLin/Documents/R/win-library/3.6'
## (as 'lib' is unspecified)## package 'ghyp' successfully unpacked and MD5 sums checked## Warning: cannot remove prior installation of package 'ghyp'## Warning in file.copy(savedcopy, lib, recursive = TRUE): problem copying C:
## \Users\FrankLin\Documents\R\win-library\3.6\00LOCK\ghyp\libs\x64\ghyp.dll to C:
## \Users\FrankLin\Documents\R\win-library\3.6\ghyp\libs\x64\ghyp.dll: Permission
## denied## Warning: restored 'ghyp'##
## The downloaded binary packages are in
## C:\Users\FrankLin\AppData\Local\Temp\RtmpiQZClE\downloaded_packageslibrary(ghyp)## Warning: package 'ghyp' was built under R version 3.6.3## Loading required package: numDeriv## Loading required package: MASS## Warning: package 'MASS' was built under R version 3.6.3##
## Attaching package: 'MASS'## The following object is masked from 'package:FRAPO':
##
## SP500library(timeSeries)
install.packages("fBasics", repos = "http://cran.us.r-project.org")## Installing package into 'C:/Users/FrankLin/Documents/R/win-library/3.6'
## (as 'lib' is unspecified)## package 'fBasics' successfully unpacked and MD5 sums checked## Warning: cannot remove prior installation of package 'fBasics'## Warning in file.copy(savedcopy, lib, recursive = TRUE): problem copying C:
## \Users\FrankLin\Documents\R\win-library\3.6\00LOCK\fBasics\libs\x64\fBasics.dll
## to C:\Users\FrankLin\Documents\R\win-library\3.6\fBasics\libs\x64\fBasics.dll:
## Permission denied## Warning: restored 'fBasics'##
## The downloaded binary packages are in
## C:\Users\FrankLin\AppData\Local\Temp\RtmpiQZClE\downloaded_packageslibrary(fBasics)## Warning: package 'fBasics' was built under R version 3.6.3##
## Attaching package: 'fBasics'## The following object is masked from 'package:TTR':
##
## volatilitydatets <- as.character(date)
HDTimeS <- timeSeries(HDRet, charvec = datets)
str(HDTimeS)## Time Series:
## Name: object
## Data Matrix:
## Dimension: 1258 1
## Column Names: TS.1
## Row Names: 2009-01-02 ... 2013-12-31
## Positions:
## Start: 2009-01-02
## End: 2013-12-31
## With:
## Format: %Y-%m-%d
## FinCenter: GMT
## Units: TS.1
## Title: Time Series Object
## Documentation: Mon May 10 23:17:08 2021head(HDTimeS)## GMT
## TS.1
## 2009-01-02 NA
## 2009-01-05 2.3752349
## 2009-01-06 2.2013896
## 2009-01-07 -2.7696121
## 2009-01-08 -0.7352627
## 2009-01-09 -2.1969851ef<-density(HDTimeS, na.rm=TRUE)
plot(ef)
ghdfit <- fit.ghypuv(HDTimeS, symmetric = FALSE, control = list(maxit = 1000), na.rm = TRUE)## Warning in .check.data(data = data, case = "uv", na.rm = na.rm): 1 NA
## observations removed## [1] "Llh: -2.33741670315059E+03; Pars: 1.000000E+00, 5.000000E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.33754570259514E+03; Pars: 1.100000E+00, 5.000000E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.33728608450772E+03; Pars: 1.000000E+00, 5.525855E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.34158775859322E+03; Pars: 1.000000E+00, 5.000000E-01, 1.783586E-01, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.29883749960255E+03; Pars: 1.000000E+00, 5.000000E-01, 7.835855E-02, 1.273345E+00, 0.000000E+00"
## [1] "Llh: -2.33997953524100E+03; Pars: 1.000000E+00, 5.000000E-01, 7.835855E-02, 1.152171E+00, 1.000000E-01"
## [1] "Llh: -2.32221899921452E+03; Pars: 1.040000E+00, 5.204054E-01, -2.164145E-02, 1.199192E+00, 4.000000E-02"
## [1] "Llh: -2.32405672745777E+03; Pars: 1.030000E+00, 5.152273E-01, 2.835855E-02, 1.187259E+00, 3.000000E-02"
## [1] "Llh: -2.32314766108921E+03; Pars: 1.056000E+00, 5.287988E-01, 3.835855E-02, 1.218533E+00, -8.400000E-02"
## [1] "Llh: -2.32247356708584E+03; Pars: 1.042000E+00, 5.214472E-01, 4.835855E-02, 1.201592E+00, -3.800000E-02"
## [1] "Llh: -2.30936183227721E+03; Pars: 9.328000E-01, 5.377577E-01, 2.635855E-02, 1.239177E+00, 8.000000E-04"
## [1] "Llh: -2.31521810382723E+03; Pars: 9.746000E-01, 5.280590E-01, 3.935855E-02, 1.216828E+00, 6.000000E-04"
## [1] "Llh: -2.30143227975550E+03; Pars: 1.005920E+00, 5.536474E-01, 5.558551E-03, 1.275793E+00, 1.120000E-03"
## [1] "Llh: -2.30824453068309E+03; Pars: 1.004440E+00, 5.397187E-01, 2.375855E-02, 1.243696E+00, 8.400000E-04"
## [1] "Llh: -2.29332418439012E+03; Pars: 1.008288E+00, 5.013458E-01, -2.356145E-02, 1.328879E+00, 1.568000E-03"
## [1] "Llh: -2.28666960924766E+03; Pars: 1.012432E+00, 4.775363E-01, -7.452145E-02, 1.427150E+00, 2.352000E-03"
## [1] "Llh: -2.28534270581671E+03; Pars: 9.544608E-01, 5.129144E-01, -4.271345E-02, 1.364992E+00, 5.570880E-02"
## [1] "Llh: -2.27790903407112E+03; Pars: 9.106912E-01, 5.087004E-01, -8.824945E-02, 1.454844E+00, 1.025632E-01"
## [1] "Llh: -2.27545168641244E+03; Pars: 9.047373E-01, 5.092858E-01, 6.433514E-04, 1.477644E+00, 2.734080E-03"
## [1] "Llh: -2.27570326362789E+03; Pars: 8.371059E-01, 5.038154E-01, 1.178575E-02, 1.640251E+00, -1.589888E-02"
## [1] "Llh: -2.27639660901425E+03; Pars: 1.000712E+00, 4.822438E-01, -5.764273E-02, 1.534262E+00, 4.270771E-02"
## [1] "Llh: -2.27782937255905E+03; Pars: 9.837341E-01, 4.955605E-01, -3.664241E-02, 1.454481E+00, 3.223078E-02"
## [1] "Llh: -2.27594085039009E+03; Pars: 9.255091E-01, 4.432343E-01, -6.212324E-02, 1.604356E+00, 5.902280E-02"
## [1] "Llh: -2.27518588647002E+03; Pars: 9.456118E-01, 4.685796E-01, -4.520279E-02, 1.515027E+00, 4.454710E-02"
## [1] "Llh: -2.29562231376693E+03; Pars: 9.096738E-01, 4.782137E-01, -1.843478E-01, 1.723154E+00, 7.796164E-02"
## [1] "Llh: -2.28209119100051E+03; Pars: 9.322553E-01, 4.835696E-01, -1.186712E-01, 1.597643E+00, 5.847123E-02"
## [1] "Llh: -2.27365283511279E+03; Pars: 8.651711E-01, 5.032303E-01, -4.912768E-02, 1.608439E+00, 9.805733E-02"
## [1] "Llh: -2.27734448716220E+03; Pars: 7.915407E-01, 5.165911E-01, -3.643079E-02, 1.707545E+00, 1.459100E-01"
## [1] "Llh: -2.27440244320156E+03; Pars: 9.185141E-01, 5.049370E-01, 2.283948E-02, 1.440671E+00, 5.777254E-02"
## [1] "Llh: -2.27356450311286E+03; Pars: 9.219494E-01, 4.995082E-01, -1.253819E-02, 1.478405E+00, 5.794721E-02"
## [1] "Llh: -2.27336506229155E+03; Pars: 9.445815E-01, 4.765037E-01, 2.270223E-02, 1.592277E+00, -4.165829E-03"
## [1] "Llh: -2.27615373180411E+03; Pars: 9.615267E-01, 4.611778E-01, 7.817807E-02, 1.665788E+00, -5.753034E-02"
## [1] "Llh: -2.27106841754615E+03; Pars: 8.321083E-01, 5.002436E-01, 2.423350E-02, 1.532439E+00, 3.694024E-02"
## [1] "Llh: -2.27027268909037E+03; Pars: 7.478063E-01, 5.094939E-01, 6.517161E-02, 1.531529E+00, 3.405651E-02"
## [1] "Llh: -2.27200882866919E+03; Pars: 8.653108E-01, 4.737592E-01, -8.241277E-03, 1.614118E+00, 8.944285E-02"
## [1] "Llh: -2.27158383999339E+03; Pars: 8.751674E-01, 4.824015E-01, -6.020120E-03, 1.578861E+00, 6.776565E-02"
## [1] "Llh: -2.27088079203931E+03; Pars: 7.962585E-01, 5.209379E-01, 5.327794E-02, 1.600493E+00, 5.691725E-02"
## [1] "Llh: -2.27078863151287E+03; Pars: 8.335968E-01, 5.073239E-01, 2.865775E-02, 1.578685E+00, 5.382471E-02"
## [1] "Llh: -2.27171306952791E+03; Pars: 8.640695E-01, 4.866425E-01, 8.831699E-02, 1.496325E+00, -1.428602E-02"
## [1] "Llh: -2.27098946180098E+03; Pars: 8.643449E-01, 4.907374E-01, 5.395582E-02, 1.523599E+00, 1.379981E-02"
## [1] "Llh: -2.27235196077317E+03; Pars: 7.842494E-01, 4.868069E-01, 7.832511E-02, 1.647662E+00, 8.165133E-03"
## [1] "Llh: -2.27110527150245E+03; Pars: 8.186744E-01, 4.899516E-01, 5.560929E-02, 1.603613E+00, 2.061065E-02"
## [1] "Llh: -2.27037568697101E+03; Pars: 7.112544E-01, 5.160215E-01, 5.624751E-02, 1.534176E+00, 8.018877E-02"
## [1] "Llh: -2.27011259174291E+03; Pars: 7.695862E-01, 5.058450E-01, 4.786119E-02, 1.548500E+00, 5.910012E-02"
## [1] "Llh: -2.27053302722858E+03; Pars: 7.384360E-01, 5.194803E-01, 1.065224E-01, 1.535242E+00, 4.791067E-03"
## [1] "Llh: -2.27023421621463E+03; Pars: 7.726189E-01, 5.099516E-01, 7.838676E-02, 1.546033E+00, 2.053471E-02"
## [1] "Llh: -2.27139029160670E+03; Pars: 7.765068E-01, 5.197264E-01, 5.400397E-02, 1.489598E+00, 5.191569E-02"
## [1] "Llh: -2.27044643657970E+03; Pars: 8.081325E-01, 4.972315E-01, 5.520796E-02, 1.574316E+00, 2.843691E-02"
## [1] "Llh: -2.27027526746344E+03; Pars: 7.083514E-01, 5.216300E-01, 5.615828E-02, 1.588497E+00, 6.458137E-02"
## [1] "Llh: -2.27012947606611E+03; Pars: 7.473498E-01, 5.137291E-01, 5.560767E-02, 1.572018E+00, 5.188598E-02"
## [1] "Llh: -2.27041509521327E+03; Pars: 7.046006E-01, 5.071144E-01, 9.223632E-02, 1.530475E+00, 2.378098E-02"
## [1] "Llh: -2.27014659880210E+03; Pars: 7.368497E-01, 5.071668E-01, 7.634168E-02, 1.542388E+00, 3.129191E-02"
## [1] "Llh: -2.27049600967366E+03; Pars: 7.015518E-01, 5.215183E-01, 7.413961E-02, 1.522196E+00, 5.031078E-02"
## [1] "Llh: -2.27018378507681E+03; Pars: 7.814873E-01, 5.031950E-01, 5.994087E-02, 1.561121E+00, 3.390538E-02"
## [1] "Llh: -2.27024127205037E+03; Pars: 7.753504E-01, 5.064401E-01, 6.208366E-02, 1.576746E+00, 4.463073E-02"
## [1] "Llh: -2.27012928658192E+03; Pars: 7.684644E-01, 5.072018E-01, 6.285564E-02, 1.565318E+00, 4.198718E-02"
## [1] "Llh: -2.27009801718890E+03; Pars: 7.488761E-01, 5.048925E-01, 4.265606E-02, 1.569718E+00, 6.673351E-02"
## [1] "Llh: -2.27027749078192E+03; Pars: 7.370047E-01, 5.023818E-01, 2.479071E-02, 1.581697E+00, 8.983291E-02"
## [1] "Llh: -2.27004708863524E+03; Pars: 7.269631E-01, 5.123616E-01, 5.418803E-02, 1.557966E+00, 6.649410E-02"
## [1] "Llh: -2.27014191274446E+03; Pars: 6.997010E-01, 5.170072E-01, 5.131161E-02, 1.556391E+00, 8.278846E-02"
## [1] "Llh: -2.27041457021433E+03; Pars: 7.676461E-01, 5.104255E-01, 2.892576E-02, 1.583240E+00, 8.318844E-02"
## [1] "Llh: -2.27005356735494E+03; Pars: 7.445488E-01, 5.079795E-01, 6.448770E-02, 1.552501E+00, 4.426605E-02"
## [1] "Llh: -2.27008470394961E+03; Pars: 7.560256E-01, 5.016419E-01, 5.321178E-02, 1.545655E+00, 5.954640E-02"
## [1] "Llh: -2.27005524788814E+03; Pars: 7.538567E-01, 5.046367E-01, 5.381075E-02, 1.552204E+00, 5.763130E-02"
## [1] "Llh: -2.27008099853787E+03; Pars: 7.290679E-01, 5.070682E-01, 4.234585E-02, 1.547056E+00, 7.570285E-02"
## [1] "Llh: -2.27004828546181E+03; Pars: 7.389170E-01, 5.071016E-01, 4.747330E-02, 1.551602E+00, 6.727393E-02"
## [1] "Llh: -2.27001040185206E+03; Pars: 7.156785E-01, 5.089331E-01, 5.718515E-02, 1.565111E+00, 6.185944E-02"
## [1] "Llh: -2.27003807221119E+03; Pars: 6.887246E-01, 5.104843E-01, 6.184712E-02, 1.573483E+00, 6.323910E-02"
## [1] "Llh: -2.27013689607207E+03; Pars: 7.231096E-01, 5.115216E-01, 6.820191E-02, 1.542140E+00, 5.227641E-02"
## [1] "Llh: -2.27004056971947E+03; Pars: 7.424344E-01, 5.065417E-01, 4.904252E-02, 1.562778E+00, 6.311924E-02"
## [1] "Llh: -2.27001267330865E+03; Pars: 7.135601E-01, 5.125528E-01, 5.513993E-02, 1.563782E+00, 6.357381E-02"
## [1] "Llh: -2.27001509599113E+03; Pars: 7.236342E-01, 5.105622E-01, 5.480763E-02, 1.560879E+00, 6.208818E-02"
## [1] "Llh: -2.27009017729805E+03; Pars: 7.104725E-01, 5.110086E-01, 4.072387E-02, 1.568017E+00, 8.466216E-02"
## [1] "Llh: -2.27002478534931E+03; Pars: 7.360297E-01, 5.087351E-01, 5.854674E-02, 1.556365E+00, 5.436508E-02"
## [1] "Llh: -2.27004980747245E+03; Pars: 7.149493E-01, 5.125522E-01, 6.216765E-02, 1.570851E+00, 5.649073E-02"
## [1] "Llh: -2.27002652284136E+03; Pars: 7.329251E-01, 5.084588E-01, 5.114689E-02, 1.556392E+00, 6.457813E-02"
## [1] "Llh: -2.27000140189311E+03; Pars: 7.292880E-01, 5.057411E-01, 5.423646E-02, 1.563802E+00, 5.650418E-02"
## [1] "Llh: -2.27000417073467E+03; Pars: 7.304505E-01, 5.024630E-01, 5.426068E-02, 1.566727E+00, 5.150921E-02"
## [1] "Llh: -2.26999801053269E+03; Pars: 7.085581E-01, 5.112282E-01, 6.145954E-02, 1.559395E+00, 5.723301E-02"
## [1] "Llh: -2.27000996036284E+03; Pars: 6.916199E-01, 5.135877E-01, 6.766805E-02, 1.557706E+00, 5.428990E-02"
## [1] "Llh: -2.26999961096343E+03; Pars: 7.083207E-01, 5.104085E-01, 6.348024E-02, 1.567001E+00, 5.283607E-02"
## [1] "Llh: -2.26999838095534E+03; Pars: 7.144718E-01, 5.099204E-01, 6.039690E-02, 1.564342E+00, 5.577159E-02"
## [1] "Llh: -2.27000147155201E+03; Pars: 6.965929E-01, 5.106063E-01, 5.682045E-02, 1.570235E+00, 6.361174E-02"
## [1] "Llh: -2.26999681070138E+03; Pars: 7.064521E-01, 5.101379E-01, 5.725202E-02, 1.566756E+00, 6.130007E-02"
## [1] "Llh: -2.26998409725805E+03; Pars: 7.162193E-01, 5.058467E-01, 6.107210E-02, 1.563977E+00, 5.349351E-02"
## [1] "Llh: -2.26997767057125E+03; Pars: 7.175489E-01, 5.025266E-01, 6.403819E-02, 1.564074E+00, 4.845336E-02"
## [1] "Llh: -2.26997614807858E+03; Pars: 7.148491E-01, 5.068694E-01, 6.176810E-02, 1.562234E+00, 4.984544E-02"
## [1] "Llh: -2.26997401112299E+03; Pars: 7.144344E-01, 5.058407E-01, 6.405957E-02, 1.560798E+00, 4.383844E-02"
## [1] "Llh: -2.26998793431863E+03; Pars: 6.952981E-01, 5.101088E-01, 6.864603E-02, 1.562340E+00, 5.013441E-02"
## [1] "Llh: -2.26998153614320E+03; Pars: 7.037956E-01, 5.090133E-01, 6.504364E-02, 1.562705E+00, 5.172685E-02"
## [1] "Llh: -2.26996375330849E+03; Pars: 7.058439E-01, 5.055678E-01, 6.434428E-02, 1.561147E+00, 4.924911E-02"
## [1] "Llh: -2.26995503809411E+03; Pars: 7.015299E-01, 5.034055E-01, 6.631798E-02, 1.559552E+00, 4.598787E-02"
## [1] "Llh: -2.26996077881237E+03; Pars: 7.089463E-01, 5.011736E-01, 6.522502E-02, 1.566163E+00, 4.328962E-02"
## [1] "Llh: -2.26996259655143E+03; Pars: 7.088492E-01, 5.036686E-01, 6.428365E-02, 1.564468E+00, 4.677547E-02"
## [1] "Llh: -2.26998851140980E+03; Pars: 7.120500E-01, 4.986958E-01, 7.262174E-02, 1.558568E+00, 3.201839E-02"
## [1] "Llh: -2.26996984900162E+03; Pars: 7.106505E-01, 5.015320E-01, 6.877931E-02, 1.560611E+00, 3.933881E-02"
## [1] "Llh: -2.26997196053182E+03; Pars: 7.174484E-01, 4.968461E-01, 6.632439E-02, 1.561769E+00, 3.663639E-02"
## [1] "Llh: -2.26996461367588E+03; Pars: 7.140352E-01, 4.998604E-01, 6.600420E-02, 1.562003E+00, 4.040900E-02"
## [1] "Llh: -2.26996011640362E+03; Pars: 7.022896E-01, 5.021897E-01, 6.811624E-02, 1.559576E+00, 3.669214E-02"
## [1] "Llh: -2.26995916831048E+03; Pars: 7.061044E-01, 5.022739E-01, 6.709673E-02, 1.560699E+00, 3.963244E-02"
## [1] "Llh: -2.26994871710677E+03; Pars: 7.020721E-01, 4.974895E-01, 6.930972E-02, 1.562811E+00, 3.962465E-02"
## [1] "Llh: -2.26994347295579E+03; Pars: 6.958910E-01, 4.933657E-01, 7.193479E-02, 1.563818E+00, 3.751776E-02"
## [1] "Llh: -2.26993692985072E+03; Pars: 6.999522E-01, 4.984793E-01, 6.585218E-02, 1.564282E+00, 4.339587E-02"
## [1] "Llh: -2.26992403247434E+03; Pars: 6.946031E-01, 4.969599E-01, 6.438861E-02, 1.566121E+00, 4.542440E-02"
## [1] "Llh: -2.26993043702639E+03; Pars: 6.887947E-01, 4.989834E-01, 6.798105E-02, 1.564534E+00, 4.433184E-02"
## [1] "Llh: -2.26993626250024E+03; Pars: 6.951048E-01, 4.992025E-01, 6.748684E-02, 1.563901E+00, 4.335113E-02"
## [1] "Llh: -2.26993431942468E+03; Pars: 6.858230E-01, 4.968048E-01, 6.986264E-02, 1.559729E+00, 4.186810E-02"
## [1] "Llh: -2.26993428278181E+03; Pars: 6.916038E-01, 4.978935E-01, 6.870324E-02, 1.561335E+00, 4.222348E-02"
## [1] "Llh: -2.26992483655267E+03; Pars: 6.828646E-01, 4.939826E-01, 6.863354E-02, 1.565444E+00, 4.656170E-02"
## [1] "Llh: -2.26992785761262E+03; Pars: 6.886746E-01, 4.960425E-01, 6.824934E-02, 1.564257E+00, 4.482938E-02"
## [1] "Llh: -2.26991582776632E+03; Pars: 6.799729E-01, 4.891610E-01, 7.033852E-02, 1.568962E+00, 4.043580E-02"
## [1] "Llh: -2.26991711571880E+03; Pars: 6.691945E-01, 4.821907E-01, 7.234879E-02, 1.573688E+00, 3.765977E-02"
## [1] "Llh: -2.26990983695611E+03; Pars: 6.792447E-01, 4.974094E-01, 6.408319E-02, 1.566738E+00, 5.007313E-02"
## [1] "Llh: -2.26990555175193E+03; Pars: 6.709215E-01, 4.994437E-01, 6.015739E-02, 1.568200E+00, 5.635081E-02"
## [1] "Llh: -2.26991180467518E+03; Pars: 6.752589E-01, 4.934994E-01, 6.389640E-02, 1.571985E+00, 5.101834E-02"
## [1] "Llh: -2.26991159254778E+03; Pars: 6.793451E-01, 4.945942E-01, 6.509811E-02, 1.569316E+00, 4.881962E-02"
## [1] "Llh: -2.26989962044424E+03; Pars: 6.742882E-01, 4.906842E-01, 6.346542E-02, 1.570688E+00, 5.070510E-02"
## [1] "Llh: -2.26989416783066E+03; Pars: 6.670350E-01, 4.865865E-01, 6.120760E-02, 1.573774E+00, 5.389173E-02"
## [1] "Llh: -2.26989839547794E+03; Pars: 6.738865E-01, 4.926697E-01, 5.984256E-02, 1.573109E+00, 5.140725E-02"
## [1] "Llh: -2.26989970286855E+03; Pars: 6.761310E-01, 4.929976E-01, 6.204030E-02, 1.571190E+00, 5.019586E-02"
## [1] "Llh: -2.26989675389898E+03; Pars: 6.538613E-01, 4.880227E-01, 6.226906E-02, 1.575233E+00, 5.493769E-02"
## [1] "Llh: -2.26989576872476E+03; Pars: 6.640468E-01, 4.902418E-01, 6.279895E-02, 1.572950E+00, 5.255936E-02"
## [1] "Llh: -2.26991118443920E+03; Pars: 6.621210E-01, 4.962413E-01, 5.330333E-02, 1.573978E+00, 6.477571E-02"
## [1] "Llh: -2.26990079577777E+03; Pars: 6.665840E-01, 4.944617E-01, 5.756212E-02, 1.572723E+00, 5.869073E-02"
## [1] "Llh: -2.26989186777879E+03; Pars: 6.576444E-01, 4.907374E-01, 5.552933E-02, 1.574989E+00, 6.034033E-02"
## [1] "Llh: -2.26989865113524E+03; Pars: 6.467940E-01, 4.888202E-01, 5.074494E-02, 1.577833E+00, 6.610069E-02"
## [1] "Llh: -2.26989734948655E+03; Pars: 6.607572E-01, 4.825659E-01, 5.861884E-02, 1.578836E+00, 5.440495E-02"
## [1] "Llh: -2.26989369746915E+03; Pars: 6.632983E-01, 4.867311E-01, 5.900347E-02, 1.576170E+00, 5.489142E-02"
## [1] "Llh: -2.26988984625459E+03; Pars: 6.637804E-01, 4.843655E-01, 6.179064E-02, 1.575675E+00, 5.054531E-02"
## [1] "Llh: -2.26989267324610E+03; Pars: 6.623785E-01, 4.793949E-01, 6.390490E-02, 1.577152E+00, 4.647259E-02"
## [1] "Llh: -2.26989541998471E+03; Pars: 6.524354E-01, 4.828329E-01, 6.028944E-02, 1.576314E+00, 5.748401E-02"
## [1] "Llh: -2.26989196136741E+03; Pars: 6.577982E-01, 4.852735E-01, 6.017772E-02, 1.575512E+00, 5.596482E-02"
## [1] "Llh: -2.26989122525592E+03; Pars: 6.597757E-01, 4.832512E-01, 5.628456E-02, 1.577501E+00, 5.769408E-02"
## [1] "Llh: -2.26989031611876E+03; Pars: 6.608435E-01, 4.849894E-01, 5.791316E-02, 1.576362E+00, 5.641040E-02"
## [1] "Llh: -2.26988961913738E+03; Pars: 6.543108E-01, 4.862416E-01, 5.655813E-02, 1.577712E+00, 5.736918E-02"
## [1] "Llh: -2.26989250474427E+03; Pars: 6.479487E-01, 4.860692E-01, 5.423339E-02, 1.579685E+00, 5.910791E-02"
## [1] "Llh: -2.26988547391039E+03; Pars: 6.544526E-01, 4.859014E-01, 5.778412E-02, 1.575929E+00, 5.736060E-02"
## [1] "Llh: -2.26988268560577E+03; Pars: 6.500298E-01, 4.854871E-01, 5.717444E-02, 1.575809E+00, 5.859519E-02"
## [1] "Llh: -2.26988610853847E+03; Pars: 6.568454E-01, 4.874467E-01, 5.540856E-02, 1.576706E+00, 5.733934E-02"
## [1] "Llh: -2.26988626171445E+03; Pars: 6.570836E-01, 4.869025E-01, 5.660085E-02, 1.576407E+00, 5.699571E-02"
## [1] "Llh: -2.26988643407527E+03; Pars: 6.566796E-01, 4.807240E-01, 6.000864E-02, 1.577917E+00, 5.176343E-02"
## [1] "Llh: -2.26988566706540E+03; Pars: 6.569208E-01, 4.832080E-01, 5.888881E-02, 1.577184E+00, 5.390766E-02"
## [1] "Llh: -2.26988086723389E+03; Pars: 6.519114E-01, 4.857059E-01, 5.801508E-02, 1.576872E+00, 5.469227E-02"
## [1] "Llh: -2.26987696995855E+03; Pars: 6.474453E-01, 4.860646E-01, 5.806604E-02, 1.577127E+00, 5.383321E-02"
## [1] "Llh: -2.26988452292355E+03; Pars: 6.424405E-01, 4.870133E-01, 5.264775E-02, 1.578141E+00, 6.187252E-02"
## [1] "Llh: -2.26988287385288E+03; Pars: 6.477755E-01, 4.863500E-01, 5.493347E-02, 1.577524E+00, 5.904072E-02"
## [1] "Llh: -2.26987499777052E+03; Pars: 6.492959E-01, 4.851775E-01, 5.723040E-02, 1.576028E+00, 5.571727E-02"
## [1] "Llh: -2.26986925096457E+03; Pars: 6.467884E-01, 4.846463E-01, 5.756654E-02, 1.575187E+00, 5.489131E-02"
## [1] "Llh: -2.26987462929515E+03; Pars: 6.427386E-01, 4.828639E-01, 5.924316E-02, 1.576427E+00, 5.476789E-02"
## [1] "Llh: -2.26987601540585E+03; Pars: 6.462653E-01, 4.840055E-01, 5.828451E-02, 1.576496E+00, 5.541075E-02"
## [1] "Llh: -2.26986931481993E+03; Pars: 6.369903E-01, 4.869607E-01, 5.590465E-02, 1.575645E+00, 5.854367E-02"
## [1] "Llh: -2.26987235444153E+03; Pars: 6.419729E-01, 4.860198E-01, 5.665069E-02, 1.576030E+00, 5.738466E-02"
## [1] "Llh: -2.26986871090209E+03; Pars: 6.418215E-01, 4.840578E-01, 6.024845E-02, 1.574555E+00, 5.321178E-02"
## [1] "Llh: -2.26986846077826E+03; Pars: 6.388445E-01, 4.829157E-01, 6.290595E-02, 1.573072E+00, 5.029731E-02"
## [1] "Llh: -2.26986243065429E+03; Pars: 6.350930E-01, 4.838892E-01, 6.030009E-02, 1.575174E+00, 5.033817E-02"
## [1] "Llh: -2.26986107633402E+03; Pars: 6.276246E-01, 4.830922E-01, 6.186291E-02, 1.574856E+00, 4.620966E-02"
## [1] "Llh: -2.26986231349326E+03; Pars: 6.297492E-01, 4.821299E-01, 6.092724E-02, 1.572950E+00, 5.205073E-02"
## [1] "Llh: -2.26986267472879E+03; Pars: 6.341732E-01, 4.831105E-01, 6.021194E-02, 1.573993E+00, 5.249635E-02"
## [1] "Llh: -2.26985569297107E+03; Pars: 6.292602E-01, 4.850304E-01, 6.042375E-02, 1.572260E+00, 5.002918E-02"
## [1] "Llh: -2.26985423566581E+03; Pars: 6.225210E-01, 4.861173E-01, 6.101405E-02, 1.570180E+00, 4.765982E-02"
## [1] "Llh: -2.26986991030554E+03; Pars: 6.292208E-01, 4.806164E-01, 6.580603E-02, 1.570855E+00, 4.189986E-02"
## [1] "Llh: -2.26986181456213E+03; Pars: 6.350479E-01, 4.853669E-01, 5.837999E-02, 1.574446E+00, 5.438272E-02"
## [1] "Llh: -2.26986432018329E+03; Pars: 6.147265E-01, 4.831987E-01, 6.446952E-02, 1.571016E+00, 4.534879E-02"
## [1] "Llh: -2.26985955568089E+03; Pars: 6.227420E-01, 4.835602E-01, 6.274377E-02, 1.572057E+00, 4.773442E-02"
## [1] "Llh: -2.26985236019552E+03; Pars: 6.162293E-01, 4.851890E-01, 5.906524E-02, 1.572722E+00, 4.891762E-02"
## [1] "Llh: -2.26985387029207E+03; Pars: 6.049218E-01, 4.863297E-01, 5.714489E-02, 1.572547E+00, 4.822777E-02"
## [1] "Llh: -2.26985403004513E+03; Pars: 6.199167E-01, 4.872110E-01, 6.029914E-02, 1.572753E+00, 4.591097E-02"
## [1] "Llh: -2.26985386758761E+03; Pars: 6.223749E-01, 4.859357E-01, 6.045617E-02, 1.572802E+00, 4.744591E-02"
## [1] "Llh: -2.26986713875932E+03; Pars: 6.095488E-01, 4.841885E-01, 6.367686E-02, 1.570602E+00, 4.080425E-02"
## [1] "Llh: -2.26985608475225E+03; Pars: 6.286731E-01, 4.850720E-01, 5.970421E-02, 1.573484E+00, 5.098810E-02"
## [1] "Llh: -2.26985217987116E+03; Pars: 6.173916E-01, 4.872648E-01, 5.933047E-02, 1.569646E+00, 5.088869E-02"
## [1] "Llh: -2.26985537127543E+03; Pars: 6.122750E-01, 4.893646E-01, 5.806425E-02, 1.567047E+00, 5.322821E-02"
## [1] "Llh: -2.26984904163297E+03; Pars: 6.201340E-01, 4.882816E-01, 5.708428E-02, 1.571475E+00, 5.062564E-02"
## [1] "Llh: -2.26984922089807E+03; Pars: 6.188300E-01, 4.906595E-01, 5.425453E-02, 1.571183E+00, 5.207125E-02"
## [1] "Llh: -2.26985340360765E+03; Pars: 6.107872E-01, 4.880455E-01, 5.907587E-02, 1.569247E+00, 4.722697E-02"
## [1] "Llh: -2.26985152206899E+03; Pars: 6.152587E-01, 4.873004E-01, 5.923296E-02, 1.570306E+00, 4.816726E-02"
## [1] "Llh: -2.26984931635582E+03; Pars: 6.140343E-01, 4.874697E-01, 5.705360E-02, 1.572600E+00, 5.075822E-02"
## [1] "Llh: -2.26984982602654E+03; Pars: 6.161560E-01, 4.871313E-01, 5.804371E-02, 1.571995E+00, 4.998362E-02"
## [1] "Llh: -2.26984937549993E+03; Pars: 6.108443E-01, 4.882671E-01, 5.625045E-02, 1.569897E+00, 5.229707E-02"
## [1] "Llh: -2.26984929260309E+03; Pars: 6.137269E-01, 4.876832E-01, 5.730188E-02, 1.570623E+00, 5.108428E-02"
## [1] "Llh: -2.26984840779602E+03; Pars: 6.159889E-01, 4.900225E-01, 5.693603E-02, 1.569139E+00, 5.169201E-02"
## [1] "Llh: -2.26984913749820E+03; Pars: 6.158686E-01, 4.924573E-01, 5.587143E-02, 1.567351E+00, 5.307921E-02"
## [1] "Llh: -2.26985007522607E+03; Pars: 6.142656E-01, 4.890378E-01, 5.571303E-02, 1.572011E+00, 5.004227E-02"
## [1] "Llh: -2.26984896360170E+03; Pars: 6.150471E-01, 4.885939E-01, 5.661739E-02, 1.571419E+00, 5.025388E-02"
## [1] "Llh: -2.26984778571840E+03; Pars: 6.163138E-01, 4.895208E-01, 5.476431E-02, 1.571796E+00, 5.359836E-02"
## [1] "Llh: -2.26984928593622E+03; Pars: 6.168414E-01, 4.906349E-01, 5.252999E-02, 1.572542E+00, 5.631391E-02"
## [1] "Llh: -2.26984827361526E+03; Pars: 6.184499E-01, 4.901734E-01, 5.602796E-02, 1.569183E+00, 5.214344E-02"
## [1] "Llh: -2.26984804513503E+03; Pars: 6.173460E-01, 4.894961E-01, 5.628437E-02, 1.570036E+00, 5.179714E-02"
## [1] "Llh: -2.26984788780714E+03; Pars: 6.202050E-01, 4.906865E-01, 5.537268E-02, 1.570923E+00, 5.210253E-02"
## [1] "Llh: -2.26984777459802E+03; Pars: 6.185855E-01, 4.899340E-01, 5.585498E-02, 1.570848E+00, 5.184797E-02"
## [1] "Llh: -2.26984731260793E+03; Pars: 6.131785E-01, 4.907480E-01, 5.509855E-02, 1.569821E+00, 5.305010E-02"
## [1] "Llh: -2.26984769616398E+03; Pars: 6.097007E-01, 4.919859E-01, 5.410569E-02, 1.568994E+00, 5.426233E-02"
## [1] "Llh: -2.26984765295916E+03; Pars: 6.175180E-01, 4.912979E-01, 5.495791E-02, 1.569237E+00, 5.454035E-02"
## [1] "Llh: -2.26984737114015E+03; Pars: 6.169003E-01, 4.906205E-01, 5.537278E-02, 1.569782E+00, 5.346873E-02"
## [1] "Llh: -2.26984812699186E+03; Pars: 6.169408E-01, 4.901046E-01, 5.401397E-02, 1.571775E+00, 5.381291E-02"
## [1] "Llh: -2.26984754642070E+03; Pars: 6.167028E-01, 4.900841E-01, 5.474448E-02, 1.571116E+00, 5.328268E-02"
## [1] "Llh: -2.26984722951157E+03; Pars: 6.153263E-01, 4.908674E-01, 5.404967E-02, 1.571309E+00, 5.430200E-02"
## [1] "Llh: -2.26984778229407E+03; Pars: 6.143164E-01, 4.915546E-01, 5.293232E-02, 1.571945E+00, 5.555443E-02"
## [1] "Llh: -2.26984733091358E+03; Pars: 6.159635E-01, 4.913823E-01, 5.528387E-02, 1.569354E+00, 5.278224E-02"
## [1] "Llh: -2.26984718361663E+03; Pars: 6.160511E-01, 4.909162E-01, 5.515398E-02, 1.569964E+00, 5.298627E-02"
## [1] "Llh: -2.26984720324635E+03; Pars: 6.126781E-01, 4.913614E-01, 5.391281E-02, 1.569949E+00, 5.498795E-02"
## [1] "Llh: -2.26984708170080E+03; Pars: 6.141549E-01, 4.910042E-01, 5.439835E-02, 1.570173E+00, 5.420295E-02"
## [1] "Llh: -2.26984714734094E+03; Pars: 6.135416E-01, 4.915795E-01, 5.488485E-02, 1.569305E+00, 5.392134E-02"
## [1] "Llh: -2.26984704265181E+03; Pars: 6.143319E-01, 4.912053E-01, 5.484976E-02, 1.569757E+00, 5.376167E-02"
## [1] "Llh: -2.26984701670248E+03; Pars: 6.123168E-01, 4.912761E-01, 5.404735E-02, 1.570627E+00, 5.385246E-02"
## [1] "Llh: -2.26984734531029E+03; Pars: 6.100251E-01, 4.916041E-01, 5.338463E-02, 1.571050E+00, 5.404433E-02"
## [1] "Llh: -2.26984708119428E+03; Pars: 6.156939E-01, 4.913598E-01, 5.390109E-02, 1.570912E+00, 5.459204E-02"
## [1] "Llh: -2.26984697268189E+03; Pars: 6.150651E-01, 4.912068E-01, 5.420046E-02, 1.570639E+00, 5.420656E-02"
## [1] "Llh: -2.26984746479855E+03; Pars: 6.134417E-01, 4.913761E-01, 5.501029E-02, 1.569156E+00, 5.330197E-02"
## [1] "Llh: -2.26984701015009E+03; Pars: 6.148551E-01, 4.909945E-01, 5.428982E-02, 1.570770E+00, 5.405199E-02"
## [1] "Llh: -2.26984702016617E+03; Pars: 6.122385E-01, 4.913586E-01, 5.356031E-02, 1.570822E+00, 5.504399E-02"
## [1] "Llh: -2.26984694799703E+03; Pars: 6.131916E-01, 4.912479E-01, 5.395873E-02, 1.570608E+00, 5.452956E-02"
## [1] "Llh: -2.26984686546594E+03; Pars: 6.137493E-01, 4.913681E-01, 5.414010E-02, 1.570787E+00, 5.395795E-02"
## [1] "Llh: -2.26984689979275E+03; Pars: 6.135465E-01, 4.915502E-01, 5.401097E-02, 1.571094E+00, 5.383544E-02"
## [1] "Llh: -2.26984749136041E+03; Pars: 6.133393E-01, 4.912321E-01, 5.340482E-02, 1.571616E+00, 5.447773E-02"
## [1] "Llh: -2.26984690778040E+03; Pars: 6.140838E-01, 4.912120E-01, 5.448852E-02, 1.570222E+00, 5.394069E-02"
## [1] "Llh: -2.26984698212992E+03; Pars: 6.160611E-01, 4.911357E-01, 5.438371E-02, 1.570583E+00, 5.442223E-02"
## [1] "Llh: -2.26984692823129E+03; Pars: 6.151251E-01, 4.911708E-01, 5.429962E-02, 1.570594E+00, 5.427979E-02"
## [1] "Llh: -2.26984681633783E+03; Pars: 6.136308E-01, 4.914878E-01, 5.414515E-02, 1.570370E+00, 5.431382E-02"
## [1] "Llh: -2.26984676330868E+03; Pars: 6.130186E-01, 4.917346E-01, 5.407281E-02, 1.570169E+00, 5.444474E-02"
## [1] "Llh: -2.26984677667722E+03; Pars: 6.126023E-01, 4.914865E-01, 5.418345E-02, 1.570313E+00, 5.425453E-02"
## [1] "Llh: -2.26984680788143E+03; Pars: 6.132180E-01, 4.914166E-01, 5.418770E-02, 1.570395E+00, 5.424254E-02"
## [1] "Llh: -2.26984672679392E+03; Pars: 6.142400E-01, 4.915408E-01, 5.451507E-02, 1.570226E+00, 5.382152E-02"
## [1] "Llh: -2.26984667960976E+03; Pars: 6.147641E-01, 4.916873E-01, 5.479324E-02, 1.570035E+00, 5.346751E-02"
## [1] "Llh: -2.26984670472379E+03; Pars: 6.121622E-01, 4.918248E-01, 5.437163E-02, 1.570017E+00, 5.374638E-02"
## [1] "Llh: -2.26984671944262E+03; Pars: 6.129029E-01, 4.916612E-01, 5.435363E-02, 1.570161E+00, 5.387973E-02"
## [1] "Llh: -2.26984654500120E+03; Pars: 6.124348E-01, 4.920289E-01, 5.413596E-02, 1.570307E+00, 5.400775E-02"
## [1] "Llh: -2.26984640106256E+03; Pars: 6.116104E-01, 4.924378E-01, 5.395968E-02, 1.570350E+00, 5.404128E-02"
## [1] "Llh: -2.26984662273693E+03; Pars: 6.119137E-01, 4.923005E-01, 5.441223E-02, 1.569567E+00, 5.402383E-02"
## [1] "Llh: -2.26984659895831E+03; Pars: 6.123726E-01, 4.920673E-01, 5.434419E-02, 1.569872E+00, 5.400736E-02"
## [1] "Llh: -2.26984646367852E+03; Pars: 6.129689E-01, 4.924145E-01, 5.443317E-02, 1.569864E+00, 5.362837E-02"
## [1] "Llh: -2.26984652332638E+03; Pars: 6.128773E-01, 4.921823E-01, 5.437074E-02, 1.569976E+00, 5.378491E-02"
## [1] "Llh: -2.26984641515971E+03; Pars: 6.125327E-01, 4.924381E-01, 5.468796E-02, 1.569886E+00, 5.311162E-02"
## [1] "Llh: -2.26984646496089E+03; Pars: 6.126542E-01, 4.922621E-01, 5.453417E-02, 1.569957E+00, 5.344490E-02"
## [1] "Llh: -2.26984628791125E+03; Pars: 6.135373E-01, 4.925933E-01, 5.451566E-02, 1.569986E+00, 5.355608E-02"
## [1] "Llh: -2.26984613834443E+03; Pars: 6.142249E-01, 4.929780E-01, 5.458768E-02, 1.569970E+00, 5.346093E-02"
## [1] "Llh: -2.26984613970813E+03; Pars: 6.107197E-01, 4.932480E-01, 5.401184E-02, 1.569941E+00, 5.383232E-02"
## [1] "Llh: -2.26984623978072E+03; Pars: 6.117308E-01, 4.928574E-01, 5.420719E-02, 1.569965E+00, 5.374112E-02"
## [1] "Llh: -2.26984609496794E+03; Pars: 6.124500E-01, 4.933399E-01, 5.432794E-02, 1.570133E+00, 5.322245E-02"
## [1] "Llh: -2.26984605752519E+03; Pars: 6.124887E-01, 4.939775E-01, 5.431981E-02, 1.570263E+00, 5.283000E-02"
## [1] "Llh: -2.26984594213927E+03; Pars: 6.116616E-01, 4.936174E-01, 5.419362E-02, 1.570300E+00, 5.328209E-02"
## [1] "Llh: -2.26984582645608E+03; Pars: 6.110079E-01, 4.942199E-01, 5.407385E-02, 1.570518E+00, 5.310894E-02"
## [1] "Llh: -2.26984586926243E+03; Pars: 6.114879E-01, 4.943073E-01, 5.369319E-02, 1.570531E+00, 5.379777E-02"
## [1] "Llh: -2.26984591725437E+03; Pars: 6.117491E-01, 4.938393E-01, 5.394188E-02, 1.570370E+00, 5.362623E-02"
## [1] "Llh: -2.26984570613932E+03; Pars: 6.123612E-01, 4.950574E-01, 5.431486E-02, 1.570140E+00, 5.277070E-02"
## [1] "Llh: -2.26984575806726E+03; Pars: 6.127367E-01, 4.963724E-01, 5.449245E-02, 1.570035E+00, 5.213541E-02"
## [1] "Llh: -2.26984605267524E+03; Pars: 6.139086E-01, 4.949686E-01, 5.438392E-02, 1.570628E+00, 5.255502E-02"
## [1] "Llh: -2.26984589155497E+03; Pars: 6.131114E-01, 4.945379E-01, 5.429090E-02, 1.570456E+00, 5.287434E-02"
## [1] "Llh: -2.26984613772047E+03; Pars: 6.099580E-01, 4.958659E-01, 5.368936E-02, 1.570793E+00, 5.269177E-02"
## [1] "Llh: -2.26984591579422E+03; Pars: 6.110247E-01, 4.951424E-01, 5.391394E-02, 1.570587E+00, 5.288406E-02"
## [1] "Llh: -2.26984562380871E+03; Pars: 6.111086E-01, 4.953291E-01, 5.379489E-02, 1.570630E+00, 5.334433E-02"
## [1] "Llh: -2.26984550825099E+03; Pars: 6.104186E-01, 4.960063E-01, 5.353242E-02, 1.570814E+00, 5.360150E-02"
## [1] "Llh: -2.26984560165649E+03; Pars: 6.123301E-01, 4.945085E-01, 5.404815E-02, 1.570396E+00, 5.357724E-02"
## [1] "Llh: -2.26984564431384E+03; Pars: 6.120038E-01, 4.946669E-01, 5.401460E-02, 1.570444E+00, 5.340395E-02"
## [1] "Llh: -2.26984548100992E+03; Pars: 6.099310E-01, 4.951012E-01, 5.357409E-02, 1.570503E+00, 5.386812E-02"
## [1] "Llh: -2.26984538202388E+03; Pars: 6.083408E-01, 4.953830E-01, 5.321569E-02, 1.570527E+00, 5.436501E-02"
## [1] "Llh: -2.26984523009319E+03; Pars: 6.102955E-01, 4.957630E-01, 5.398080E-02, 1.570427E+00, 5.317158E-02"
## [1] "Llh: -2.26984496178184E+03; Pars: 6.096993E-01, 4.964925E-01, 5.412461E-02, 1.570375E+00, 5.285849E-02"
## [1] "Llh: -2.26984500677237E+03; Pars: 6.102521E-01, 4.967615E-01, 5.362044E-02, 1.570383E+00, 5.376023E-02"
## [1] "Llh: -2.26984516774755E+03; Pars: 6.104410E-01, 4.961249E-01, 5.373380E-02, 1.570416E+00, 5.359741E-02"
## [1] "Llh: -2.26984496754453E+03; Pars: 6.080551E-01, 4.966032E-01, 5.310166E-02, 1.570858E+00, 5.449429E-02"
## [1] "Llh: -2.26984506960830E+03; Pars: 6.091316E-01, 4.962163E-01, 5.340496E-02, 1.570678E+00, 5.406339E-02"
## [1] "Llh: -2.26984482872678E+03; Pars: 6.063762E-01, 4.979957E-01, 5.298978E-02, 1.570786E+00, 5.405456E-02"
## [1] "Llh: -2.26984478778154E+03; Pars: 6.033992E-01, 4.997486E-01, 5.246060E-02, 1.570981E+00, 5.429323E-02"
## [1] "Llh: -2.26984434577850E+03; Pars: 6.054800E-01, 4.979869E-01, 5.307678E-02, 1.570435E+00, 5.430700E-02"
## [1] "Llh: -2.26984388585019E+03; Pars: 6.030107E-01, 4.989801E-01, 5.284896E-02, 1.570246E+00, 5.465975E-02"
## [1] "Llh: -2.26984410856223E+03; Pars: 6.054258E-01, 5.000585E-01, 5.324682E-02, 1.570610E+00, 5.366139E-02"
## [1] "Llh: -2.26984430768876E+03; Pars: 6.061546E-01, 4.988855E-01, 5.323904E-02, 1.570589E+00, 5.383729E-02"
## [1] "Llh: -2.26984397178189E+03; Pars: 6.015840E-01, 4.999922E-01, 5.269262E-02, 1.570845E+00, 5.422662E-02"
## [1] "Llh: -2.26984416053374E+03; Pars: 6.037510E-01, 4.991826E-01, 5.292458E-02, 1.570729E+00, 5.411003E-02"
## [1] "Llh: -2.26984373982931E+03; Pars: 6.011926E-01, 5.015140E-01, 5.304778E-02, 1.570365E+00, 5.338550E-02"
## [1] "Llh: -2.26984357259401E+03; Pars: 5.977613E-01, 5.039876E-01, 5.302084E-02, 1.570118E+00, 5.283111E-02"
## [1] "Llh: -2.26984378401555E+03; Pars: 5.947732E-01, 5.046413E-01, 5.158333E-02, 1.570745E+00, 5.501035E-02"
## [1] "Llh: -2.26984372727454E+03; Pars: 5.985047E-01, 5.025917E-01, 5.221865E-02, 1.570653E+00, 5.447238E-02"
## [1] "Llh: -2.26984274624552E+03; Pars: 5.991154E-01, 5.024923E-01, 5.315055E-02, 1.570008E+00, 5.364728E-02"
## [1] "Llh: -2.26984196346267E+03; Pars: 5.969735E-01, 5.038699E-01, 5.349553E-02, 1.569521E+00, 5.332430E-02"
## [1] "Llh: -2.26984272767334E+03; Pars: 5.937079E-01, 5.037084E-01, 5.246381E-02, 1.569944E+00, 5.414428E-02"
## [1] "Llh: -2.26984293685578E+03; Pars: 5.966373E-01, 5.027934E-01, 5.265956E-02, 1.570110E+00, 5.402356E-02"
## [1] "Llh: -2.26984236219963E+03; Pars: 5.943992E-01, 5.052695E-01, 5.292649E-02, 1.569348E+00, 5.354611E-02"
## [1] "Llh: -2.26984262972720E+03; Pars: 5.961954E-01, 5.039450E-01, 5.286802E-02, 1.569722E+00, 5.371624E-02"
## [1] "Llh: -2.26984259563302E+03; Pars: 5.895278E-01, 5.088375E-01, 5.280117E-02, 1.569587E+00, 5.266752E-02"
## [1] "Llh: -2.26984253773337E+03; Pars: 5.928986E-01, 5.063550E-01, 5.281312E-02, 1.569752E+00, 5.316558E-02"
## [1] "Llh: -2.26984196792282E+03; Pars: 5.917915E-01, 5.066908E-01, 5.366927E-02, 1.568821E+00, 5.233217E-02"
## [1] "Llh: -2.26984216948624E+03; Pars: 5.934698E-01, 5.056629E-01, 5.330661E-02, 1.569279E+00, 5.286722E-02"
## [1] "Llh: -2.26984134285719E+03; Pars: 5.901469E-01, 5.063696E-01, 5.312645E-02, 1.568836E+00, 5.377386E-02"
## [1] "Llh: -2.26984100126589E+03; Pars: 5.863397E-01, 5.075649E-01, 5.317925E-02, 1.568196E+00, 5.424524E-02"
## [1] "Llh: -2.26984121147502E+03; Pars: 5.912531E-01, 5.081984E-01, 5.396965E-02, 1.568312E+00, 5.250108E-02"
## [1] "Llh: -2.26984140562475E+03; Pars: 5.918668E-01, 5.070722E-01, 5.359319E-02, 1.568720E+00, 5.291188E-02"
## [1] "Llh: -2.26984132557550E+03; Pars: 5.914042E-01, 5.062775E-01, 5.408296E-02, 1.567928E+00, 5.321398E-02"
## [1] "Llh: -2.26984129752598E+03; Pars: 5.917778E-01, 5.062969E-01, 5.376550E-02, 1.568384E+00, 5.320188E-02"
## [1] "Llh: -2.26984069415329E+03; Pars: 5.888550E-01, 5.077775E-01, 5.430518E-02, 1.567946E+00, 5.269576E-02"
## [1] "Llh: -2.26984055071041E+03; Pars: 5.860830E-01, 5.090362E-01, 5.499453E-02, 1.567246E+00, 5.227058E-02"
## [1] "Llh: -2.26984055941823E+03; Pars: 5.891794E-01, 5.072895E-01, 5.409251E-02, 1.567842E+00, 5.388506E-02"
## [1] "Llh: -2.26984065235877E+03; Pars: 5.898324E-01, 5.071397E-01, 5.398670E-02, 1.568087E+00, 5.349684E-02"
## [1] "Llh: -2.26984114283436E+03; Pars: 5.808797E-01, 5.115116E-01, 5.450505E-02, 1.566472E+00, 5.311724E-02"
## [1] "Llh: -2.26984072443227E+03; Pars: 5.849032E-01, 5.095904E-01, 5.425267E-02, 1.567234E+00, 5.316900E-02"
## [1] "Llh: -2.26984025865653E+03; Pars: 5.833255E-01, 5.103816E-01, 5.442995E-02, 1.567148E+00, 5.322651E-02"
## [1] "Llh: -2.26984027066531E+03; Pars: 5.790994E-01, 5.124362E-01, 5.476217E-02, 1.566531E+00, 5.323882E-02"
## [1] "Llh: -2.26984136685673E+03; Pars: 5.806792E-01, 5.093445E-01, 5.440991E-02, 1.566755E+00, 5.421748E-02"
## [1] "Llh: -2.26984066569516E+03; Pars: 5.886096E-01, 5.084847E-01, 5.407972E-02, 1.567923E+00, 5.293018E-02"
## [1] "Llh: -2.26984055290969E+03; Pars: 5.865005E-01, 5.103497E-01, 5.556050E-02, 1.566761E+00, 5.194730E-02"
## [1] "Llh: -2.26984042928719E+03; Pars: 5.864603E-01, 5.096521E-01, 5.496519E-02, 1.567120E+00, 5.252178E-02"
## [1] "Llh: -2.26984014614632E+03; Pars: 5.885600E-01, 5.083458E-01, 5.477209E-02, 1.567678E+00, 5.276464E-02"
## [1] "Llh: -2.26983997221694E+03; Pars: 5.903884E-01, 5.077247E-01, 5.503180E-02, 1.567900E+00, 5.256247E-02"
## [1] "Llh: -2.26984033419597E+03; Pars: 5.855650E-01, 5.091465E-01, 5.532588E-02, 1.566980E+00, 5.285638E-02"
## [1] "Llh: -2.26984024523345E+03; Pars: 5.863261E-01, 5.089810E-01, 5.501434E-02, 1.567215E+00, 5.287483E-02"
## [1] "Llh: -2.26984027931067E+03; Pars: 5.838540E-01, 5.110260E-01, 5.568181E-02, 1.566809E+00, 5.149741E-02"
## [1] "Llh: -2.26984020731567E+03; Pars: 5.851853E-01, 5.100893E-01, 5.528448E-02, 1.567068E+00, 5.209432E-02"
## [1] "Llh: -2.26983984809497E+03; Pars: 5.865913E-01, 5.096937E-01, 5.489577E-02, 1.567335E+00, 5.304138E-02"
## [1] "Llh: -2.26983961493457E+03; Pars: 5.868455E-01, 5.100228E-01, 5.484639E-02, 1.567379E+00, 5.342678E-02"
## [1] "Llh: -2.26983958096359E+03; Pars: 5.863680E-01, 5.092259E-01, 5.487759E-02, 1.567564E+00, 5.315218E-02"
## [1] "Llh: -2.26983924479674E+03; Pars: 5.863218E-01, 5.090129E-01, 5.483380E-02, 1.567787E+00, 5.346738E-02"
## [1] "Llh: -2.26983958318198E+03; Pars: 5.907013E-01, 5.079522E-01, 5.557438E-02, 1.567791E+00, 5.254380E-02"
## [1] "Llh: -2.26983964426286E+03; Pars: 5.888574E-01, 5.085584E-01, 5.528827E-02, 1.567631E+00, 5.271448E-02"
## [1] "Llh: -2.26983919417754E+03; Pars: 5.894508E-01, 5.089379E-01, 5.521400E-02, 1.567954E+00, 5.276307E-02"
## [1] "Llh: -2.26983892265143E+03; Pars: 5.910131E-01, 5.089164E-01, 5.531384E-02, 1.568324E+00, 5.270718E-02"
## [1] "Llh: -2.26983923878202E+03; Pars: 5.929228E-01, 5.073646E-01, 5.495560E-02, 1.568605E+00, 5.378872E-02"
## [1] "Llh: -2.26983922266931E+03; Pars: 5.909884E-01, 5.080444E-01, 5.503782E-02, 1.568221E+00, 5.336512E-02"
## [1] "Llh: -2.26983872652767E+03; Pars: 5.879597E-01, 5.098559E-01, 5.521069E-02, 1.567901E+00, 5.364164E-02"
## [1] "Llh: -2.26983845041451E+03; Pars: 5.867454E-01, 5.109248E-01, 5.530014E-02, 1.567902E+00, 5.418123E-02"
## [1] "Llh: -2.26983858848489E+03; Pars: 5.914625E-01, 5.079174E-01, 5.557759E-02, 1.568631E+00, 5.307911E-02"
## [1] "Llh: -2.26983873074867E+03; Pars: 5.903083E-01, 5.084430E-01, 5.539479E-02, 1.568318E+00, 5.316602E-02"
## [1] "Llh: -2.26983818567854E+03; Pars: 5.879112E-01, 5.099739E-01, 5.485090E-02, 1.568554E+00, 5.417621E-02"
## [1] "Llh: -2.26983789737912E+03; Pars: 5.865161E-01, 5.109878E-01, 5.448916E-02, 1.568936E+00, 5.499241E-02"
## [1] "Llh: -2.26983846518337E+03; Pars: 5.923684E-01, 5.097001E-01, 5.545362E-02, 1.569019E+00, 5.386264E-02"
## [1] "Llh: -2.26983835936263E+03; Pars: 5.908567E-01, 5.095282E-01, 5.529866E-02, 1.568711E+00, 5.376382E-02"
## [1] "Llh: -2.26983765496740E+03; Pars: 5.876492E-01, 5.112678E-01, 5.535394E-02, 1.568780E+00, 5.412438E-02"
## [1] "Llh: -2.26983732106318E+03; Pars: 5.859796E-01, 5.128872E-01, 5.551200E-02, 1.569060E+00, 5.450401E-02"
## [1] "Llh: -2.26983796863253E+03; Pars: 5.856110E-01, 5.119810E-01, 5.515718E-02, 1.568972E+00, 5.550105E-02"
## [1] "Llh: -2.26983783644133E+03; Pars: 5.869616E-01, 5.112131E-01, 5.519635E-02, 1.568810E+00, 5.480258E-02"
## [1] "Llh: -2.26983791293561E+03; Pars: 5.833612E-01, 5.143169E-01, 5.474093E-02, 1.568737E+00, 5.581851E-02"
## [1] "Llh: -2.26983777869769E+03; Pars: 5.853865E-01, 5.127095E-01, 5.495009E-02, 1.568710E+00, 5.513366E-02"
## [1] "Llh: -2.26983774603486E+03; Pars: 5.875348E-01, 5.120031E-01, 5.487837E-02, 1.569790E+00, 5.509736E-02"
## [1] "Llh: -2.26983766367383E+03; Pars: 5.873375E-01, 5.117333E-01, 5.498381E-02, 1.569317E+00, 5.486833E-02"
## [1] "Llh: -2.26983744684956E+03; Pars: 5.820157E-01, 5.142941E-01, 5.475390E-02, 1.569223E+00, 5.595657E-02"
## [1] "Llh: -2.26983747987267E+03; Pars: 5.842260E-01, 5.130985E-01, 5.489009E-02, 1.569095E+00, 5.540838E-02"
## [1] "Llh: -2.26983760933512E+03; Pars: 5.845562E-01, 5.141497E-01, 5.566930E-02, 1.569112E+00, 5.511365E-02"
## [1] "Llh: -2.26983752467095E+03; Pars: 5.850462E-01, 5.133574E-01, 5.537426E-02, 1.569068E+00, 5.508334E-02"
## [1] "Llh: -2.26983735448722E+03; Pars: 5.833447E-01, 5.147843E-01, 5.503328E-02, 1.569342E+00, 5.541578E-02"
## [1] "Llh: -2.26983739635855E+03; Pars: 5.842489E-01, 5.138892E-01, 5.507405E-02, 1.569209E+00, 5.526248E-02"
## [1] "Llh: -2.26983723497554E+03; Pars: 5.841029E-01, 5.141118E-01, 5.531280E-02, 1.569694E+00, 5.519755E-02"
## [1] "Llh: -2.26983726773052E+03; Pars: 5.834611E-01, 5.148144E-01, 5.549416E-02, 1.570187E+00, 5.522949E-02"
## [1] "Llh: -2.26983724044106E+03; Pars: 5.808582E-01, 5.160488E-01, 5.541069E-02, 1.569238E+00, 5.559457E-02"
## [1] "Llh: -2.26983725004245E+03; Pars: 5.824780E-01, 5.149665E-01, 5.530397E-02, 1.569257E+00, 5.541301E-02"
## [1] "Llh: -2.26983705633659E+03; Pars: 5.814742E-01, 5.154932E-01, 5.503480E-02, 1.569555E+00, 5.558405E-02"
## [1] "Llh: -2.26983693796842E+03; Pars: 5.796882E-01, 5.165644E-01, 5.486507E-02, 1.569798E+00, 5.583441E-02"
## [1] "Llh: -2.26983703198760E+03; Pars: 5.835737E-01, 5.154617E-01, 5.569964E-02, 1.569630E+00, 5.466195E-02"
## [1] "Llh: -2.26983706190711E+03; Pars: 5.831842E-01, 5.151696E-01, 5.546320E-02, 1.569528E+00, 5.498561E-02"
## [1] "Llh: -2.26983681069883E+03; Pars: 5.823364E-01, 5.152418E-01, 5.568680E-02, 1.569626E+00, 5.490121E-02"
## [1] "Llh: -2.26983659665773E+03; Pars: 5.818322E-01, 5.154708E-01, 5.601356E-02, 1.569768E+00, 5.464393E-02"
## [1] "Llh: -2.26983716299065E+03; Pars: 5.780425E-01, 5.181881E-01, 5.540871E-02, 1.570191E+00, 5.586896E-02"
## [1] "Llh: -2.26983698550601E+03; Pars: 5.800268E-01, 5.168578E-01, 5.543453E-02, 1.569908E+00, 5.552772E-02"
## [1] "Llh: -2.26983681632330E+03; Pars: 5.828314E-01, 5.153362E-01, 5.551955E-02, 1.570282E+00, 5.475165E-02"
## [1] "Llh: -2.26983683362762E+03; Pars: 5.823381E-01, 5.155143E-01, 5.549234E-02, 1.570021E+00, 5.496238E-02"
## [1] "Llh: -2.26983650856964E+03; Pars: 5.790780E-01, 5.177703E-01, 5.570014E-02, 1.570061E+00, 5.497031E-02"
## [1] "Llh: -2.26983631540355E+03; Pars: 5.765655E-01, 5.196093E-01, 5.589380E-02, 1.570244E+00, 5.485670E-02"
## [1] "Llh: -2.26983641052949E+03; Pars: 5.768040E-01, 5.180724E-01, 5.539097E-02, 1.570370E+00, 5.558381E-02"
## [1] "Llh: -2.26983651019068E+03; Pars: 5.784964E-01, 5.174185E-01, 5.546814E-02, 1.570185E+00, 5.535334E-02"
## [1] "Llh: -2.26983612601770E+03; Pars: 5.790617E-01, 5.171584E-01, 5.563865E-02, 1.570277E+00, 5.474048E-02"
## [1] "Llh: -2.26983578908970E+03; Pars: 5.785792E-01, 5.173087E-01, 5.574072E-02, 1.570461E+00, 5.434686E-02"
## [1] "Llh: -2.26983599715289E+03; Pars: 5.789567E-01, 5.177502E-01, 5.655837E-02, 1.570652E+00, 5.383877E-02"
## [1] "Llh: -2.26983610821474E+03; Pars: 5.791396E-01, 5.174535E-01, 5.613504E-02, 1.570439E+00, 5.433768E-02"
## [1] "Llh: -2.26983558166775E+03; Pars: 5.742637E-01, 5.199552E-01, 5.631941E-02, 1.570316E+00, 5.455637E-02"
## [1] "Llh: -2.26983519414952E+03; Pars: 5.699799E-01, 5.222802E-01, 5.671934E-02, 1.570333E+00, 5.445872E-02"
## [1] "Llh: -2.26983561572950E+03; Pars: 5.705219E-01, 5.225554E-01, 5.610772E-02, 1.571056E+00, 5.459001E-02"
## [1] "Llh: -2.26983565698935E+03; Pars: 5.733495E-01, 5.207752E-01, 5.608418E-02, 1.570734E+00, 5.460349E-02"
## [1] "Llh: -2.26983525734637E+03; Pars: 5.730373E-01, 5.217263E-01, 5.701701E-02, 1.570728E+00, 5.325262E-02"
## [1] "Llh: -2.26983541520855E+03; Pars: 5.739790E-01, 5.208104E-01, 5.661050E-02, 1.570638E+00, 5.383542E-02"
## [1] "Llh: -2.26983470960695E+03; Pars: 5.718645E-01, 5.210298E-01, 5.696346E-02, 1.571048E+00, 5.333810E-02"
## [1] "Llh: -2.26983414578445E+03; Pars: 5.695140E-01, 5.217415E-01, 5.749829E-02, 1.571451E+00, 5.257880E-02"
## [1] "Llh: -2.26983427238381E+03; Pars: 5.656962E-01, 5.245094E-01, 5.667486E-02, 1.570959E+00, 5.385203E-02"
## [1] "Llh: -2.26983462229666E+03; Pars: 5.690113E-01, 5.228113E-01, 5.664574E-02, 1.570882E+00, 5.384872E-02"
## [1] "Llh: -2.26983469884395E+03; Pars: 5.609205E-01, 5.278677E-01, 5.786617E-02, 1.571350E+00, 5.314602E-02"
## [1] "Llh: -2.26983462130384E+03; Pars: 5.653352E-01, 5.252079E-01, 5.733481E-02, 1.571127E+00, 5.344623E-02"
## [1] "Llh: -2.26983382052316E+03; Pars: 5.669031E-01, 5.236271E-01, 5.799000E-02, 1.570783E+00, 5.244535E-02"
## [1] "Llh: -2.26983332227818E+03; Pars: 5.650937E-01, 5.241637E-01, 5.893114E-02, 1.570647E+00, 5.137301E-02"
## [1] "Llh: -2.26983318150762E+03; Pars: 5.612103E-01, 5.254380E-01, 5.784637E-02, 1.571079E+00, 5.303090E-02"
## [1] "Llh: -2.26983242399540E+03; Pars: 5.552967E-01, 5.273037E-01, 5.826104E-02, 1.571254E+00, 5.292004E-02"
## [1] "Llh: -2.26983265554044E+03; Pars: 5.583944E-01, 5.268944E-01, 5.876072E-02, 1.571843E+00, 5.120932E-02"
## [1] "Llh: -2.26983302242353E+03; Pars: 5.612908E-01, 5.257370E-01, 5.825037E-02, 1.571465E+00, 5.202167E-02"
## [1] "Llh: -2.26983199204663E+03; Pars: 5.602629E-01, 5.246295E-01, 5.871561E-02, 1.571334E+00, 5.132706E-02"
## [1] "Llh: -2.26983116995934E+03; Pars: 5.577267E-01, 5.243406E-01, 5.940602E-02, 1.571437E+00, 5.026747E-02"
## [1] "Llh: -2.26983186629989E+03; Pars: 5.567140E-01, 5.252606E-01, 6.046802E-02, 1.571694E+00, 4.948743E-02"
## [1] "Llh: -2.26983204562117E+03; Pars: 5.589596E-01, 5.250727E-01, 5.951973E-02, 1.571510E+00, 5.057858E-02"
## [1] "Llh: -2.26983090895374E+03; Pars: 5.477763E-01, 5.294690E-01, 6.083249E-02, 1.571299E+00, 4.952411E-02"
## [1] "Llh: -2.26983073695027E+03; Pars: 5.369075E-01, 5.333755E-01, 6.249959E-02, 1.571223E+00, 4.799677E-02"
## [1] "Llh: -2.26982994352646E+03; Pars: 5.409220E-01, 5.307076E-01, 6.082702E-02, 1.572334E+00, 4.937940E-02"
## [1] "Llh: -2.26982934426356E+03; Pars: 5.288362E-01, 5.340102E-01, 6.177495E-02, 1.573179E+00, 4.838259E-02"
## [1] "Llh: -2.26983002457693E+03; Pars: 5.357980E-01, 5.307979E-01, 6.220314E-02, 1.571672E+00, 4.841239E-02"
## [1] "Llh: -2.26983005492584E+03; Pars: 5.414471E-01, 5.298193E-01, 6.134253E-02, 1.571714E+00, 4.911163E-02"
## [1] "Llh: -2.26983127366296E+03; Pars: 5.310962E-01, 5.317888E-01, 6.427964E-02, 1.572427E+00, 4.489862E-02"
## [1] "Llh: -2.26983029913706E+03; Pars: 5.371463E-01, 5.306640E-01, 6.277499E-02, 1.572134E+00, 4.690397E-02"
## [1] "Llh: -2.26982912244208E+03; Pars: 5.218519E-01, 5.360473E-01, 6.299545E-02, 1.572164E+00, 4.729785E-02"
## [1] "Llh: -2.26982950546061E+03; Pars: 5.044208E-01, 5.415234E-01, 6.425916E-02, 1.572399E+00, 4.620307E-02"
## [1] "Llh: -2.26983182773751E+03; Pars: 5.064892E-01, 5.417517E-01, 6.549323E-02, 1.572711E+00, 4.532996E-02"
## [1] "Llh: -2.26982996343508E+03; Pars: 5.449173E-01, 5.286402E-01, 6.092782E-02, 1.571756E+00, 4.903309E-02"
## [1] "Llh: -2.26982893614454E+03; Pars: 5.305124E-01, 5.306786E-01, 6.177095E-02, 1.573139E+00, 4.801519E-02"
## [1] "Llh: -2.26982902411127E+03; Pars: 5.273149E-01, 5.293352E-01, 6.140663E-02, 1.574097E+00, 4.802440E-02"
## [1] "Llh: -2.26982858607834E+03; Pars: 5.276200E-01, 5.333961E-01, 6.109393E-02, 1.572629E+00, 4.955248E-02"
## [1] "Llh: -2.26982840533131E+03; Pars: 5.228568E-01, 5.347674E-01, 6.025340E-02, 1.572877E+00, 5.087673E-02"
## [1] "Llh: -2.26982829414501E+03; Pars: 5.237918E-01, 5.348532E-01, 6.088589E-02, 1.573574E+00, 4.902979E-02"
## [1] "Llh: -2.26982850421686E+03; Pars: 5.177888E-01, 5.368924E-01, 6.022727E-02, 1.574526E+00, 4.933848E-02"
## [1] "Llh: -2.26982853414710E+03; Pars: 5.062223E-01, 5.395520E-01, 6.214444E-02, 1.574218E+00, 4.840777E-02"
## [1] "Llh: -2.26982835639933E+03; Pars: 5.158961E-01, 5.368031E-01, 6.184028E-02, 1.573602E+00, 4.856410E-02"
## [1] "Llh: -2.26982774394464E+03; Pars: 5.171274E-01, 5.352420E-01, 6.132344E-02, 1.572963E+00, 4.913087E-02"
## [1] "Llh: -2.26982755802155E+03; Pars: 5.112730E-01, 5.358589E-01, 6.109768E-02, 1.572856E+00, 4.950502E-02"
## [1] "Llh: -2.26982827764651E+03; Pars: 5.198802E-01, 5.331329E-01, 5.934383E-02, 1.574256E+00, 5.109848E-02"
## [1] "Llh: -2.26982796693660E+03; Pars: 5.203731E-01, 5.338600E-01, 6.025673E-02, 1.573733E+00, 5.014832E-02"
## [1] "Llh: -2.26982787698435E+03; Pars: 5.071639E-01, 5.398155E-01, 5.996264E-02, 1.573518E+00, 5.123439E-02"
## [1] "Llh: -2.26982773963969E+03; Pars: 5.130010E-01, 5.375167E-01, 6.041472E-02, 1.573423E+00, 5.042959E-02"
## [1] "Llh: -2.26982765057149E+03; Pars: 5.108772E-01, 5.367880E-01, 6.154473E-02, 1.573998E+00, 4.819400E-02"
## [1] "Llh: -2.26982763862559E+03; Pars: 5.138721E-01, 5.362822E-01, 6.122189E-02, 1.573718E+00, 4.886468E-02"
## [1] "Llh: -2.26982729818413E+03; Pars: 5.170284E-01, 5.345448E-01, 5.971049E-02, 1.573319E+00, 5.062686E-02"
## [1] "Llh: -2.26982734031917E+03; Pars: 5.175946E-01, 5.334191E-01, 5.864559E-02, 1.573178E+00, 5.165824E-02"
## [1] "Llh: -2.26982779832964E+03; Pars: 5.064273E-01, 5.363698E-01, 6.019472E-02, 1.573245E+00, 5.080000E-02"
## [1] "Llh: -2.26982751191515E+03; Pars: 5.107684E-01, 5.359902E-01, 6.036751E-02, 1.573327E+00, 5.035745E-02"
## [1] "Llh: -2.26982722922715E+03; Pars: 5.060041E-01, 5.382243E-01, 6.086818E-02, 1.572925E+00, 4.976512E-02"
## [1] "Llh: -2.26982739947713E+03; Pars: 4.988196E-01, 5.404197E-01, 6.117391E-02, 1.572521E+00, 4.957351E-02"
## [1] "Llh: -2.26982736041344E+03; Pars: 5.105774E-01, 5.348442E-01, 6.089158E-02, 1.573035E+00, 4.921806E-02"
## [1] "Llh: -2.26982728326523E+03; Pars: 5.111833E-01, 5.355111E-01, 6.077237E-02, 1.573132E+00, 4.952094E-02"
## [1] "Llh: -2.26982769754837E+03; Pars: 5.086308E-01, 5.357669E-01, 5.990460E-02, 1.572506E+00, 5.104547E-02"
## [1] "Llh: -2.26982737839004E+03; Pars: 5.125618E-01, 5.361533E-01, 6.089257E-02, 1.573415E+00, 4.940988E-02"
## [1] "Llh: -2.26982719335413E+03; Pars: 5.117454E-01, 5.363079E-01, 5.994677E-02, 1.573592E+00, 5.036708E-02"
## [1] "Llh: -2.26982733371552E+03; Pars: 5.119815E-01, 5.365325E-01, 5.937131E-02, 1.573960E+00, 5.079811E-02"
## [1] "Llh: -2.26982695861918E+03; Pars: 5.126408E-01, 5.363036E-01, 6.050864E-02, 1.573226E+00, 4.951850E-02"
## [1] "Llh: -2.26982680753287E+03; Pars: 5.135770E-01, 5.364603E-01, 6.057920E-02, 1.573175E+00, 4.909903E-02"
## [1] "Llh: -2.26982684891856E+03; Pars: 5.112535E-01, 5.362633E-01, 5.985823E-02, 1.573042E+00, 5.034173E-02"
## [1] "Llh: -2.26982694122026E+03; Pars: 5.115806E-01, 5.362358E-01, 6.011682E-02, 1.573135E+00, 5.010877E-02"
## [1] "Llh: -2.26982691417804E+03; Pars: 5.044769E-01, 5.385680E-01, 6.109941E-02, 1.573027E+00, 4.901070E-02"
## [1] "Llh: -2.26982691755247E+03; Pars: 5.076148E-01, 5.375594E-01, 6.075218E-02, 1.573100E+00, 4.941474E-02"
## [1] "Llh: -2.26982666087406E+03; Pars: 5.076394E-01, 5.388217E-01, 6.016835E-02, 1.573172E+00, 4.991252E-02"
## [1] "Llh: -2.26982656727206E+03; Pars: 5.058675E-01, 5.404846E-01, 5.986635E-02, 1.573192E+00, 5.010831E-02"
## [1] "Llh: -2.26982670054901E+03; Pars: 5.127640E-01, 5.370049E-01, 5.967180E-02, 1.573486E+00, 4.980563E-02"
## [1] "Llh: -2.26982667465840E+03; Pars: 5.110740E-01, 5.373095E-01, 5.997090E-02, 1.573346E+00, 4.979550E-02"
## [1] "Llh: -2.26982632339341E+03; Pars: 5.067542E-01, 5.393261E-01, 6.060287E-02, 1.572721E+00, 4.897503E-02"
## [1] "Llh: -2.26982618847514E+03; Pars: 5.042585E-01, 5.408416E-01, 6.093092E-02, 1.572286E+00, 4.827901E-02"
## [1] "Llh: -2.26982654008812E+03; Pars: 5.139353E-01, 5.379685E-01, 5.938283E-02, 1.572989E+00, 5.003873E-02"
## [1] "Llh: -2.26982645310695E+03; Pars: 5.115707E-01, 5.381183E-01, 5.981197E-02, 1.572999E+00, 4.978172E-02"
## [1] "Llh: -2.26982641131685E+03; Pars: 5.072856E-01, 5.410274E-01, 6.060551E-02, 1.572956E+00, 4.848370E-02"
## [1] "Llh: -2.26982637083196E+03; Pars: 5.082776E-01, 5.398324E-01, 6.041869E-02, 1.572978E+00, 4.894821E-02"
## [1] "Llh: -2.26982614591272E+03; Pars: 5.028423E-01, 5.421859E-01, 5.982033E-02, 1.572745E+00, 4.966607E-02"
## [1] "Llh: -2.26982618350617E+03; Pars: 4.974750E-01, 5.450716E-01, 5.944089E-02, 1.572530E+00, 4.994959E-02"
## [1] "Llh: -2.26982610686658E+03; Pars: 5.020527E-01, 5.432889E-01, 6.036841E-02, 1.572334E+00, 4.891783E-02"
## [1] "Llh: -2.26982624300983E+03; Pars: 4.975420E-01, 5.463036E-01, 6.056716E-02, 1.571829E+00, 4.847900E-02"
## [1] "Llh: -2.26982609927969E+03; Pars: 5.057333E-01, 5.412165E-01, 6.067378E-02, 1.572145E+00, 4.812882E-02"
## [1] "Llh: -2.26982637163892E+03; Pars: 5.056662E-01, 5.415828E-01, 6.107750E-02, 1.571622E+00, 4.713908E-02"
## [1] "Llh: -2.26982617093426E+03; Pars: 4.976950E-01, 5.448461E-01, 6.107288E-02, 1.571997E+00, 4.779425E-02"
## [1] "Llh: -2.26982608323277E+03; Pars: 5.011640E-01, 5.431563E-01, 6.075765E-02, 1.572247E+00, 4.829112E-02"
## [1] "Llh: -2.26982607112989E+03; Pars: 4.981427E-01, 5.444513E-01, 6.060175E-02, 1.571725E+00, 4.836493E-02"
## [1] "Llh: -2.26982639728811E+03; Pars: 4.930753E-01, 5.467755E-01, 6.069328E-02, 1.571099E+00, 4.807330E-02"
## [1] "Llh: -2.26982609921777E+03; Pars: 4.997154E-01, 5.448834E-01, 5.995784E-02, 1.572193E+00, 4.906850E-02"
## [1] "Llh: -2.26982603623128E+03; Pars: 5.008512E-01, 5.438701E-01, 6.020111E-02, 1.572216E+00, 4.887113E-02"
## [1] "Llh: -2.26982632755256E+03; Pars: 5.003352E-01, 5.442070E-01, 6.122075E-02, 1.571522E+00, 4.736346E-02"
## [1] "Llh: -2.26982604145640E+03; Pars: 5.022155E-01, 5.426905E-01, 6.017043E-02, 1.572439E+00, 4.909042E-02"
## [1] "Llh: -2.26982601768237E+03; Pars: 5.011900E-01, 5.428628E-01, 6.059349E-02, 1.571975E+00, 4.818074E-02"
## [1] "Llh: -2.26982609347626E+03; Pars: 5.007587E-01, 5.426498E-01, 6.070602E-02, 1.571795E+00, 4.781219E-02"
## [1] "Llh: -2.26982606662706E+03; Pars: 4.956921E-01, 5.456040E-01, 6.025599E-02, 1.572096E+00, 4.899051E-02"
## [1] "Llh: -2.26982602083297E+03; Pars: 4.982024E-01, 5.445038E-01, 6.036044E-02, 1.572108E+00, 4.877509E-02"
## [1] "Llh: -2.26982597810341E+03; Pars: 4.990768E-01, 5.441945E-01, 6.001324E-02, 1.571938E+00, 4.902180E-02"
## [1] "Llh: -2.26982601562230E+03; Pars: 4.980332E-01, 5.447143E-01, 5.964103E-02, 1.571784E+00, 4.938715E-02"
## [1] "Llh: -2.26982604974660E+03; Pars: 5.024717E-01, 5.427976E-01, 5.993373E-02, 1.572546E+00, 4.921074E-02"
## [1] "Llh: -2.26982600142572E+03; Pars: 5.013894E-01, 5.432106E-01, 6.010074E-02, 1.572340E+00, 4.899929E-02"
## [1] "Llh: -2.26982602679665E+03; Pars: 4.980684E-01, 5.447675E-01, 6.033717E-02, 1.571792E+00, 4.844880E-02"
## [1] "Llh: -2.26982599435182E+03; Pars: 4.991052E-01, 5.442475E-01, 6.029549E-02, 1.571954E+00, 4.860921E-02"
## [1] "Llh: -2.26982597811975E+03; Pars: 4.987343E-01, 5.437368E-01, 6.034424E-02, 1.571910E+00, 4.856332E-02"
## [1] "Llh: -2.26982597233896E+03; Pars: 4.992636E-01, 5.437701E-01, 6.030846E-02, 1.571987E+00, 4.864027E-02"
## [1] "Llh: -2.26982595491516E+03; Pars: 5.018076E-01, 5.428112E-01, 6.016412E-02, 1.571969E+00, 4.860543E-02"
## [1] "Llh: -2.26982596961222E+03; Pars: 5.036102E-01, 5.419668E-01, 6.006597E-02, 1.571900E+00, 4.852061E-02"
## [1] "Llh: -2.26982600948405E+03; Pars: 4.990670E-01, 5.444313E-01, 5.975933E-02, 1.572100E+00, 4.936966E-02"
## [1] "Llh: -2.26982597497684E+03; Pars: 4.995978E-01, 5.440387E-01, 5.996787E-02, 1.572069E+00, 4.907243E-02"
## [1] "Llh: -2.26982602674997E+03; Pars: 4.981509E-01, 5.444144E-01, 6.019893E-02, 1.571626E+00, 4.858037E-02"
## [1] "Llh: -2.26982597138746E+03; Pars: 5.005798E-01, 5.435113E-01, 6.012529E-02, 1.572162E+00, 4.889456E-02"
## [1] "Llh: -2.26982594142713E+03; Pars: 5.010250E-01, 5.430830E-01, 5.993610E-02, 1.572096E+00, 4.908460E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982594471585E+03; Pars: 5.022167E-01, 5.424595E-01, 6.011562E-02, 1.572204E+00, 4.879219E-02"
## [1] "Llh: -2.26982594389106E+03; Pars: 5.014317E-01, 5.428927E-01, 6.009003E-02, 1.572137E+00, 4.884959E-02"
## [1] "Llh: -2.26982594739004E+03; Pars: 5.024293E-01, 5.421573E-01, 6.020985E-02, 1.572100E+00, 4.865243E-02"
## [1] "Llh: -2.26982594226217E+03; Pars: 5.017214E-01, 5.426270E-01, 6.014936E-02, 1.572092E+00, 4.875743E-02"
## [1] "Llh: -2.26982600596593E+03; Pars: 5.039849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982601703091E+03; Pars: 5.029849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593692232E+03; Pars: 5.029849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982648904168E+03; Pars: 5.029849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982622427318E+03; Pars: 5.029849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982673380035E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982650276366E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982634975799E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982608743866E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982592837482E+03; Pars: 4.999849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982594427892E+03; Pars: 5.009849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593219231E+03; Pars: 5.009849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982644559247E+03; Pars: 5.009849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982619037159E+03; Pars: 5.009849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982655664369E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982660203476E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982631201624E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982604763454E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982601703091E+03; Pars: 5.029849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982594427892E+03; Pars: 5.009849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982603378376E+03; Pars: 5.019849E-01, 5.435878E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982648802453E+03; Pars: 5.019849E-01, 5.430445E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982624282128E+03; Pars: 5.019849E-01, 5.430445E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982675530377E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982649888600E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982635851745E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982609626304E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982593692232E+03; Pars: 5.029849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593219231E+03; Pars: 5.009849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982594153788E+03; Pars: 5.019849E-01, 5.414178E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982645687492E+03; Pars: 5.019849E-01, 5.419594E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982618204865E+03; Pars: 5.019849E-01, 5.419594E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982654539461E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982661614842E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982631350231E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982604905506E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982648904168E+03; Pars: 5.029849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982644559247E+03; Pars: 5.009849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982648802453E+03; Pars: 5.019849E-01, 5.430445E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982645687492E+03; Pars: 5.019849E-01, 5.419594E-01, 6.075641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982776645592E+03; Pars: 5.019849E-01, 5.425017E-01, 6.175641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982712554715E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982710185477E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982735401364E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982607407936E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982622427318E+03; Pars: 5.029849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982619037159E+03; Pars: 5.009849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982624282128E+03; Pars: 5.019849E-01, 5.430445E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982618204865E+03; Pars: 5.019849E-01, 5.419594E-01, 5.875641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982724657402E+03; Pars: 5.019849E-01, 5.425017E-01, 5.775641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982693366787E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982677351837E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982607744046E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982683065911E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982673380035E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982655664369E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982675530377E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982654539461E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982712554715E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982693366787E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.573741E+00, 4.932229E-02"
## [1] "Llh: -2.26982864415765E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.575315E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982700991037E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 5.032229E-02"
## [1] "Llh: -2.26982677272208E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 4.832229E-02"
## [1] "Llh: -2.26982650276366E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982660203476E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982649888600E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982661614842E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982710185477E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982677351837E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.570596E+00, 4.932229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982846080091E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.569027E+00, 4.932229E-02"
## [1] "Llh: -2.26982694458615E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 5.032229E-02"
## [1] "Llh: -2.26982665494056E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 4.832229E-02"
## [1] "Llh: -2.26982634975799E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982631201624E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982635851745E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982631350231E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982735401364E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982607744046E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 5.032229E-02"
## [1] "Llh: -2.26982700991037E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 5.032229E-02"
## [1] "Llh: -2.26982694458615E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 5.032229E-02"
## [1] "Llh: -2.26982721740931E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 5.132229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982608743866E+03; Pars: 5.029849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982604763454E+03; Pars: 5.009849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982609626304E+03; Pars: 5.019849E-01, 5.430445E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982604905506E+03; Pars: 5.019849E-01, 5.419594E-01, 5.975641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982607407936E+03; Pars: 5.019849E-01, 5.425017E-01, 6.075641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982683065911E+03; Pars: 5.019849E-01, 5.425017E-01, 5.875641E-02, 1.572168E+00, 4.832229E-02"
## [1] "Llh: -2.26982677272208E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.573741E+00, 4.832229E-02"
## [1] "Llh: -2.26982665494056E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.570596E+00, 4.832229E-02"
## [1] "Llh: -2.26982593779657E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.932229E-02"
## [1] "Llh: -2.26982669075794E+03; Pars: 5.019849E-01, 5.425017E-01, 5.975641E-02, 1.572168E+00, 4.732229E-02"hypfit <- fit.hypuv(HDTimeS, symmetric = FALSE, control = list(maxit = 1000), na.rm = TRUE)## Warning in .check.data(data = data, case = "uv", na.rm = na.rm): 1 NA
## observations removed## [1] "Llh: -2.33741670315059E+03; Pars: 5.000000E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.33731546514863E+03; Pars: 5.358867E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.33918920416620E+03; Pars: 5.000000E-01, 1.476733E-01, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.30877625938294E+03; Pars: 5.000000E-01, 7.835855E-02, 1.234866E+00, 0.000000E+00"
## [1] "Llh: -2.33803609975631E+03; Pars: 5.000000E-01, 7.835855E-02, 1.152171E+00, 6.931472E-02"
## [1] "Llh: -2.32272931777261E+03; Pars: 5.176325E-01, 9.043833E-03, 1.192802E+00, 3.465736E-02"
## [1] "Llh: -2.32543578674348E+03; Pars: 5.131669E-01, 4.370119E-02, 1.182512E+00, 2.599302E-02"
## [1] "Llh: -2.32021465288350E+03; Pars: 5.266805E-01, 4.370119E-02, 1.213652E+00, -5.198604E-02"
## [1] "Llh: -2.32210557252292E+03; Pars: 5.198798E-01, 5.236553E-02, 1.197980E+00, -2.166085E-02"
## [1] "Llh: -2.30822629100700E+03; Pars: 5.405500E-01, 2.637251E-02, 1.245612E+00, -8.664340E-03"
## [1] "Llh: -2.29880356234911E+03; Pars: 5.620421E-01, 3.794936E-04, 1.295137E+00, -1.299651E-02"
## [1] "Llh: -2.29527238745195E+03; Pars: 5.165124E-01, -1.261702E-02, 1.320633E+00, -1.516259E-02"
## [1] "Llh: -2.28859681623889E+03; Pars: 5.070896E-01, -5.810480E-02, 1.413887E+00, -2.274389E-02"
## [1] "Llh: -2.28771031111283E+03; Pars: 5.292537E-01, 2.312339E-02, 1.388845E+00, -7.852058E-02"
## [1] "Llh: -2.28464195736404E+03; Pars: 5.351617E-01, 3.016316E-02, 1.498637E+00, -1.351095E-01"
## [1] "Llh: -2.27893088637109E+03; Pars: 5.243328E-01, -1.830299E-02, 1.516803E+00, -3.343894E-02"
## [1] "Llh: -2.28652094801119E+03; Pars: 5.231629E-01, -4.930508E-02, 1.695692E+00, -2.416539E-02"
## [1] "Llh: -2.30075633474120E+03; Pars: 5.655914E-01, -1.012911E-01, 1.652183E+00, -1.021444E-01"
## [1] "Llh: -2.28831457950941E+03; Pars: 5.484280E-01, -5.637870E-02, 1.536203E+00, -7.660833E-02"
## [1] "Llh: -2.29735031312417E+03; Pars: 4.970263E-01, -5.169116E-02, 1.715651E+00, -1.209538E-01"
## [1] "Llh: -2.28753931814660E+03; Pars: 5.125388E-01, -3.867349E-02, 1.599192E+00, -9.396451E-02"
## [1] "Llh: -2.28939131083366E+03; Pars: 5.538402E-01, 1.650879E-02, 1.671371E+00, -1.468168E-01"
## [1] "Llh: -2.28508871668243E+03; Pars: 5.183936E-01, -3.945140E-02, 1.474277E+00, -5.376211E-02"
## [1] "Llh: -2.27918085111231E+03; Pars: 4.978742E-01, 2.324634E-02, 1.506969E+00, -8.152922E-02"
## [1] "Llh: -2.28102180341072E+03; Pars: 5.100581E-01, 3.340078E-03, 1.514225E+00, -8.029900E-02"
## [1] "Llh: -2.28258129465431E+03; Pars: 5.250600E-01, 3.650105E-02, 1.405251E+00, -5.795540E-02"
## [1] "Llh: -2.28118000062715E+03; Pars: 5.219013E-01, 1.770741E-02, 1.451412E+00, -6.695768E-02"
## [1] "Llh: -2.27742995916953E+03; Pars: 5.208833E-01, 6.585836E-02, 1.512450E+00, -1.047556E-01"
## [1] "Llh: -2.27557444332925E+03; Pars: 5.221326E-01, 1.185132E-01, 1.531906E+00, -1.302523E-01"
## [1] "Llh: -2.27364332559690E+03; Pars: 4.983826E-01, 4.041884E-02, 1.504290E+00, -2.097953E-02"
## [1] "Llh: -2.27147610963971E+03; Pars: 4.809520E-01, 4.554668E-02, 1.507124E+00, 3.608549E-02"
## [1] "Llh: -2.27308443575350E+03; Pars: 4.905872E-01, 6.679423E-02, 1.582766E+00, -3.760982E-02"
## [1] "Llh: -2.27362935965881E+03; Pars: 4.982350E-01, 5.452252E-02, 1.548853E+00, -4.494678E-02"
## [1] "Llh: -2.27154143938611E+03; Pars: 5.104853E-01, 8.302925E-02, 1.562281E+00, -1.078568E-03"
## [1] "Llh: -2.27214011664031E+03; Pars: 5.073029E-01, 6.808352E-02, 1.548266E+00, -2.119123E-02"
## [1] "Llh: -2.27314589256076E+03; Pars: 4.782820E-01, 1.752447E-01, 1.575249E+00, -3.298867E-02"
## [1] "Llh: -2.27173235809256E+03; Pars: 4.894009E-01, 1.268578E-01, 1.560430E+00, -3.310124E-02"
## [1] "Llh: -2.27162807766339E+03; Pars: 4.649992E-01, 4.260072E-02, 1.574173E+00, 1.124002E-01"
## [1] "Llh: -2.27104508826473E+03; Pars: 4.786680E-01, 6.157885E-02, 1.563498E+00, 5.173711E-02"
## [1] "Llh: -2.27199073690629E+03; Pars: 4.888503E-01, 9.171205E-02, 1.514287E+00, 6.443121E-02"
## [1] "Llh: -2.27118728135853E+03; Pars: 4.892839E-01, 8.548259E-02, 1.531124E+00, 3.892095E-02"
## [1] "Llh: -2.27127778898239E+03; Pars: 4.899767E-01, 1.096092E-02, 1.521472E+00, 9.593373E-02"
## [1] "Llh: -2.27102346013799E+03; Pars: 4.898327E-01, 3.993513E-02, 1.531119E+00, 6.367498E-02"
## [1] "Llh: -2.27140997205190E+03; Pars: 4.601391E-01, 3.324238E-02, 1.504436E+00, 9.628783E-02"
## [1] "Llh: -2.27105418043301E+03; Pars: 4.722399E-01, 4.568910E-02, 1.518693E+00, 7.194623E-02"
## [1] "Llh: -2.27151513233552E+03; Pars: 4.839507E-01, 7.079615E-02, 1.565469E+00, 7.705415E-02"
## [1] "Llh: -2.27106803899880E+03; Pars: 4.817000E-01, 5.185905E-02, 1.521503E+00, 4.632765E-02"
## [1] "Llh: -2.27110715753799E+03; Pars: 4.720084E-01, 1.404847E-02, 1.536081E+00, 7.792203E-02"
## [1] "Llh: -2.27096401951549E+03; Pars: 4.762693E-01, 3.190700E-02, 1.534840E+00, 6.817176E-02"
## [1] "Llh: -2.27097136252657E+03; Pars: 4.767295E-01, 3.769599E-02, 1.552555E+00, 8.143739E-02"
## [1] "Llh: -2.27091723106353E+03; Pars: 4.779673E-01, 4.123676E-02, 1.544733E+00, 7.265996E-02"
## [1] "Llh: -2.27123996145114E+03; Pars: 4.892195E-01, 4.163977E-02, 1.568706E+00, 5.617567E-02"
## [1] "Llh: -2.27093035906706E+03; Pars: 4.764287E-01, 4.467677E-02, 1.531045E+00, 6.800359E-02"
## [1] "Llh: -2.27140332390003E+03; Pars: 4.815195E-01, 1.729898E-02, 1.507854E+00, 8.451804E-02"
## [1] "Llh: -2.27092042681997E+03; Pars: 4.793793E-01, 5.050888E-02, 1.549398E+00, 5.993234E-02"
## [1] "Llh: -2.27084350986892E+03; Pars: 4.654963E-01, 4.422957E-02, 1.548905E+00, 7.070884E-02"
## [1] "Llh: -2.27084170674628E+03; Pars: 4.537853E-01, 4.637679E-02, 1.557875E+00, 7.422577E-02"
## [1] "Llh: -2.27098926265992E+03; Pars: 4.673157E-01, 5.949260E-02, 1.556701E+00, 6.923907E-02"
## [1] "Llh: -2.27089084139161E+03; Pars: 4.740149E-01, 3.880340E-02, 1.540277E+00, 6.843859E-02"
## [1] "Llh: -2.27097722627148E+03; Pars: 4.659743E-01, 4.378615E-02, 1.565258E+00, 6.962474E-02"
## [1] "Llh: -2.27088075615561E+03; Pars: 4.737934E-01, 4.445411E-02, 1.539527E+00, 6.840888E-02"
## [1] "Llh: -2.27083677608152E+03; Pars: 4.604003E-01, 3.492665E-02, 1.541783E+00, 8.193426E-02"
## [1] "Llh: -2.27084631376701E+03; Pars: 4.511944E-01, 2.713553E-02, 1.537990E+00, 9.293522E-02"
## [1] "Llh: -2.27076696249692E+03; Pars: 4.531952E-01, 4.104372E-02, 1.544961E+00, 7.384379E-02"
## [1] "Llh: -2.27070165627614E+03; Pars: 4.412948E-01, 4.094721E-02, 1.545076E+00, 7.443571E-02"
## [1] "Llh: -2.27076748212732E+03; Pars: 4.409201E-01, 4.454898E-02, 1.551843E+00, 8.106372E-02"
## [1] "Llh: -2.27077324740188E+03; Pars: 4.489706E-01, 4.311259E-02, 1.548943E+00, 7.790744E-02"
## [1] "Llh: -2.27073531634487E+03; Pars: 4.255479E-01, 3.894570E-02, 1.558796E+00, 8.742085E-02"
## [1] "Llh: -2.27074228721773E+03; Pars: 4.371279E-01, 4.032281E-02, 1.553957E+00, 8.266786E-02"
## [1] "Llh: -2.27070659100857E+03; Pars: 4.302649E-01, 3.330748E-02, 1.540893E+00, 8.820150E-02"
## [1] "Llh: -2.27071545550390E+03; Pars: 4.360282E-01, 3.657481E-02, 1.545121E+00, 8.470757E-02"
## [1] "Llh: -2.27064027363006E+03; Pars: 4.099689E-01, 4.394803E-02, 1.556526E+00, 8.362663E-02"
## [1] "Llh: -2.27063220132912E+03; Pars: 3.868643E-01, 4.845873E-02, 1.563951E+00, 8.447282E-02"
## [1] "Llh: -2.27055947137921E+03; Pars: 4.009830E-01, 3.628058E-02, 1.552457E+00, 8.620172E-02"
## [1] "Llh: -2.27048764005294E+03; Pars: 3.823922E-01, 3.214637E-02, 1.552763E+00, 8.877072E-02"
## [1] "Llh: -2.27051216067091E+03; Pars: 3.938382E-01, 3.848419E-02, 1.542538E+00, 8.051952E-02"
## [1] "Llh: -2.27054193524812E+03; Pars: 4.015369E-01, 3.859957E-02, 1.546587E+00, 8.224485E-02"
## [1] "Llh: -2.27045904798037E+03; Pars: 3.726673E-01, 4.671077E-02, 1.561293E+00, 7.589789E-02"
## [1] "Llh: -2.27045629886348E+03; Pars: 3.468278E-01, 5.341242E-02, 1.571595E+00, 6.974608E-02"
## [1] "Llh: -2.27046866063755E+03; Pars: 3.221223E-01, 4.530365E-02, 1.570372E+00, 8.731886E-02"
## [1] "Llh: -2.27044432786898E+03; Pars: 3.484960E-01, 4.421454E-02, 1.564009E+00, 8.409807E-02"
## [1] "Llh: -2.27033945828977E+03; Pars: 3.487471E-01, 3.567003E-02, 1.551449E+00, 7.709438E-02"
## [1] "Llh: -2.27035707899710E+03; Pars: 3.311210E-01, 2.927569E-02, 1.545236E+00, 7.340516E-02"
## [1] "Llh: -2.27044925935580E+03; Pars: 3.223650E-01, 4.423749E-02, 1.577522E+00, 7.933510E-02"
## [1] "Llh: -2.27040065606033E+03; Pars: 3.389147E-01, 4.279917E-02, 1.568703E+00, 7.963121E-02"
## [1] "Llh: -2.27041121371803E+03; Pars: 3.125701E-01, 5.590171E-02, 1.575157E+00, 6.651415E-02"
## [1] "Llh: -2.27038000769675E+03; Pars: 3.287288E-01, 4.996287E-02, 1.569528E+00, 7.207829E-02"
## [1] "Llh: -2.27034483773064E+03; Pars: 3.355092E-01, 3.291088E-02, 1.555259E+00, 8.670490E-02"
## [1] "Llh: -2.27034701568932E+03; Pars: 3.383037E-01, 3.803627E-02, 1.559327E+00, 8.246519E-02"
## [1] "Llh: -2.27029391921204E+03; Pars: 3.276225E-01, 3.645694E-02, 1.558424E+00, 7.365632E-02"
## [1] "Llh: -2.27030507402156E+03; Pars: 3.176594E-01, 3.257814E-02, 1.555639E+00, 6.843544E-02"
## [1] "Llh: -2.27031399534474E+03; Pars: 3.312249E-01, 3.470120E-02, 1.548663E+00, 7.513573E-02"
## [1] "Llh: -2.27030471353397E+03; Pars: 3.331308E-01, 3.672569E-02, 1.553649E+00, 7.625960E-02"
## [1] "Llh: -2.27044592293375E+03; Pars: 3.437669E-01, 2.091890E-02, 1.539998E+00, 8.477930E-02"
## [1] "Llh: -2.27032298544888E+03; Pars: 3.324255E-01, 4.270188E-02, 1.562093E+00, 7.525355E-02"
## [1] "Llh: -2.27029626082563E+03; Pars: 3.352682E-01, 4.286639E-02, 1.557538E+00, 6.442703E-02"
## [1] "Llh: -2.27029834785658E+03; Pars: 3.353284E-01, 4.037751E-02, 1.556968E+00, 6.999649E-02"
## [1] "Llh: -2.27028545115486E+03; Pars: 3.162474E-01, 4.370542E-02, 1.564424E+00, 6.770387E-02"
## [1] "Llh: -2.27030470363036E+03; Pars: 3.011516E-01, 4.772311E-02, 1.570952E+00, 6.300861E-02"
## [1] "Llh: -2.27029338493317E+03; Pars: 3.236008E-01, 3.717534E-02, 1.554923E+00, 6.576986E-02"
## [1] "Llh: -2.27028583764032E+03; Pars: 3.257848E-01, 3.855697E-02, 1.556713E+00, 6.814078E-02"
## [1] "Llh: -2.27029476263207E+03; Pars: 3.193355E-01, 4.406717E-02, 1.564915E+00, 6.070439E-02"
## [1] "Llh: -2.27028655367435E+03; Pars: 3.227298E-01, 4.223180E-02, 1.562091E+00, 6.459319E-02"
## [1] "Llh: -2.27028099864993E+03; Pars: 3.113099E-01, 3.760918E-02, 1.563287E+00, 7.262005E-02"
## [1] "Llh: -2.27028878872665E+03; Pars: 2.999807E-01, 3.498057E-02, 1.566169E+00, 7.671657E-02"
## [1] "Llh: -2.27027926290818E+03; Pars: 3.105426E-01, 4.459474E-02, 1.564834E+00, 6.287263E-02"
## [1] "Llh: -2.27029056055493E+03; Pars: 3.023394E-01, 4.866364E-02, 1.568049E+00, 5.748079E-02"
## [1] "Llh: -2.27027374673761E+03; Pars: 3.092409E-01, 4.000136E-02, 1.562531E+00, 7.107547E-02"
## [1] "Llh: -2.27027439713648E+03; Pars: 3.027093E-01, 3.888613E-02, 1.562751E+00, 7.431661E-02"
## [1] "Llh: -2.27030453163345E+03; Pars: 2.984614E-01, 4.439837E-02, 1.570857E+00, 6.899523E-02"
## [1] "Llh: -2.27027741731368E+03; Pars: 3.187279E-01, 4.001732E-02, 1.560237E+00, 6.835439E-02"
## [1] "Llh: -2.27027329261481E+03; Pars: 3.086658E-01, 3.740588E-02, 1.561021E+00, 6.975741E-02"
## [1] "Llh: -2.27027840159154E+03; Pars: 3.049435E-01, 3.425612E-02, 1.559322E+00, 7.078418E-02"
## [1] "Llh: -2.27027018514271E+03; Pars: 3.122270E-01, 4.340047E-02, 1.561023E+00, 6.340990E-02"
## [1] "Llh: -2.27027134372113E+03; Pars: 3.126865E-01, 4.629612E-02, 1.559893E+00, 5.880482E-02"
## [1] "Llh: -2.27027659701036E+03; Pars: 3.138462E-01, 3.581778E-02, 1.557580E+00, 7.342595E-02"
## [1] "Llh: -2.27027274131002E+03; Pars: 3.130170E-01, 3.801202E-02, 1.559390E+00, 7.078762E-02"
## [1] "Llh: -2.27026697996216E+03; Pars: 3.030343E-01, 3.939254E-02, 1.561745E+00, 6.916081E-02"
## [1] "Llh: -2.27026601812264E+03; Pars: 2.954796E-01, 3.908015E-02, 1.562500E+00, 6.956401E-02"
## [1] "Llh: -2.27027016400170E+03; Pars: 3.053025E-01, 3.894791E-02, 1.559437E+00, 6.568400E-02"
## [1] "Llh: -2.27026838027787E+03; Pars: 3.062824E-01, 3.921127E-02, 1.560210E+00, 6.703187E-02"
## [1] "Llh: -2.27026338559031E+03; Pars: 3.046880E-01, 4.244607E-02, 1.560540E+00, 6.563929E-02"
## [1] "Llh: -2.27026183635168E+03; Pars: 3.027183E-01, 4.496617E-02, 1.560300E+00, 6.358024E-02"
## [1] "Llh: -2.27026542581721E+03; Pars: 2.954686E-01, 4.531702E-02, 1.562628E+00, 6.100539E-02"
## [1] "Llh: -2.27026382115610E+03; Pars: 2.997613E-01, 4.349077E-02, 1.561818E+00, 6.345095E-02"
## [1] "Llh: -2.27026341178140E+03; Pars: 2.902428E-01, 3.997370E-02, 1.561390E+00, 6.840363E-02"
## [1] "Llh: -2.27026256747741E+03; Pars: 2.955893E-01, 4.083040E-02, 1.561298E+00, 6.715520E-02"
## [1] "Llh: -2.27026416295405E+03; Pars: 2.906654E-01, 4.497247E-02, 1.562749E+00, 6.484333E-02"
## [1] "Llh: -2.27026240456560E+03; Pars: 2.944933E-01, 4.353217E-02, 1.562114E+00, 6.539047E-02"
## [1] "Llh: -2.27026301705289E+03; Pars: 3.007888E-01, 4.732960E-02, 1.560265E+00, 6.022441E-02"
## [1] "Llh: -2.27026183422679E+03; Pars: 2.994526E-01, 4.526724E-02, 1.560824E+00, 6.255931E-02"
## [1] "Llh: -2.27026068813951E+03; Pars: 2.963398E-01, 4.380722E-02, 1.560451E+00, 6.589166E-02"
## [1] "Llh: -2.27026147755843E+03; Pars: 2.946437E-01, 4.396545E-02, 1.559767E+00, 6.711202E-02"
## [1] "Llh: -2.27026361030590E+03; Pars: 3.009036E-01, 4.795600E-02, 1.560546E+00, 6.155564E-02"
## [1] "Llh: -2.27026142341861E+03; Pars: 2.969090E-01, 4.261180E-02, 1.561110E+00, 6.575531E-02"
## [1] "Llh: -2.27026119963202E+03; Pars: 3.032597E-01, 4.479404E-02, 1.559230E+00, 6.350280E-02"
## [1] "Llh: -2.27026092192229E+03; Pars: 3.010439E-01, 4.447857E-02, 1.559950E+00, 6.397471E-02"
## [1] "Llh: -2.27026088165840E+03; Pars: 2.942029E-01, 4.311625E-02, 1.560867E+00, 6.551026E-02"
## [1] "Llh: -2.27026075093232E+03; Pars: 2.963090E-01, 4.357873E-02, 1.560725E+00, 6.502775E-02"
## [1] "Llh: -2.27026093234765E+03; Pars: 2.958462E-01, 4.197092E-02, 1.560294E+00, 6.776541E-02"
## [1] "Llh: -2.27026064874040E+03; Pars: 2.967437E-01, 4.279500E-02, 1.560427E+00, 6.646388E-02"
## [1] "Llh: -2.27026045019252E+03; Pars: 2.982976E-01, 4.471796E-02, 1.559666E+00, 6.492370E-02"
## [1] "Llh: -2.27026111290850E+03; Pars: 2.989943E-01, 4.577105E-02, 1.558945E+00, 6.450789E-02"
## [1] "Llh: -2.27026127609550E+03; Pars: 2.928554E-01, 4.297088E-02, 1.560684E+00, 6.717879E-02"
## [1] "Llh: -2.27026054722621E+03; Pars: 2.989756E-01, 4.410165E-02, 1.560134E+00, 6.477573E-02"
## [1] "Llh: -2.27026064841907E+03; Pars: 2.988708E-01, 4.413219E-02, 1.559613E+00, 6.599973E-02"
## [1] "Llh: -2.27026048750723E+03; Pars: 2.982283E-01, 4.399383E-02, 1.559891E+00, 6.575674E-02"
## [1] "Llh: -2.27026035595086E+03; Pars: 2.997905E-01, 4.399700E-02, 1.559609E+00, 6.506836E-02"
## [1] "Llh: -2.27026052547876E+03; Pars: 3.015309E-01, 4.409189E-02, 1.559188E+00, 6.465671E-02"
## [1] "Llh: -2.27026087616147E+03; Pars: 3.009155E-01, 4.561022E-02, 1.559224E+00, 6.379838E-02"
## [1] "Llh: -2.27026042600254E+03; Pars: 2.977812E-01, 4.349881E-02, 1.560126E+00, 6.579751E-02"
## [1] "Llh: -2.27026043004688E+03; Pars: 2.980720E-01, 4.400215E-02, 1.559512E+00, 6.599742E-02"
## [1] "Llh: -2.27026036690370E+03; Pars: 2.982976E-01, 4.402702E-02, 1.559668E+00, 6.569200E-02"
## [1] "Llh: -2.27026024527217E+03; Pars: 2.988536E-01, 4.412657E-02, 1.559643E+00, 6.498405E-02"
## [1] "Llh: -2.27026021745346E+03; Pars: 2.991667E-01, 4.419294E-02, 1.559519E+00, 6.459770E-02"
## [1] "Llh: -2.27026040220075E+03; Pars: 2.992192E-01, 4.313992E-02, 1.559794E+00, 6.565409E-02"
## [1] "Llh: -2.27026030548139E+03; Pars: 2.989885E-01, 4.353443E-02, 1.559762E+00, 6.547149E-02"
## [1] "Llh: -2.27026029755477E+03; Pars: 3.003451E-01, 4.437689E-02, 1.559153E+00, 6.461727E-02"
## [1] "Llh: -2.27026025593069E+03; Pars: 2.997020E-01, 4.415737E-02, 1.559396E+00, 6.491233E-02"
## [1] "Llh: -2.27026050555051E+03; Pars: 3.005300E-01, 4.391385E-02, 1.559475E+00, 6.433294E-02"
## [1] "Llh: -2.27026027321802E+03; Pars: 2.988542E-01, 4.399873E-02, 1.559620E+00, 6.535223E-02"
## [1] "Llh: -2.27026014322479E+03; Pars: 2.985661E-01, 4.394474E-02, 1.559540E+00, 6.509851E-02"
## [1] "Llh: -2.27026007120420E+03; Pars: 2.979558E-01, 4.391861E-02, 1.559505E+00, 6.511359E-02"
## [1] "Llh: -2.27026017888294E+03; Pars: 2.988495E-01, 4.459939E-02, 1.559258E+00, 6.451643E-02"
## [1] "Llh: -2.27026016565997E+03; Pars: 2.988843E-01, 4.433315E-02, 1.559384E+00, 6.475520E-02"
## [1] "Llh: -2.27026013394482E+03; Pars: 2.989989E-01, 4.430231E-02, 1.559282E+00, 6.433717E-02"
## [1] "Llh: -2.27026013681400E+03; Pars: 2.989627E-01, 4.422641E-02, 1.559367E+00, 6.459094E-02"
## [1] "Llh: -2.27026006495956E+03; Pars: 2.978030E-01, 4.421613E-02, 1.559449E+00, 6.448950E-02"
## [1] "Llh: -2.27026006532191E+03; Pars: 2.968581E-01, 4.424552E-02, 1.559475E+00, 6.427809E-02"
## [1] "Llh: -2.27026001154158E+03; Pars: 2.976552E-01, 4.419216E-02, 1.559291E+00, 6.475003E-02"
## [1] "Llh: -2.27025998493285E+03; Pars: 2.969023E-01, 4.419177E-02, 1.559178E+00, 6.482620E-02"
## [1] "Llh: -2.27025998429111E+03; Pars: 2.969470E-01, 4.398125E-02, 1.559323E+00, 6.462803E-02"
## [1] "Llh: -2.27026002204697E+03; Pars: 2.959831E-01, 4.380530E-02, 1.559293E+00, 6.456445E-02"
## [1] "Llh: -2.27026001651251E+03; Pars: 2.958129E-01, 4.385157E-02, 1.559445E+00, 6.519149E-02"
## [1] "Llh: -2.27025999026360E+03; Pars: 2.966062E-01, 4.396426E-02, 1.559404E+00, 6.497791E-02"
## [1] "Llh: -2.27026001616331E+03; Pars: 2.961755E-01, 4.425810E-02, 1.559172E+00, 6.434723E-02"
## [1] "Llh: -2.27025998829172E+03; Pars: 2.966196E-01, 4.417323E-02, 1.559255E+00, 6.453882E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026008906597E+03; Pars: 2.947109E-01, 4.380061E-02, 1.558973E+00, 6.524921E-02"
## [1] "Llh: -2.27025995441786E+03; Pars: 2.963301E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026051778638E+03; Pars: 2.960339E-01, 4.493912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026019938758E+03; Pars: 2.960339E-01, 4.293912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026070162106E+03; Pars: 2.960339E-01, 4.393912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026057684623E+03; Pars: 2.960339E-01, 4.393912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026028409744E+03; Pars: 2.960339E-01, 4.393912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026016621765E+03; Pars: 2.960339E-01, 4.393912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026003268662E+03; Pars: 2.951471E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026056749024E+03; Pars: 2.954424E-01, 4.493912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026022808802E+03; Pars: 2.954424E-01, 4.293912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026068760175E+03; Pars: 2.954424E-01, 4.393912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026066899547E+03; Pars: 2.954424E-01, 4.393912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026032376221E+03; Pars: 2.954424E-01, 4.393912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026020484932E+03; Pars: 2.954424E-01, 4.393912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026051778638E+03; Pars: 2.960339E-01, 4.493912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026056749024E+03; Pars: 2.954424E-01, 4.493912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026187970732E+03; Pars: 2.957380E-01, 4.593912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026121034564E+03; Pars: 2.957380E-01, 4.493912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026121487364E+03; Pars: 2.957380E-01, 4.493912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026137439504E+03; Pars: 2.957380E-01, 4.493912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026022263808E+03; Pars: 2.957380E-01, 4.493912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026019938758E+03; Pars: 2.960339E-01, 4.293912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026022808802E+03; Pars: 2.954424E-01, 4.293912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026122209707E+03; Pars: 2.957380E-01, 4.193912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026095664603E+03; Pars: 2.957380E-01, 4.293912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026081055028E+03; Pars: 2.957380E-01, 4.293912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026001214973E+03; Pars: 2.957380E-01, 4.293912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026092709421E+03; Pars: 2.957380E-01, 4.293912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026070162106E+03; Pars: 2.960339E-01, 4.393912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026068760175E+03; Pars: 2.954424E-01, 4.393912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026121034564E+03; Pars: 2.957380E-01, 4.493912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026095664603E+03; Pars: 2.957380E-01, 4.293912E-02, 1.560691E+00, 6.499598E-02"
## [1] "Llh: -2.27026274355285E+03; Pars: 2.957380E-01, 4.393912E-02, 1.562253E+00, 6.499598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026099322550E+03; Pars: 2.957380E-01, 4.393912E-02, 1.560691E+00, 6.599598E-02"
## [1] "Llh: -2.27026090646277E+03; Pars: 2.957380E-01, 4.393912E-02, 1.560691E+00, 6.399598E-02"
## [1] "Llh: -2.27026057684623E+03; Pars: 2.960339E-01, 4.393912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026066899547E+03; Pars: 2.954424E-01, 4.393912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026121487364E+03; Pars: 2.957380E-01, 4.493912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27026081055028E+03; Pars: 2.957380E-01, 4.293912E-02, 1.557573E+00, 6.499598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026260287782E+03; Pars: 2.957380E-01, 4.393912E-02, 1.556016E+00, 6.499598E-02"
## [1] "Llh: -2.27026095454705E+03; Pars: 2.957380E-01, 4.393912E-02, 1.557573E+00, 6.599598E-02"
## [1] "Llh: -2.27026080435701E+03; Pars: 2.957380E-01, 4.393912E-02, 1.557573E+00, 6.399598E-02"
## [1] "Llh: -2.27026028409744E+03; Pars: 2.960339E-01, 4.393912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026032376221E+03; Pars: 2.954424E-01, 4.393912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026137439504E+03; Pars: 2.957380E-01, 4.493912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026001214973E+03; Pars: 2.957380E-01, 4.293912E-02, 1.559131E+00, 6.599598E-02"
## [1] "Llh: -2.27026099322550E+03; Pars: 2.957380E-01, 4.393912E-02, 1.560691E+00, 6.599598E-02"
## [1] "Llh: -2.27026095454705E+03; Pars: 2.957380E-01, 4.393912E-02, 1.557573E+00, 6.599598E-02"
## [1] "Llh: -2.27026113544454E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.699598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026016621765E+03; Pars: 2.960339E-01, 4.393912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026020484932E+03; Pars: 2.954424E-01, 4.393912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026022263808E+03; Pars: 2.957380E-01, 4.493912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026092709421E+03; Pars: 2.957380E-01, 4.293912E-02, 1.559131E+00, 6.399598E-02"
## [1] "Llh: -2.27026090646277E+03; Pars: 2.957380E-01, 4.393912E-02, 1.560691E+00, 6.399598E-02"
## [1] "Llh: -2.27026080435701E+03; Pars: 2.957380E-01, 4.393912E-02, 1.557573E+00, 6.399598E-02"
## [1] "Llh: -2.27025998414876E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.499598E-02"
## [1] "Llh: -2.27026089870873E+03; Pars: 2.957380E-01, 4.393912E-02, 1.559131E+00, 6.299598E-02"nigfit <- fit.NIGuv(HDTimeS, symmetric = FALSE, control = list(maxit = 1000), na.rm = TRUE)## Warning in .check.data(data = data, case = "uv", na.rm = na.rm): 1 NA
## observations removed## [1] "Llh: -2.35705793847316E+03; Pars: 5.000000E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.35326688078383E+03; Pars: 5.358867E-01, 7.835855E-02, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.35907524964687E+03; Pars: 5.000000E-01, 1.476733E-01, 1.152171E+00, 0.000000E+00"
## [1] "Llh: -2.32717417163025E+03; Pars: 5.000000E-01, 7.835855E-02, 1.234866E+00, 0.000000E+00"
## [1] "Llh: -2.35799132150721E+03; Pars: 5.000000E-01, 7.835855E-02, 1.152171E+00, 6.931472E-02"
## [1] "Llh: -2.34062327359152E+03; Pars: 5.176325E-01, 9.043833E-03, 1.192802E+00, 3.465736E-02"
## [1] "Llh: -2.34360494418594E+03; Pars: 5.131669E-01, 4.370119E-02, 1.182512E+00, 2.599302E-02"
## [1] "Llh: -2.33675726937094E+03; Pars: 5.266805E-01, 4.370119E-02, 1.213652E+00, -5.198604E-02"
## [1] "Llh: -2.33937199424744E+03; Pars: 5.198798E-01, 5.236553E-02, 1.197980E+00, -2.166085E-02"
## [1] "Llh: -2.32296705393143E+03; Pars: 5.405500E-01, 2.637251E-02, 1.245612E+00, -8.664340E-03"
## [1] "Llh: -2.31098436788160E+03; Pars: 5.620421E-01, 3.794936E-04, 1.295137E+00, -1.299651E-02"
## [1] "Llh: -2.31019502626412E+03; Pars: 5.165124E-01, -1.261702E-02, 1.320633E+00, -1.516259E-02"
## [1] "Llh: -2.30099683510631E+03; Pars: 5.070896E-01, -5.810480E-02, 1.413887E+00, -2.274389E-02"
## [1] "Llh: -2.29867603122070E+03; Pars: 5.292537E-01, 2.312339E-02, 1.388845E+00, -7.852058E-02"
## [1] "Llh: -2.29013650741856E+03; Pars: 5.351617E-01, 3.016316E-02, 1.498637E+00, -1.351095E-01"
## [1] "Llh: -2.28441123023391E+03; Pars: 5.243328E-01, -1.830299E-02, 1.516803E+00, -3.343894E-02"
## [1] "Llh: -2.28268368109417E+03; Pars: 5.231629E-01, -4.930508E-02, 1.695692E+00, -2.416539E-02"
## [1] "Llh: -2.30221262328985E+03; Pars: 5.649600E-01, -1.167922E-01, 1.746896E+00, -9.750767E-02"
## [1] "Llh: -2.29063527285927E+03; Pars: 5.479688E-01, -6.800448E-02, 1.601789E+00, -7.313075E-02"
## [1] "Llh: -2.29980649391654E+03; Pars: 4.962636E-01, -7.300510E-02, 1.852320E+00, -1.145783E-01"
## [1] "Llh: -2.28967449713162E+03; Pars: 5.119488E-01, -5.465895E-02, 1.693815E+00, -8.918284E-02"
## [1] "Llh: -2.29508143572601E+03; Pars: 5.526720E-01, -1.279788E-02, 1.857131E+00, -1.380504E-01"
## [1] "Llh: -2.28778407481140E+03; Pars: 5.409059E-01, -2.412461E-02, 1.734745E+00, -1.092237E-01"
## [1] "Llh: -2.28174076179290E+03; Pars: 5.081352E-01, 1.904175E-02, 1.705954E+00, -1.057100E-01"
## [1] "Llh: -2.28026991581559E+03; Pars: 4.893179E-01, 6.256487E-02, 1.760550E+00, -1.219996E-01"
## [1] "Llh: -2.30101316146820E+03; Pars: 4.975114E-01, -6.292505E-02, 1.976295E+00, -3.717626E-02"
## [1] "Llh: -2.28417844436993E+03; Pars: 5.254901E-01, 6.891110E-03, 1.605962E+00, -1.106262E-01"
## [1] "Llh: -2.27768533978277E+03; Pars: 5.269030E-01, 5.267209E-02, 1.702617E+00, -9.382466E-02"
## [1] "Llh: -2.27495160506424E+03; Pars: 5.345431E-01, 1.063376E-01, 1.707036E+00, -9.614558E-02"
## [1] "Llh: -2.27438148613383E+03; Pars: 4.957536E-01, 8.736886E-02, 1.649115E+00, -6.724466E-02"
## [1] "Llh: -2.27468960750551E+03; Pars: 4.746111E-01, 1.431156E-01, 1.607898E+00, -4.625512E-02"
## [1] "Llh: -2.27733624961961E+03; Pars: 4.956484E-01, 9.659202E-02, 1.805132E+00, -4.415141E-02"
## [1] "Llh: -2.27619765311148E+03; Pars: 5.029461E-01, 7.416679E-02, 1.753136E+00, -6.077011E-02"
## [1] "Llh: -2.27704479513048E+03; Pars: 4.881392E-01, 2.145241E-01, 1.738328E+00, -1.489146E-01"
## [1] "Llh: -2.27542508000889E+03; Pars: 4.966689E-01, 1.485668E-01, 1.727570E+00, -1.177273E-01"
## [1] "Llh: -2.27278165052314E+03; Pars: 5.258083E-01, 1.456552E-01, 1.658533E+00, -4.894419E-02"
## [1] "Llh: -2.27482541054354E+03; Pars: 5.450616E-01, 1.872003E-01, 1.609764E+00, -1.241646E-02"
## [1] "Llh: -2.27401310058933E+03; Pars: 5.230582E-01, 1.697975E-01, 1.619985E+00, -1.042608E-01"
## [1] "Llh: -2.27356591164127E+03; Pars: 5.179560E-01, 1.458898E-01, 1.652294E+00, -9.338809E-02"
## [1] "Llh: -2.27262538439643E+03; Pars: 5.409006E-01, 9.405889E-02, 1.607744E+00, -3.513396E-02"
## [1] "Llh: -2.27326399324263E+03; Pars: 5.644725E-01, 6.680491E-02, 1.550985E+00, 6.162715E-03"
## [1] "Llh: -2.27414223569792E+03; Pars: 5.055558E-01, 1.301487E-01, 1.579054E+00, -2.620987E-02"
## [1] "Llh: -2.27314647656303E+03; Pars: 5.126518E-01, 1.241960E-01, 1.610120E+00, -4.369380E-02"
## [1] "Llh: -2.27296924723779E+03; Pars: 5.543248E-01, 1.675311E-01, 1.615074E+00, -4.333536E-02"
## [1] "Llh: -2.27263913573461E+03; Pars: 5.390632E-01, 1.474905E-01, 1.623518E+00, -4.931268E-02"
## [1] "Llh: -2.27230669377316E+03; Pars: 5.412659E-01, 1.098105E-01, 1.597868E+00, 4.845779E-03"
## [1] "Llh: -2.27289045018426E+03; Pars: 5.533113E-01, 9.177083E-02, 1.571332E+00, 5.396272E-02"
## [1] "Llh: -2.27183786799344E+03; Pars: 5.619208E-01, 1.243116E-01, 1.633470E+00, -2.057872E-02"
## [1] "Llh: -2.27172083174993E+03; Pars: 5.883035E-01, 1.243694E-01, 1.645271E+00, -9.021187E-03"
## [1] "Llh: -2.27183358053937E+03; Pars: 5.795050E-01, 9.220943E-02, 1.579436E+00, 4.633164E-03"
## [1] "Llh: -2.27181415742392E+03; Pars: 5.655874E-01, 1.055709E-01, 1.598849E+00, -8.761174E-03"
## [1] "Llh: -2.27149971973496E+03; Pars: 5.789952E-01, 6.941430E-02, 1.601196E+00, 2.527741E-02"
## [1] "Llh: -2.27170915106028E+03; Pars: 6.000571E-01, 3.037619E-02, 1.590151E+00, 6.257246E-02"
## [1] "Llh: -2.27177439266740E+03; Pars: 5.970022E-01, 1.105236E-01, 1.613609E+00, 4.130437E-02"
## [1] "Llh: -2.27153057056389E+03; Pars: 5.824536E-01, 1.064074E-01, 1.612141E+00, 2.219479E-02"
## [1] "Llh: -2.27119637272321E+03; Pars: 6.188823E-01, 9.307050E-02, 1.630817E+00, 9.999143E-03"
## [1] "Llh: -2.27157531461427E+03; Pars: 6.617695E-01, 8.470051E-02, 1.647546E+00, 1.257582E-02"
## [1] "Llh: -2.27144803950302E+03; Pars: 6.195454E-01, 9.105993E-02, 1.646030E+00, 3.298625E-02"
## [1] "Llh: -2.27128017128445E+03; Pars: 6.055915E-01, 9.468766E-02, 1.634106E+00, 2.254940E-02"
## [1] "Llh: -2.27128297631794E+03; Pars: 6.043112E-01, 5.742058E-02, 1.594149E+00, 4.903156E-02"
## [1] "Llh: -2.27117541646852E+03; Pars: 6.002689E-01, 7.415778E-02, 1.606779E+00, 3.451837E-02"
## [1] "Llh: -2.27132060941195E+03; Pars: 6.196451E-01, 5.925768E-02, 1.624202E+00, 2.397737E-02"
## [1] "Llh: -2.27117753372813E+03; Pars: 6.101303E-01, 7.104512E-02, 1.621178E+00, 2.353173E-02"
## [1] "Llh: -2.27145664872772E+03; Pars: 6.398868E-01, 9.706624E-02, 1.645476E+00, 2.002191E-02"
## [1] "Llh: -2.27121380883395E+03; Pars: 6.240884E-01, 9.015325E-02, 1.634293E+00, 2.133578E-02"
## [1] "Llh: -2.27111993201571E+03; Pars: 6.210571E-01, 6.952566E-02, 1.612430E+00, 2.214312E-02"
## [1] "Llh: -2.27128053588645E+03; Pars: 6.289374E-01, 5.694466E-02, 1.601700E+00, 2.193997E-02"
## [1] "Llh: -2.27133496617252E+03; Pars: 6.011847E-01, 6.374628E-02, 1.601425E+00, 2.376040E-02"
## [1] "Llh: -2.27112360432431E+03; Pars: 6.182819E-01, 8.355151E-02, 1.626013E+00, 2.194194E-02"
## [1] "Llh: -2.27125271173681E+03; Pars: 6.059478E-01, 5.606954E-02, 1.602472E+00, 4.106843E-02"
## [1] "Llh: -2.27112464001453E+03; Pars: 6.156230E-01, 8.382026E-02, 1.623685E+00, 1.776647E-02"
## [1] "Llh: -2.27105645064529E+03; Pars: 6.174002E-01, 8.448248E-02, 1.613246E+00, 2.465322E-02"
## [1] "Llh: -2.27106161681673E+03; Pars: 6.210676E-01, 9.120117E-02, 1.609295E+00, 2.521396E-02"
## [1] "Llh: -2.27120041480061E+03; Pars: 6.364349E-01, 8.653218E-02, 1.630976E+00, 8.733995E-03"
## [1] "Llh: -2.27110218853438E+03; Pars: 6.091130E-01, 7.725138E-02, 1.612794E+00, 2.807228E-02"
## [1] "Llh: -2.27105797220705E+03; Pars: 6.172719E-01, 7.358526E-02, 1.608572E+00, 3.063881E-02"
## [1] "Llh: -2.27106220389424E+03; Pars: 6.168593E-01, 7.614401E-02, 1.612337E+00, 2.742072E-02"
## [1] "Llh: -2.27113629034034E+03; Pars: 6.141151E-01, 6.887088E-02, 1.597631E+00, 3.081177E-02"
## [1] "Llh: -2.27107997253790E+03; Pars: 6.172376E-01, 7.988135E-02, 1.618871E+00, 2.415940E-02"
## [1] "Llh: -2.27113668607700E+03; Pars: 6.094881E-01, 8.807457E-02, 1.614304E+00, 3.161873E-02"
## [1] "Llh: -2.27107791817750E+03; Pars: 6.181444E-01, 7.416289E-02, 1.612898E+00, 2.451202E-02"
## [1] "Llh: -2.27102724347402E+03; Pars: 6.260297E-01, 7.880461E-02, 1.613992E+00, 2.390944E-02"
## [1] "Llh: -2.27100754367969E+03; Pars: 6.346634E-01, 7.958123E-02, 1.614591E+00, 2.182803E-02"
## [1] "Llh: -2.27100975040601E+03; Pars: 6.264494E-01, 7.602458E-02, 1.605806E+00, 2.665664E-02"
## [1] "Llh: -2.27102169121931E+03; Pars: 6.241336E-01, 7.698877E-02, 1.609062E+00, 2.603233E-02"
## [1] "Llh: -2.27098802866576E+03; Pars: 6.297185E-01, 8.267389E-02, 1.608205E+00, 2.737633E-02"
## [1] "Llh: -2.27097844146434E+03; Pars: 6.355865E-01, 8.692938E-02, 1.605863E+00, 2.880848E-02"
## [1] "Llh: -2.27097806761779E+03; Pars: 6.398953E-01, 8.992358E-02, 1.611171E+00, 2.033437E-02"
## [1] "Llh: -2.27100218695689E+03; Pars: 6.515161E-01, 9.809274E-02, 1.612473E+00, 1.518215E-02"
## [1] "Llh: -2.27091655862498E+03; Pars: 6.513130E-01, 8.174690E-02, 1.605471E+00, 2.416054E-02"
## [1] "Llh: -2.27088146909527E+03; Pars: 6.689617E-01, 8.037911E-02, 1.601597E+00, 2.391420E-02"
## [1] "Llh: -2.27095985071666E+03; Pars: 6.633278E-01, 9.238207E-02, 1.610793E+00, 2.078590E-02"
## [1] "Llh: -2.27094043814539E+03; Pars: 6.539095E-01, 8.829270E-02, 1.609545E+00, 2.225358E-02"
## [1] "Llh: -2.27092352307997E+03; Pars: 6.645961E-01, 9.318116E-02, 1.599525E+00, 2.582729E-02"
## [1] "Llh: -2.27091796876974E+03; Pars: 6.569831E-01, 8.978118E-02, 1.603278E+00, 2.482748E-02"
## [1] "Llh: -2.27092175293084E+03; Pars: 6.747087E-01, 8.725890E-02, 1.606922E+00, 1.685634E-02"
## [1] "Llh: -2.27091306209247E+03; Pars: 6.647080E-01, 8.717652E-02, 1.606657E+00, 1.984437E-02"
## [1] "Llh: -2.27086682509025E+03; Pars: 6.830352E-01, 8.289118E-02, 1.599383E+00, 2.508545E-02"
## [1] "Llh: -2.27086237565551E+03; Pars: 7.056838E-01, 7.937497E-02, 1.593521E+00, 2.746098E-02"
## [1] "Llh: -2.27084713297786E+03; Pars: 6.943535E-01, 8.006319E-02, 1.593010E+00, 2.576993E-02"
## [1] "Llh: -2.27085518006072E+03; Pars: 7.155041E-01, 7.594844E-02, 1.584806E+00, 2.752811E-02"
## [1] "Llh: -2.27091215696754E+03; Pars: 7.104888E-01, 7.371572E-02, 1.594107E+00, 2.366727E-02"
## [1] "Llh: -2.27087588677190E+03; Pars: 6.967171E-01, 7.773209E-02, 1.596395E+00, 2.395732E-02"
## [1] "Llh: -2.27087327859312E+03; Pars: 7.189411E-01, 7.159816E-02, 1.585666E+00, 3.070684E-02"
## [1] "Llh: -2.27085796540213E+03; Pars: 7.049815E-01, 7.549275E-02, 1.590888E+00, 2.799123E-02"
## [1] "Llh: -2.27089461956217E+03; Pars: 7.333499E-01, 7.595239E-02, 1.585349E+00, 2.867553E-02"
## [1] "Llh: -2.27086214092764E+03; Pars: 6.845083E-01, 7.927243E-02, 1.597519E+00, 2.510453E-02"
## [1] "Llh: -2.27084135589626E+03; Pars: 6.979384E-01, 7.936959E-02, 1.591075E+00, 2.920602E-02"
## [1] "Llh: -2.27083815509281E+03; Pars: 6.985499E-01, 8.018834E-02, 1.588421E+00, 3.183036E-02"
## [1] "Llh: -2.27084125208953E+03; Pars: 6.855784E-01, 7.813338E-02, 1.591392E+00, 2.788705E-02"
## [1] "Llh: -2.27084172607331E+03; Pars: 6.905504E-01, 7.844378E-02, 1.591924E+00, 2.778053E-02"
## [1] "Llh: -2.27083687794774E+03; Pars: 7.073392E-01, 7.766640E-02, 1.584361E+00, 3.163475E-02"
## [1] "Llh: -2.27085124223263E+03; Pars: 7.190386E-01, 7.686338E-02, 1.577823E+00, 3.489986E-02"
## [1] "Llh: -2.27084122187887E+03; Pars: 6.879450E-01, 8.253290E-02, 1.587699E+00, 3.056982E-02"
## [1] "Llh: -2.27083538680448E+03; Pars: 6.921651E-01, 8.077287E-02, 1.588495E+00, 2.992517E-02"
## [1] "Llh: -2.27083903569221E+03; Pars: 6.973734E-01, 7.831730E-02, 1.583336E+00, 3.486874E-02"
## [1] "Llh: -2.27083452143695E+03; Pars: 6.966172E-01, 7.875377E-02, 1.585749E+00, 3.259403E-02"
## [1] "Llh: -2.27085121924225E+03; Pars: 7.119629E-01, 8.055730E-02, 1.582133E+00, 3.510512E-02"
## [1] "Llh: -2.27083502421861E+03; Pars: 6.920814E-01, 7.873936E-02, 1.589072E+00, 2.969156E-02"
## [1] "Llh: -2.27083359131171E+03; Pars: 6.954998E-01, 7.777786E-02, 1.585417E+00, 3.009240E-02"
## [1] "Llh: -2.27083795866057E+03; Pars: 6.939797E-01, 7.657263E-02, 1.583916E+00, 2.922341E-02"
## [1] "Llh: -2.27084105006576E+03; Pars: 6.810850E-01, 8.035553E-02, 1.590008E+00, 2.951683E-02"
## [1] "Llh: -2.27083382178229E+03; Pars: 7.006822E-01, 7.833868E-02, 1.585771E+00, 3.110527E-02"
## [1] "Llh: -2.27083494067415E+03; Pars: 7.002853E-01, 7.603197E-02, 1.584510E+00, 3.181646E-02"
## [1] "Llh: -2.27083358941568E+03; Pars: 6.982464E-01, 7.721720E-02, 1.585505E+00, 3.134364E-02"
## [1] "Llh: -2.27083595973503E+03; Pars: 7.034825E-01, 7.730439E-02, 1.582157E+00, 3.287611E-02"
## [1] "Llh: -2.27083355077447E+03; Pars: 6.949143E-01, 7.838062E-02, 1.587340E+00, 3.048770E-02"
## [1] "Llh: -2.27083582682430E+03; Pars: 6.980472E-01, 7.710341E-02, 1.586268E+00, 2.892047E-02"
## [1] "Llh: -2.27083339776385E+03; Pars: 6.969744E-01, 7.834118E-02, 1.585878E+00, 3.167564E-02"
## [1] "Llh: -2.27083422561824E+03; Pars: 6.921588E-01, 7.751975E-02, 1.586299E+00, 3.069442E-02"
## [1] "Llh: -2.27083330733011E+03; Pars: 6.985416E-01, 7.813395E-02, 1.585903E+00, 3.100256E-02"
## [1] "Llh: -2.27083459309273E+03; Pars: 6.988396E-01, 7.825861E-02, 1.586897E+00, 3.216237E-02"
## [1] "Llh: -2.27083321545078E+03; Pars: 6.963332E-01, 7.789805E-02, 1.585787E+00, 3.060989E-02"
## [1] "Llh: -2.27083337932565E+03; Pars: 6.951364E-01, 7.915970E-02, 1.586949E+00, 3.054426E-02"
## [1] "Llh: -2.27083320545602E+03; Pars: 6.959126E-01, 7.867408E-02, 1.586588E+00, 3.074410E-02"
## [1] "Llh: -2.27083339475006E+03; Pars: 6.989711E-01, 7.814301E-02, 1.584739E+00, 3.152840E-02"
## [1] "Llh: -2.27083318870854E+03; Pars: 6.979547E-01, 7.820241E-02, 1.585389E+00, 3.126822E-02"
## [1] "Llh: -2.27083323893680E+03; Pars: 6.973950E-01, 7.811306E-02, 1.585955E+00, 3.013674E-02"
## [1] "Llh: -2.27083312823853E+03; Pars: 6.972898E-01, 7.817009E-02, 1.585936E+00, 3.052147E-02"
## [1] "Llh: -2.27083312639941E+03; Pars: 6.952067E-01, 7.833837E-02, 1.585946E+00, 3.056928E-02"
## [1] "Llh: -2.27083342133196E+03; Pars: 6.935452E-01, 7.844057E-02, 1.585968E+00, 3.035265E-02"
## [1] "Llh: -2.27083319308045E+03; Pars: 6.968471E-01, 7.879442E-02, 1.586142E+00, 3.094165E-02"
## [1] "Llh: -2.27083309645179E+03; Pars: 6.967186E-01, 7.857033E-02, 1.586053E+00, 3.085871E-02"
## [1] "Llh: -2.27083317982760E+03; Pars: 6.976719E-01, 7.796652E-02, 1.585074E+00, 3.086474E-02"
## [1] "Llh: -2.27083308697778E+03; Pars: 6.972317E-01, 7.814341E-02, 1.585453E+00, 3.083458E-02"
## [1] "Llh: -2.27083314189318E+03; Pars: 6.952702E-01, 7.840869E-02, 1.586305E+00, 3.012380E-02"
## [1] "Llh: -2.27083306798886E+03; Pars: 6.959404E-01, 7.835712E-02, 1.586076E+00, 3.040991E-02"
## [1] "Llh: -2.27083314270251E+03; Pars: 6.952594E-01, 7.853452E-02, 1.585829E+00, 3.081477E-02"
## [1] "Llh: -2.27083306616990E+03; Pars: 6.967817E-01, 7.826120E-02, 1.585909E+00, 3.059479E-02"
## [1] "Llh: -2.27083314261661E+03; Pars: 6.981322E-01, 7.832766E-02, 1.585799E+00, 3.077971E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083307201279E+03; Pars: 6.973302E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083319761197E+03; Pars: 6.966332E-01, 7.933569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083365834863E+03; Pars: 6.966332E-01, 7.733569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083390368515E+03; Pars: 6.966332E-01, 7.833569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083348649613E+03; Pars: 6.966332E-01, 7.833569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083336290147E+03; Pars: 6.966332E-01, 7.833569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083322808316E+03; Pars: 6.966332E-01, 7.833569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083323525009E+03; Pars: 6.945464E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083328354192E+03; Pars: 6.952413E-01, 7.933569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083373598408E+03; Pars: 6.952413E-01, 7.733569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083369705280E+03; Pars: 6.952413E-01, 7.833569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083385579516E+03; Pars: 6.952413E-01, 7.833569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083344543677E+03; Pars: 6.952413E-01, 7.833569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083330871799E+03; Pars: 6.952413E-01, 7.833569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083319761197E+03; Pars: 6.966332E-01, 7.933569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083328354192E+03; Pars: 6.952413E-01, 7.933569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083413695022E+03; Pars: 6.959369E-01, 8.033569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083390955123E+03; Pars: 6.959369E-01, 7.933569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083382542187E+03; Pars: 6.959369E-01, 7.933569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083403552339E+03; Pars: 6.959369E-01, 7.933569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083290041296E+03; Pars: 6.959369E-01, 7.933569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083365834863E+03; Pars: 6.966332E-01, 7.733569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083373598408E+03; Pars: 6.952413E-01, 7.733569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083505018306E+03; Pars: 6.959369E-01, 7.633569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083441029920E+03; Pars: 6.959369E-01, 7.733569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083423769333E+03; Pars: 6.959369E-01, 7.733569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083349278633E+03; Pars: 6.959369E-01, 7.733569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083435635235E+03; Pars: 6.959369E-01, 7.733569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083390368515E+03; Pars: 6.966332E-01, 7.833569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083369705280E+03; Pars: 6.952413E-01, 7.833569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083390955123E+03; Pars: 6.959369E-01, 7.933569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083441029920E+03; Pars: 6.959369E-01, 7.733569E-02, 1.587496E+00, 3.062189E-02"
## [1] "Llh: -2.27083578950335E+03; Pars: 6.959369E-01, 7.833569E-02, 1.589084E+00, 3.062189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083408610720E+03; Pars: 6.959369E-01, 7.833569E-02, 1.587496E+00, 3.162189E-02"
## [1] "Llh: -2.27083396823238E+03; Pars: 6.959369E-01, 7.833569E-02, 1.587496E+00, 2.962189E-02"
## [1] "Llh: -2.27083348649613E+03; Pars: 6.966332E-01, 7.833569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083385579516E+03; Pars: 6.952413E-01, 7.833569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083382542187E+03; Pars: 6.959369E-01, 7.933569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083423769333E+03; Pars: 6.959369E-01, 7.733569E-02, 1.584324E+00, 3.062189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083553271644E+03; Pars: 6.959369E-01, 7.833569E-02, 1.582741E+00, 3.062189E-02"
## [1] "Llh: -2.27083397601882E+03; Pars: 6.959369E-01, 7.833569E-02, 1.584324E+00, 3.162189E-02"
## [1] "Llh: -2.27083382226277E+03; Pars: 6.959369E-01, 7.833569E-02, 1.584324E+00, 2.962189E-02"
## [1] "Llh: -2.27083336290147E+03; Pars: 6.966332E-01, 7.833569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083344543677E+03; Pars: 6.952413E-01, 7.833569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083403552339E+03; Pars: 6.959369E-01, 7.933569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083349278633E+03; Pars: 6.959369E-01, 7.733569E-02, 1.585909E+00, 3.162189E-02"
## [1] "Llh: -2.27083408610720E+03; Pars: 6.959369E-01, 7.833569E-02, 1.587496E+00, 3.162189E-02"
## [1] "Llh: -2.27083397601882E+03; Pars: 6.959369E-01, 7.833569E-02, 1.584324E+00, 3.162189E-02"
## [1] "Llh: -2.27083419896770E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.262189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083322808316E+03; Pars: 6.966332E-01, 7.833569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083330871799E+03; Pars: 6.952413E-01, 7.833569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083290041296E+03; Pars: 6.959369E-01, 7.933569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083435635235E+03; Pars: 6.959369E-01, 7.733569E-02, 1.585909E+00, 2.962189E-02"
## [1] "Llh: -2.27083396823238E+03; Pars: 6.959369E-01, 7.833569E-02, 1.587496E+00, 2.962189E-02"
## [1] "Llh: -2.27083382226277E+03; Pars: 6.959369E-01, 7.833569E-02, 1.584324E+00, 2.962189E-02"
## [1] "Llh: -2.27083306414405E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 3.062189E-02"
## [1] "Llh: -2.27083392747133E+03; Pars: 6.959369E-01, 7.833569E-02, 1.585909E+00, 2.862189E-02"## 2. Plot the combined density functions
ghddens <- dghyp(ef$x, ghdfit)
hypdens <- dghyp(ef$x, hypfit)
nigdens <- dghyp(ef$x, nigfit)
nordens <- dnorm(ef$x, mean = mean(HDTimeS, na.rm=TRUE), sd = sd(c(HDTimeS[, 1]), na.rm=TRUE))
col.def <- c("black", "red", "blue", "green", "orange")
plot(ef, xlab = "", ylab = expression(f(x)), ylim = c(0, 0.25))
lines(ef$x, ghddens, col = "red")
lines(ef$x, hypdens, col = "blue")
lines(ef$x, nigdens, col = "green")
lines(ef$x, nordens, col = "orange")
legend("topleft",
legend = c("empirical", "GHD", "HYP", "NIG", "NORM"),
col = col.def, lty = 1)
## 3. Create a Q-Q plot
qqghyp(ghdfit, line = TRUE, ghyp.col = "red", plot.legend = FALSE,
gaussian = FALSE, main = "", cex = 0.8)
qqghyp(hypfit, add = TRUE, ghyp.pch = 2, ghyp.col = "blue",
gaussian = FALSE, line = FALSE, cex = 0.8)
qqghyp(nigfit, add = TRUE, ghyp.pch = 3, ghyp.col = "green",
gaussian = FALSE, line = FALSE, cex = 0.8)
legend("topleft", legend = c("GHD", "HYP", "NIG"),
col = col.def[-c(1,5)], pch = 1:3)
## 4. Make a model recommendation using lik.ratio.test
LRghdnig <- lik.ratio.test(ghdfit, nigfit)
LRghdnig## $statistic
## L
## 0.3652671
##
## $p.value
## [1] 0.155828
##
## $df
## [1] 1
##
## $H0
## [1] TRUELRghdhyp <- lik.ratio.test(ghdfit, hypfit)
LRghdhyp## $statistic
## L
## 0.6478822
##
## $p.value
## [1] 0.3514839
##
## $df
## [1] 1
##
## $H0
## [1] TRUE## Both tests have no signficant result to say there is a relationship between data sets## 5. Calculate and plot the VaR (using all models)
p <- seq(0.001, 0.05, 0.001)
ghd.VaR <- abs(qghyp(p, ghdfit))
hyp.VaR <- abs(qghyp(p, hypfit))
nig.VaR <- abs(qghyp(p, nigfit))
nor.VaR <- abs(qnorm(p, mean = mean(HDTimeS, na.rm=TRUE), sd = sd(c(HDTimeS[, 1]), na.rm = TRUE)))
emp.VaR <- abs(quantile(x = HDTimeS, probs = p, na.rm=TRUE)) plot(emp.VaR, type = "l", xlab = "", ylab = "VaR", axes = FALSE,
ylim = range(c(hyp.VaR, nig.VaR, ghd.VaR, nor.VaR, emp.VaR)))
box()
axis(1, at = seq(along = p), labels = names(emp.VaR), tick = FALSE)
axis(2, at = pretty(range(emp.VaR, ghd.VaR, hyp.VaR,
nig.VaR, nor.VaR)))
lines(seq(along = p), ghd.VaR, col = "red")
lines(seq(along = p), hyp.VaR, col = "blue")
lines(seq(along = p), nig.VaR, col = "green")
lines(seq(along = p), nor.VaR, col = "orange")
legend("topright",
legend = c("Empirical", "GHD", "HYP", "NIG", "Normal"),
col = col.def, lty = 1)
## 6. Calculate and plot the ES (using all models)
ghd.ES <- abs(ESghyp(p, ghdfit))
hyp.ES <- abs(ESghyp(p, hypfit))
nig.ES <- abs(ESghyp(p, nigfit))
nor.ES <- abs(mean(HDTimeS, na.rm=TRUE) - sd(c(HDTimeS[, 1]), na.rm = TRUE) *
dnorm(qnorm(1 - p)) / p)
obs.p <- ceiling(p * length(HDTimeS))
emp.ES <- sapply(obs.p, function(x) abs(mean(sort(c(HDTimeS))[1:x])))plot(emp.ES, type = "l", xlab = "", ylab = "ES", axes = FALSE,
ylim = range(c(hyp.ES, nig.ES, ghd.ES, nor.ES, emp.ES), na.rm = TRUE))
box()
axis(1, at = 1:length(p), labels = names(emp.VaR), tick = FALSE)
axis(2, at = pretty(range(emp.ES, ghd.ES, hyp.ES, nig.ES, nor.ES)))
lines(1:length(p), ghd.ES, col = "red")
lines(1:length(p), hyp.ES, col = "blue")
lines(1:length(p), nig.ES, col = "green")
lines(1:length(p), nor.ES, col = "orange")
legend("topright",
legend = c("Empirical", "GHD", "HYP", "NIG", "Normal"),
col = col.def, lty = 1)