Hello,
I have some parameters from Mclust function. The parameters are in the form
*parametersDf *
* mu_1 mu_2 var_mc1 var_mc2 c1
c2 *
*2 1.357283 2.962736 0.466154 0.1320129 0.5258975
0.4741025 *
*21 8.357283 9.962736 0.466154 0.1320129 0.5258975
0.4741025 *
Each row in the above data frame represent the segment and each segment
contain two-mixture component. (Latter I will increase the segment and the
mixture component)
So for every segment with n mixture component I have to find the
*max f(x) with –infinity < x < infinity*
*f(x)=sum_i c_i N(O,mu_i,sigma_i). *
*Since i need to calculate the derivate f’ and set f’(x)=0 I thought of
using Newton method *
*In R the function I selected for my problem is optim.*
*Below is my code *
*l = 2 # represent the number of mixture component*
*i=1# represent the segment 1 if i=2 represent the 2 segment*
*parametersDf #which I have posted above*
*xnorm= function(x)*
*{ *
* #observation= append(segment1, segment2)*
* mu.vector = parametersDf[1:l] *
* var.vector =parametersDf[(l+1):((2*l))]*
* coeff.vector=parametersDf[(2*l+1):(3*l)]*
* val_df=data.frame()*
* sum_prob=0*
* for (j in 1:l){*
*
val=(coeff.vector[i,j]*(dnorm(x,mu.vector[i,j],sqrt(var.vector[i,j]))))*
* sum_prob=sum_prob+val*
* } *
* return(sum_prob)*
*}*
*#segment1*
*val1=optim(parametersDf$mean_mc_1[1],** fn=xnorm,*
* method = "BFGS")*
*I am not getting right result from optim function. Is there something I
can do?*
*Thanks in Advance for the help *
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
