Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- log.pr.data<-function(lambda.s,n1,n0,lambda.star,psi0,tobs,g=T.chan){
- print("Function log.pr.data")
- print(g)
- psi.s<-boundary(lambda.s,g,psi0,tobs,n1,n0)
- -my.dbinom(n0*lambda.s,n0,lambda.star,log=TRUE)
- }
- nlm(log.pr.data,p=0.6,n1=n1,n0=n0,lambda.star=lambda.star,psi0=psi0,tobs=tobs)
- log.pr.data(l,n1,n0,lambda.star,psi0,tobs,T.chan)
- nlm(log.pr.data,p=0.6,n1=n1,n0=n0,lambda.star=lambda.star,psi0=psi0,tobs=tobs,g=T.chan)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement