similar to: Constraint maximum (likelihood) using nlm

Displaying 20 results from an estimated 500 matches similar to: "Constraint maximum (likelihood) using nlm"

2005 Sep 29
1
Error using a data frame as the "start" parameter in mle()
Dear R-Users, I am trying to use mle() to optimize two (or more) parameters, but I want to specify those parmeters in a data frame rather than having to spell them out separately in the "start" variable of mle(). My call is > mle(negll, start=list(aps=init), fixed=list(measphot=newphot, formod=formod, Nbands=Nbands), method="BFGS") where negll is a function I have
2006 Aug 26
1
problems with loop
Dear all, I am trying to evaluate the optimisation behaviour of a function. Originally I have optimised a model with real data and got a set of parameters. Now I am creating simulated data sets based on these estimates. With these simulations I am estimating the parameters again to see how variable the estimation is. To this end I have written a loop which should generate a new simulated data
2011 Aug 17
2
An example of very slow computation
This message is about a curious difference in timing between two ways of computing the same function. One uses expm, so is expected to be a bit slower, but "a bit" turned out to be a factor of >1000. The code is below. We would be grateful if anyone can point out any egregious bad practice in our code, or enlighten us on why one approach is so much slower than the other. The problem
2012 Mar 03
3
[RESEND][RFC GIT PATCHES] acpioff: COM32 module to shut off machine using ACPI
[resending, cutting out git pull-request output to make message smaller] Hi all, I have written a COM32 module called "acpioff", that, not surprisingly, powers off a machine using ACPI. I have tested it only using pxelinux.0 with a SeaBIOS/Qemu virtual machine. Since the changeset pulls in a "Linux-ized" version of the open source ACPI Component Architecture, it was too big
2007 Dec 11
1
R computing speed
Dear helpers, I am using R version 2.5.1 to estimate a multinomial logit model using my own maximum likelihood function (I work with share data and the default function of R cannot deal with that). However, the computer (I have an Athlon XP 3200+ with 512 GB ram) takes quite a while to estimate the model. With 3 categories, 5 explanatory variables and roughly 5000 observations it takes 2-3 min.
2017 Nov 20
1
Intel Bluetooth problems - "Invalid request code (56)" - on CentOS 6.9
Hi, I have Lenovo laptop with an Intel wireless and Bluetooth adapter, running CentOS 6.9. The wireless works just fine, but I can't seem to get Bluetooth to work. When I try to force it up from the command line I get # hciconfig hci0 up Can't init device hci0: Invalid request code (56) There is also a similar message message in /var/log/messages from start-up on boot. The device
2007 Jan 08
4
Export dataframe to txt
Hi all, Is there a function to export a dataframe to a text file? I want to store a large set of data which I have saved in a dataframe in my workspace and copy and past doesn't cut it. Thank you, Benjamin
2011 Mar 01
1
How to prove the MLE estimators are normal distributed?
Dear List, I'm now working on MLE and OSL estimators.I just noticed that the textbook argues they are joint normal distributed.But how to prove the conclusion? Thanks for your time in advance! Best, Ning
2014 Oct 08
0
Open Software License v. 3.0
Dear All, I would like to propose adding the OSL-3.0 license to the list of "standard" licenses bundled with R: Index: share/licenses/license.db =================================================================== --- share/licenses/license.db (revision 66733) +++ share/licenses/license.db (working copy) @@ -317,3 +317,12 @@ URL:
2006 Oct 26
4
Header of dataframe
Hi I am fairly new to R and I would appreciate some help to hopefully a trivial problem. I created a function: summary.aggregate <- function(y, ...) { temp.mean <- aggregate(y, FUN=mean, ...) temp.sd <- aggregate(y, FUN=sd, ...) temp.length <- aggregate(y, FUN=length, ...) temp <- data.frame(cbind(mean=temp.mean$x,stdev=temp.sd$x,n=temp.length$x)) } this outputs e.g.:
2023 Jun 24
2
Mirror problems with elfutils-debuginfod-client
The package elfutils-debuginfod-client is needed for even a minimal install, but it is not available on most mirrors. I suspect some are excluding mirroring debuginfo packages with just a *debuginfo* pattern to rsync, where they should do something like *-debuginfo-*.rpm (which should be good for now as I don't see any package with just "debuginfo" in the name, even in Fedora). The
2007 Mar 14
4
abs(U) > 0 where U is a vector?
Hi, I am looking for a way to compare if every element of a vector is > 0. i.e. while(abs(U) > 0) { .. } is there a function for this or do I have to write one? I'd appreciate your help! Benjamin
2007 Jan 08
3
Speeding things up
Hi, is it possible to do this operation faster? I am going over 35k data entries and this takes quite some time. for(cnt in 2:length(sdata$date)) { if(sdata$value[cnt] < sdata$value[cnt - 1]) { sdata$ddtd[cnt] <- sdata$ddtd[cnt - 1] + sdata$value[cnt - 1] - sdata$value[cnt] } else sdata$ddtd[cnt] <- 0 } return(sdata) Thank you, Benjamin
2007 Jan 15
4
How to format R code in LaTex documents
Hi, I am planning on putting some R script in an appendix of a LaTex document. Can anyone recommend me a way of how to format it? Is there a way to keep all line breaks without having to insert \\ in every single line? Thank you! Benjamin
2007 Jan 09
2
Logical operations or selecting data from data.frames
Hi all, why doesn't something like this does not work? speedy <- (sdata$VaR < sdata$DdtdAbs) && sdata$DdtdDuration >= qpois(pct,lambda) && sdata$Ddtd > MinDD or sdata$Ddtd[sdata$Ddtd > 0 && sdata$VaR < sdata$DdtdAbs] sdata looks like this: dataId date value Ddtd VaR DdtdAbs DdtdDuration 18948 79637 2004-07-27 10085.10
2007 Nov 06
1
[RFC] Enable _OSI(Linux) for Xen HVM domains
I''d like to send this patch to upstream lkml/linux-acpi, but I want to make sure that the x86 side wants this too. The _OSI method is a way for platform ACPI firmware to determine the features and capabilities of the OS running on the system. This is usually done by calling the _OSI method with various strings in the _INI method of an ACPI device. The strings are things like
2006 Nov 07
2
Boxplot
Hi, I am new to R and am still trying to get a grip of it. I have data in this format: Run Lab Batch Y 1 1 1 1 608.781 2 2 1 2 569.670 3 3 1 1 689.556 4 4 1 2 747.541 5 5 1 1 618.134 6 6 1 2 612.182 7 7 1 1 680.203 8 8 1 2 607.766 9 9 1 1 726.232 and I want to make a boxplot of the Y values for each
2018 Oct 04
0
Copy to smb share fails with "invalid argument" on CentOS 7
Hi, I've had problems copying files to Windows shares from my CentOS 7 machine lately. I originally got this in the desktop file manager, but find that I can also reproduce using gvfs-copy. "cp" to the directory mounted by gvfs works just fine, on the other hand. Also, the problem does not occur with small files - I think anything below 64k is OK. The following command sequence
2012 Mar 03
0
[RFC GIT PATCHES] acpioff: COM32 module to shut off machine using ACPI
Hi all, I have written a COM32 module called "acpioff", that, not surprisingly, powers off a machine using ACPI. I have tested it only using pxelinux.0 with a SeaBIOS/Qemu virtual machine. Since the changeset pulls in a "Linux-ized" version of the open source ACPI Component Architecture, it was too big to post as a patches to the list. See the git pull-request output below.
2007 Feb 18
1
Help with pair plot
Hi, I'm trying to create a plot using pair. Currently I'm doing the following pairs(~mpg + hp + wt , data=cars, labels = c("Miles per Gallon","Horsepower","Weight"),pch = c(24,25)[unclass(cars$tr + 1)], bg = c("red", "green3")[unclass(cars$tr + 1)],panel=panel.smooth) for the attached dataset. However,instead of using panel.smooth I