similar to: Problems with rank()

Dear Sirs, Could you please be so kind as to send us some information on residuals in multinomial logistic models? Is it possible to use R software? We thank you in advance. Sincerely yours Luciana Alves,MSc Beatriz Leimann, MD -- Luciana Correia Alves Doutoranda em Sa??de P??blica ENSP - Fiocruz

Hello, I've run into a problem which is both maddening and rather hard to replicate, therefore I wondered if someone might know of a plausible explanation for it. I couldn't find anything in the archives, though maybe I'm searching for the wrong thing. I'm calling some C code using .C, and get the vector I'm interested in back as the 7th location in the returned list.

This is a follow-up question for PR#15052 <http://bugs.r-project.org/bugzilla3/show_bug.cgi?id=15052> There is another thing I would like to discuss wrt how nlminb() should proceed with NAs. The question is: What would be a successful way to deal with an evaluation point of the objective function where the gradient and the hessian are not well defined? If the gradient and the hessian both

Hi R-users, Recently, I use maxNR function to find maximizer. I have error appears as follows Error in maxNRCompute(fn = fn, grad = grad, hess = hess, start = start, : NA in the initial gradient My code is mu=2 s=1 n=300 library(maxLik) set.seed(1004) x<-rcauchy(n,mu,s) loglik<-function(mu) { log(prod(dcauchy(x,mu,s))) } maxNR(loglik,start=median(x))$estimate Does anyone know how

Hi! I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb I use the function optmin optmin(start,fn,gr,method="L-BFGS-B", lower, upper, hess,...) But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn Then I used optmin(start,fn,gr,method="BFGS", hess, ...) But then I get the Error in optmin ...: initial value in vmmin is not

Greetings, I was wondering if somebody is working on a solution to the NAT/SIP-issues? It seems to me that the problem has been identified, is that correct? Just hoping that someone with more skills will provide us with a solution sooner or later... Regards, Stig -------------- next part -------------- An HTML attachment was scrubbed... URL:

To see what I mean, try x11() par(mfrow=c(2,1)) x <- 1:100 plot(x, sin(pi*x/100)) options(digits= 20) plot(x, sin(pi*x/100)) The problem really is in axis(.) which uses format(.) for the axis labels. I think the result of axis(.) should NOT be influenced by options()$digits. Should I change it? -- Martin

Haven't seen email since the 29th.. just testing. -------------- next part -------------- A non-text attachment was scrubbed... Name: mhess.vcf Type: text/x-vcard Size: 288 bytes Desc: not available Url : http://lists.digium.com/pipermail/asterisk-users/attachments/20050801/96555713/mhess.vcf

Full_Name: Lara Jamieson Version: 1.1.1 OS: Linux Submission from: (NULL) (131.111.20.203) I've checked for other reports on this and can't find anything. It seems that the .Random.seed isn't updated properly after using rm(.Random.seed). Mu understanding from the help file indicates that a new 'random' seed will be generated. It appears that sometimes the old one is just

In that case, I can't work out why the first model fails but not the second. I would start looking at "Data" to see what it contains. if: object2 <- polr(Inc ~ Training ,Data,Hess = T,method = "logistic" ) works, the problem may be with the "Adopt" variable. Jim On Tue, Mar 20, 2018 at 10:55 AM, Willy Byamungu <wmulimbi at email.uark.edu> wrote: >

Hi, I'm getting an error message using polr(): Error in optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) : initial value in 'vmmin' is not finite The outcome variable is ordinal and factored, and the independant variable is continuous. I've checked the source code for both polr() and optim() and can't find any variable called

Hello R users, I am a researcher at the University of Michigan looking for a solution to an R problem. I have loaded my data in from a mysql database and it looks like this > data ds c1 c2 1 2010-04-03 100 0 2 2010-04-30 11141 15 3 2010-05-01 3 16 4 2010-05-02 7615 14 5 2010-05-03 6910 17 6 2010-05-04

hello everybody solved the problem with summary, now I have another one eg I estimate > try.op <- polr( > as.ordered(sod.sit.ec.fam) ~ > log(y) + > log(1 + nfiglimin) + > log(1 + nfiglimagg) + > log(ncomp - nfiglitot) + > eta + > I(eta^2) + >

Dear all, I have 5 plots that I would like to include in a single figure, spread over two rows. If I use mfrow=c(2,3), and produce my plots one after the other, I will end up with three plots in the first row, and 2 in the second row, which is what I want. However, I would like the two plots in the second row to be moved to the centre of that row, rather than occupying the two left-most cells.

I found this great search source("http://www.statslab.cam.ac.uk/~djw1005/Stats/Interests/search.R") helpHTML() Has (or will) this become the standard search method? Is R 'Free Software'? The dependence on Java seems a bit of a pain for 'freeness'. Did the above make it into CRAN? Cheers, Dan.

Hello Mark, Did you ever solve this problem? I saw your post but no replies. I'm in the exact same situation. If you have any pointers on where to get this I'd appreciate it. - Paul Hess >From: "Mark Mestdagh" <mestdagh@dma.be> >To: "samba digest" <samba@samba.anu.edu.au> >Subject: TCP/IP >Date: Thu, 26 Mar 1998 21:40:59

There have been a number of posts to this list by people having trouble with the HTML search engine. Often these troubles are caused by incorrect setups (user hasn't installed Java properly, or Java is disabled, or Javascript is disabled). Sometimes the trouble persists even when Java is installed properly. I have written an alternative HTML search engine, which is based on Javascript rather

Dear Sir, I get an error message when I use polr() in MASS package. My data is "ord.dat". I made "y" a factor. y y1 y2 x lx 1 0 0 0 3.2e-02 -1.49485 2 0 0 0 3.2e-02 -1.49485 3 0 0 0 1.0e-01 -1.00000 4 0 0 0 1.0e-01 -1.00000 5 0 0 0 3.2e-01 -0.49485 6 0 0 0 3.2e-01 -0.49485 7 1 1 0 1.0e+00 0.00000 8 0 0 0 1.0e+00 0.00000 9 1 1 0

Hi, I have a dataset that I am trying to analyze and plot as an ordered logistic regression (y = ordinal categories 1-3, x = continuous variable with values 3-9). First is a problem with cdplot: Produces a beautiful plot, with the "right" trend, but my independent factor values are transformed. The factor has values from 3-9, but the plot produces an x-axis with values from 20-140.

I am trying to forecast volatility 2 periods forward using a ARCH(1) model: predict(garch(fit2,order=c(0,1),n.ahead=2)) ***** ESTIMATION WITH ANALYTICAL GRADIENT ***** Error in qr(com.hess$hess, ...) : unused argument(s) (n.ahead ...) What did I do wrong? Thank you. Best regards, Peter Arnold, CFA President PRA Investment Counsel, Inc. 704-341-8193 www.prainvestment.com