similar to: EM algorithm

Dear mailing list, I would like to use the optim() command in order to maximize the logged likelihood of the following function, where p is the parameter of interest and should be constrained between 0 and positive infinity. y = 1/2 * ((te - x)/(te - tc))^p x and y are given by x <- c(5.18, 6.28, 7.00, 7.08, 7.54, 7.90, 8.24, 8.64, 12.17, 12.89, 14.27, 15.38, 15.80, 16.46, 20.41, 21.27,

Hi everyone, I am looking to do some manual maximum likelihood estimation in R. I have done a lot of work in Stata and so I have been using output comparisons to get a handle on what is happening. I estimated a simple linear model in R with lm() and also my own maximum likelihood program. I then compared the output with Stata. Two things jumped out at me. Firstly, in Stata my coefficient

I've converted the bdsmatrix package (used by coxme) to the newer style of S4, and uploaded it. It is in the pkg directory of survival, on r-forge. I'm stuck on something that is almost certainly a namespace issue. I've borrowed liberally from Matrix (a big help) and read the documentation, and just upgraded to 2.8.1 - but still stuck. All the tests work when I run them "by

Jean-Mark, Tim, Could either of you expound on the following comment in cwrs.c? /*If _k==0, the following do-while loop will overflow the buffer.*/ ---------------------------------------------------------------- ...because the following do-loop does overflow the buffer when k=126 k=2; do _u[k]=(k<<1)-1; while(++k<len); Thanks, MikeH -------------- next

Hi all: I'm testing a new installation of *, bringing up some ATA186. In * environment, all stuff works greats. The only thing that don't work is a Call Transfer, but the 3Party works ok. Some time ago I read that somebody had proven this functionality successfully. If somebody knows what I missing, please let me know. Thanks in advance, Gus -------------- next part -------------- An

On 13 January 2015 at 08:21, Dan Tenenbaum wrote: | Where should the package source be downloaded from? I see it in CRAN (but presumably the latest version that causes the issue is not yet downloadable) and in github. The "presumable" assumption is incorrect AFAIK. The error should presumably came up in both versions as annoylib.h did not change there. Feel free to prove me wrong :)

Hi, I needed a function to check if a file is existing on the Debian box I am configuring. So I used this one module Puppet::Parser::Functions newfunction(:file_exists, :type => :rvalue) do |args| if File.exists?(args[0]) return 1 else return 0 end end end It work greats for some service where apache2 doesn''t work for me. I am using it this way

Hello, I notice that these entities are defined like this: def emdash opts "&#8212;" end def endash opts "&#8211;" end def arrow opts "&#8594;" end However, any space from the original markup is stripped, so "hello -- world" becomes "hello&#8212;world". I''ve overridden the methods to deal with

Hello, I want to have a layout which is resizable by the user. So I am useing em for widths, heights, etc of the elements. However if I resize the screen (via browser) the draggables do not change their position according to the resizeing. Is it possible to do that? TIA, Alex -- Posted via http://www.ruby-forum.com/.

hello can somebody help me hackin the moveby effect so i can use em''s instead of px? any feedback very much appriciated. thanx alot & have a nice time, daniel :)

I have 2 questions concerning the EM algorithm. Is it true that the EM algorithm gives unique answers for the means and variances of a mixture of 2 normals? I am using mixtools and I am surprised that it works better than a Bayesian program I wrote. If so can someone say why the mixing probabilities are so good? -- Thanks, Jim. [[alternative HTML version deleted]]

Atencao: openssh-unix-dev at mindrot.org Um virus foi encontrado numa mensagem de Email que acabou de ser enviada por voce. Este scanner de Email a interceptou e impediu a mensagem de chegar no seu destino. O virus foi reportado como sendo: Worm.Somefool.AR Por favor atualize seu antivirus ou contate o seu suporte tecnico o mais rapido possivel pois voce tem um virus no seu computador.

When I call EM::stop from within a Deferrable callback, I''m getting the following: terminate called after throwing an instance of ''std::runtime_error'' what(): already initialized Abort trap Is there any way to trace this a little further? -- Duane Johnson

Hi, I am looking for a package to do the clustering analysis using the expectation maximization algorithm. Thanks in advance. Ming

What repository should I pull from to get the latest snapshot of the EM java code for jruby? cr

I couldn't find a direct answer in CRAN to this question, so I'm asking with some trepidation. I have a multivariate dataset (data.frame) with columns that can be expressed as a set of mixed normals (at least I think) and need to impute values that have constraints (truncated mixture of normals where the values cannot be below zero). If there isn't a package that can do this, is there

Atencao: samba@samba.org. Um Virus foi encontrado numa mensagem de Email que acabou de ser enviada por voce. Este scanner de Email a interceptou e impediu a mensagem de chegar no seu destino. O Virus foi reportado como sendo: the W32/Netsky.p@MM virus !!! Por favor atualize seu antivirus ou contate o seu suporte tecnico o mais rapido possivel pois voce tem um virus no seu computador.

Hi, I am trying to understand how the functions em() and me() from the mclust package work. I cannot make sense of what the algorithm returns. Here is a basic, simple example: ######################################### # two bivariate normals, centered at (-5,0) and (5,0), with Id covariance x1 = cbind(rep(-5, 100), rep(0, 100)) + matrix(rnorm(100*2), 100, 2) x2 = cbind(rep(5, 100),

hi, I want to use EM for MLE estimation.WHAT would be the relevant package or code?? Regards, Subhajit. [[alternative HTML version deleted]]

Hi, I woul like to know if it is possible to have a "R code" to generate EM Algorithm for a normal bivariate mixture. Best regard, S.F.