Verildi

"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: R Programlama Dili, İstatistiksel Analiz, İstatistik

Daha fazlasını görün: GINI Index as goodness function , Using GINI Index as goodness function i.e. CART Algorithm., vipcollection space index php r ad271e0895b056de97db868c4751642c, http onlinevission info index php r shahadat, forbidden you don t have permission to access index php on this server additionally a 404 not found error was encountered while, a href http www kayads com index php r spgo img src http www kayads com img banner2 png border 0 width 468 height 60 a, gini index international trade, calculate correlation gini index openness excel, calculate correlation openness gini index, normal facebook log, statistics php normal distribution, normal distribution curve, normal distribution curve php, normal distribution, pligg beta site index

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

Proje NO: #15594267

1 freelancer is bidding on average $20 for this job

karanrajput1112

Your work will be handled with confidentiality. Professionalism is my core value. • descriptive statistics • chi square and t test • Analysis of variance • correlation and regression. • and lot more. Relevant Skills a Daha fazlası

in 0 gün içinde20$ USD
(42 Değerlendirme)
5.4