Displaying 9 results from an estimated 9 matches for "ll3".
Did you mean:
ll
2007 Dec 06
1
suggested modification to the 'mle' documentation?
..., 12, 20, 5, 9, 8, 5, 4, 8)
# Pass data via function arguments rather than global variables
ll.5 <- function(ymax=15, xhalf=6, x., y.)
-sum(stats::dpois(y., lambda=ymax/(1+x./xhalf), log=TRUE))
(fit.5 <- mle(ll.5, start=list(ymax=15, xhalf=6),
fixed=list(x.=x, y.=y)))
ll3 <- function(lymax=log(15), lxhalf=log(6), x., y.)
-sum(stats::dpois(y.,
lambda=exp(lymax)/(1+x./exp(lxhalf)), log=TRUE))
(fit3 <- mle(ll3, start=list(lymax=0, lxhalf=0),
fixed=list(x.=x, y.=y)))
2012 May 13
2
Discrete choice model maximum likelihood estimation
...L1 <- function(b1,b2,b,t) {
P11 <- P1(b1,b,t)
P22 <- P2(b2,b,t)
L1l <- (P11=1)*(P22=0)
return(L1l)
}
# Log-likelihood function
llfn <- function(par,a,t) {
b1 <- par[1]
b2 <- par[2]
b3 <- par[3]
b <- par[4]
lL1 <- log(L1(b1,b2,b,t))
lL2 <- log(L2(b1,b2,b3,b,t))
lL3 <- log(L3(b1,b2,b3,b,t))
llfn <- (a==1)*lL1+(a==2)*lL2+(a==3)*lL3
}
est <- optim(par,llfn, method = c("CG"),control=list(trace=2,maxit=2000),
hessian=TRUE)
And when I run this code I get "cannot coerce type 'closure' to vector of
type 'double'" error....
2012 Sep 11
1
Strange result from GAMLSS
Hi Folks! Just started using the gamlss package and I tried a simple code
example (see below). Why the negative sigma?
John
> y <- rt(100, df=1)> m1<-fitDist(y, type="realline")Warning messages:1: In MLE(ll3, start = list(eta.mu = eta.mu, eta.sigma = eta.sigma, :
possible convergence problem: optim gave code=1 false convergence
(8)2: In MLE(ll4, start = list(eta.mu = eta.mu, eta.sigma = eta.sigma,
:
possible convergence problem: optim gave code=1 false convergence
(8)3: In MLE(ll4, start = list(e...
2006 Jan 25
0
Log-Likelihood 3d-plot and contourplot / optim() starting values
...tion(N,t) {
i<-seq(along=t)
n<-length(t)
s1<-sum(1/(N-(i-1)))
s2<-n/(N-(1/sum(t[i]))*(sum((i-1)*t[i])))
(s1-s2)
}
# you get this function as usual, if you set the loglikelihood equal
zero and differentiate for N
# i take the squares of ll2 in order to get the minimum easier
ll3<-function(N,t) {
ll2(N,t)^2
}
# then i do an iteration to get the mle of N
xmin<-optim(10,ll3,method="BFGS",control=list(reltol=.Machine$double.eps),t=t)
# the problem is that this method only works well, if the starting
values for optim() are very close to the real value!
# s...
2009 Jan 07
3
[LLVMdev] LLVM optmization
...0 ; 00000064H
jl SHORT $LL6 at main
; 61 :
; 62 : //WaitForMultipleObjects( num_th, (const HANDLE* )mainThread, TRUE, INFINITE);
; 63 : for( i=0; i < num_th; i++)
mov ebx, DWORD PTR __imp__WaitForSingleObject at 8
mov ebp, DWORD PTR __imp__CloseHandle at 4
xor esi, esi
$LL3 at main:
; 64 : {
; 65 : WaitForSingleObject( mainThread[i], INFINITE );
mov edi, DWORD PTR _mainThread$[esp+esi*4+420]
push -1
push edi
call ebx
; 66 : CloseHandle(mainThread[i]);
push edi
call ebp
add esi, 1
cmp esi, 100 ; 00000064H
jl SHORT $LL3 at main
; 67 : }...
2020 Feb 27
2
[PATCH] Update the 5 year logo to 10 year logo
...4WWu!<z%mUWG}=Zh{sRhO+}@bIsaGAAm{1a_;Ekd8mz=>vK64KLqhCJUX+KPCH_
zUVe9WiT=TQj1<+rzcr6<|MrDFbXUK}$z+<`Tf0A8*k6KZ97SOB)JDSuC)}?POLLJF
zMJ!do=#?UY^5WcLC{)>wZ8|vyCVa4C?`^%Bl?)MVfsB!?xjZQyaxv!oJgBEQ<yN$f
zF;leugbCb#hXMK~hA*g0j_l>GkOSAiC(8+^yW+6y=rmIPh|2ocAKODRc!!i~tHLE}
zWkrCt42*^!LL3(#FV<2$rT~=WVzPZ!C}(%qSdJ2{wY3$I04OG3J}5QHlk-TItXkP@
zEaCPt_zo!ZfrCUXX8RAd(CQkUovLI+L*5&nFcAw@(s1VRLK0_}U4Af{B~(QlSq%TA
zcc~@#sn|VWKIaW3zjgTC0bVH~akO_!fY^8CsR9oFU&b>3FY|JP+UZ<#)}0G{-moYi
zt7RgaY8@{cnU!F_8In*1O3ojGzl9i5V`X9B(OPzxzbPyG=Q}nxKzw&~K1ZK(eGg$~
z`-<SumvM!){?wFwEA}NW&...
2008 May 21
0
Errors in using gdb (PR#11496)
...S/HEi7u1q16Q7DrosYbmZE7BhOHx+hscc6/p5Pf1e
4NNN6XYZmJ+dB8MmCCU6V+9Rl7GTgdbqZHiMlJjdG+527ahFvjqClBX79EG7J+TIBVkh15l3
YZ5rYHUduR2BCqSUx525nxg+CY9x9R7FtFvQ7hjtiGdwjdaaTT2765Nd7z1qYKtr0zHaEQ92
GIlpV9mTgS8hymI/6puaGrSvBdr/e/svh+uqz2J2Qs3qENKQ0x2R89n92dnurMNjPi6BxL/e
zsiZB53jw8+zk4xYqn/p47gVx+08ns9SkCwR85Ill3mwxw8JaAcg5NovZzO+O47d8ikDk2QM
P2YyD6ctQszny7OOQA5VbMrcArRHbY+GAI42OWf3LMSVkorMZ6uoSx4V0ZvH68AAfayQxxKD
dnpPlvD7LeZTxLQfl91n8alSul0mwj/HGI/n8yxfrrioy1jQbq1OgvZX8evPXcnbRC2y6khS
1kG7PTnKCjlyQVbIdeZ93Bt4xaKO3I4A0E5vwGoO2l2910G7Ndohz1AE7fXs/kSK571HDWx1
bbqdq+4ZIpbg1nIdRZ42p2OWuUfsNjXdALQTluN/t+F8qKIItD...
2008 Jun 30
4
Rebuild of kernel 2.6.9-67.0.20.EL failure
Hello list.
I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without
modifications.
How did I try it?
Created a (non-root) build environment (not a mock )
Installed the kernel.scr.rpm and did a
rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee
prep-out.log
The build failed at the end:
Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL
Checking
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...Qiy|d&(cpq#2nUiQxw;u`jk$rEk~&<
zWIsZ|WE_AikF&ul|3osoho%sT;Omp24tsbo!G0N2=}982v9pb_ZXok6U{!e*99?Y}
z$;xsf$U at A)R5!j1<n9JgfpSNl^9Vtk-J1I~uaBRLB(P0a#g`4*i%79u9LT=%0T2Ew
zL>zhG7w4{^;b-S=ZuqrnLQSm53tmbtPi3mN%VQ5_y`4|S)@1Zkd&g#<QmBI5Wp|?^
zt8E+8Ak0V*mhVztX!r`jnWMq65m%FoHuWiw+lL34in90=Zm9GRoebXVk#O3i!{j<(
zxP%6JZqH0TRg*0qtY4MdJ1m0lAsNh8+&zGT7`X!|ETlz?X&;$`YMRkPoID4~J<8a5
z!vc2OM|gkLQ<;;X4g4%kqM|q}a>_lM3Pa<|5M5vt7tyrLjxC)UCT at k;Xh81WVO(36
zlT2h+>&ZBN<a`39(lUTWL7$G~$m0PW{!rvt%-Zv))#VKrzjaUPPBHX^ls@&eVFoZc
zj)R1%a_i|FU}E$zDa&mi#dDh-&a...