search for: p0l

...> >>>On Sat, Mar 9, 2013 at 8:57 PM, arun <smartpink111 at yahoo.com> wrote: >>> >>>After running this: >>>>I am getting: >>>> >>>> >>>>search(0.01,0.02,0.15,0.20) >>>> >>>>? m1 n1 m n cterm2_P0L cterm2_P0H N EN???????? BH??????? BL??????? AH??????? AL >>>>1? 2? 2 4 4? 0.8145062? 0.6634204 8? 8 0.10011292 0.3164062 0.3365796 0.1854938 >>>>2? 2? 2 5 4? 0.7737809? 0.7961045 9? 9 0.20022583 0.2373047 0.2038955 0.2262191 >>>>3? 3? 2 5 4? 0.7737809? 0.7961045...

...ue ese error no aparece en las versiones anteriores. #Doble Pareto con p=1/2 y mismo lambda (simétrica) probando<-function(data=NULL,p0=c(1,1),method=1,iters=10000,print.level=0){ x<-data n<-length(x) x.pos<-x[x>=0] n.pos<-length(x.pos) x.neg<-x[x<0] n.neg<-length(x.neg) p0l<-p0[1] p0s<-p0[2] logL1<-function(p){ l1<-exp(p[1]) l2<-l1 s<-exp(p[2]) -(n.neg*log(l1)+n.pos*log(l2) -l1*sum(log(1-x.neg/s))-l2*sum(log(1+x.pos/s)) -sum(log(1+abs(x)/s)) -n*log(s)+n*log(0.5) ) } if (method==1){ p0<-c(log(p0l),log(p0s)) fit<-nlm(logL1,p=p0,hes...

Thank you very much for your reply. Your code work well with this example. I modified a little to fit my real data, I got an error massage. Error in split.default(x = seq_len(nrow(x)), f = f, drop = drop, ...) : Group length is 0 but data length > 0 On Thu, Jan 31, 2013 at 12:21 PM, arun kirshna [via R] < ml-node+s789695n4657196h87@n4.nabble.com> wrote: > Hi, > Try this: >

depends on "autoconf: regenerate configure scripts with 4.4 version" This series removes some of the really old deadwood from the tools build and makes some other things which are on their way out configurable at build time with a default depending on how far down the slope I judge them to be. * nuke in tree copy of libaio * nuke obsolete tools: xsview, miniterm, lomount & sv *