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

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 cannot figure out how to draw a certain plot: could someone help me out? I have this data.frame from a survey my.data that looks like something like this: col1 col2 col3 col4 1 5 5 4 5 2 3 5 3 1 3 2 3 4 5 4 3 1 1 2 5 5 5 4 5 6 4 2 5 5 .... Each row represents a single questionnaire

Hi,guys, HTML help is not right for math formula: For example: I got the following formula in help page of ca.jo in urca package: <p align="center"><i><B>X</B>_t = <B>Pi</B>_1 <B>X</B>_{t-1} + ... + <B>Pi</B>_k <B>X</B>_{t-k} + <B>&mu;</B> + <B>Phi D</B>_t +

Hello, I've ran into problem with Dovecot and dsync replication. Everything works perfectly, including replication of sieve scripts, except fact that if user activates the 'managesieve' ruleset (I'm using currently Roundcubemail) on "mail1" host, it wouldn't be activated on "mail2" host, by creating symlink ".dovecot.sieve ->

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

Hello Everybody, We develop a Windows based radio broadcast system for professional radios and we would like to implant an interface in this system to make it able to automatically update metadata of an OGG (!) stream (recorded from the radio) sent to an Icecast 2.0.2 server. We read about the curl --user admin:hackme

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, in the last days I switched from Xen 4.2.1 to Xen 4.2.2 (source from git repository tag 4.2.2) commit a125ec6a4a2f028f97f0bdd8946d347a36dbe76b Author: Jan Beulich <jbeulich@suse.com> Date: Tue Apr 23 18:42:55 2013 +0200 update Xen version to 4.2.2 My "base system" is a debian squeeze. For Xen I use the xl toolstack. Since updating to 4.2.2 the autosave and

I''m actually looking to create my own markup language, after a look to the greats markdown and textilize markup languages they seem to be more complexes that what I need, so not feating with my application. Any help is welcome ;-) -- Grosjean S?bastien

Hello I am new to this post. My current setup is this. 2003 R2 with Identity managment installed - I have statically mapped Unique Unix attributes (UID and GID)to each user. Multiple CentOS 4.4 servers with Samba 3.0.25a-32 installed. Everything works greats wbinfo -u -g -t, getent passwd group. But when I access any shares I get this error message. [2007/06/11 11:27:35, 1]

I have a new Sun Fire AMD64 server running Solaris 10 12/06 and Samba 3.0.21b. I have migrated a configuration file from an older RHEL machine running Samba 2.2.x. On this new machine, if one tries to copy a large number of files from a Windows XP client to the server, the copy starts and runs for a few minutes, but then stops. The error message on the XP machine is "Cannot copy

Hello. I have a problem with theora when using 2 video screens. The second is connected to the dvi output of the video card. When I watch the video, I move the renderer to the other screen and video still remains black, or freezes, for a while. Considerer this test : construct a graph DirectShow using GraphEdit like this : - Webcam filter - Theora Encoder - Theora Decoder - Video Renderer Run

Am 12.02.2015 um 15:47 schrieb Martin ?tefany: > Hello, > > I've ran into problem with Dovecot and dsync replication. Everything > works perfectly, including replication of sieve scripts, except fact > that if user activates the 'managesieve' ruleset (I'm using currently > Roundcubemail) on "mail1" host, it wouldn't be activated on

Hello all. Again I call upon the greats of this mailing list to help me with an issue related to MEMDISK - notably when paired with various Dell OptiPlex and Latitude systems. The issue presented itself sometime ago and I best dealt with it by remaining with version 3.83 of the SYSLINUX package. But valuing the improved changes fostered by newer versions, I would prefer to move forward with

Hello, I just build my new Rails application, i placed the website like that : at http://www.***.net/ = /var/www/devkedare/public/ (with mod_rails) http://media.***.net/ = /var/www/media/ I want to use the javascript :cache features, the "timed" files to prevent caching when file change (ex: x.js?1249517571), and all others greats files features, I configured the asset_host to use the

Is it possible to send logs to syslog and files set by log_path and log_info_path on config? Just two logs at the same time (not realtime). -- Antonio P?rez-Aranda Alcaide aperezaranda at yaco.es Yaco Sistemas S.L. http://www.yaco.es/ C/ Rioja 5, 41001 Sevilla Tel?fono +34 954 50 00 57 Fax ? ? ?+34 954 50 09 29