search for: likelyhood

...t;attempt to find suitable starting values failed", for example with : require(MASS) data(iris) polr(Species~Sepal.Length+Sepal.Width+Petal.Length+Petal.Width,iris) (I know the response variable Species should be nominal but I do as levels were ordered for the example). I think this is a likelyhood maximization problem ; I tried to solve this by setting the "start" option of polr to a null or a random vector by it doesn't garantee to find "a good" solution at the end. Does anyone have a clue ? Thanks a lot ! -- View this message in context: http://r.789695.n4.nabb...

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

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

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

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

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 jon at jgcomp.com 11226 South Shore Rd. (703) 787-0688 (H) Resto...

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

...E 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 explicit tutorial about a poisson regression with MLE or, maybe better, about maximum likelyhood estimation? Best wishes Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 Mobile: 0044 (0)79 64925523 www.lshtm.ac.uk/pehru/ antonio.gasparrini at lshtm.ac.uk

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

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 -- "The best way to predict the future is to invent it." - XEROX PARC slogan, circa 1978 -- This message has been scanned for viru...

..., 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 Ph.D. Département de Biologie Université du Québec à Rimouski 300 Allée des Ursulines Rimouski, Québec, G5L 3A1 Tél. : 723-1986 poste 1653 E-mail : arnaud_mosnier@uqar.qc.ca [[alternative HTML version deleted]]

...ted 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 weight, so that after all 730,000 [x,y] are plotted, a maximum likelyhood becomes apparent through the size of the point that are made up by accumulated weights. Thank you in advance for any thoughts! Best, Akko Eleveld

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

...;Study"]&logit[,"Tr"]==1,] (number of row in each group (1-3) is difference but equal in colume) fnTr1 <- function (p) sum( n/2*log(2*pi)+log(1/logitTr1$sd)+1/2*(logitTr1$logitp*logitTr1$logitp-2*logitTr1$logitp*p+p^2) *1/logitTr1$sd*logitTr1$sd ) (maximum likelyhood function) outTr1<- nlm (fnTr1,p=c(10),hessian=TRUE) minimumTr1 <- outTr1$minimum valueTr1 <- outTr1$estimate (estimate the value) The problem is the program couldn't work by using logitTr[i],fnTr[i] outTr[i],minimumTr[i],value[i] in a loop. The function logitTr[i] is the data matrix...

...a(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 (checked with SAS 9.1) except nu (degrees of freedom of the student) and the log-likelyhood. I've really checked everything and can't find the estimated series sigma (volatility) and eta, such that eps = sigma * eta and eta is centered and reduced... I've tryed combinations of all s$x,s$h,s$z and nothing looks looks correct. Also, is it possible to fit EGARCH and TGARCH with...

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

...rrently doesn't fill in a frame number, so MSC will always be zero for now, and this makes mplayer choose target_msc = 0 as 'swap at next interval', which appears to be broken by the commit. The following patch fixes frame numbers in the nouveau kernel module, but it only decreases the likelyhood of occuring, it still happens when seeking, for example. ~Maarten --- diff --git a/drivers/gpu/drm/nouveau/nouveau_display.c b/drivers/gpu/drm/nouveau/nouveau_display.c index 4610c3a..fd3cb62 100644 --- a/drivers/gpu/drm/nouveau/nouveau_display.c +++ b/drivers/gpu/drm/nouveau/nouveau_display.c @@...

...ergeometric" show this term in the vignettes of package's ?combinat? and ?forward?. But i'm unable to make sense of the these pachakege functions in the context of my aforementioned apllication. Can some one suggest a function, script or method to achieve my goal of estimating the likelyhood of overlap between 3 lists of genes, ideally using the multivariate hypergeometric, or anything else for that matter? cheers in advance, Karl #example attempt with two gene lists m & n N <- 45101 # total number balls in urn m <- 720 # number of 'white' or 'special'...