similar to: Bug or Feature? LogLik.nls and non-central F distribution.

Hello, I have been using R to fit my data using non-linear least squares. I have used the optimize routine to minimize the sum of squared errors (using Nealder Meade optimization routine), but couldn't get the non-linear model in R to converge to the estimates acheived in a convergent Neadler Meade routine. It tells me about problems with the gradient. I was wondering if there is any way

As the original coder (in mid 1970s) of BFGS, CG and Nelder-Mead in optim(), I've been pushing for some time for their deprecation. They aren't "bad", but we have better tools, and they are in CRAN packages. Similarly, I believe other optimization tools in the core (optim::L-BFGS-B, nlm, nlminb) can and should be moved to packages (there are already 2 versions at least of LBFGS

The question is general for all functions, but here is a specific example - # For the logistic regression of the following correlated variables: C <- c(457, 1371, 4113, 12339, 37017, 111051, 333153, 999459) E <- c(0.003858377, 0.014334578, 0.014092836, 0.737950754, 0.996371828, 0.997482379, 1.005569257, 0.994382856) # The nls function: A = nls(E~(Em*C^p)/(C50^p + C^p),

Hello list, my name is Matteo, i'm new on that list. I'm working on a multitenancy platform with linux containers through libvirt on a production system with Red Hat 6.4. Every container run a separate instance of OpenSSH and Apache HTTPd and I need to give root privileges to the developers and I try to configure SELinux using svirt and MCS. I try the secmodel type dynamic and static in

Colleagues, I am fitting an Emax model using nls. The code is: START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL <- list(maxiter=1000, warnOnly=T) #FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula FORMULA <- as.formula(YVAR ~ EMAX - EFFECT / (1 + (C50/XVAR)^GAMMA)) FIT <-

On Thu, Oct 31, 2013 at 04:32:45PM +0100, Matteo Piccinini wrote: > Hello list, > > my name is Matteo, i'm new on that list. > I'm working on a multitenancy platform with linux containers through libvirt on a production system with Red Hat 6.4. > Every container run a separate instance of OpenSSH and Apache HTTPd and I need to give root privileges to the developers and I

I was working with the rmetrics portfolioBacktesting function and dug into the code to try to find why my formula with 113 items, i.e. A1 thru A113, was being truncated and I only get 85 items, not 113. Is it due to a string length limitation in R or is it a bug in the strsplit or gsub functions, or in my string? I'd very much appreciate any suggestions ============Input script:

Hi, I have 22 files (A1, A2, ..., A22) with different number of columns, totaling 10,000 columns: c1, c2, c3, ..., c10000 I have another file with a list of 100 columns that I need to extract. These 100 columns are distributed in 22 files. How to extract the 100 columns of the 22 files? I have done it "manually" with the following commands, for example: cromo1 = read.table ("~

Hi All, I am trying to use the C50 package to build classification trees in R. Unfortunately there is not enought documentation around its use. Can anyone explain to me - how to prune the decision trees? Regards, Indrajit [[alternative HTML version deleted]]

I have a Asterisk server (PBX in a Flash) with Aastra 57i phones. When there is an incoming call the phone will display two buttons "answer" and "ignore". If you press "ignore" the call is dropped instead of sent to voice mail. The following is the log: -- Called 111 -- SIP/111-00001c14 is ringing -- Got SIP response 486 "Busy Here" back from

R users, I have df like this a <- data.frame( indx = 1:20, var1 = rep(c("I20", "I40", "A50", "B60"), each=5), var1_lab= rep(c("cat", "dog", "mouse", "horse"), each=5), var2 = rep(c("B20", "X40", "D50",

R users, My intention is to take factors out of DF, create list of tables and export these tables separately using write.table or sink function. write.table writes tables out as DF:s, should I use sink instead? Here is my example: a <- data.frame( indx = 1:20, var1 = rep(c("I20", "I40", "A50", "B60"), each=5),

Dear R users, I am a beginner in R. I have 506 text files (data frame) in one folder namely DATA. The files are called A1 to A253 (253 files) and B1 to B253 (another 253 files). Each file has two columns; V1 (row number) and V2 (the value for each row name). Now I would like to add the values of V2 in each A-file with its relative value in B-file and save it as a new data frame named as C (e.g. C1

Hi all, I'm having problems with the installation of CentOS 5 for i386 to my Dell Percision T4500 Workstation. The installation run just fine. The installation DVD had been tested before installation. After installation, the installer ask me to reboot. after reboot, when the msg for LVM saying detection of my LVMs, 1 or 2 lines after that my monitor turned off, but i believed the cpu is still

hi, i need to do a demo of asterisk in a office that doesn't have an ISDN line, can someone tell me if asterisk works also if i use a modem card mapped on /dev/ttyl0 insted of use a isdn card? did someone try to do something similar before? thanks for your help!! Angelo

hi there ;-) can you help ? ouer network has gone unreasonly slow newerdays, i?m nearly shure its no virus or shit.... i?d used tcpdump to see whats going on and found these errormessages ouer network is quite small, the linuxserver(samba), some windoofclients and the plotter_pc ( this one runs as a MSdosclient and searches a pathstructure for new files) what is this SNAP length ??? and

On Mon, Mar 4, 2019 at 5:01 PM J C Nash <profjcnash at gmail.com> wrote: > As the original coder (in mid 1970s) of BFGS, CG and Nelder-Mead in > optim(), I've > been pushing for some time for their deprecation. They aren't "bad", but > we have > better tools, and they are in CRAN packages. Similarly, I believe other > optimization > tools in the core

Hi, Is there a way to turn off the automatic inclusion of "./" at the beginning of a path specified in an \SweaveInput{} instruction? I'd like to create some reusable "template modules" of Sweave code and put them in a standard directory like /Resources/Affymetrix Then the corresponding file that uses one of these would include a command like

Some time ago I (and some other CRAN package maintainers) got an email from Brian Ripley (copied below) regarding compiler warnings under gcc 8.1. In my case this concerns package PropClust (https://CRAN.R-project.org/package=PropClust). The check with gcc 8.1 (https://www.stats.ox.ac.uk/pub/bdr/gcc8/PropClust.out) produced a warning about an uninitialized variable. I cannot reproduce the warning

I concur with Avraham that capabilities need to be ensured e.g., in recommended packages. I should have mentioned that. My concern is that the core should be focused on the programming language aspects. The computational math and some of the more intricate data management could better be handled by folk outside the core. JN On 2019-03-04 9:12 a.m., Avraham Adler wrote: > On Mon, Mar 4, 2019