Kapalı

"maximum likelihood estimation of beta-normal" & "simulation the gini index of beta-normal distribution" in R

Regarding "beta-normal distribution formula" which has been uploaded, I want to debug my scripts to estimate the parameters of the beta-normal distribution.

I've used the maxLik package in r:

------------------------------------------------------------------------------------------------------------------

library(VGAM)

library(maxLik)

alfa=4;beta=3;mu=0;sigma=1

n=100

x=rbetanorm(n,alfa,beta,mu,sigma)

logLikFun=function(w){

alfa=w[1]

beta=w[2]

mu=w[3]

sigma=w[4]

z=(x-mu)/sigma

ll={-1*(n*gamma(alfa+beta)-n*log(gamma(alfa))-n*log(gamma(beta))-(n/2)*log(2*pi)-n*log(sigma)+sum((alfa-1)*(pnorm(z,mean=0,sd=1))+(beta-1)*(log(1-pnorm(z,mean=0,sd=1)))-((z^2)/2)))}

ll

}

mle=maxLik(logLikFun,start=c(alfa=5,beta=2.5,mu=1,sigma=1))

summary(mle)

------------------------------------------------------------------------------------------------------------------

the results are quite wrong as shown below:

------------------------------------------------------------------------------------------------------------------

------------------------------------------—

Maximum Likelihood estimation

Newton-Raphson maximisation, 3 iterations

Return code 3: Last step could not find a value above the current.

Boundary of parameter space?

Consider switching to a more robust optimisation method temporarily.

Log-Likelihood: -131157.7

4 free parameters

Estimates:

Estimate Std. error t value Pr(> t)

alfa 5.1037852 NA NA NA

beta 2.3962128 NA NA NA

mu 1.5954536 0.0099872 159.8 <2e-16 *

sigma 0.0444610 0.0001726 257.6 <2e-16 *

—-

Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘’ 0.05 ‘.’ 0.1 ‘ ’ 1

—------------------------------------------

Warning messages:

1: In sqrt(diag(vc)) : NaNs produced

2: In sqrt(diag(vc)) : NaNs produced

------------------------------------------------------------------------------------------------------------------

the problem is the NA value for the errors which is not acceptable.

I would be pleased if someone could solve this problem and send me the current scripts which I can edit or add it to my scripts.

More over, regarding "gini index of beta-normal distribution" formula which has been uploaded, I want to simulate this formula by the above method

please inform me about the duration of the project.

Beceriler: Veri Bilimi, Makine Öğrenimi, R Programlama Dili, İstatistiksel Analiz, İstatistik

Daha fazlasını gör: maximum likelihood excel vba, fminsearch matlab maximum likelihood, fminsearch maximum likelihood, maximum likelihood matlab code fminsearch, maximum likelihood matlab fminsearch, estimation parameters using maximum likelihood via algorithm, estimation parameters using maximum likelihood via, matlab code maximum likelihood estimate, fminsearch maximum likelihood matlab, matlab fminsearch maximum likelihood, maximum likelihood matlab, fminsearch function matlab maximum likelihood, maximum likelihood matlab image, maximum likelihood image processing matlab, maximum likelihood expectation maximization matlab

İşveren Hakkında:
( 0 değerlendirme ) Turkey

Proje NO: #15451253

Bu iş için 5 freelancer ortalamada $33 teklif veriyor

aki003iitr

Hi, I have more than 4 years of experience in working on R as a data scientist and can finish the task and give you the correct scripts in a few hours . I have completed 23 projects in R and your requirements overlap w Daha Fazla

1 gün içinde %bids___i_sum_sub_32%%project_currencyDetails_sign_sub_33% USD
(29 Değerlendirme)
5.5
mohamedrdait

Hello sir I can achieve this project perfectly using R Relevant Skills and Experience I Master R and statistics analysis please check my work history and my portfolio for more details Proposed Milestones $30 USD - mi Daha Fazla

1 gün içinde %bids___i_sum_sub_32%%project_currencyDetails_sign_sub_33% USD
(26 Değerlendirme)
4.4
in %bids___i_period_sub_35% gün içinde35%project_currencyDetails_sign_sub_37% %project_currencyDetails_code_sub_38%
(23 Değerlendirme)
4.1
1 gün içinde %bids___i_sum_sub_32%%project_currencyDetails_sign_sub_33% USD
(0 Değerlendirme)
0.0
kanyekicharles

i got enough knowledge in statistcs to hundle this

1 gün içinde %bids___i_sum_sub_32%%project_currencyDetails_sign_sub_33% USD
(0 Değerlendirme)
0.0