similar to: Warning: number of items to replace is not a multiple of replacement length

Dear all, I ran into problems with the function "optim" when I tried to do an mle estimation of a simple lognormal regression. Some warning message poped up saying NANs have been produced in the optimization process. But I could not figure out which part of my code has caused this. I wonder if anybody would help. The code is in the following and the data is in the attachment. da <-

Dear List, I am not sure how should i optimize a log-likelihood numerically: Here is a Text book example from Statistical Inference by George Casella, 2nd Edition Casella and Berger, Roger L. Berger (2002, pp. 355, ex. 7.4 # 7.2.b): data = x = c(20.0, 23.9, 20.9, 23.8, 25.0, 24.0, 21.7, 23.8, 22.8, 23.1, 23.1, 23.5, 23.0, 23.0) n <- length(x) # likelihood from a 2 parameter Gamma(alpha,

Hello, I have a problem with Hausman test. I am performing my analysis with these commands: > library(plm) > data<-read.csv2("paolo.csv",header=TRUE) > data< pdata.frame(data,index=c("FIRM","YEAR"),drop.index=TRUE,row.names=TRUE) > RECEIV~LSIZE+LAGE+LAGE2+CFLOW+STLEV+FCOST+PGROWTH+NGROWTH+TURN+GPROF+GPROF2 >

Dear list, here is an example of stepAIC that I do not understand. The data is n=42, Lage is the only factor and there are four other variables treated as continuous. First you see the stepAIC-forward solution (fs7). The strange thing here is that apparently not all interactions are tried for inclusion, but only WQ:Lage. In particular, I think that WFL:Lage should be tried in the last two

Note the typo in your 3rd line: data < Don't know if this means anything... Bert On Feb 11, 2018 7:33 AM, "PAOLO PILI" <paolo.pili at student.unife.it> wrote: > Hello, > > I have a problem with Hausman test. I am performing my analysis with these > commands: > > > library(plm) > > data<-read.csv2("paolo.csv",header=TRUE) >

Hi: I am using R 1.7.0 on Windows. I am having trouble getting "outer" to work on one of my functions. Here is a simple example illustrating my problem: > b1 <- c(1.2,2.3) > b2 <- c(0.5,0.6) > x <- c(3e+01, 1e+02, 3e+02, 5e+02, 1e+03, 1e+04, 1e+05, 1e+06) > y <- c(2,4,2,5,2,3,1,1) > n <- c(5,8,3,6,2,3,1,1) > outer(b1,b2,FUN=bpllkd,x,y,n)

Hi R users, I have a couple of questions about some problems that I am facing with regard to numerical integration and optimization of likelihood functions. Let me provide a little background information: I am trying to do maximum likelihood estimation of an econometric model that I have developed recently. I estimate the parameters of the model using the monthly US unemployment rate series

you are right about the 3rd line but it doesn't help me for my problem. I remove the 3rd line but there is still the same problem: Error in solve.default (dvcov): the system is numerically unique: reciprocity condition value = 1.63418e-19 Paolo 2018-02-11 16:54 GMT+01:00 Bert Gunter <bgunter.4567 at gmail.com>: > Note the typo in your 3rd line: data < > > Don't

Full_Name: Lyle W. Konigsberg Version: 1.5.1 OS: Windows Submission from: (NULL) (160.36.64.99) I apparently found an error in how "deriv" puts a derivative to the screen (though the internal representation is correct). Here's what I was doing for a class example: > lnlk<-expression(15*log(p)+11*log(1-p)) > p<-15/26 > D(D(lnlk,"p"),"p") -15 *

Hello, First time poster here so let me know if you need any more information. I am trying to run an ordered probit with maximum likelihood model in R with a very simple model (model <- econ3 ~ partyid). Everything looks ok until i try to run the optim() command and that's when I get " Error in rep(1, nrow(x)) : invalid 'times' argument". I had to adapt the code from a 4

Hi all, Sorry if this has been answered already, but I couldn't find it in the archives or general internet. Is it possible to implement the gompertz distribution as survreg.distribution to use with survreg of the survival library? I haven't found anything and recent attempts from my side weren't succefull so far. I know that other packages like 'eha' and

I have a user who has several folders in his mailbox, which we can not delete, neither via IMAP nor via doveadm: root at ka-nfs-mail:~# doveadm mailbox list -u <username> | grep hbereiche | cat -v INBOX/[Fachbereiche ^M INBOX/Fachbereiche ^M INBOX/hbereiche^M INBOX/hbereiche/LAGen]^M INBOX/hbereiche/LAG^M INBOX/[Fachbereiche^M INBOX/[Fachbereiche/LAGen]^M INBOX/[Fachbereiche]^M

One more update, I tried renaming that "[Test " folder to just "Test", which creates an also not deletable folder "Test". Good news is, if I delete these folders on both replication nodes from mdbox/mailboxes and mdbox/subscriptions they seem to be really gone. I can event create the folder "Test" again and it behaves normal aka. I can delete it afterwards

Hello. I am trying to write an R function to plot the survival function (and associated hazard and density) for a Siler competing hazards model. This model is similar to the Gompertz-Makeham, with the addition of a juvenile component that includes two parameters---one that describes the initial infant mortality rate, and a negative exponential that describes typical mortality decline over the

I'm looking at what would be involved in converting from MeetMe to ConfBridge and there seems to be a lot of missing administrative things, but I hope I'm just missing it. We all know about the missing realtime linkage. That's a major nuisance, but can be worked around. More serious is that the CLI command to display users in a ConfBridge don't show the caller ID information, so

Hi, all. I'm working with published paleodemographic data (counts of skeletons that have been assigned into an age-range category, based upon observed morphological characteristics). For example, the following is the age distribution from a prehistoric cemetery in Egypt: naga <-

Dear R users, Is there a simple way to use the effect function (library(effects)) with a GEE estimated model? library(gee) library(MASS) library(effects) attach(epil) b = gee(y ~ lbase*trt + lage + V4, family=poisson, id=subject, corstr="exchangeable") plot(effect("lbase*trt", b)) # Errore in effect("lbase*trt", b) : nessun metodo applicabile per

Hello listmates, Have any of you used SquashFS? See: http://en.wikipedia.org/wiki/SquashFS We have tried it on small file sets, seems good but if you have experience using it, especially to store/archive lage volumes of files, please share your experience. Thanks. Boris.

On Thu, 3 May 2018 19:18:45 +0100 Rowland Penny via samba <samba at lists.samba.org> wrote: > On Thu, 3 May 2018 14:59:18 -0300 > "Ethy H. Brito via samba" <samba at lists.samba.org> wrote: > > > I run the pdbedit command. > > I got a lage amount of users (and groups). > > The admin of the S3 server deleted (userdel) 75 users and these are >

Hi all. I?m using the binary distribution of SAMBA 1.9.16 on SCO3.2v4.2. Last weekend we got a new server. Installed all applications and all worked fine (as expected) ... but: the win/nt clients could map the shares, but they could not browse the SAMBA-box! I checked the "browseable=" options, but everything looked ok. Increasing the debug level and stepping through the logfiles