search for: mll

...best way of fitting the same negative log-likelihood function to more than one set of data, using mle() from the stats4 package. Here's what I would have thought would work: -------------- library(stats4) ## simulate values r = rnorm(1000,mean=2) ## very basic neg. log likelihood function mll <- function(mu,logsigma) { -sum(dnorm(r,mean=mu,sd=exp(logsigma),log=TRUE)) } mle(minuslogl=mll,start=list(mu=1,logsigma=0)) r2 = rnorm(1000,mean=3) ## second "data set" with(list(r=r2), mle(minuslogl=mll,start=list(mu=1,logsigma=0)) ) ------------- but this doesn'...

...ate the # samples, I use BIC/AIC with MLE (computed from dgamma function). Here is the code I have. __BEGIN__ mlogl <- function( x_func, theta_func, samp) { # computing log_likelihood return( - sum(dgamma(samp, shape = x_func, scale=theta_func, log = TRUE))) } find_bic <- function(mll,smpl,k) { bic <- (-2 * mll) + (k * log(length(smpl))) bic } find_aic <- function(mll,smpl,k) { aic <- (-2 * mll) + (k * 2) aic } mlogl_process <- function(smpl,error,start ) { # EM algorithm to estimate max loglikelihood thetalast <- 0 thetac...

Dear R list, I have a function specifying my log-likelihood, and now I need to set the constraint that *alpha > kappa*, could anyone help me with setting this in my function? the function is defined as follows: mll <- function(param){ n <- length(x) psi <- numeric(n) psi[1] <- 1.0 a0 <- exp(param[1]); a1 <-exp(param[2]); b1 <- exp(param[3]); *alpha *<- exp(param[4]); *kappa *<- exp(param[5]); for (i in 2:n) {psi[i] <- a0 + a1*x[i-1] + b1*psi[i-1]} lam <- gamma(kap...

...xt._(text) end bindtextdomain should only be called once on initialization and than cached for every language. I want to select the language based on the current request. Is there anyone out there already using Gettext for localization? If yes, how do you do it? Sascha Ebach PS: I know of the MLL project, but I want to be able to use gettext tools to work with the language files (MLL uses YAML to store the strings not the native PO/MO format). MLL extends the the Base classes of ActionView for example, but I want to be able to translate everything in the application, so also error messa...

...TO DO: torture-test on some real problems! better documentation? e.g. ?profile.mle-class doesn't give details on arguments -- have to look at profile.nls (allow "which" to be a character vector -- match names)? HARDER: fancy formula interface [cf. svymle in survey package] e.g. mll <- mLL(type="independent",distrib="normal",resp=y,mean=~a+b*x,sd=~s, param=~a+b+s) allow for fitting of transformed parameters (exp/log, tanh/atanh = logistic/logit) 2D profiles (quadratic or thin-plate spline???) EASIER but breaking backward compatibility: merge ab...

...e got .... RootDev <- 600 splits <- NULL for (cutpoint in cutpoints) { LRGdata <- split(gdata, gdata[,4] <= cutpoint) vgrmL <- variog(as.geodata(LRGdata$"TRUE"), covar.col=covcol)) vgrmR <- variog(as.geodata(LRGdata$"FALSE"), covar.col=covcol)) mlL<- likfit(as.geodata(LRGdata$"TRUE"), ini=expand.grid(seq(0, max(vgrmL$v), l=10), seq(0,max(vgrmL$u), l=10))) mlR<- likfit(as.geodata(LRGdata$"FALSE"), ini=expand.grid(seq(0, max(vgrmR$v), l=10), seq(0,max(vgrmR$u), l=10))) LeftDev <- summary(mlL)[[8]]$log.L...

Greetings, I have a meta-analysis problem in which I have fixed effects regression coefficients (and estimated standard errors) from identical models fit to different data sets. I would like to use these results to create pooled estimated regression coefficients and estimated standard errors for these pooled coefficients. In particular, I would like to estimate the model \beta_{i} = \mu +

...3fF+9Y5fDcn338RsbeScqHv+DXIt2O5RcSWEkze4XQyzx6W6fjSXvy6KHW/7IgNf SMS0uVhBBreQkT4ywHC14hKd3c7gpjPSwyjYGtDNktigLjaaGLpQm4ua6NFEOppo4GtMuWH0tXXM maJCUZiUJ2UxPfokbqqhOZ97B3Qk7dycndhM7ZkV64vyoHa+NMa9g5znRhLtBSVtg3nnxcqk4V/0 H2a0uy/QnxDv/MpPq3+q5KuXEXd67Ywej61DyLLTz+j1XruT0Rm+euFFLkZ69Vjc1yf2oEGv7b6D wRA+mlL/QQi0PkuLcg6NOL2madEHc5Hcv0EvFQTdIl4m/Xg0t1kBv4dzeJ+8ps1RLwHjDlotwwUI jZp2+Q/f6f7hULB4nW63333lZ5ofpyw786ZZs28aO/3mwLmdnmkddJumrXu6pq3jZq1pvH08actn tDno+WzyNq0t+3ZOp2VTbdv3cVqX4p4knExcAtl+bm+w+fL8DpIl5nl8AzorW6dRyufptTs/MzMl XPIHUEsDBAoAAAAGAEt5bSsmmHjeiQMAAOkEAAANAAAAUjJkaXN0aWxsLnBkZg8HAQIzFAJMAgQc...

Hello: I wish to again express my appreciation to all who have contributed to making R what it is today. At this moment, I'm particularly grateful for whoever modified the 'mle' code so data no longer need be passed via global variables. I remember struggling with this a couple of years ago, and I only today discovered that it is no longer the case. I'd

Still don't recommend you use 'make install', but here we go. Rich.

Hi, I want to create a data frame including a column containing some special characters, like "ø". when I print that data frame out, the content change to <U+00F8>, and when save the data frame to a txt file, the content keep in that style, but I need it in its original form, anybody can explain? > x <- data.frame(part = c("målløs", "ny")) > x

Hi All I have a symmetric matrix of genes ( 100x100 matrix). I also have a matrix (100x2) of two columns where column 1 has the gene names and column 2 has the cluster it belongs to (they are sorted and grouped based on the cluster no). I would like to order the rows and columns of the 100x 100 matrix such that the first n genes correspond to cluster 1 and next n genes correspond to cluster 2

...cSoft$ dir PhotoImpression\ 5 dstock76 at dstock76-desktop:~/.wine/drive_c/Program Files/ArcSoft$ cd PhotoImpression\ 5 dstock76 at dstock76-desktop:~/.wine/drive_c/Program Files/ArcSoft/PhotoImpression 5$ dir Albums creative_content EXIF.dll Modules Pfc1033.dll Pfc.dll Res_Frame.mll UI_Framework AlignSplit.dll DGUI.dll EzDll.dll msvcrt.dll Pfc1034.dll Pfc.inf Share Background.jpg Director.dll FPXLIB.DLL PadusDll.dll Pfc1036.dll photoimpression.exe ShareObjs.mll CheckUpdate dtype.inf gdiplus.dll Pcdlib32.dll Pfc1040.dll PHOTOIMPRESSION.hlp SUPPORT...

This is something I've been working on: Reworking inspection so it's not a big mess of ad hoc C code, but instead uses a well-defined domain-specific language to describe how we inspect guests. The best introduction to this is the manual page, which I include below (it's also included in patch 2/3). Rich. ---------------------------------------------------------------------- NAME

...$(builddir) $^ | \ $(SED) 's/ *$$//' | \ $(SED) -e :a -e '/ *\\$$/N; s/ *\\\n */ /; ta' | \ - $(SED) 's/ :/:/' | \ + $(SED) -e 's,$(abs_srcdir)/,$(builddir)/,g' | \ sort > $@-t mv $@-t $@ @@ -167,7 +167,7 @@ SUFFIXES = .cmo .cmi .cmx .ml .mli .mll .mly if HAVE_PERLDOC -supermin.1: supermin.pod +supermin.1: $(srcdir)/supermin.pod pod2man \ -u \ --section 1 \ @@ -178,12 +178,12 @@ supermin.1: supermin.pod noinst_DATA = \ ../html/supermin.1.html -../html/supermin.1.html: supermin.pod +../html/supermin.1.html: $(srcdir)/super...

- run bindtests.opt only if available - use ocamlc -custom - Don't try to install native code if it hasn't been built (Thanks to Olaf Hering) --- mllib/Makefile.am | 2 +- ocaml/Makefile.am | 7 ++++++- ocaml/examples/Makefile.am | 11 ++++------- ocaml/run-bindtests | 2 ++ 4 files changed, 13 insertions(+), 9 deletions(-) diff --git a/mllib/Makefile.am b/mllib/Makefile.am index 2e35d1e..0e317b2 100644 --- a/mllib/M...

Previously posted: https://www.redhat.com/archives/libguestfs/2015-December/msg00038.html Inspection now really succeeds on a small number of simple guests. To test it out: $ ./run guestfish -v -x -a /tmp/centos-6.img ><fs> run ><fs> debug sh "guestfs-inspection --verbose" Rich.

This is a more working version. Inspection (partially) succeeds on a real guest this time :-) You can test it out on a real guest (in this case, a CentOS disk image located at /tmp/centos-6.img) by doing: $ ./run guestfish -v -x -a /tmp/centos-6.img ><fs> run ><fs> debug sh "guestfs-inspection --verbose" which will print lots of debugging, and at the end the

For background on this change, see: https://rwmj.wordpress.com/2015/12/06/inspection-now-with-added-prolog/ v2 was previously posted here: https://www.redhat.com/archives/libguestfs/2015-December/msg00038.html To test this patch series on a real guest, you can do: $ ./run guestfish -v -x -a /var/tmp/centos-6.img ><fs> run ><fs> debug sh "guestfs-inspection

...2 +- inspector/Makefile.am | 4 ---- java/Makefile.am | 3 +-- java/examples/Makefile.am | 6 ++--- lua/Makefile.am | 4 +--- lua/examples/Makefile.am | 4 ---- make-fs/Makefile.am | 4 ---- mllib/Makefile.am | 6 ----- ocaml/Makefile.am | 4 +--- ocaml/examples/Makefile.am | 7 ++---- p2v/Makefile.am | 12 ++-------- perl/Makefile.am | 6 ++--- perl/examples/Makefile.am | 4 ----...