search for: niharikasinghal1990

...states) prop.table(table(head(data, -1), tail(data, -1)), 1)[states,] checkdf=data.frame(clusterNum=c(3,2,3,1,1,3,1,3,2,1,1,3,2,1,3,2)) states=c(2,3) transition_matrix= generate_transition_matrix(checkdf$clusterNum, states) transition_matrix 2017-08-28 18:09 GMT+02:00 niharika singhal <niharikasinghal1990 at gmail.com>: > yes or S_i= number of state > > M[S_i][S_j]= transition from state S_i to state S_j / number of > transition from state S_i where S_i is in the set S > For example 1 row of 3 states will be > x= s1->s1/no of transition in s1 > y= s1->s...

Please keep the list in cc. Sorry, it didn't work as expected. Maybe someone else have an appropriate solution. Best, Ulrik On Sa., 26. Aug. 2017, 12:57 niharika singhal <niharikasinghal1990 at gmail.com> wrote: > Hi > > Thanks for you mail, > I really appreciate your time on my problem > > I have posted this problem on > > > https://stats.stackexchange.com/questions/299590/to-find-maxima-for-gaussian-mixture-model > > > The plot I am getting usi...

...tion from state i but from what you just answered it looks like you want to compute a matrix M for a set of states S such that: M[S_i][S_j]= transition from state S_i to state S_j / number of transition from state S_i where S_i is in the set S ? 2017-08-28 17:50 GMT+02:00 niharika singhal <niharikasinghal1990 at gmail.com>: > Sorry I have typed wrong value for output > (2/6, 1/6, 3/6) > > On Mon, Aug 28, 2017 at 5:49 PM, niharika singhal < > niharikasinghal1990 at gmail.com> wrote: > >> Hi Elie, >> >> Thanks for the mail >> >> If the number of st...

Hi keep your emails to R help, I do not offer private consultance and others could have different opinion how to solve your problem. Did you even try my suggestion? If not, why not? If yes in what respect it does not comply with your expectations. Cheers Petr From: niharika singhal [mailto:niharikasinghal1990 at gmail.com] Sent: Tuesday, October 10, 2017 2:34 PM To: PIKAL Petr <petr.pikal at precheza.cz> Subject: Re: [R] "Time Series Plotting" Hello, My Energy column looks like in the image in the attachment, and my sequence column correspond to state. I wanted to look at the graph and...

...ns[k,] <- transitions[k,] / sum(transitions[k,]) transitions } checkdf=data.frame(clusterNum=c(3,2,3,1,1,3,4,3,2,1,1,3,2,1,3,2)) no_of_state=4 transition_matrix= generate_transition_matrix(checkdf$clusterNum, no_of_state) transition_matrix 2017-08-28 16:37 GMT+02:00 niharika singhal <niharikasinghal1990 at gmail.com>: > Hello, > > I am trying to implement a formula > > aij= transition from state S_i to S_j/no of transition at state S_i > > > > Code I have written is working with three state {1,2,3 }, but if the number > of states become={1,2,3,4,......n} then the...

Hello, I am trying to implement a formula aij= transition from state S_i to S_j/no of transition at state S_i Code I have written is working with three state {1,2,3 }, but if the number of states become={1,2,3,4,......n} then the code will not work, so can some help me with this. For and some rows of my data frame look like

> On 26 Aug 2017, at 16:39, niharika singhal <niharikasinghal1990 at gmail.com> wrote: > > Hi, > > Thanks for your mail, and time > > It is not working for some arguments, when mean value is like >6. > > > case > > mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. > 1134463,0.2074027) > > rv &...

Hi, Thanks for your mail, and time It is not working for some arguments, when mean value is like >6. case mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. 1134463,0.2074027) rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133), Norm(664.0713, 3.674773), Norm(669.0484, 4.101381),

---------- Forwarded message ---------- From: niharika singhal <niharikasinghal1990 at gmail.com> Date: Sun, Aug 27, 2017 at 11:57 AM Subject: Re: Find maxima of a function To: "David Winsemius [via R]" <ml+s789695n4745009h56 at n4.nabble.com>, "Ismail SEZEN [via R]" <ml+s789695n4744993h60 at n4.nabble.com>, Ulrik Stervbo <ulrik.stervbo at gma...

I have a Gaussian mixture model with some parameters mean=(506.8644,672.8448,829.902) sigma=(61.02859,9.149168,74.84682) c=(0.1241933, 0.6329082, 0.2428986) And the plot look something like below.[image: enter image description here] <https://i.stack.imgur.com/4uUQ9.png> Also, if I change my parameters to mean=(2.15,2.0,2.9) sigma=(0.1,0.1,0.1) c=(1/3,1/3,1/3) Then plot would change

...to investigate your failed case. At this point, you can use ?dput? function to extract, copy and paste your x-y vectors. So, we can copy the statement, run in our environment without requiring distr or any other package. Best wishes, isezen > On 27 Aug 2017, at 12:59, niharika singhal <niharikasinghal1990 at gmail.com> wrote: > > Dear David and Ismail, > > The actual problem is I am getting the parameters from the Kmeans cluster > on the data set obtained from the mclust package. > > In mclust method I am changing the value of G according to user input, so > the number...

Hello, I have initial parameters for HMM model and I want to use depmixS4 package. The parameters are in the form intial_prob_matrix=matrix(c(0.07614213, 0.45177665, 0.47208122), nrow=1, ncol=3, byrow = TRUE) transition_matrix=matrix(c(0.46666667,0.46666667,0.06666667, 0.06741573,0.5617978,0.37078652, 0.02173913,0.3478261,0.63043478), nrow = 3, ncol =

...At this point, you can use ?dput? function to extract, copy and paste your x-y vectors. So, we can copy the statement, run in our environment without requiring distr or any other package. > > Best wishes, > > isezen > > > > On 27 Aug 2017, at 12:59, niharika singhal <niharikasinghal1990 at gmail.com> wrote: > > > > Dear David and Ismail, > > > > The actual problem is I am getting the parameters from the Kmeans cluster > > on the data set obtained from the mclust package. > > > > In mclust method I am changing the value of G according...

Hello, I am having some problem in retrieving the list elements Code Mean_of_mc is a list as described below mean_mc_1 mean_mc_2 mean_mc_3 mean_mc 0.9577246 2.000516 2.962521 mean_mc 7.957725 9.000516 9.962521 optim(par=c(lVal=2, mu=Mean_of_mc), fn=xnorm, method = "BFGS", control=list(fnscale=-1)) xnorm= function(par) { print(par)

Hello I have a vector v=c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.1134463,0.2074027) when i do sum(v) or 0.08844446+0.1744455+0.1379778+0.1209769+0.1573065+0.1134463+0.2074027 i am getting output as 1 But if i add them manually i get 1.00000026 I do not want to round of my value since it effect my code further Can anyone suggest how can i avoid this. Thanks & Regards Niharika

Hello, I need some help in plotting time series. I have dataframe Data with two column and thousands of row, I want wherever the sequence corresponding to Energy column is changed the color change should be reflected in Time Series plot, some rows of dataframe are below Energy sequence 13536 3 14335 3 14638 3 25363 3 18511 2 18371 2 14555 3 I am able

Hello, I am facing a problem with optimization in R from 2-3 weeks. I have some Gaussian mixtures parameters and I want to find the maximum in that *Parameters are in the form * mean1 mean2 mean3 sigma1 sigma2 sigma3 c1 c2 c3 506.8644 672.8448 829.902 61.02859 9.149168 74.84682 0.1241933 0.6329082 0.2428986 I have used optima and optimx to find the

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