similar to: (no subject)

Chuck (Is it fine to call you Chuck?) has far more R jutsu than I do obviously. I don't know much about pmin and factor but it might worth looking into if you want to manipulate states by names (I assume this is why one might want to use it?) generate_transition_matrix <- function(data, states) prop.table(table(head(data, -1), tail(data, -1)), 1)[states,]

Ok, I assumed you wanted to compute a matrix M for all states such that M[i][j]= transition from state i to state j / number of transition 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

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

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 > > >

Hi, I think you overthought this one a little bit, I don't know if this is the kind of code you are expecting but I came up with something like that: generate_transition_matrix <- function(data, n_states) { #To be sure I imagine you should check n_states is right at this point transitions <- matrix(0, n_states, n_states) #we could improve a little bit here because at

I have not followed the history of this thread, but I am quite flummoxed as to why the OP is rewriting code to estimate parameters from an univariate Gaussian mixture model when alternatives such as EMCluster (which generally appears to handle initialization better than MClust) exist. Or perhaps there is more to it in which case I apologize. But I thought that I would make the OP aware of the

---------- 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 gmail.com>

Dear Niharika, As I said before, the problem is basically an optimization issue. You should isolate the problematic part from the rest of your study. Sometimes, more information does not help to solution. All the answers from us (Ulrik, David, me) are more or less are correct to find a maximum point. Newton?s method is also correct. But after answers, you only say, it didn?t find the right

... and following up on Spencer's answer, see the "digits" argument of ?options and ?print.default. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Aug 22, 2017 at 7:30 AM, Spencer Graves <spencer.graves

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),

On 2017-08-22 9:26 AM, niharika singhal wrote: > 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 No:? That's only the display: > sum(v)-1 [1] 1.6e-07 ????? hope this helps.? Spencer

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

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 =

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 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

webroot. On Windows and Centos that's easy as that's just ONE directory. On Debian and it's ancestors, Ubuntu being one of them, there are TWO directories and he put the files in the wrong one of both and doesn't symlink them from the right one. I tried to explain this to him on IRC and I summarized this here on the mailing list. I'd suggest anyone with web serving related

> 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

Greetings! I just signed up to be a part of the CentOS docs Wiki contributors. At the moment, I do not have a proposed subject or location. I presume that once I have a content of some sort, that's when I propose the subject and location of my Wiki contribution(s). In the meantime, I'd like to request a personal homepage. My username is AdrianCho. If you need any further information

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

--- Ryan Lane <ryanlane@gmail.com> wrote: > On Fri, 11 Mar 2005 07:21:39 -0800 (PST), Steven > Vishoot > <sir_funzone@yahoo.com> wrote: > > ok now i am frustrated. this is the third day in a > row > > that my server has crashed running centos4, i can > not > > ping or access the console but the hardware is > still > > running. it does not leave