similar to: likelyhood maximization problem with polr

Dear friends i am trying to fit an Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD method. i found these formulas on http://www.sitmo.com/article/calibrating-the-ornstein-uhlenbeck-model/ this is the mean-reverting process http://r.789695.n4.nabble.com/file/n4637271/process.txt process.txt and this is the script that i am using....... ouFit.ML=function(spread) { n=length(spread)

Hello. I have a wireless system (link bandwidth = 700 kbit/s) that transmits realtime udp video data and other various tcp protocols. I don''t care if I lose the occassional udp packet but i''d like all of the tcp traffic to get through. I''m thinking that a Generalized RED will give me the best results: 2 queues, one for TCP (very low probability of dropping) and one

I am trying to optimize custom likelyhood with nlminb() Arguments h and f are meant to be fixed. example.R: compute.hyper.log.likelyhood <- function(a, h, f) { a1 <- a[1] a2 <- a[2] l <- 0.0 for (j in 1:length(f)) { l <- l + lbeta(a1 + f[j], a2 + h - f[j]) - lbeta(a1, a2) } return(l) } compute.optimal.hyper.params <- function(start, limits, h_, f_) { result

What is the likelyhood that major inexpensive webhosts like godaddy, bluehost, etc. will upgrade to RoR 1.1? Is this going to be like PHP 5 where it has to percolate for a year or more before it becomes widly available? Your thoughts? Along the same lines... is it possible to adopt some of the new improved Ajax / javascript capabilities without actually upgrading the ruby installation?

Hi everyone, I am wondering what the likelyhood would be to get a "delete machine script" implemented in the samba core that would delete a machine upon successful un-join of a samba domain. Alternatively, why is this not a good idea? tom.

A friend runs 6.8 but would like to add postscreen to his email server. It is provided with postfix version 2.8 and up. The newest postfix for 6.8 is 2.6.x. Any likelyhood it will be further updated? Or is there another source for newer postfix versions for CentOS 6.8? Similar query for the clamav db, now at 0.99.1 with 0.99.2 out for quite some time. Jon -- Jon H. LaBadie

Begging indulgence of the list - asking here because due to the higher likelyhood of getting an answer... I have an application that often loses connection from one machine to another when one (or more) of the machines has particular brands of 1gig nics, but which runs rock solid when on 10meg nics and some 1 gig nics. The application senses (falsely) that the connection has been shut down

Dear all, I'm trying to estimate the parameters of a special case of a poisson model, where the specified equation has an integral and several fixed parameters. I think that the MLE command in STATS4 package could be a good choice, but it's a little complicated. I've got some problems with the offset and I don't understand some of the functions. Do you know where can I find some

I am testing the likelyhood of updating a couple RH9 systems to centos 5 386. I found at document that talked about updating to centos 4. steps being: 1.) backup 2.) rpm --import http://pubmirrors.reflected.net/centos/5/os/i386/RPM-GPG-KEY-centOS-5 rpm -Uvh --nodeps http://pubmirrors.reflected.net/centos/5/os/i386/CentOS/yum-3.0.5-1.el5.centos.5.noarch.rpm

Dear R Users, May be this message should be directy send to Douglas Bates ... I just want to know if I can use the AIC value given in the output of an lmer model to classify my logistic models. I heard that the AIC value given in GLIMMIX output (SAS) is false because it come from a calculation based on pseudo-likelyhood. Is it the same for lmer ??? thanks, Arnaud Arnaud MOSNIER Biologiste

Does anyone here on this list have experience with HA clustering? I'm previewing drbd as a potential tool, and wanted to know if anyone here has experiemented with it at all... How stable is it? Does the additional likelyhood of failure given the additional complexity actually get compensated by a better overall system? http://www.drbd.org/ Any feedback is welcome... -Ben --

Anyone who can help me with the following question? How can I add weight to [x,y] coordinates on a graph/scatterplot? Background: Monte Carlo simulation generated 730,000 [x,y] coordinates with a weight attached (from 0-0.5). Both x and y are rounded and fit on a raster with x-axis 0-170 months (smalles unit = 1 month) and y-axis 0-6 (smallest unit=0.1). I would like every [x,y] to add its

Begging indulgence of the list - asking here because due to the higher likelyhood of getting an answer... I have an application that often loses connection from one machine to another when one (or more) of the machines has particular brands of 1gig nics, but which runs rock solid when on 10meg nics and some 1 gig nics. The application senses (falsely) that the connection has been shut

Dear everyone I am a new user,would like to combine these code together by using a loop,each function has three value as Tr = 1 - 3,how can i combine together? logitTr1 <-logit[logit[,"Study"]&logit[,"Tr"]==1,] (number of row in each group (1-3) is difference but equal in colume) fnTr1 <- function (p) sum(

I'm trying to use garchFit from fSeries, with Student or Skewed Student conditionnal distribution. Let's say that eps (vector) is my series of daily log-returns: data(EuStockMarkets) eps = diff(log(EuStockMarkets[,"CAC"])) library(fSeries) g = garchFit(series = eps, formula.var = ~garch(2,2), cond.dist = "dstd") s = g at fit$series All the coefficients are ok

So there is no way around having to re-encode them all? --Pritpaul --- Geoff Shang <gshang@uq.net.au> wrote: > Hi: > > My guess is that ices can't handle 16000 or 32000hz > sampling rates. The > 32kbps file is 22050hz and the 80kbps file is > 32000hz. I couldn't get > mpg123 to play test2 remotely, but lame defaults to > 32000hz with 40kbps > mono and

There's one vulnerability that's bugged me for some time. It applies to nearly all crypto software, including ssh. That's the swapping of sensitive info (such as keys and key equivalents) onto hard drives where they could possibly be recovered later. The Linux kernel provides a system call, mlock(), that inhibits swapping of a specified region of virtual memory. It locks it into real

Hey Mario, It seems that your ddx commit b4231dd715a8 is causing a regression when I use mplayer -vo vdpau rendering with mesa 9.1. It fails to start drawing here, leaving the mplayer screen black, I can make it recover usually by seeking or moving another window in front. When I revert the commit it behaves normally. Changing can_sync_to_vblank(draw) to 0 in nouveau_dri2_schedule_swap also fixes

Dear R Users, I employed the phyper() function to estimate the likelihood that the number of genes overlapping between 2 different lists of genes is due to chance. This appears to work appropriately. Now i want to try this with 3 lists of genes which phyper() does not appear to support. Some googling suggests i can utilize the Multivariate hypergeometric distribution to achieve this. eg.:

Full_Name: Greg Kochanski Version: 2.2.1 OS: Debian Linux (testing) Submission from: (NULL) (212.159.16.190) This is really a feature request. When you do mosaicplot() on a data set where the probability of several nearby rows is small, then the labels for those rows are plotted overlapping each other. This situation can be improved by calling mosaicplot() with a large value of