similar to: Optim Problem

Hello, When i run the code below from Weibull distribution with 30% censoring by using optim i get an error form R, which states that Error in optim(start, fn = z, data = q, hessian = T) :? ? objective function in optim evaluates to length 25 not 1 can somebody?help me remove this error. Is my censoring approach correct. n=25;rr=1000 p=1.5;b=1.2 for (i in 1:rr){ q<-c(t,cen)

Dear All, Can you help me, with the code below how do I obtain the fisher information from it. Is my q<-replicate(1000,x) the right way to do simulation. thank you. x<-rweibull(100,0.8,1.5) q<-replicate(1000,x) z<-function(p){ beta<-p[1] eta<-p[2] log1<-(n*log(beta)-n*beta*log(eta)+(beta-1)*sum(log(x))-sum((x/eta)^beta)) return(-log1) } zz<-optim(c(0.5,0.5),z) zz Chris

Your function is giving NaN's during the optimization. The R-forge version of optimx() has functionality specifically intended to deal with this. NOTE: the CRAN version does not, and the R-forge version still has some glitches! However, I easily ran the code you supplied by changing optim to optimx in the penultimate line. Here's the final output. KKT condition testing Number of

Dear all, I am a newcomer to R programming, I met the problem: I have a lot of .txt files in my directory. Firstly, I check whether the file satisfies the conditions: 1.empty 2.the "Rep" column of the file has no "useractivity_idle" or "useractivity_act" 3.even The "rep" has both of them, numbers of "useractivity_idle"==numbers of

Dear All, I have many files in my directory. I want to transfer each data into one which is readable. They have so many possibilities, i have collected(manually and visually) all possibilities and represent them as different numbers. Rep[grep('context_log',log1$Remain[1:length(log1$Date)]),]<-"2" Rep[grep('gs',log1$Remain[1:length(log1$Date)]),]<-"5"

Dear all, I am a newcomer to R programming, I met the problem: I have a lot of .txt files in my directory. Firstly, I check whether the file satisfies the conditions: 1.empty 2.the "Rep" column of the file has no "useractivity_idle" or "useractivity_act" 3.even The "rep" has both of them, numbers of "useractivity_idle"==numbers of

hello list, I'm trying to write a script that will search through a directory of trace logs for an oracle database. From what I understand new files are always being created in the directory and it's not possible to know the exact names of the files before they are created. The purpose of this is to create service checks in nagios. Because you don't know the names of the files ahead

Hello, I have being trying to estimate the parameters of the?generalized?exponential distribution. The random number generation for the GE distribution is?x<-(-log(1-U^(1/p1))/b), where U stands for uniform dist. The data i have generated to estimate the parameters is right censored and the code is given below; The problem is that, the newton-Raphson approach isnt working and i do not know what

Hi, R-users, I have a problem: Because there are few files which can't be readed into R completely, so on the following subsequence programme, I use write.table, which creates the "NA" files for those incomplete files autimatically. I don't want those NA files. My programes formats looks like: name_c<-Sys.glob("C:/Documents and Settings/lma/My

Hi, Wine documentation says that it can be installed anywhere, locate its libraries and just run. This is certainly true for self-compiled sources which could e.g. be run from their compiled location. I wonder if this is also possible with pre-built Debian/Ubuntu binary packages. Normally, Debian installs exactly one version of a package in /usr/. Is there some possibility to use multiple

I'm on freeBSD 7.0, wine 0.9.53. When I run Orcad, it can load projects, but schematics never open (unlike other users' experiences in the apps db from Orcad v9.0). You open schematics by expanding them in the 'project' treeview window, and under wine, the little "+" box next to the treeview just dissapears when you click on it instead of expanding. Steve [steve at

Hi, everyone This proposal covers the OpenProjects idea 'Port the Bigloo Scheme compiler, from Manuel Serrano at INRIA Sophia-Antipolis, to output LLVM bytecode.' There is not a popular functional language frontend yet on LLVM. To my knowledge, this is mostly because LLVM is lack of runtime support for the functional language needs, and C foreign function interface calls are relatively

Hi all, I was hoping to get some advice about a problem that I realize will be difficult to reproduce for some people. I'm running R 1.7.1 on an openMosix (Linux) cluster and have been experiencing some odd slow-downs. If anyone has experience with such a setup (or a similar one) I'd appreciate any help. Here's a simplified version of the problem. I'm trying to run the

Hello! I got something to ask..whether you can help me with the R program...i got this for example 5x4 matrix..and i want to find: ?i) mean for each row of the matrix ii) median for each column of the matrix and i need to do this using a loop function...below is my program..u try to check it for me as the output that i got is not what i desired...thanks.. data<-rnorm(20,0,1)

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20080304/9de37092/attachment.pl

Hello, ?can i implement this as 10% censored data where t gives me failure and x censored. Thank you p=2;b=120 n=50 set.seed(132); r<-sample(1:50,45) t<-rweibull(r,shape=p,scale=b) t set.seed(123);? cens <- sample(1:50, 5)? x<-runif(cens,shape=p,scale=b)? x Chris Guure Researcher, Institute for Mathematical Research UPM

https://bugzilla.netfilter.org/show_bug.cgi?id=977 --- Comment #6 from Netbug <b1b30ee4 at opayq.com> --- I upgraded to ulogd 2.0.4 & kernel 3.17.2 but these did not make a difference, here's how the output looks now; Sat Nov 1 17:07:03 2014 <5> ulogd.c:843 building new pluginstance stack:

https://bugzilla.netfilter.org/show_bug.cgi?id=1218 Bug ID: 1218 Summary: ULOGD PCAP Plugin Missing Ethernet Headers Product: ulogd Version: SVN (please provide timestamp) Hardware: All OS: All Status: NEW Severity: blocker Priority: P5 Component: ulogd Assignee:

> good morging, > > I need to install R-2.x on my host: > > Linux, 2.6.5-7.308-smp #1 SMP Mon Dec 10 11:36:40 UTC 2007 x86_64 > x86_64 x86_64 GNU/Linux > > now I have check the packages on your document: > http://cran.r-project.org/bin/linux/suse/ReadMe.html > but I have a problem whit xorg-x11-lib, in my host I have installed > package:

HI,Barry. Thanks a lot for your details reply. I have already rewritten the programs. Now it totally works. Your comment: " Note the use of comments and breaking the code up into small independent, testable functions" is really Valuable! Thanks again. Hi,David, You are right, My previous code was really horrible. Its not good to make multiple codes at the same time, sometime I feel